Off topic:DEXSeqDataSet example code
0
0
Entering edit mode
pensona ▴ 10
@pensona-9299
Last seen 8.4 years ago
United States

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
DEXSeq • 666 views
ADD COMMENT
This thread is not open. No new answers may be added
Traffic: 907 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