data <- c(AirPassengers)
level <- polynomial_block(rate = 1, order = 2, D = 0.95)
season <- harmonic_block(rate = 1, order = 2, period = 12, D = 0.975)
outcome <- Poisson(lambda = "rate", data)
fitted.data <- fit_model(level, season,
AirPassengers = outcome
)
model.coef <- coef(fitted.data)
plot(model.coef)$plot
Run the code above in your browser using DataLab