R/Figure2.R

In mikemahoney218/BeaverForaging:

#'
#'
#' @import dplyr betareg ggplot2 cowplot
#'
Figure2 <- function(){
  CCBetaRegression()
  BetaModelPred <- tibble(PlotCode = rep(seq(0, 64, 0.1), 2),
                          LakeRiver = c(rep("Lake", 641),
                                        rep("River", 641)),
                          LogCode = log(PlotCode))
  BetaModelPred$Predict <- predict(LightBetaReg, BetaModelPred)

  Fig2 <<- ggplot(data = BetaModelPred, aes(PlotCode, linetype = LakeRiver, shape = LakeRiver)) +
    geom_line(data = BetaModelPred, aes(y = Predict * 100), color = "black", size = 0.8) +
    geom_jitter(data = LightData, aes(y = CanopyProp * 100), alpha = 1,
                size = 2,
                height = 0,
                width = 2) +
    background_grid() +
    scale_shape_manual(name = "Site Type",
                       labels = c("Lake",
                                  "Stream"),
                       values = c(2, 1)) +
    scale_linetype(name = "Site Type",
                   labels = c("Lake",
                              "Stream")) +
    theme(text = element_text(size = 10),
          axis.text = element_text(size = 10),
          aspect.ratio = (9/16)) +
    guides(shape = guide_legend(override.aes = list(alpha = 1))) +
    scale_y_continuous(name = "% Canopy cover") +
    scale_x_continuous(name = "Distance from impoundment (m)",
                       breaks = c(0, 4, 14, 24, 34, 44, 54, 64)) +
    theme_cowplot()

  return(Fig2)

}