envelope {ssym}

Building of normal probability plots with simulated envelope of the deviance-type residuals.

Description

envelope is used to calculate and display normal probability plots with simulated envelope of the deviance-type residuals.

Usage

envelope(object, reps, conf, xlab.mu, ylab.mu, main.mu, xlab.phi, ylab.phi, main.phi)

Arguments

object

an object of the class ssym. This object returned from the call to ssym.l() or ssym.nl().

reps

a positive integer representing the number of iterations in which the simulated envelopes are based. Default is reps=25.

conf

value within the interval (0,1) that represents the confidence level of the simulated envelopes. Default is conf=0.95.

xlab.mu

character. An optional label for the x axis for the graph of the deviance-type residuals for the median submodel.

ylab.mu

character. An optional label for the y axis for the graph of the deviance-type residuals for the median submodel.

main.mu

character. An optional overall title for the plot for the graph of the deviance-type residuals for the median submodel.

xlab.phi

character. An optional label for the x axis for the graph of the deviance-type residuals for the skewness submodel.

ylab.phi

character. An optional label for the y axis for the graph of the deviance-type residuals for the skewness submodel.

main.phi

character. An optional overall title for the plot for the graph of the deviance-type residuals for the skewness submodel.

Details

Objects of the class ssym obtained from the application of ssym.l2() are not supported. The smoothing parameters are assumed to be known.

Author(s)

Luis Hernando Vanegas <[e-mail]> and Gilberto A. Paula

References

Atkinson, A. C. (1985) Plots, transformations and regression: an introduction to graphical methods of diagnostic regression analysis. Oxford Science Publications, Oxford.

Examples

###################################################################################
################# Blood flow Data - a log-power-exponential model #################
###################################################################################
#data("la", package="gamlss.nl")
#fit <- ssym.nl(log(PET60) ~ log(bflow) + log(1+b1*exp(-b2/bflow)) | bflow,
#       data=la, start=c(b1=-0.6,b2=98), family="Powerexp", xi=-0.45)
#summary(fit)
#
################## Simulated envelopes ##################
#envelope(fit,reps=50,conf=0.99)