gamlss.add (version 5.1-13)

centilesTwo: Centiles contour plots in GAMLSS

Description

This function centilesTwo() plots two dimensional centiles contour plots for GAMLSS models.

Usage

centilesTwo(object, grid.x1, grid.x2, x1.name, x2.name, 
       cent = 0.05, dist = 0.01, points = TRUE, 
       other = list(), point.col = 1, point.pch = ".", 
       image = FALSE, image.col = heat.colors(12), ...)

Value

Produce a contour plot.

Arguments

object

an gamlss object

grid.x1

grid values for x-variable one

grid.x2

grid values for x-variable two

x1.name

the name of x-variable on

x2.name

the name of x-variable two

cent

the required centiles

dist

the distance

points

whether to plot the data points

other

a list having other explanatory variables at fixed values

point.col

the colour of the data points

point.pch

the type of the data point

image

whether to plot using the image9 function

image.col

the colour scheme

...

for extra arguments for the contour() function

Author

Mikis Stasinopoulos , Bob Rigby, Fernanda De Bastiani

Details

The function uses the function exclude.too.far() of the package mgcv.

References

Rigby, R. A. and Stasinopoulos D. M. (2005). Generalized additive models for location, scale and shape,(with discussion), Appl. Statist., 54, part 3, pp 507-554.

Rigby R.A., Stasinopoulos D. M., Heller G., and De Bastiani F., (2019) Distributions for Modeling Location, Scale and Shape: Using GAMLSS in R, Chapman and Hall/CRC.

Stasinopoulos D. M. Rigby R.A. (2007) Generalized additive models for location scale and shape (GAMLSS) in R. Journal of Statistical Software, 23(7), 1--46, tools:::Rd_expr_doi("10.18637/jss.v023.i07")

Stasinopoulos D. M., Rigby R.A., Heller G., Voudouris V., and De Bastiani F., (2017) Flexible Regression and Smoothing: Using GAMLSS in R, Chapman and Hall/CRC.

(see also https://www.gamlss.com/).

Wood S.N. (2006) Generalized Additive Models: An Introduction with R. Chapman and Hall/CRC Press.

See Also

centiles

Examples

if (FALSE) { 
data(plasma)
m1 <- gamlss(betadiet ~ ga(~te(age, fiber)), sigma.formula = ~1,  
    nu.formula = ~ga(~te(age, fiber)), tau.formula = ~1,  
    family = BCTo, data = plasma)
centilesTwo(m1, 18:90, seq(2.5,38, 0.5), age, fiber, cent=0.05, dist=.1,
            xlab="age", ylab='fiber') 
centilesTwo(m1, 18:90, seq(2.5,38, 0.5), age, fiber, cent=0.95, dist=.1)
}