2022-09-29 summarise group_by
2022-09-29 本文已影响0人
学习生信的小兔子
#先看看基因是否有重复
table(exprSet$symbol)
#可看出基因有重复
#########接下来根据最大值去重#########
exp=exprSet %>%
mutate(rowmean=rowMeans(.[,-1])) %>%
arrange(rowmean) %>%
distinct(symbol,.keep_all = TRUE) %>%
arrange(-rowmean)
#########探针平均值#########
exp1=exprSet %>%
group_by(symbol) %>%
summarise_all(mean)
### 3.重复的探针都不要(比较极端)
data4 <- exprSet%>%
group_by(symbol) %>%
mutate(n=n()) %>%
filter(n==1) %>%
dplyr::select(-n)