rsem-run-ebseq --ngvector failed with "subscript out of bounds".

deweylab / RSEM

RSEM: accurate quantification of gene and isoform expression from RNA-Seq data

http://deweylab.biostat.wisc.edu/rsem/

GNU General Public License v3.0

404 stars 117 forks source link

rsem-run-ebseq --ngvector failed with "subscript out of bounds". #70

Open wangyugui opened 7 years ago

wangyugui commented 7 years ago

Hi.

rsem-run-ebseq --ngvector failed with "subscript out of bounds".

rsem version 1.3.0

#rsem-run-ebseq --ngvector /usr/bio-ref/GRCh38.p10/rsem.bowtie2.ngvec /biowrk/TCM.bio.RNA/rsem.spec/trans.matrix.txt 12,11,11 /biowrk/TCM.bio.RNA/rsem.spec/trans.matrix.ebseq

Loading required package: testthat
Removing transcripts with 100 th quantile < = 0
109861 transcripts will be tested
Error in Data[NameList[[i]], ] : subscript out of bounds
Calls: EBMultiTest -> sapply -> lapply -> FUN

wangyugui commented 7 years ago

I have re-run rsem-generate-ngvector.

and there is no relationship to rsem-generate-ngvector.

wangyugui commented 7 years ago

I can split the trans.matrix.txt to small files(with heading line) and then run rsem-generate-ngvector on every file.

Will the result be same?