inst/doc/vetresearchlmm-plotting.R

In VetResearchLMM: Linear Mixed Models - An Introduction with Applications in Veterinary Research

## -----------------------------------------------------------------------------
#| include: false
library(methods)
if (requireNamespace("lme4", quietly = TRUE)) library(lme4)
if (requireNamespace("lmerTest", quietly = TRUE)) library(lmerTest)
load(file.path("..", "data", "ex124.RData"))
load(file.path("..", "data", "ex125.RData"))
load(file.path("..", "data", "ex33.RData"))
load(file.path("..", "data", "ex31.RData"))


## -----------------------------------------------------------------------------
if (requireNamespace("collapse", quietly = TRUE) && requireNamespace("ggplot2", quietly = TRUE)) {
  ex124_plot <-
    ex124 |>
    collapse::fmutate(herd = factor(herd))

  ggplot2::ggplot(
    ex124_plot,
    ggplot2::aes(x = herd, y = PCVdif, shape = dose, color = drug)
  ) +
    ggplot2::geom_point(size = 2.4, position = ggplot2::position_jitter(width = 0.08)) +
    ggplot2::facet_wrap(~ drug, nrow = 1) +
    ggplot2::labs(
      x = "Herd",
      y = "Change in packed cell volume",
      shape = "Dose",
      color = "Drug",
      title = "Packed cell volume response by herd, drug, and dose"
    ) +
    ggplot2::theme_minimal()
}


## -----------------------------------------------------------------------------
if (requireNamespace("ggplot2", quietly = TRUE)) {
  ggplot2::ggplot(
    ex125,
    ggplot2::aes(x = dose, y = Pcv, group = Drug, color = Drug)
  ) +
    ggplot2::stat_summary(fun = mean, geom = "line") +
    ggplot2::stat_summary(fun = mean, geom = "point", size = 2.6) +
    ggplot2::facet_wrap(~ Region) +
    ggplot2::labs(
      x = "Dose",
      y = "Mean packed cell volume",
      title = "Mean response by region, drug, and dose"
    ) +
    ggplot2::theme_minimal()
}


## -----------------------------------------------------------------------------
if (requireNamespace("ggplot2", quietly = TRUE)) {
  ggplot2::ggplot(
    ex33,
    ggplot2::aes(x = time, y = PCV, group = animal_id, color = breed)
  ) +
    ggplot2::geom_line(linewidth = 0.5) +
    ggplot2::geom_point(size = 1.4) +
    ggplot2::labs(
      x = "Time",
      y = "Packed cell volume",
      color = "Breed",
      title = "Longitudinal packed cell volume by animal"
    ) +
    ggplot2::theme_minimal()
}


## -----------------------------------------------------------------------------
if (requireNamespace("collapse", quietly = TRUE) && requireNamespace("ggplot2", quietly = TRUE)) {
  ex31_plot <-
    ex31 |>
    collapse::fmutate(
      treatment = ifelse(drug == "BERENIL", paste(drug, dose), as.character(drug))
    )

  ggplot2::ggplot(
    ex31_plot,
    ggplot2::aes(x = PCV1, y = PCV2, color = treatment)
  ) +
    ggplot2::geom_point(size = 2) +
    ggplot2::geom_smooth(method = "lm", se = FALSE, linewidth = 0.8) +
    ggplot2::labs(
      x = "PCV at treatment",
      y = "PCV one month later",
      color = "Treatment",
      title = "Relationship between baseline and follow-up PCV"
    ) +
    ggplot2::theme_minimal()
}