我有以下代码返回具有最低 AIC 的模型,但我希望所有模型的 AIC 以升序或降序排列,而不使用 R 中的内置排序函数
sp <- rnorm(100) ## just some toy data to make code work!
spfinal.aic <- Inf
spfinal.order <- c(0,0,0)
for (p in 1:4) for (d in 0:1) for (q in 1:4) {
spcurrent.aic <- AIC(arima(sp, order=c(p, d, q)))
if (spcurrent.aic < spfinal.aic) {
spfinal.aic <- spcurrent.aic
spfinal.order <- c(p, d, q)
spfinal.arima <- arima(sp, order=spfinal.order)
}
}
我想spfinal.order<-c(p,d,p)
成为 AIC 升序或降序的所有模型的列表。我怎样才能做到这一点?