R/plot_ts.R

In Enpeterson/outputsim: Output Simulation Results

plot_ts <- function(df, Nyears) {

senspars <- paste("sens.ct[1,", seq(1:Nyears), "]", sep="")
include <- which(df$parameter %in% c(senspars) & df$sim ==1)
dfsens <- df[include,]
year <- rep(NA, length(dfsens$parameter))
for(i in 1:length(dfsens$parameter)){
  year[i] <- as.numeric(gsub("]", "", strsplit(dfsens$parameter[i], ",")[[1]][2]))
}
dfsens$year <- year
dfsens <- dfsens[order(dfsens$year, dfsens$runsetting),]


g <-ggplot( dfsens, aes(x = year, y = truth, color="black")) +
  geom_line()+ geom_point() +  facet_grid(runsetting~.) +
  labs(x = "Year", y = "Sensitivity") +
  geom_point(aes(x = year, y= median, col="darkblue"))+
geom_errorbar(aes(ymin=lower , ymax=upper), colour="darkblue")
g1 <- g + theme(legend.position="none")






specpars <- paste("spec.ct[1,", seq(1:Nyears), "]", sep="")
include <- which(df$parameter %in% c(specpars) & df$sim ==1)
dfspec <- df[include,]
year <- rep(NA, length(dfspec$parameter))
for(i in 1:length(dfspec$parameter)){
  year[i] <- as.numeric(gsub("]", "", strsplit(dfspec$parameter[i], ",")[[1]][2]))
}
dfspec$year <- year
dfspec <- dfspec[order(dfspec$year, dfspec$runsetting),]


g2 <-ggplot( dfspec, aes(x = year, y = truth, color="black")) +
  geom_line()+ geom_point() +  facet_grid(runsetting~.) +
  labs(x = "Year", y = "specitivity") +
  geom_point(aes(x = year, y= median, col="darkblue"))+
  geom_errorbar(aes(ymin=lower , ymax=upper), colour="darkblue")
g3 <- g2 + theme(legend.position="none")



out <-multiplot(g1, g3, cols=1)
print(out)

}