level <- polynomial_block(rate = 1, order = 2, D = 0.95)
season <- harmonic_block(rate = 1, order = 2, period = 12, D = 0.975)
# Only first 100 observations (for the sake of the example)
outcome <- Poisson(lambda = "rate", data = c(AirPassengers)[1:100])
fitted.data <- fit_model(level, season,
AirPassengers = outcome
)
updated.fit <- update(fitted.data, AirPassengers = list(data = c(AirPassengers)[101:144]))
# If a offset was present, the user should pass its value when updating
# updated.fit=update(fitted.data,
# AirPassengers=list(
# data=c(AirPassengers)[101:144],
# offset= ... ))
Run the code above in your browser using DataLab