Plot a multi-panel plot of a parametric predictor across multiple QGAMs in a MQGAM

facet_parametric_plot creates a multi-panel plot of point plots with confidence interval ranges using ggplot2.

facet_parametric_plot(
  qgam,
  pred,
  cond = NULL,
  print.summary = FALSE,
  order = NULL,
  ncol = 1,
  xlab = "fit",
  ylab = NULL,
  scales = "free",
  size = 0.5,
  color = NULL,
  alpha = 1)

Arguments

qgam

An mqgam object created with qgam::mqgam.

pred

The predictor term to plot. Note: This is no longer identical to the pred argument specified for itsadug::plot_parametric.

cond

A named list of the values to use for the other predictor terms (not in view). Used for choosing between smooths that share the same view predictors.

print.summary

Logical: whether or not to print summary.

order

Specify the order with which the levels given in pred should be plotted.

ncol

The number of columns of the multi-panel plot.

xlab

The x-axis label.

ylab

The y-axis label.

scales

Should scales be free ("free", the default) or fixed ("fixed")?

size

Size argument for the ggplot object; specifies the size of points and lines.

color

Color argument for the ggplot object; specifies the color of points and lines.

alpha

Alpha argument for the ggplot object; specifies the transparency of points and lines.

Value

A ggplot object.

References

Fasiolo M., Goude Y., Nedellec R., & Wood S. N. (2017). Fast calibrated additive quantile regression. URL: https://arxiv.org/abs/1707.03307

van Rij J, Wieling M, Baayen R, & van Rijn H (2020). itsadug: Interpreting Time Series and Autocorrelated Data Using GAMMs. R package version 2.4.

Wickham, H. (2016). ggplot2: Elegant Graphics for Data Analysis. Springer-Verlag New York.

Author

D. Schmitz

Examples

# basic example
facet_parametric_plot(qgam = mtqgam_mqgam,
  pred = "factor_3")
#> i Plotting with default order of predictor levels.


# combining facet_parametric_plot with ggplot2
facet_parametric_plot(qgam = mtqgam_mqgam,
  pred = "factor_3") +
  theme_void() +
  labs(subtitle = "This is a subtitle")
#> i Plotting with default order of predictor levels.