Hello,
I am trying to use DEXSeq for the first time and I am encountering this error.  I am closely following this paper: "Inferring differential exon usage in RNA-Seq data with the DEXSeq package "  I am using my own data.

The files correspond to:

countFiles = 37 samples with counts generated as suggested in the paper.

flattenedfile = GTF annotation file as provided by Illumina in iGenomes.

Could the problem be related to the very large GTF file?

Thanks

> dxd = DEXSeqDataSetFromHTSeq(
+     countFiles,
+     sampleData=sampleTable,
+     design= ~ sample + exon + condition:exon,
+     flattenedfile=flattenedFile )

Error in FUN(X[[i]], ...) : subscript out of bounds

I flattened the GTF file via that py script as specified in that paper.  I used the same flattened file when running HTSeq.

This is from the paper:

python /path/to/library/DEXSeq/python_scripts/dexseq_count.py
             Dmel.BDGP5.25.62.DEXSeq.chr.gff untreated1.sam untreated1fb.txt

Hi Alejandro,

I'll proceed to check that.  I'll let you know.  Thanks

1. All the count files have exactly the same number of lines as reported by "wc -l"  A total of 66028.
2. This is the command I used to create the gff file: python /data5/bsi/adhoc/s200555.R-infrastructure/install/R-3.2.3/new/DEXSeq/python_scripts/dexseq_prepare_annotation.py genes.gtf genes.gff
3. I did not see this number "66028" in the gff file.  Here are just the top lines of the file:

   KN196472.1    dexseq_prepare_annotation.py    aggregate_gene    14223    14535    .    -    .    gene_id "ENSG00000281764.1"

   KN196472.1    dexseq_prepare_annotation.py    exonic_part    14223    14535    .    -    .    transcripts "ENST00000630776.1"; exonic_part_number "001"; gene_id "ENSG00000281764.1"

   KN196472.1    dexseq_prepare_annotation.py    aggregate_gene    81784    81917    .    +    .    gene_id "ENSG00000281299.1"

   KN196472.1    dexseq_prepare_annotation.py    exonic_part    81784    81917    .    +    .    transcripts "ENST00000626225.1"; exonic_part_number "001"; gene_id "ENSG00000281299.1"

   KN196472.1    dexseq_prepare_annotation.py    aggregate_gene    111604    133487    .    +    .    gene_id "ENSG00000281518.2"

4. Finally, the original gtf file is the annotation file provided in iGenomes and it is located here in the iGenomes directory ( Homo_sapiens/NCBI/GRCh38/Annotation/Genes.gencode )

How can I make files available?  Is there an ftp location where I can put the files?  I believe that I can ftp but I need to check with firewall restrictions.

Hi Alejandro,

I have uploaded the files in dropbox.  I'll need your e-mail address to share that folder with you.

This is the command I used to create the sampleTable.

sampleTable = data.frame(
  row.names = c(  "ME26501", "ME26502", "ME26503", "ME26504", "ME26505", "ME26506", "ME26507", "ME26508", "ME26509", "ME26510", "ME26511", "ME26512", "ME26513", "ME26514", "ME26515", "ME26516", "ME26517", "ME26518", "ME26519", "ME26520", "ME26521", "ME26522", "ME26523", "ME26524", "ME26525", "ME26526", "ME26527", "ME26528", "ME26529", "ME26530", "ME26531", "ME26532", "ME26533", "ME26534", "ME26535", "ME26536", "ME26537" ),
  condition = c(  "yes", "yes", "yes", "yes", "yes", "yes", "yes", "yes", "yes", "yes", "yes", "yes", "yes", "yes", "yes", "no", "no", "no", "no", "no", "no", "no", "no", "no", "no", "no", "no", "no", "no", "no", "no", "no", "no", "no", "no", "no", "no" ),
  libType = c(  "paired-end", "paired-end", "paired-end", "paired-end", "paired-end", "paired-end", "paired-end", "paired-end", "paired-end", "paired-end", "paired-end", "paired-end", "paired-end", "paired-end", "paired-end", "paired-end", "paired-end", "paired-end", "paired-end", "paired-end", "paired-end", "paired-end", "paired-end", "paired-end", "paired-end", "paired-end", "paired-end", "paired-end", "paired-end", "paired-end", "paired-end", "paired-end", "paired-end", "paired-end", "paired-end", "paired-end", "paired-end" ) )

and this is for dxd:

dxd = DEXSeqDataSetFromHTSeq(
  countFiles,
  sampleData=sampleTable,
  design= ~ sample + exon + condition:exon,
  flattenedfile=flattenedFile )