How to create better separation in PCA plot
1
0
Entering edit mode
@rajeshparmar4-22077
Last seen 4.1 years ago
UCLA

Hello, I am using dds object for my RNA-seq data and I am creating PCA with my gene expression data, I am not getting clear visualization of my treatment or control. There is any other exploratory analysis or way to visualize clear separation between control or treated groups. I am using these commands in R studio:

dds <- DESeqDataSetFromMatrix(countData = new_count, colData = sample.data, design = ~ outcome_txt)
vsddds <- vst(dds, blind = FALSE)
plotPCA(vsddds, "outcome_txt")
PCA plot RNAseq • 881 views
ADD COMMENT
1
Entering edit mode
Kevin Blighe ★ 3.9k
@kevin
Last seen 1 day ago
Republic of Ireland

This is not quite a Bioconductor question and would be more suited for Biostars.

One does not require visual separation of groups via a PCA bi-plot in order for an effect / difference to exist between groups. Also, the plotPCA function will [by default] only plot PC1 versus PC2, whereas the differences between your groups may be reflected by PC6, PC8, or PC25.

It is important to first understand what is PCA, what it does to your data, and what it is measuring. See my post on Biostars in order to get a quick intro on this:

Also see our vignette for PCAtools, a Bioconductor package, which helps an individual to explore their data via PCA.

Kevin

ADD COMMENT
0
Entering edit mode

Thank you so much. I will look into it

ADD REPLY

Login before adding your answer.

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