refit: Refit a GAMLSS model

Description Usage Arguments Details Value Note Author(s) References See Also Examples

Description

This function refits a GAMLSS model. It is useful when the algorithm has not converged after 20 outer iteration (the default value)

Usage

1
refit(object, ...)

Arguments

object

a GAMLSS fitted model which has not converged

...

for extra arguments

Details

This function is useful when the iterations have reach the maximum value set by the code(n.cyc) of the gamlss.control function and the model has not converged yet

Value

Returns a GAMLSS fitted model

Note

The function update does a very similar job

Author(s)

Mikis Stasinopoulos d.stasinopoulos@londonmet.ac.uk, Bob Rigby

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. Z., and De Bastiani, F. (2019) Distributions for modeling location, scale, and shape: Using GAMLSS in R, Chapman and Hall/CRC. An older version can be found in https://www.gamlss.com/.

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, https://www.jstatsoft.org/v23/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/).

See Also

gamlss, update.gamlss

Examples

1
2
3
4
data(aids)
h<-gamlss(y~poly(x,3)+qrt, family=PO, data=aids) # 
refit(h)
rm(h)

Example output

Loading required package: splines
Loading required package: gamlss.data
Loading required package: gamlss.dist
Loading required package: MASS
Loading required package: nlme
Loading required package: parallel
 **********   GAMLSS Version 5.1-2  ********** 
For more on GAMLSS look at http://www.gamlss.org/
Type gamlssNews() to see new features/changes/bug fixes.

GAMLSS-RS iteration 1: Global Deviance = 416.8014 
GAMLSS-RS iteration 2: Global Deviance = 416.8014 
GAMLSS-RS iteration 3: Global Deviance = 416.8014 

Family:  c("PO", "Poisson") 
Fitting method: RS() 

Call:  gamlss(formula = y ~ poly(x, 3) + qrt, family = PO,  
    data = aids, start.from = h, iter = 2, n.cyc = 40) 

Mu Coefficients:
(Intercept)  poly(x, 3)1  poly(x, 3)2  poly(x, 3)3         qrt2         qrt3  
   4.814473     8.256092    -3.353310     0.938816    -0.156381     0.009986  
       qrt4  
  -0.115705  

 Degrees of Freedom for the fit: 7 Residual Deg. of Freedom   38 
Global Deviance:     416.801 
            AIC:     430.801 
            SBC:     443.448 

gamlss documentation built on March 31, 2021, 5:10 p.m.