forked from wchnicholas/site_B_landscape
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Plot_TrackPref.R
35 lines (33 loc) · 1.26 KB
/
Plot_TrackPref.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
#R code
library(ggplot2)
library(scales)
library(RColorBrewer)
library(readr)
library(tidyr)
library(reshape)
library(stringr)
library(dplyr)
library(ggrepel)
library(gridExtra)
require(cowplot)
strain_levels <- c('HK68','Bk79','Bei89','Mos99','Bris07','NDako16')
YearPrefTable <- read_tsv('result/Prefs_ByYear.tsv') %>%
filter(Background!='Bris07P194') %>%
mutate(Background=factor(Background,levels=strain_levels))
palette <- c(brewer.pal(9,"Set3"))
palette <- c(palette[3:12])
textsize <- 7
p <- ggplot(YearPrefTable, aes(x=Year,y=MeanPref,group=Background,color=Background)) +
geom_line(size=0.4) +
geom_point(size=0.4) +
geom_errorbar(aes(ymin=MeanPref-StdPref, ymax=MeanPref+StdPref), width=0) +
scale_color_manual(values=palette) +
theme_cowplot(12) +
theme(plot.title=element_text(size=textsize,face="bold"),
axis.title=element_text(size=textsize,face="bold"),
axis.text=element_text(size=textsize,face="bold"),
legend.title=element_blank(),
legend.text=element_text(size=textsize,face="bold"),
legend.position='right') +
ylab(bquote(bold(Mean~preference)))
ggsave('graph/PrefsByYear.png',p,height=1.8,width=4.5)