Hi all,

I've watched this video and wants to visualize the PCA scree plot to check my PCA plot that was generated in DESeq2. 

Is there any way I can do it in DESeq2 or in other Bioconductor packages?

I'd recommend just using your own code to do the PCA.

You can copy out the first lines of the plotPCA function to get started, see ?plotPCA.

Then you'll just want to plot the amounts in 'percentVar'.

Hi Micheal,

Thank you so much!!!

I manage to generate the PCA scree plot by following your instruction.

Here is my script:

     library(DESeq2)
     library(matrixStats)
     
     #How to get PCA plot?
     ##how to obtain d.deseq was described in DESeq2 manual
     
     cds=estimateDispersions(dds)
     vsd=varianceStabilizingTransformation(cds)
     plotPCA(vsd,intgroup=c("condition","timepoints"))
     
     #How to get PCA scree plot?
     
     ## calculate the variance for each gene
     rv <- rowVars(assay(vsd))

     ## select the ntop genes by variance
     select <- order(rv, decreasing=TRUE)[seq_len(min(500, length(rv)))]
     
     ## perform a PCA on the data in assay(x) for the selected genes
     pca <- prcomp(t(assay(vsd)[select,]))
     
     ## the contribution to the total variance for each component
     percentVar <- pca$sdev^2 / sum( pca$sdev^2 )
     
     ##plot the "percentVar"
     scree_plot=data.frame(percentVar)
     scree_plot[,2]<- c(1:24)

     colnames(scree_plot)<-c("variance","component_number")
     ggplot(scree_plot, mapping=aes(x=component_number, y=variance))+geom_bar(stat="identity")