Entering edit mode
liushang
•
0
@liushang-16396
Last seen 5.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.