一、导入数据
num<-read.csv("D:\\nhtemp.csv",header = T)
num

二、画出时序图
number<-ts(num[,2],start = 1912)
number
plot(number)

三、平稳性检验和白噪声检验
ADF检验
library(aTSA)
adf.test(number, nlag = 2)
从ADF检验结果上看,在95%的显著性水平下,P<0.05,拒绝原假设,该序列为平稳性序列。
白噪声检验(纯随机性检验)
lag=c()
Q.value=