plot.Mort1Dsmooth: Plot of the Mort1Dsmooth outcomes

In MortalitySmooth: Smoothing and Forecasting Poisson Counts with P-Splines

Description

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

Usage

## 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. Not in used.

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 are given and cannot be changed.

Value

None. Function produces a plot

Author(s)

Carlo G Camarda

See Also

Mort1Dsmooth for computing Mort1Dsmooth.object.

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)
plot(fit, "deaths")
par(mfrow=c(1,1))

Example output

Loading required package: svcm
Loading required package: Matrix
Loading required package: splines
Loading required package: lattice

MortalitySmooth documentation built on May 2, 2019, 6:07 a.m.