R语言中实现统计基因的缺失率
001、
gene1 <- c(0,0,1,0) gene2 <- c(1,1,0,1) gene3 <- c(0,1,0,1) dat <- data.frame(gene1, gene2, gene3) ## 生成测试数据, 0表示基因缺失, 1表示基因存在 rate_0 <- function(x) ## 生成统计函数 { count = sum(x == 0) return(count/length(x)) } dat apply(dat, 2, rate_0) ## 借助apply函数输出所有基因的缺失率