假设表达矩阵data,行为基因,列为样本,则代码如下:
data=data[which(apply(data,1,function(x){return(sum(x>10))})>ncol(data)*0.25),]
如果觉得我的回答对您有用,请随意打赏。你的支持将鼓励我继续创作!