> #輸出RPKM值
> geneRT=read.table("geneLength.txt",sep="\t",header=T,check.names=F)
> geneRT=geneRT[,c(1,2,2)]
> geneRT=as.matrix(geneRT)
> rownames(geneRT)=geneRT[,1]
> exp1=geneRT[,2:ncol(geneRT)]
> dimnames1=list(rownames(exp1),colnames(exp1))
> data1=matrix(as.numeric(as.matrix(exp1)),nrow=nrow(exp1),dimnames=dimnames1)
> data1=avereps(data1)
> data1=data1[row.names(y$counts),]
Error in data1[row.names(y$counts), ] : 下標出界
> geneLen=data1[,1]
> rpkm=rpkm(y,geneLen)
Warning message:
In y/gene.length.kb : 長的物件長度不是短的物件長度的整倍數
> rpkm=rbind(id=colnames(rpkm),rpkm)
> write.table(rpkm,file="TCGA_rpkm.txt",sep="\t",quote=F,col.names=F)
大神們，我該怎么修改啊。
轉載請註明出處，本文鏈接：https://www.uj5u.com/qita/23084.html
標籤：非技術區