DESeq2 setting significant p-value
1
0
Entering edit mode
Merlin ▴ 10
@merlin-15723
Last seen 4.4 years ago
Vancouver

Hello everyone, I have a question about the p-value that comes out from DESeq2.

Here is the command line that I use:

res <- results(dds, alpha = 0.05)
resultsNames(dds)
resLFC <- lfcShrink(dds, coef=2, res=res, type="apeglm")
resFile <-na.omit(resLFC)
resFile <-resFile[order(resFile$padj),]
write.csv(as.data.frame(resFile),file="~/RNAseq/salmon_folder/file.csv", quote=F)

The p-value and padj have lots of values greated then 0.05 and I’m not sure if I have to remove this values before plotting this data like when I do an heatmap for example.

Thank you

DESeq2 • 808 views
ADD COMMENT
2
Entering edit mode
@mikelove
Last seen 12 hours ago
United States

This is totally up to you, but in the workflow we simply plot the top genes by variance. Others will subset to the genes in the FDR set.

ADD COMMENT
0
Entering edit mode

after running DESeq2 does it sort the value on the basis of pvalue, or anything, or does it give the data back in the input format

I have passed an input in DESeq2 where my data has sorted gene names by A-Z, now in results it is not showing the gene names in result. now I want to know if it DESeq2 is sorting the result or not on the basis of pvalue or anything, or i just can use my previous expression file's list and then manually take out the DEGs.

ADD REPLY
0
Entering edit mode

It does not resort. Also the dataset and results will have the same row names from object creation throughout the analysis.

ADD REPLY
0
Entering edit mode

Sir but when I put in the expression data via the csv file, it was showing the gene names then but now when i set it in the form factor it is not showing the result.

ADD REPLY
0
Entering edit mode
> counts <- matrix(1:16, ncol=4, dimnames=list(letters[1:4],NULL))
> coldata <- data.frame(x=1:4)
> DESeqDataSetFromMatrix(counts, coldata, ~1)
class: DESeqDataSet
dim: 4 4
metadata(1): version
assays(1): counts
rownames(4): a b c d
rowData names(0):
colnames: NULL
colData names(1): x

Here the rownames are present in the dataset and will be propagated to results.

ADD REPLY

Login before adding your answer.

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