# 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)
sortScore(aic,score = "aic")
sortScore(bic,score = "bic")
# }
Run the code above in your browser using DataLab