Description Usage Arguments Value Author(s) See Also Examples
Re-fit a "hhh4"
model with a modified control
list.
1 2 3 |
object |
a fitted |
... |
components modifying the original control list for
|
S |
a named list of numeric vectors serving as argument for
|
subset.upper |
if a scalar value, restrict the new fit to those epochs of
|
use.estimates |
logical specifying if |
evaluate |
logical indicating if the updated model should be fitted directly
(defaults to |
If evaluate = TRUE
the re-fitted object, otherwise the updated
control
list for hhh4
.
Sebastian Meyer
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 | data("salmonella.agona")
## convert to sts class
salmonella <- disProg2sts(salmonella.agona)
## fit a basic model
fit0 <- hhh4(salmonella,
list(ar = list(f = ~1), end = list(f = addSeason2formula(~t))))
## update: Poisson -> NegBin1, component seasonality
fit1 <- update(fit0, family = "NegBin1", S = list(end=2, ar=2))
## compare fits
AIC(fit0, fit1)
opar <- par(mfrow=c(2,2))
plot(fit0, type="fitted", names="fit0", par.settings=NULL)
plot(fit1, type="fitted", names="fit1", par.settings=NULL)
plot(fit0, fit1, type="season", components=c("end", "ar"), par.settings=NULL)
par(opar)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.