time series analysis

ARIMA模型

读取数据

data <-read.table(file='D:/python/anaconda/ARIMA/airPassenger.csv',header=T,sep=',',row.names='time')

数据格式如下:

num

1949/1/1 112

1949/2/1 118

1949/3/1 132

1949/4/1 129

1949/5/1 121

1949/6/1 135

row.names(data) ---返回索引值(时间)

data$num ---返回变量值

****将数据转为ts对象 ****

****series <-ts(data,start=c(1949,1),end=c(1959,12),frequency=”periodic”) ****

****绘制时序图 ****

****plot.ts(series,xlab='时间',ylab='乘客数/个') ****

****boxplot(series) ---箱线图