Filter the columns

单细胞测序count table QC要filter 掉low quality cells(columns)
比如:
cells with fewer than 500 detected genes --excluded

count <- function(x){ length((which(x != 0))) }
obj= which(apply(data,2,count)>500)
data=data[,obj]

你可能感兴趣的:(Filter the columns)