R: calculate prevalence


统计丰度>0.0001的细菌的出现率

target="China"
tmp=df[which(df$Country==target),]
n_row=nrow(tmp)
rate=c()
for(i in 3:529)
{
    rate=c(rate, nrow(tmp[tmp[,i]>0.0001,])/n_row)
}
out_china=data.frame(Country=rep("China",527),
                     Genome=colnames(tmp)[3:529],
                     Prevalence=rate)
©著作权归作者所有,转载或内容合作请联系作者
平台声明：文章内容（如有图片或视频亦包括在内）由作者上传并发布，文章内容仅代表作者本人观点，简书系信息发布平台，仅提供信息存储服务。

推荐阅读更多精彩内容