envelope: 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.

Author

Luis Hernando Vanegas <hvanegasp@gmail.com> and Gilberto A. Paula

Details

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

References

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

Examples

Run this code

###################################################################################
################# 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)

Run the code above in your browser using DataCamp Workspace