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