Adding everything to make structural plots.

wilkelab · Mar 6, 2015 · dbc0fc6 · dbc0fc6
1 parent c197da2
commit dbc0fc6
Show file tree

Hide file tree

Showing 12 changed files with 540 additions and 0 deletions.
diff --git a/HA/dNdS/distances/correlations_violin_HA.idraw b/HA/dNdS/distances/correlations_violin_HA.idraw
diff --git a/HA/dNdS/distances/correlations_violin_HA.pdf b/HA/dNdS/distances/correlations_violin_HA.pdf
diff --git a/HA/dNdS/distances/correlations_violin_structures_HA.pdf b/HA/dNdS/distances/correlations_violin_structures_HA.pdf
diff --git a/HA/dNdS/distances/scale.idraw b/HA/dNdS/distances/scale.idraw
diff --git a/HA/dNdS/distances/scale.png b/HA/dNdS/distances/scale.png
diff --git a/NA/dNdS/distances/correlations/.Rhistory b/NA/dNdS/distances/correlations/.Rhistory
diff --git a/NA/dNdS/distances/correlations_violin_NA.idraw b/NA/dNdS/distances/correlations_violin_NA.idraw
diff --git a/NA/dNdS/distances/correlations_violin_NA.pdf b/NA/dNdS/distances/correlations_violin_NA.pdf
diff --git a/NA/dNdS/distances/correlations_violin_structures_NA.pdf b/NA/dNdS/distances/correlations_violin_structures_NA.pdf
diff --git a/NA/dNdS/distances/scale.idraw b/NA/dNdS/distances/scale.idraw
diff --git a/NA/dNdS/distances/scale.png b/NA/dNdS/distances/scale.png
diff --git a/cor_violin.R b/cor_violin.R
@@ -0,0 +1,28 @@
+rm(list = ls())
+
+library(ggplot2)
+library(grid)
+library(cowplot)
+
+make.plot <- function(protein) {
+  correlations <- data.frame(ID=numeric(), correlation=numeric())
+
+  count <- 0
+  for(file in list.files()) {
+    count <- count + 1
+    tmp.cor <- read.table(file, head=F)$V1
+    correlations <- rbind(correlations, data.frame(ID=rep(count, length(tmp.cor)), correlation=tmp.cor))
+  }
+
+  p <- ggplot(correlations, aes(factor(ID), correlation)) + geom_violin(fill = 'gray') +
+    xlab('Time (Months)') +
+    ylab('Distribution of Correlations')
+
+  ggsave(filename=paste('../correlations_violin_', protein, '.pdf', sep=''), height=5, width=15, useDingbats=F)
+}
+
+setwd('~/Google Drive/Data/influenza_pH1N1_timecourse/HA/dNdS/distances/correlations/')
+make.plot('HA')
+
+setwd('~/Google Drive/Data/influenza_pH1N1_timecourse/NA/dNdS/distances/correlations/')
+make.plot('NA')