我有一个4列多行的大数据框架(附件中有一个例子)。
#what I have
Arm <- c("5prime","3prime","5prime","CoMature","3prime","5prime","3prime","3prime")
Family <- c("LET-7","LET-7","LET-7","MIR-10","MIR-103","MIR-124","MIR-124","MIR-124")
Sequence <- c("ATCGGCA","ATGCTAC","ATCGGCA","ATCGTTT","TGAGGAG","TGATCAG","AATTCAG","AATTCAG")
Star_seq <- c("TTCAGGT","TATACTG","TTCAGGT","GAGATCA","CAAAAGC","CACATGC","AATATGC","AATATGC")
my_data_frame <- data.frame(Arm,Family,Sequence,Star_seq)我想要做的基本上是,对于Family列中的每个i,计算Arm列中出现'5prime‘、'3prime’或'CoMature‘的次数。然后,对于最频繁的一个('5prime','3prime‘或'CoMature'),取第三和第四列。总而言之,我需要一个最终文件,它在Family列中显示每个i的最频繁的arm (在第一行中),并在第三和第四列中显示它们的相对序列。
#what I want as output
five_prime_counts <- c("2","0","0","1")
three_prime_counts <- c("1","0","1","2")
CoMature_counts <- c("0","1","0","0")
Arm_new <- c("5prime","CoMature","3prime","3prime")
Family_new <- c("LET-7","MIR-10","MIR-103","MIR-124")
Sequence_new <- c("ATCGGCA","ATCGTTT","TGAGGAG","AATTCAG")
Star_seq_new <- c("TTCAGGT","GAGATCA","CAAAAGC","AATATGC")
my_data_frame_new <- data.frame(five_prime_counts,three_prime_counts,CoMature_counts,Arm_new,Family_new,Sequence_new,Star_seq_new)发布于 2020-01-28 20:45:40
我们可以为每个Family和Arm添加一个计数变量,为每个Family中的最大计数获得相应的Sequence、Star_seq和Arm值,并获得宽格式的数据。
library(dplyr)
my_data_frame %>%
add_count(Family, Arm) %>%
group_by(Family) %>%
mutate(Sequence = Sequence[which.max(n)],
Star_seq = Star_seq[which.max(n)],
Arm_new = Arm[which.max(n)]) %>%
distinct() %>%
tidyr::pivot_wider(names_from = Arm, values_from = n, values_fill = list(n = 0))
# Family Sequence Star_seq Arm_new `5prime` `3prime` CoMature
# <fct> <fct> <fct> <fct> <int> <int> <int>
#1 LET-7 ATCGGCA TTCAGGT 5prime 2 1 0
#2 MIR-10 ATCGTTT GAGATCA CoMature 0 0 1
#3 MIR-103 TGAGGAG CAAAAGC 3prime 0 1 0
#4 MIR-124 AATTCAG AATATGC 3prime 1 2 0https://stackoverflow.com/questions/59948871
复制相似问题