# NOT RUN {
data(warming)
model.poly1 = polyDlm(x = warming$NoMotorVehicles , y = warming$Warming ,
q = 2 , k = 2 , show.beta = TRUE , show.summary = TRUE)
model.poly2 = polyDlm(x = warming$NoMotorVehicles , y = warming$Warming ,
q = 3 , k = 2 , show.beta = TRUE , show.summary = TRUE)
model.poly3 = polyDlm(x = warming$NoMotorVehicles , y = warming$Warming ,
q = 4 , k = 2 , show.beta = TRUE , show.summary = TRUE)
aic = AIC(model.poly1$model, model.poly2$model, model.poly3$model)
bic = BIC(model.poly1$model, model.poly2$model, model.poly3$model)
mase = MASE(model.poly1$model, model.poly2$model, model.poly3$model)
sortScore(aic , score = "aic")
sortScore(bic , score = "bic")
sortScore(mase , score = "mase")
# }
Run the code above in your browser using DataLab