--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/04_pca_plot.Rmd	Sun Dec 30 12:45:56 2018 -0500
@@ -0,0 +1,20 @@
+---
+title: 'Principle component analysis plot'
+output:
+    html_document:
+      highlight: pygments
+---
+
+```{r setup, include=FALSE, warning=FALSE, message=FALSE}
+knitr::opts_chunk$set(error = TRUE, echo = FALSE)
+```
+
+
+```{r}
+vsd <- vst(dds, blind=FALSE)
+p = plotPCA(vsd, intgroup=c(opt$X_D)) + 
+  scale_color_discrete(name = 'Group') +
+  theme_classic()
+ggplotly(p)
+```
+