- library(readxl)
- library(tseries)
- hx <- read_xlsx("时序.xlsx")
- hx <- ts(hx$词频)
- plot(hx,col=4,lwd=2,pch=8,type="o")
- #ADF检验和白噪声LB检验【在这一步报错】
- adf.test(hx)
- for(i in 1:4) print(Box.test(hx,type="Ljung-Box",lag=3*i))
- #绘制自相关和偏自相关图
- opar<-par(mfrow=c(2,1))
- acf(hx,col=4,lwd=2)
- pacf(hx,col=4,lwd=2)
- par(opar)