#' Plotting the difference between infected and and uninfected hibernating
#' bats.
#'
#' \code{dangerZone} plots the difference in hibernation length potentials
#' between infected and uninfected bats accross the environmental space.
#'
#' @return returns an environmental surface accross which the difference
#' between infected and uninfected hibernation times is plotted.
#'
#' @param mod.df results data from \code{\link{hibernationModel}}
#' @param title title for the figure
#' @param save.name a name and file path relative to the working directory
#' where the figure will be saved.
#' @param ... used to pass arguments to \code{\link[ggplot2]{ggsave}}
#'
#' @details Function will be updated in the future to include more plotting
#' parameter accessability, and potentially a raster method as well.
#' \url{https://www.youtube.com/watch?v=kyAn3fSs8_A}
#' @author L. Cain, K. Loggins
#' @family Plot Functions
#' @seealso \code{\link{survivalMultiplot}}; \code{\link{survivalPlotter}}; \code{\link{survivalHistogram}}
#' @export
dangerZone <-function(mod.df, title, save.name=NULL, ...){
mod.dif <- mod.df %>%
group_by_(~Ta, ~pct.rh) %>%
summarise_(max.null = ~max(time*surv.null),max.inf = ~max(time*surv.inf)) %>%
mutate_(diff = ~hour.to.month(max.inf - max.null))
dz <- ggplot(mod.dif, aes_(~Ta, ~pct.rh, z = ~diff)) +
scale_fill_gradientn("Difference\n(months)",
colors = c("#e66101", "#fdb863","#ffffff", "#b2abd2", "#5e3c99"), #purp low orange hi
limits = c(-12,0)
) +
geom_raster(aes_(fill = ~diff), interpolate = T) +
scale_x_continuous(expand = c(0,0))+
scale_y_continuous(expand = c(0,0))+
geom_contour(binwidth = 1,
colour = "grey15") +
ggtitle(title) +
xlab("Temperature (C)") +
ylab("Relative Humidity (%)")+
theme_minimal()+
theme(plot.title = element_text(size = 18, family="serif"),
axis.title = element_text(size = 16, family="serif"),
axis.text = element_text(size = 16, family="serif"),
aspect.ratio = 1,
legend.key.size = unit(42, "points"),
legend.title = element_text(size = 16, family="serif"),
legend.text = element_text(size = 16, family="serif"))
if(!is.null(save.name)){
ggsave( filename = save.name, ...)
return(invisible(NULL))
}
return(dz)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.