5 如何过滤掉表达丰度低的基因?比如原始表达值小于10,超过25%样本就把它剔除掉;老师能不能给个R语言代码,非常感谢!

请先 登录 后评论

1 个回答

祝让飞 - 生物信息工程师

假设表达矩阵data,行为基因,列为样本,则代码如下:

data=data[which(apply(data,1,function(x){return(sum(x>10))})>ncol(data)*0.25),]

请先 登录 后评论