compareML: Function for comparing two GAMM models.

View source: R/test.R

compareMLR Documentation

Function for comparing two GAMM models.


Function for comparing two GAMM models.


  signif.stars = TRUE, = FALSE,
  print.output = TRUE



First model.


Second model.


Logical (default = TRUE). Whether or not to display stars indicating the level of significance on 95% confidence level.

Logical (default = FALSE). Whether or not to present a suggestion on how one could report the information. If print.output is set to FALSE, will set to FALSE too. Please inspect yourself whether the label between square bracket fits your own standards. Note: the X2 should be replaced by a proper Chi-Square symbol χ^2.


Logical: whether or not to print the output. By default set to true, even if the the messages are not allowed by a global package option using the function infoMessages.


As an Chi-Square test is performed on two times the difference in minimized smoothing parameter selection score (GCV, fREML, REML, ML), and the difference in degrees of freedom specified in the model. The degrees of freedom of the model terms are the sum of 1) the number of estimated smoothing parameters for the model, 2) number of parametric (non-smooth) model terms including the intercept, and 3) the sum of the penalty null space dimensions of each smooth object.

This method is preferred over other functions such as AIC for models that include an AR1 model or random effects (especially nonlinear random smooths using bs='fs'). CompareML also reports the AIC difference, but that value should be treated with care.

Note that the Chi-Square test will result in a very low p-value when the difference in degrees of freedom approaches zero. Use common sense to determine if the difference between the two models is meaningful. A warning is presented when the difference in score is smaller than 5.

The order of the two models is not important. Model comparison is only implemented for the methods GCV, fREML, REML, and ML.


Optionally returns the Chi-Square test table.


For suppressing the output and all warnings, set infoMessages to FALSE (infoMessages('off') ), set the argument print.output to FALSE, and use the function suppressWarnings to suppress warning messages.


Jacolien van Rij. With many thanks to Simon N. Wood for his feedback.

See Also

For models without AR1 model or random effects AIC can be used.

Other Testing for significance: plot_diff2(), plot_diff(), report_stats(), wald_gam()



## Not run: 
# some arbitrary models:
m1 <- bam(Y~Group + s(Time, by=Group), method='fREML', data=simdat)
m2 <- bam(Y~Group + s(Time), method='fREML', data=simdat)

compareML(m1, m2)

# exclude significance stars:
compareML(m1, m2, signif.stars=FALSE)

m3 <- bam(Y~Group + s(Time, by=Group, k=25), method='fREML', 
compareML(m1, m3)

# do not print output, but save table for later use:
cml <- compareML(m1, m2, print.output=FALSE)$table

# Use suppressWarnings to also suppress warnings:
suppressWarnings(cml <- compareML(m1, m2, print.output=FALSE)$table)

## End(Not run)

itsadug documentation built on June 17, 2022, 5:05 p.m.