plot.resStMoMo {StMoMo} | R Documentation |
Plot the residuals of a Stochastic Mortality Model
Description
Plots the deviance residuals of a Stochastic Mortality Model which are
of class "resStMoMo"
. Three types of plots
are available: scatter plot of residuals by age, period and cohort,
colour map (heatmap) of the residuals, and a black and white signplot
of the residuals.
Usage
## S3 method for class 'resStMoMo'
plot(x, type = c("scatter", "colourmap",
"signplot"),
reslim = NULL, plotAge = TRUE,
plotYear = TRUE, plotCohort = TRUE,
pch = 20, col = NULL, ...)
Arguments
x |
an object of class |
type |
the type of the plot. The alternatives are
|
reslim |
optional numeric vector of length 2, giving the range of the residuals. |
plotAge |
logical value indicating if the age scatter plot should be
produced. This is only used when |
plotYear |
logical value indicating if the calendar year scatter plot
should be produced. This is only used when |
plotCohort |
logical value indicating if the cohort scatter plot
should be produced. This is only used when |
pch |
optional symbol to use for the points in a scatterplot.
This is only used when |
col |
optional colours to use in plotting. If
|
... |
other plotting parameters to be passed to the plotting functions. This can be used to control the appearance of the plots. |
Details
When type = "scatter"
scatter plots of the residuals against age,
calendar year and cohort (year of birth) are produced.
When type = "colourmap"
a two dimensional colour map of the
residuals is plotted. This is produced using function
image.plot
. See image.plot
for further parameters that can be passed to this type of plots.
When type = "signplot"
a two dimensional black and white map of the
residuals is plotted with dark grey representing negative residuals and
light grey representing positive residuals. This is produced using
function image.default
.
@seealso residuals.fitStMoMo
Examples
CBDfit <- fit(cbd(), data = central2initial(EWMaleData), ages.fit = 55:89)
CBDres <- residuals(CBDfit)
plot(CBDres)
plot(CBDres, type = "signplot")
plot(CBDres, type = "colourmap")