Question: Is it required to specify Outliers function in DESeq2
0
gravatar for angajalaanusha
8 weeks ago by
angajalaanusha10 wrote:

Please see following code.I am wondering DESeq() automatically removes the outliers or should I include or specify the outliers function? When I did the analysis in around 100 samples the results looks fine. But in case of higher numbers of samples, the results looks different. Please advice how can I visualize the normalized data in box plots.

## QNBCvsTNBC ARP
dds <- DESeqDataSetFromMatrix(countData = cts,
                              colData = coldata,
                              design = ~ TNBCstatus)
dds <- DESeq(dds)
##results
colData(dds)[, "TNBCstatus"]
resultsNames(dds)
res <- results(dds, contrast=c("TNBCstatus", "QNBC", "TNBCARP"))
res
mcols(res, use.names = TRUE)
resOrdered <- res[order(res$pvalue),]
write.csv(as.data.frame(resOrdered), 
          file="QNBC_vs_TNBC.csv")
resSig <- subset(resOrdered, padj < 0.05)
resSig
write.csv(resSig,"QNBC_VS_TNBC_significant")
plotMA(res, ylim=c(-2,2))
vsd <- varianceStabilizingTransformation(dds,blind = TRUE)
n <- assay(vsd)
write.csv(n,"normalized data_log2values_QNBC_VS_TNBCdds")
plotMA(resSig, ylim=c(-2,2))
deseq2 • 66 views
ADD COMMENTlink modified 8 weeks ago by Michael Love23k • written 8 weeks ago by angajalaanusha10
Answer: Is it required to specify Outliers function in DESeq2
0
gravatar for Michael Love
8 weeks ago by
Michael Love23k
United States
Michael Love23k wrote:

Please take a look at the DESeq2 paper and vignette, which both have sections about how outliers are treated.

ADD COMMENTlink written 8 weeks ago by Michael Love23k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 290 users visited in the last hour