residuals.dlmFiltered {dlm}R Documentation

One-step forecast errors

Description

The function computes one-step forecast errors for a filtered dynamic linear model.

Usage

## S3 method for class 'dlmFiltered':
residuals(object, ..., type = c("standardized", "raw"), sd = TRUE)

Arguments

object an object of class "dlmFiltered", such as the output from dlmFilter
... unused additional arguments.
type should standardized or raw forecast errors be produced?
sd when sd = TRUE, standard deviations are returned as well.

Value

A vector or matrix (in the multivariate case) of one-step forecast errors, standardized if type = "standardized". Time series attributes of the original observation vector (matrix) are retained by the one-step forecast errors.
If sd = TRUE then the returned value is a list with the one-step forecast errors in component res and the corresponding standard deviations in component sd.

Note

The object argument must include a component y containing the data. This component will not be present if object was obtained by calling dlmFilter with simplify = TRUE.

Author(s)

Giovanni Petris GPetris@uark.edu

References

Harrison and West, Bayesian forecasting and dynamic models (2nd ed.), Springer (1997).

See Also

dlmFilter

Examples

## diagnostic plots 
nileMod <- dlmModPoly(1, dV = 15100, dW = 1468)
nileFilt <- dlmFilter(Nile, nileMod)
res <- residuals(nileFilt, sd=FALSE)
qqnorm(res)
tsdiag(nileFilt)

[Package dlm version 0.99-0 Index]