lca.dev.res: Miscellaneous utility functions for 'lca' and 'lca.rh' type...

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

Description

A simple utility function to replace the original residuals (e.g. logrates, rates, deaths) of a LC fit with deviance residuals without the need to re-estimate the regression parameters. We note that the estimation of the parameters can be particularly slow in the case of the APC model.

Usage

1
lca.dev.res(lca.obj, pop, clip = 0)

Arguments

lca.obj

an object of class lca

pop

matrix of population data corresponding to the fitted mortality rates

clip

number of years to clip from start and end of cohort years

Details

The Lee-Carter regression object contains the type of residuals specified in the original function call, which might need to be changed for further analysis, but without actually re-running the entire iterative estimation process.

Value

An identical regression object as lca.obj containing the corresponding deviance residuals

Author(s)

Zoltan Butt, Steven Haberman and Han Lin Shang

See Also

lca.rh

Examples

1
2
3
4
# original model object with 'logrates' residuals
mod6 <- lca.rh(dd.cmi.pens, mod="lc", error="gauss", max=110, interpolate=TRUE)
# adjusted model object with 'deviance' residuals: 
dev6 <- lca.dev.res(mod6, insp.dd(dd.cmi.pens, "pop"))

ilc documentation built on May 2, 2019, 5:07 a.m.