A question about DEXSeq of estimateDispersions
2
0
Entering edit mode
liushang • 0
@liushang-16396
Last seen 3.3 years ago

HI, I meet an error when I run a 119311* 749 dataframe:

       error: matrix multiplication: incompatible matrix dimensions: 1498x751 and 0x1
       Error in fitBeta(ySEXP = ySEXP, xSEXP = xSEXP, nfSEXP = nfSEXP, alpha_hatSEXP = alpha_hatSEXP,  :
       matrix multiplication: incompatible matrix dimensions: 1498x751 and 0x1
       Calls: estimateDispersions ... tryCatch -> tryCatchList -> tryCatchOne -> <Anonymous>
       Execution halted

using the command:

      matall=cbind(mat_filter,mat_used) ## the 119311 * 749 datatframe
      sampleData <- data.frame(condition = c(rep('filter',length(colnames(mat_filter))),rep('used',length(colnames(mat_used)))))
      groupID <- as.character(rowann$V3)
      featureID <- as.character(rowann$V1)
      design <- formula( ~ sample + exon + condition:exon  )
      dim(matall)
      dim(sampleData)
      length(groupID)
      all <- DEXSeqDataSet( matall, sampleData, design,featureID, groupID)
      dxd <- all
      dxd = estimateSizeFactors( dxd )
      dxd = estimateDispersions( dxd )

I guess that the error results from the large dataframe, so I creat a 119311 * 10 dataframe from the 119311 * 749 datatframe, and then it works and didnot produce the error.

Any help is greatly appreciated.

dexseq • 473 views
ADD COMMENT
0
Entering edit mode
glmazzo • 0
@glmazzo-12778
Last seen 3.2 years ago

HI liushang, I am having the same issue. Any luck with it?

Thank you. Gianluca

ADD COMMENT
0
Entering edit mode
glmazzo • 0
@glmazzo-12778
Last seen 3.2 years ago

HI liushang, I am having the same issue. Any luck with it?

Thank you. Gianluca

ADD COMMENT

Login before adding your answer.

Traffic: 441 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6