Manually feed DESeq with dispersions
1
0
Entering edit mode
@ruhollahmb-20544
Last seen 5.5 years ago

I am dealing with a very special problem: I have to perform DESeq on the count data that belong to a single gene (only ONE gene, i.e. the count data has one row & 100 columns). The dispersion value is already available for the gene. (Actually, the dispersions were calculated by DESeq over a master count data with 19k genes )

Question: How can I use DESeq on this special count data?! I have to somehow feed the DESeq with the already-known dispersion value, becauseDESeq cannot calculate dispersion when there's only one gene in the count data.

deseq2 • 741 views
ADD COMMENT
0
Entering edit mode
@mikelove
Last seen 1 day ago
United States

You can just use dispersions(dds) <- ... to supply your own dispersion estimates rather than estimating them with DESeq2. Instead of running DESeq() you would run the individual steps, see ?DESeq.

ADD COMMENT

Login before adding your answer.

Traffic: 911 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