forecast: Run YourCast forecasts given selected weights.
In IQSS/AutoCast: Automated Bayesian Forecasting with YourCast

Description Usage Arguments Value Examples

forecast runs final YourCast forecasts given a "gridcast" object and a "weightcast" object as inputs.

1 2	forecast(obj, weights, point.estimate = median, lower.bound = min, upper.bound = max)

`obj`	Object of class `"gridcast"`.
`weights`	Object of class `"weightcast"`.
`point.estimate`	Function applied to calculate the point estimate of the final forecast. Default: `median`.
`lower.bound`	Function applied to calculate the lower bound the final forecast. Default: `min`.
`upper.bound`	Function applied to calculate the upper bound of the final forecast. Default: `max`.

forecast returns an object of class "forecast", which contains the following components:

yhat A matrix of forecast point estimates (age groups by time periods).
yhat.lower A matrix of forecast lower bound estimates (age groups by time periods).
yhat.upper A matrix of forecast upper bound estimates (age groups by time periods).
y A matrix of observed data (age groups by time periods).
aux List. A list of inputs to yourcast.
validation List. A list containing all sigma combinations and the associated diagnostics.
weights The "weightcast" object provided as an input to forecast.

data(netherlands)
ff <- log(brst3/popu3) ~ log(hc) + log(gdp) + log(tobacco3) + log(fat) + time
out <- gridcast(netherlands_data, formula=ff, model="map",
                           sample.frame=c(1950,2000,2001,2030), verbose=FALSE)
netherlands_forecast <- forecast(out, netherlands_weights)

ageplot(netherlands_forecast)
timeplot(netherlands_forecast)