# output RPKM value
GeneRT=read. Table (" geneLength. TXT ", sep="\ t", the header=t, check the names=F)
,2,2 geneRT=geneRT [, c (1)]
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),]
The Error in data1 [row names (y $counts),] : subscript out
GeneLen=data1 [1]
RPKM=RPKM (y, geneLen)
Warning message:
In y/gene. Length. KB: short length of the object is not a multiple of that of the object length
RPKM=rbind (id=colnames (RPKM), RPKM)
Write. The table (RPKM file="TCGA_rpkm. TXT", sep="\ t", quote=F, col. Names=F)

The great god, how can I modify,

CodePudding user response:

I also met the same problem, can you tell me you solved?