您好,登錄后才能下訂單哦!
這篇文章給大家分享的是有關怎么使用R語言篩選基因的內容。小編覺得挺實用的,因此分享給大家做個參考,一起跟隨小編過來看看吧。
篩選基因
myfpkm<-read.table("All_gene_fpkm.txt",header=TRUE,comment.char="",sep = "\t",check.names=FALSE,row.names=1) head(myfpkm) myfpkm[order(rowSums(myfpkm),decreasing=T)[1:5000],] #篩選表達量高的前5000個基因 myfpkm[rowSums(myfpkm)>1,] #篩選掉表達量低的基因 minRowFPKM=rowMeans(myfpkm)>2 #按平均數篩選 minNumFPKM=rowSums(myfpkm>0)>10 #表達量不為0的樣品個數篩選 myfpkm=myfpkm[minRowFPKM & minNumFPKM,] #聯合一下
感謝各位的閱讀!關于“怎么使用R語言篩選基因”這篇文章就分享到這里了,希望以上內容可以對大家有一定的幫助,讓大家可以學到更多知識,如果覺得文章不錯,可以把它分享出去讓更多的人看到吧!
免責聲明:本站發布的內容(圖片、視頻和文字)以原創、轉載和分享為主,文章觀點不代表本網站立場,如果涉及侵權請聯系站長郵箱:is@yisu.com進行舉報,并提供相關證據,一經查實,將立刻刪除涉嫌侵權內容。