gaulss: Gaussian location-scale model family

Description Usage Arguments Details Value References Examples

View source: R/gamlss.r

Description

The gaulss family implements Gaussian location scale additive models in which the mean and the logb of the standard deviation (see details) can depend on additive smooth predictors. Useable only with gam, the linear predictors are specified via a list of formulae.

Usage

1
gaulss(link=list("identity","logb"),b=0.01)

Arguments

link

two item list specifying the link for the mean and the standard deviation. See details.

b

The minumum standard deviation, for the "logb" link.

Details

Used with gam to fit Gaussian location - scale models. gam is called with a list containing 2 formulae, the first specifies the response on the left hand side and the structure of the linear predictor for the mean on the right hand side. The second is one sided, specifying the linear predictor for the standard deviation on the right hand side.

Link functions "identity", "inverse", "log" and "sqrt" are available for the mean. For the standard deviation only the "logb" link is implemented: eta = log(sigma-b) and sigma = b + exp(eta). This link is designed to avoid singularities in the likelihood caused by the standard deviation tending to zero. Note that internally the family is parameterized in terms of the tau=1/sigma, i.e. the standard deviation of the precision, so the link and inverse link are coded to reflect this, however the reltaionships between the linear predictor and the standard deviation are as given above.

The fitted values for this family will be a two column matrix. The first column is the mean, and the second column is the inverse of the standard deviation. Predictions using predict.gam will also produce 2 column matrices for type "link" and "response". The second column when type="response" is again on the reciprocal standard deviation scale (i.e. the square root precision scale). The second column when type="link" is log(sigma-b). Also plot.gam will plot smooths relating to sigma on the log(sigma-b) scale (so high values correspond to high standard deviation and low values to low standard deviation). Similarly the smoothing penalties are applied on the (log) standard deviation scale, not the log precision scale.

The null deviance reported for this family is the sum of squares of the difference between the response and the mean response divided by the standard deviation of the response according to the model. The deviance is the sum of squares of residuals divided by model standard deviations.

Value

An object inheriting from class general.family.

References

Wood, S.N., N. Pya and B. Saefken (2016), Smoothing parameter and model selection for general smooth models. Journal of the American Statistical Association 111, 1548-1575 doi: 10.1080/01621459.2016.1180986

Examples

1
2
3
4
5
library(mgcv);library(MASS)
b <- gam(list(accel~s(times,k=20,bs="ad"),~s(times)),
            data=mcycle,family=gaulss())
summary(b) 
plot(b,pages=1,scale=0)

Example output

Loading required package: nlme
This is mgcv 1.8-28. For overview type 'help("mgcv-package")'.

Family: gaulss 
Link function: identity logb 

Formula:
accel ~ s(times, k = 20, bs = "ad")
~s(times)

Parametric coefficients:
              Estimate Std. Error z value Pr(>|z|)    
(Intercept)   -25.3249     1.8530  -13.67   <2e-16 ***
(Intercept).1   2.6131     0.0632   41.35   <2e-16 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Approximate significance of smooth terms:
              edf Ref.df Chi.sq p-value    
s(times)   10.584 11.951  548.1  <2e-16 ***
s.1(times)  7.265  8.233  272.7  <2e-16 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Deviance explained = 98.7%
-REML = 580.04  Scale est. = 1         n = 133

mgcv documentation built on June 1, 2021, 5:11 p.m.