fix_family_cdf: Fix $cdf() slot of a mgcv family

In dill/mgcvUtils: Some extra utilities for generalized additive modelling

Fix $cdf() slot of a mgcv family

Description

Fix $cdf() slot of a mgcv family.

Usage

fix_family_cdf(fam)

Arguments

fam	an mgcv family

Examples

library(mgcvUtils); library(MASS)
b <- gam(list(accel~s(times,k=20,bs="ad"),~s(times)),
         data=mcycle,family=gaulss())

b$family <- fix_family_cdf( b$family )

# Quantile residuals should look uniformly distributed
qres <- b$family$cdf(mcycle$accel, b$fitted.values,
                     wt = rep(1, nrow(mcycle)), scale = 1)
hist(qres)
# Not too bad

dill/mgcvUtils documentation built on Oct. 14, 2024, 7:30 p.m.