A question about NgList=GetNg(IsoNames, IsosGeneNames)

13 views
Skip to first unread message

Aimin Yan

unread,
Dec 27, 2020, 10:33:41 AM12/27/20
to EBSeq users
I have a question about isoform DE analysis when using EBseq

From the output of RSEM,  I calculated the number of isoform for each gene, and find the max number is 239, so I set 239 in GetNg like this,  then do the following: 

NgList=GetNg(IsoNames, IsosGeneNames,TrunThre = 239)
Sizes=MedianNorm(IsoMat)  EBOut=EBTest(Data=IsoMat,NgVector=NgList$IsoformNgTrun,Conditions=as.factor(rep(c('trt','ck'),each=3)),sizeFactors=Sizes, maxround=5)

but I got this :

Error in DataList[[i]][GoodData[GoodData %in% rownames(DataList[[i]])],  : 
  incorrect number of dimensions

Does anyone know the reason for this error?

What I want is to use the exact number of isoform in each gene to set uncertainty groups.

Aimin




Reply all
Reply to author
Forward
0 new messages