How to filter genes inside a normalized count data file based on p-adjust value?
0
0
Entering edit mode
16 months ago
Amr ▴ 160

How to filter genes inside a normalized count data file based on p-adjust value?

I want to add a filter to only take the genes under specific p-adjust value?

Here are the normalization steps:

dds <- DESeqDataSetFromMatrix(countData = countdata,colData = coldata,design= ~ condition)
dds <- DESeq(dds)
count_data_normalized <- counts(dds, normalized=TRUE)
normalization p-value DESeq2 • 878 views
ADD COMMENT
0
Entering edit mode

Based on which p-values ? I can't see any differential expression analysis here

ADD REPLY
0
Entering edit mode

Do I have to extract the results table from results function: results(dds), then I make filtration based on p-adjust value? But If I do that, then how can I extract again the genes under specific p-value?

ADD REPLY
0
Entering edit mode

If you want to select genes from a differential analysis results, yes you should use results() but pay attention to your study design to extract the contrast you need. Then you can use dplyr and the filter function to retain the genes which meet your specific criteria

ADD REPLY
0
Entering edit mode

So, use results command to extract the p-adjust, then I need to select the filtered genes to pick them only from the count data normalized file, right? If so then how? by left_join? Thanks a lot @Basti

ADD REPLY
0
Entering edit mode

When you have the filtered genes, just select from counts the corresponding rownames

ADD REPLY
0
Entering edit mode

But the genes will gonna be several as the data already have 24,000 genes

ADD REPLY
0
Entering edit mode

Please clarify your problem, and explain the desired output

ADD REPLY

Login before adding your answer.

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