Update method for cumulative link models

Update method for cumulative link models fitted with clm2. This makes it possible to use e.g. update(obj, location = ~ . - var1, scale = ~ . + var2)

Usage

# S3 method for class 'clm2'
update(object, formula., location, scale, nominal,...,
        evaluate = TRUE)
# S3 method for class 'clmm2'
update(object, formula., location, scale, nominal,...,
        evaluate = TRUE)

Arguments

object

a clm2 object.

formula.

not used—unfortunately this argument is part of the default method.

location

an optional new formula for the location; see update.formula for details.

scale

an optional new formula for the scale; see update.formula for details.

nominal

an optional new formula for nominal effects; see update.formula for details.

...

additional arguments to the call, or arguments with changed values.

evaluate

if true evaluate the new call else return the call.

Value

If evaluate = TRUE the fitted object is returned, otherwise the updated call.

Author

Rune Haubo B Christensen

Examples

options(contrasts = c("contr.treatment", "contr.poly"))

m1 <-  clm2(SURENESS ~ PROD, scale = ~PROD, data = soup,
           link = "logistic")

m2 <- update(m1, link = "probit")
m3 <- update(m1, link = "cloglog")
m4 <- update(m1, link = "loglog")
anova(m1, update(m1, scale = ~.-PROD))
#> Likelihood ratio tests of cumulative link models
#> 
#> Response: SURENESS
#>            Model Resid. df -2logLik   Test    Df LR stat.  Pr(Chi)
#> 1    PROD | 1 |       1841 5380.664                               
#> 2 PROD | PROD |       1840 5375.489 1 vs 2     1 5.174937 0.022915
mT1 <- update(m1, threshold = "symmetric")