(1,1,1)把。

你可以自己写个程序来得到AIC/BIC值最小时的p,d值。。。

myfunction<-function(p,d,q,data){
AIC=99999999999
P=NULL
Q=NULL
for(i in 1:p){
for(j in 1:q){
fit<-arima(data,order=c(i,d,j),method="ML")
if(AIC>AIC(fit)){
P<-i
Q<-j
AIC<-AIC(fit)
}}}
L=cbind(P,Q,AIC)
return(L)
}
##使用这个函数语法 myfunction(3,1,3,ap)##
如果你要显示出所有的AIC值，你可以这样：
myfunction<-function(p,d,q,data){
AIC=NULL
aic=99999999999
P=NULL
Q=NULL
for(i in 0:p){
for(j in 0:q){
fit<-arima(data,order=c(i,d,j),method="ML")
AIC=cbind(AIC,AIC(fit))
if(aic>AIC(fit)){
P<-i
Q<-j
aic<-AIC(fit)
}}}
L=cbind(aic,P,Q,AIC)
return(L)
}
其中aic为最小AIC值，P,Q为最小AIC时对应的P,Q值，AIC存储了所有的值