Question: A question about DEXSeq of estimateDispersions
0
gravatar for liushang
3 months ago by
liushang0
liushang0 wrote:

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 • 102 views
ADD COMMENTlink modified 6 weeks ago by glmazzo0 • written 3 months ago by liushang0
Answer: A question about DEXSeq of estimateDispersions
0
gravatar for glmazzo
6 weeks ago by
glmazzo0
glmazzo0 wrote:

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

Thank you. Gianluca

ADD COMMENTlink written 6 weeks ago by glmazzo0
Answer: A question about DEXSeq of estimateDispersions
0
gravatar for glmazzo
6 weeks ago by
glmazzo0
glmazzo0 wrote:

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

Thank you. Gianluca

ADD COMMENTlink written 6 weeks ago by glmazzo0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 454 users visited in the last hour