plot.Mort1Dsmooth {MortalitySmooth}R Documentation

Plot of the Mort1Dsmooth outcomes

Description

It produces a plot where the x-axis is the predictor of the given Mort1Dsmooth object and the y-axis is either the actual and fitted log-rates, or the actual and fitted counts.

Usage

## S3 method for class 'Mort1Dsmooth'
## S3 method for class 'Mort1Dsmooth':
plot(x, type = c("logrates", "deaths"), ...)

Arguments

x An object of class "Mort1Dsmooth", usually, a result of a call to Mort1Dsmooth.
type The type of plot which should be returned. The alternatives are: logrates (default) and death.
... Other plotting parameters passed to plot.

Details

The function provides a simple tool for checking the outcomes of a Mort1Dsmooth object. In a Poisson setting, one would commonly look at rates in log-scale, but actual and fitted counts can be plotted too.

Plotting parameters can be change to enhance the output. Only exceptions are the x values, y values, and the ylim of the main plot, and the arguments for the fitted lines.

Value

None. Function produces a plot

Author(s)

Carlo G Camarda

See Also

Mort1Dsmooth for computing Mort1Dsmooth.object, plot, par.

Examples

# selected data
years <- 1930:2006
death <- selectHMDdata("Denmark", "Deaths", "Females",
                       ages = 60, years = years)
exposure <- selectHMDdata("Denmark", "Exposures", "Females",
                          ages = 60, years = years)
# fit
fit <- Mort1Dsmooth(x=years, y=death, offset=log(exposure),
                    method=3, lambda=100)
# plotting actual and fitted data
par(mfrow=c(1,2))
plot(fit, "deaths")
plot(fit, "logrates")
par(mfrow=c(1,1))

[Package MortalitySmooth version 1.0 Index]