提取半数以上样本中CPM大于2的基因

CPM<-read.table("~/Desktop/OneDrive/GTEx/cpm.tsv",header = T,row.names = 1)#读取文档

TPM<-read.table("~/Desktop/OneDrive/GTEx/pm.tsv",header = T,row.names = 1)

#设置一个计数用的函数

MyFunction<-function(x){

  ifelse(x>2,1,0)

}

#设置一个空的变量储存基因

keep_genes<-vector()

row_names<-row.names(CPM)

#Run一个循环，这里样本总数为130也可以用ncol(CPM)来实现

for (i in row_names){

  cpm<-CPM[i,]

  count<-apply(cpm,2,MyFunction)

  a<-sum(count)

  if (a>=65)

  keep_genes<-c(keep_genes,i)

}

#最终还是要回归到TPM中，提取TPM中符合条件的基因

keep_TPM<-TPM[keep_genes,]

write.csv(keep_TPM,"tpm_cpm>2.csv")