summary.gamlss
is the GAMLSS specific method for the generic function summary
which summarize
objects returned by modelling functions.## S3 method for class 'gamlss':
summary(object, type = c("vcov", "qr"),
robust=FALSE, save = FALSE,
hessian.fun = c("R", "PB"), ...)
vcov
uses the vcov()
method for gamlss to get the
variance-covariance matrix of the estimated beta coefficients, see details below.
The alternative qr
is the original method used in gamlss tooptimHess()
or a function based on Pinheiro and Bates nlme
package, "PB".type="vcov"
, the vcov()
method for gamlss is used to get
the variance covariance matrix (and consequently the standard errors)of the beta parameters.
The variance covariance matrix is calculated using the inverse of the numerical second derivatives
of the observed information matrix. This is a more reliable method since it take into the account the
inter-correlation between the all the parameters. The type="qr"
assumes that the parameters are fixed
at the estimated values. Note that both methods are not appropriate and should be used with caution if smoothing
terms are used in the fitting.Stasinopoulos D. M., Rigby R.A. and Akantziliotou C. (2003) Instructions on how to use the GAMLSS package in R.
Accompanying documentation in the current GAMLSS help files, (see also
Stasinopoulos D. M. Rigby R.A. (2007) Generalized additive models for location scale and shape (GAMLSS) in R.
Journal of Statistical Software, Vol. 23, Issue 7, Dec 2007,
gamlss
, deviance.gamlss
, fitted.gamlss
data(aids)
h<-gamlss(y~poly(x,3)+qrt, family=PO, data=aids) #
summary(h)
rm(h)
Run the code above in your browser using DataLab