library(uotm)
ts <- arma.sim(ars = c(-0.9, -1.4, -0.7, -0.6), mas = c(0.5, -0.4), nobs = 100)
tmcb1 <- arma.mcb(ts, method = "aic", max.p = 4, max.q = 4,
bsamples = 50, stepwise = TRUE, blength = 12)
arma.muc(c("aic"), tmcb1)
# compare the uncertainty of different model selection methods using muc
tmcb2 <- arma.mcb(ts, method = "bic", max.p = 4, max.q = 4,
bsamples = 50, stepwise = TRUE, blength = 12)
arma.muc(c("aic", "bic"), tmcb1, tmcb2)
Run the code above in your browser using DataLab