fix_family_qf: Fix $qf() slot of a mgcv family

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

Fix $qf() slot of a mgcv family

Description

Fix $qf() slot of a mgcv family.

Usage

fix_family_qf(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_qf( b$family )
qfit <- sapply(c(0.1, 0.5, 0.9),
               function(.q)
                 b$family$qf(.q, b$fitted.values))

plot(mcycle, ylim = c(-150, 100))
for(ii in 1:ncol(qfit)){
  lines(mcycle$times, qfit[ , ii], col = 2)
}

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