Generic checking function

Generic function for checking R objects which produces, for instance, convergence tests or diagnostic plots. For qgam objects check.qgam() will be used.

check(obj, ...)

Arguments

obj	the object to be checked.
...	extra arguments, mainly used by graphic functions.

Value

Reports the results of convergence tests and/or produces diagnostic plots.

Examples

#######
# Using check.qgam
#######
library(qgam)
set.seed(0)
dat <- gamSim(1, n=200)
#> Gu & Wahba 4 term additive model
b<-qgam(y~s(x0)+s(x1)+s(x2)+s(x3), data=dat, qu = 0.5)
#> Estimating learning rate. Each dot corresponds to a loss evaluation. 
#> qu = 0.5...............done 
plot(b, pages=1)
check(b, pch=19, cex=.3)
#> Theor. proportion of neg. resid.: 0.5   Actual proportion: 0.52
#> Integrated absolute bias |F(mu) - F(mu0)| = 0.04500293
#> Method: REML   Optimizer: outer newton
#> full convergence after 6 iterations.
#> Gradient range [-0.0001503744,2.710312e-05]
#> (score 454.5555 & scale 1).
#> Hessian positive definite, eigenvalue range [0.02110521,2.627179].
#> Model rank =  37 / 37 
#> 
#> Basis dimension (k) check: if edf is close too k' (maximum possible edf) 
#> it might be worth increasing k. 
#> 
#>       k'  edf
#> s(x0)  9 2.75
#> s(x1)  9 2.58
#> s(x2)  9 7.44
#> s(x3)  9 1.21