mqgamV: Fit multiple QGAM models and get a mgamViz object

In mfasiolo/mgcViz: Visualisations for Generalized Additive Models

Fit multiple QGAM models and get a mgamViz object

Description

These are wrapper that fits multple QGAM models using qgam::mqgam and converts it to a mgamViz object using the getViz function. It is essentially a shortcut.

Usage

mqgamV(form, data, qu, lsig = NULL, err = NULL, aQgam = list(), aViz = list())

Arguments

form, data, qu, lsig, err	same arguments as in qgam::mqgam.
aQgam	list of further arguments to be passed to qgam::mqgam.
aViz	list of arguments to be passed to getViz.

Value

An object of class "mgamViz" which can, for instance, be plotted using plot.mgamViz.

Examples

library(mgcViz)
set.seed(2) ## simulate some data...
dat <- gamSim(2,n=500,dist="normal",scale=0.25)$data

# Fit GAM and get gamViz object
b <- mqgamV(y~s(x) + s(z) + I(x*z), data = dat, qu = c(0.25, 0.5, 0.75),
            aQgam = list(argGam = list(select = TRUE)), aViz = list("nsim" = 0))

# This is equivalent to doing
# 1. Fit QGAM
# b <- mqgam(y~s(x) + s(z) + I(x*z), data=dat, 
#            qu = c(0.25, 0.5, 0.75), argGam = list(select = TRUE))
# 2. Convert to gamViz object
# b <- getViz(b, nsim = 0)

# Either way, we all effects by doing
print(plot(b, allTerms = TRUE), pages = 1)

mfasiolo/mgcViz documentation built on April 19, 2024, 8:16 a.m.