R语言利用arima模型的时间序列分析-代码

代码如下

library(readxl)
library(forecast)
dat <- read_excel("rawdatacompet.xlsx",sheet=1,na="NA")
#对数据进行处理，取出待分析的数据并转化
series1 <- dat[,1]
series1 <- series1[1:60,]
series1 <- ts(series1)
series1 <- as.double(series1)
#自动分析适合该数据的参数值，并进行预测
fore <- forecast(auto.arima(series1,trace = T),h=5)
fore
plot(fore)

结果如图
数据如下（蓝色部分表示待预测值）：

代码系编程小白所写，如有错漏，敬请指正。