R-时间序列分析

                                      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=12)

 

绘制时序图

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

 

 

boxplot(series) ---箱线图

boxplot.stats(series)

检查含有的异常值:

value = boxplot.stats(series)$out

index = which(data$num %in% value)