[BioC] question about DeSeq

[Martin Morgan]

On 12/22/2010 05:24 AM, Simon Anders wrote:
> Hi
> 
> On 12/22/2010 08:23 AM, vasu punj wrote:
>> What kind of input DEseq can take Can it take Cufflinks out put or we
>> have to use FKPM conversion.
> 
> DESeq needs integer counts, i.e., just the raw number of reads that
> align to a feature, without any normalization for sequencing depth or
> transcript length. Hence, the FPKM values reported by cufflinks are not
> suitable. the easiest is to use HTSeq-count to get a table of raw counts.
> 
> http://www-huber.embl.de/users/anders/HTSeq/
> http://www-huber.embl.de/users/anders/HTSeq/doc/count.html

It is as straight-forward and flexible to count overlaps with
Bioconductor IRanges countOverlaps, with the added advantages of
integration with annotation resources (e.g., GenomicFeatures
TranscriptDb) and familiar (to those on the Bioconductor mailing list)
programming environment.

Martin

> 
>   Simon
> 
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at r-project.org
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives:
> http://news.gmane.org/gmane.science.biology.informatics.conductor


-- 
Computational Biology
Fred Hutchinson Cancer Research Center
1100 Fairview Ave. N. PO Box 19024 Seattle, WA 98109

Location: M1-B861
Telephone: 206 667-2793