Minor changes to figures

03_Figures_Simulations.R

@@ -98,7 +98,7 @@ ggplot(data=outp.summary %>%
   facet_grid(corr~inter) + 
   theme_bw() + 
   scale_x_continuous(name="Effect of factor e1 -> trait t1") + 
-  scale_y_continuous(name="Prop. of significant tests") + 
+  scale_y_continuous(name="Prop. of significant r(XY)") + 
   theme(panel.grid = element_blank())
 
 ggsave(filename="_pics/R1/Fig2_CorrInte_02March.png", width=6, height=5, device="png", dpi = 300, last_plot())
@@ -170,7 +170,7 @@ fig3 <- ggplot(data=outp.summary2) +
   geom_line(aes(x=main, y=power, group=trait, col=trait.col, lty=linetype)) + 
   guides(linetype = FALSE) + 
   scale_x_continuous(name="Effect of factor e1 -> trait t1", n.breaks = 4) + 
-  scale_y_continuous(name="Prop. of significant tests") + 
+  scale_y_continuous(name="Prop. of significant r(XY)") + 
   facet_grid(sel.ntraits.lab~ntraits.lab) +
   scale_color_identity(guide = "legend", 
                        labels= hugepalette$trait) +
@@ -240,7 +240,7 @@ ggplot(data=outp.summary %>%
   facet_grid(corr~inter) + 
   theme_bw() + 
   scale_x_continuous(name="Effect of factor e1 -> trait t1") + 
-  scale_y_continuous(name="Prop. of significant tests") + 
+  scale_y_continuous(name="Prop. of significant r(XY)") + 
   theme(panel.grid = element_blank())
 
 ggsave(filename="_pics/R1/FigS4_Extra_CorrInte_08Jul20.png", width=6, height=5, device="png", dpi = 300, last_plot())
@@ -285,7 +285,7 @@ ggplot(data=outp.summary %>%
   facet_grid(corr~inter) + 
   theme_bw() + 
   scale_x_continuous(name="Effect of factor e1 -> trait t1") + 
-  scale_y_continuous(name="Prop. of significant tests") + 
+  scale_y_continuous(name="Prop. of significant r(XW)") + 
   theme(panel.grid = element_blank())
 
 ggsave(filename="_pics/R1/FigSXXX_CorrInte_30November_XW.png", width=6, height=5, device="png", dpi = 300, last_plot())
@@ -332,7 +332,50 @@ ggplot(data=outp.summary %>%
   scale_y_continuous(name="Prop. of significant tests") + 
   theme(panel.grid = element_blank())
 
-ggsave(filename="_pics/R1/FigSXXX_CorrInte_30November_SbM.png", width=6, height=5, device="png", dpi = 300, last_plot())
+ggsave(filename="_pics/R1/FigSXXY_CorrInte_30November_SbM.png", width=6, height=5, device="png", dpi = 300, last_plot())
+
+
+#### FIGURE SXXO - XE ####
+mypath <- "_data/Experiment_02Mar2020_FactorInteraction&TraitCorr"
+myfiles <- list.files(path=mypath, pattern = "Summary.txt", recursive = T)
+outp.summary <- FormatData(myfiles) 
+
+mypalette <- palette(c("#e41a1c",  #1 - red)
+                       "#ff7f00", #12 - orange
+                       "#984ea3", #13 - violet
+                       "##ffed6f", #2 - yellow
+                       "#4daf4a", #23 - green
+                       "#377eb8")) #3 - blue
+
+ggplot(data=outp.summary %>% 
+         mutate(main=main/100) %>% 
+         mutate(corr=factor(corr/10, levels=c(0, 0.4, 0.8), labels=paste0("Correlation = ", c(0, 0.4, 0.8)))) %>% 
+         mutate(inter=factor(inter/10, levels=c(0, 0.3, 0.5), labels=paste0("Interaction = ", c(0, 0.3, 0.5)))) %>% 
+         ungroup() %>% 
+         dplyr::filter(stat.type=="XE") %>% 
+         dplyr::filter(envir==1) %>% 
+         dplyr::filter(trait %in% c("1", "2", "1 2", "3", "1 3", "2 3")) %>% 
+         mutate(trait=factor(trait, levels=c("1", "2", "3", "1 2", "1 3", "2 3"), 
+                             labels=c("t1", "t2", "tn", "t1 t2", "t1 tn", "t2 tn"))) %>% 
+         mutate(linetype=as.factor(ifelse(grepl(pattern="tn", x = trait), 2, 1)))) + 
+  geom_line(aes(x=main, y=power, group=trait, col=trait, lty=linetype)) + 
+  guides(linetype = FALSE) + 
+  scale_color_manual("Trait\ncomb.",
+                     values=c("#e41a1c",  #1 - red)
+                              "#e6ab02", #2 - yellow
+                              "#377eb8", #3 - blue
+                              "#d95f02", #12 - orange
+                              "#984ea3", #13 - violet
+                              "#4daf4a" #23 - green
+                     ) 
+  ) +
+  facet_grid(corr~inter) + 
+  theme_bw() + 
+  scale_x_continuous(name="Effect of factor e1 -> trait t1") + 
+  scale_y_continuous(name="Prop. of significant r(XE)") + 
+  theme(panel.grid = element_blank())
+
+ggsave(filename="_pics/R1/FigSXXO_CorrInte_02March_XE.png", width=6, height=5, device="png", dpi = 300, last_plot())