model <- SSModel(log(drivers) ~ SSMtrend(1, Q = list(1)) +
SSMseasonal(period = 12, sea.type = "trigonometric") +
log(PetrolPrice) + law, data = Seatbelts, H = 1)
out <- KFS(model)
confint(out, parm = "regression")
Run the code above in your browser using DataLab