This is a test version of Biostars. For the public version, visit https://www.biostars.org.
DEXSeq for exon counting & package name for human?

Hi,

I want to get exon count so I tried HTSeq-count first (with exon_id option) but I found that DEXSeq is better for this purpose.

I followed the manual and got a result of reads counting as a *.txt file by running the script python_count.py.

So my questions are,

  1. If I want to get exon counts only, is it a right approach? (is the txt file from python_count.py the final output for this case?)
  2. In the manual, they used the pasilla dataset so they used the package pasilla (p6 of the manual). However, my samples are about human and I want to compare test (knocked down) with control. So I wonder what is the right package name for this case.
rna-seq dexseq
• 3,222 views • link
updated by Ram Ram • written by tinmad • tinmad

1 answer

  1. Yes, the DEXseq script is a wrapper around htseq-count that produces appropriate numbers. Use it rather than htseq-count
  2. Your data is not in a package. See section 2.6 of the vignette for examples of loading data.
• 0 views • link
written by Devon Ryan Devon Ryan

Thank you for your help!
For the second question, then is it okay if I follow the same package with the vignette?
or Do I have to process any other alternative way?

• 0 views • link
written by tinmad • tinmad

I'm not sure what you mean by "follow the same package". If you mean follow DEXSeq's vignette, then yes you should do that. If you mean, "should I still load the pasilla package?", then no, that's just an example dataset. More specifically, you'll need to use the DEXSeqDataSetFromHTSeqCount() function.

• 0 views • link
written by Devon Ryan Devon Ryan

Sorry for the confusion, yes, it was about pasilla package!

About your specific comment, do you mean that I have to use DEXSeqDataSetFromHTSeq() function without loading any package? I'm just confused that what should I type instead of inDir = system.file("extdata", package="pasilla") from the vignette. If I don't need to load any package, do I have to type something like inDir = system.file("mysample")?

• 0 views • link
updated by Ram Ram • written by tinmad • tinmad

Yes, you MUST use DEXSeqDataSetFromHTSeq() to load your data. You will not use any system.file() commands, only DEXSeqDataSetFromHTSeq().

• 0 views • link
updated by Ram Ram • written by Devon Ryan Devon Ryan

Log in to answer this question.

More posts like this