Entering edit mode
The example DEXSeqDataSet code gives and error in version 1.17.4 and 1.16.3.
Thanks
Alex
countData <- matrix( rpois(10000, 100), nrow=1000 ) sampleData <- data.frame( condition=rep( c("untreated", "treated"), each=5 ) ) design <- formula( ~ sample + exon + condition:exon ) groupID <- rep( paste0("gene", 1:10), each= 100 ) featureID <- rep( paste0("exon", 1:10), times= 100 ) DEXSeqDataSet( countData, sampleData, design, featureID, groupID )
Error:
converting counts to integer mode Error in `$<-.data.frame`(`*tmp*`, "dispersion", value = NA) : replacement has 1 row, data has 0