R：应用时间序列分析--基于R（1）第一章 时间序列分析简介

注：教材为《时间序列分析——基于R》 王燕 编著 中国人民大学出版社

1.1 引言

1.2时间序列的定义

1.3 时间序列的分析方法

1.3.1 描述性时序分析

1.3.2 统计时序分析

频域分析方法（谱分析法）

时域分析方法

1.4 R简介

安装程序包

install.packages("tseries")##安装程序包
library(tseries)###调用程序包

赋值

x<-3
y<-2*x+1
z<-y^2

x
y
z

输入时间序列

price<-c(101,82,66,35,31,7)
price<-ts(price,start=c(2005,1),frequency = 12)
price
price<-scan()  ###使用列输入
price<-ts(price,start=c(2005,1),frequency = 12)

x<-read.table("C:/Users/Daisy/Desktop/时间序列/data/file1.csv",sep=",",header=T)
x
yield
x$yield

yield<-ts(x$yield,start=1884,frequency=1)
yield

x$yield<-ts(x$yield,start=1884,frequency=1)
x$yield

#取对数
y<-log(x$yield)
y

取子集

z<-subset(x,year>1925,select=yield)
z

缺失值插值

install.packages("zoo")
library(zoo)
a<-c(1:7)
a[4]<-NA
a

线性插值

y1<-na.approx(a)
y1

样条插值

y2<-na.spline(a)
y2

绘制时序图

plot(x$yield,type="o")

绘制各种类型时序图

par(mfrow=c(2,3))
plot(x$yield,type="p",main='tpye="p"')
plot(x$yield,type="l",main='tpye="l"')
plot(x$yield,type="b",main='tpye="b"')
plot(x$yield,type="o",main='tpye="o"')
plot(x$yield,type="h",main='tpye="h"')
plot(x$yield,type="s",main='tpye="s"')

绘制各种符号类型的时序图

par(mfrow=c(2,2))
plot(x$yield,type="o",pch=2,main='pch=2')
plot(x$yield,type="o",pch=9,main='pch=9')
plot(x$yield,type="o",pch=17,main='pch=17')
plot(x$yield,type="o",pch=25,main='pch=25')

绘制各种连线类型的时序图

par(mfrow=c(1,2))
plot(x$yield,lty=2,main='lty=2')
plot(x$yield,lty=6,main='lty=6')

绘制不同线宽的时序图

par(mfrow=c(1,3))
plot(x$yield,lwd=1,main='lwd=1')
plot(x$yield,lwd=2,main='lwd=2')
plot(x$yield,lwd=-2,main='lwd=-2')

绘制各种符号类型和颜色的时序图

par(mfrow=c(2,2))
plot(x$yield,type="o",pch=2, col=1, main='pch=2,col=1')
plot(x$yield,type="o",pch=9,col=2, main='pch=9,col=2')
plot(x$yield,type="o",pch=17,col=3,main='pch=17,col=3')
plot(x$yield,type="o",pch=25,col=4,bg="yellow",main='pch=25,col=4,bg=yellow')

绘制添加文本时序图

par(mfrow=c(1,1))
plot(x$yield,type="o",main='yield时序图', xlab="年份", ylab="yield")

分别制定横坐标和纵坐标范围

par(mfrow=c(1,2))
plot(x$yield,xlim=c(1910,1930),main='横轴范围1910-1930')
plot(x$yield,ylim=c(15,17),main='纵轴范围15-17')

添加参照线

par(mfrow=c(1,1))
plot(x$yield,type="o")
abline(v=1915,h=15,lty=2)

对报错进行修正

install.packages("readxl")
library(readxl)
x<-read_excel("F:/应用时间序列分析/时间序列分析——基于R（第2版）案例数据/A1_1.xlsx",sheet="Sheet1",Header=T)
x
yield
x$yield

yield<-ts(x$yield,start=1884,frequency=1)
yield

x$yield<-ts(x$yield,start=1884,frequency=1)
x$yield

#取对数
y<-log(x$yield)
y

#取子集
z<-subset(x,Year>1925,select=yield)
z

#绘制时序图
plot(x$yield,type="o")

#绘制各种类型时序图
par(mfrow=c(2,3))
win.graph(width=4.875, height=2.5,pointsize=8)
plot(x$yield,type="p",main='tpye="p"')
plot(x$yield,type="l",main='tpye="l"')
plot(x$yield,type="b",main='tpye="b"')
plot(x$yield,type="o",main='tpye="o"')
plot(x$yield,type="h",main='tpye="h"')
plot(x$yield,type="s",main='tpye="s"')

#绘制各种符号类型的时序图
par(mfrow=c(2,2))
win.graph(width=4.875, height=2.5,pointsize=8)
plot(x$yield,type="o",pch=2,main='pch=2')
plot(x$yield,type="o",pch=9,main='pch=9')
plot(x$yield,type="o",pch=17,main='pch=17')
plot(x$yield,type="o",pch=25,main='pch=25')

#绘制各种连线类型的时序图
par(mfrow=c(1,2))
win.graph(width=4.875, height=2.5,pointsize=8)
plot(x$yield,lty=2,main='lty=2')
plot(x$yield,lty=6,main='lty=6')

#绘制不同线宽的时序图
par(mfrow=c(1,3))
win.graph(width=4.875, height=2.5,pointsize=8)
plot(x$yield,lwd=1,main='lwd=1')
plot(x$yield,lwd=2,main='lwd=2')
plot(x$yield,lwd=-2,main='lwd=-2')

#绘制各种符号类型和颜色的时序图
par(mfrow=c(2,2))
win.graph(width=4.875, height=2.5,pointsize=8)
plot(x$yield,type="o",pch=2, col=1, main='pch=2,col=1')
plot(x$yield,type="o",pch=9,col=2, main='pch=9,col=2')
plot(x$yield,type="o",pch=17,col=3,main='pch=17,col=3')
plot(x$yield,type="o",pch=25,col=4,bg="yellow",main='pch=25,col=4,bg=yellow')

#绘制添加文本时序图
par(mfrow=c(1,1))
win.graph(width=4.875, height=2.5,pointsize=8)
plot(x$yield,type="o",main='yield时序图', xlab="年份", ylab="yield")

#分别制定横坐标和纵坐标范围
par(mfrow=c(1,2))
plot(x$yield,xlim=c(1910,1930),main='横轴范围1910-1930')
plot(x$yield,ylim=c(15,17),main='纵轴范围15-17')

#添加参照线
par(mfrow=c(1,1))
win.graph(width=4.875, height=2.5,pointsize=8)
plot(x$yield,type="o")
abline(v=1915,h=15,lty=2)