Search
Question: DESeq2 for single gene
0
14 months ago by
asd0
asd0 wrote:

I have 950 cases and ~60k genes in a matrix and I would like to test differential expression for one gene.  The conditions are true only for that gene.

If the countData parameter of DESeqDataSetFromMatrix  function has a matrix only with one row, it doesnt work.

How can I do it?

modified 14 months ago by Michael Love18k • written 14 months ago by asd0
0
14 months ago by
Michael Love18k
United States
Michael Love18k wrote:

Can you add more detail? The whole point of DESeq2 is to perform more rigorous testing across many genes by sharing information across those genes (take a look at the DESeq2 paper, and the methods section).

If you only want to test the counts of one gene, you shouldn't use DESeq2.

If I have ~60k gene and filter the matrix to contain just 10 gene and after that I run DESeq(dds). Does the filtering affect the result of the method?

To clarify the original question: I assume that one gene in my count dataset has different distribution if I consider given cases. I would have liked to know how can I test for this with DESeq2, instead of U-test.