How to add dots to the PCA (in addition to labels)?

I made a PCA plot after RUVseq normalization with the command

plotPCA(threeset, col=colors[threex], cex=1.2)

The resulting plot is a bit unreadable. How can I add the dots (and remain labels), or at least do the plot in dots, and put labels into legend? It looks like this now: https://ibb.co/0nqfmc9 Thanks a lot in advance!

PCA

rna-seq edger

• 1,853 views • link

updated 6.2 years ago by Kevin Blighe

• written 6.2 years ago by tanya_fiskur •

1 answer

Hey Tanya,

Looks like you're using prcomp() to perform PCA?

You can re-use some of the code that I posted here: A: PCA plot from read count matrix from RNA-Seq

Or, there is also a dedicated package that I developed: PCAtools: everything Principal Component Analysis

Kevin

• 0 views • link

written 6.2 years ago by Kevin Blighe

More posts like this

Help with adding variables to PCA
written 4.0 years ago by Marina •

Hi there, I'm currently running a PCA on an RNAseq dataset. I'm rather new to both R and PCA but so far I have a …
Network plot using NetCoMi and igraph
written 4.9 years ago by Bioinfonext

Hi, I am trying to plot the Network plot as suggested here (https://github.com/stefpeschel/NetCoMi#single-association-network-on-genus-level) by using igraph and NetCoMi. But I am not getting the network …
WGCNA soft thresholding problem
written 5.5 years ago by lilingjoyo •

Hi there, I'm using WGCNA to construct co-expression network. My input is z-score normalized expression from RNAseq. Here is the distribution of normalized expression cross …
plot principal component analysis: how to improve the graphics
written 6.8 years ago by SeaStar •

I have created this PCA plot using BioGenerics with my data: boxplotPCA= plotPCA(table, labels =TRUE, isLog= FALSE, main= "PCA") obtaining this plot: ![enter image description …
Atypical Volcano plot of RNAseq data
written 6.9 years ago by saamar.rajput

I performed differential gene expression analysis using EgdeR on RNAseq data and using the DE i generated volcano plot. The volcano plot does not look …
DESeq2 unpaired sample RNAseq analysis
written 6.9 years ago by Adeler001 •

Hello I'm new to coding on Rstudio. I'm doing a RNA seq analysis to test for differential gene expression using deseq2, by using unpaired samples. …
Need help correcting overlappinge genes on valcano plot
written 6.9 years ago by Adeler001 •

I'm doing an RNA seq analysis and I'm trying to show my results using a volcano plot on R studio. I used the the follow …
WGCNA: Error in plot.window(...) : need finite 'ylim' values
written 7.2 years ago by ashish

I am trying to learn WGCNA through this [tutorial][1]. During the network construction step I am getting this error `Error in plot.window(...) : need finite …
DESeq2 Heatmap issue
written 8.9 years ago by dazhudou1122

Dear Community, I stole the codes from Stephen Turner to analyze RNAseq data using DESeq2 (out put of featureCount). Everything worked well, except that I …
SNPRelate: how to give specific color to a population in PCA plot
written 11.4 years ago by kumar.vinod81

I am using SNPRelate for PCA analysis. Its using default color for different populations but I want to color them according to me. Plotting commands …