Results analysis of differential gene expression using DESeq2 package
1
0
Entering edit mode
N • 0
@c1d5a219
Last seen 2.0 years ago
United Arab Emirates

I wanted to do differential expression analysis of about 10,340 genes collected from smoking vs non smoking mothers and the upregulated and downregulated genes results just does not make sense with both being 0. Is there something I don't understand or something I'm missing? Moreover, the data is pure counts data without any normalization so could that also be the issue? Plus do I need to adjust the alpha value because most of the adjusted p-value are really high.

Code should be placed in three backticks as shown below

dseq <- DESeqDataSetFromMatrix(countData=exp30, 
                              colData=meta30, 
                              design= ~ smoking_status)
dseq
dseq<- DESeq(dseq)
result <-results(dseq)
result
summary(result)

out of 10340 with nonzero total read count adjusted p-value < 0.1 LFC > 0 (up) : 0, 0% LFC < 0 (down) : 1, 0.0097% outliers [1] : 0, 0% low counts [2] : 0, 0% (mean count < 6) [1] see 'cooksCutoff' argument of ?results [2] see 'independentFiltering' argument of ?results

RnaSeqSampleSize statistics DESeq2 • 972 views
ADD COMMENT
0
Entering edit mode

You need to add any details for help. No DEGs can be biological or technical (confounders/noisy data, wrong input data type) etc.

ADD REPLY
0
Entering edit mode
swbarnes2 ★ 1.4k
@swbarnes2-14086
Last seen 4 hours ago
San Diego

Pure count data is an acceptable input for DESeq. Follow the vignette to generate a PCA plot. If you see no separation between your sample types in the first PCAs, it might be that your condition really doesn't do anything big enough to see in your experiment.

ADD COMMENT

Login before adding your answer.

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