R-4.0.2 forecast预测


library("forecast")

skirts <- scan("http://robjhyndman.com/tsdldata/roberts/skirts.dat",skip=5)
skirtsts<- ts(skirts,start = c(1866))
plot.ts(skirtsts)

skirtsarima<-arima(skirtsts,order=c(1,2,5))
skirtsarima2<-arima(skirtsts,order=c(2,2,0))
skirtsarima3<-arima(skirtsts,order=c(2,2,1))
skirtsarima

skirtsarimadict.fore<-forecast(skirtsarima,h=10)

plot(skirtsarimadict.fore)

R-4.0.2 forecast预测_第1张图片

你可能感兴趣的:(R语言)