应用时间序列分析上机实验-

应用时间序列分析上机实验教案
第一次实验：（2014.9.22）
Exercise 1.1
要求：在R软件中命名、输入时间序列数据，并图示；
案例：在R软件中将样本容量为10的数据：10，20，30，40，50，60，70，80，90，100赋值给变量X；
命令：X<-c(10,20,30,40,50,60,70,80,90,100)
Plot(X)
运行结果：
图1.1
Exercise 1.2
要求：在Exercise 1.1的图1.1中对X轴与Y轴重新标注；
案例：在图1.1中将X轴标注为：“序号”，将Y轴标注为：“10-100的整10倍数”
命令：X<-c(10,20,30,40,50,60,70,80,90,100) （略，若没有上前的命令，要写）plot(X,xlab='序号',ylab='10-100的整十倍数')
注解：注意：1、xlab与ylab的标注是单引号；2、汉字不能用五笔输入法；
运行结果：
图1.2
Exercise 1.3
要求：将Exercise 1.2的图1.2转换为拆线图；
案例：在图1.2中的plot命令中加入选项“type=’o’；
命令：plot(X,xlab='序号',ylab='10-100的整十倍数',type=’o’)
注解：没有“type=’o’”时，图形为点图，加入后将点连接为线，即拆线；
运行结果：
图1.3
Exercise 1.4
要求：改变Exercise 1.3的图1.3的高与宽及图中点的大小；
案例：图1.4的高度为2.5，宽度为4.875，点大小为10；
命令：win.graph(width=4.875, height=2.5,pointsize=10)
plot(X,xlab='序号',ylab='10-100的整十倍数',type=’o’)
注解：“width=4.875”--宽度为4.875，“height=2.5”--高度为2.5，“pointsize=10”--点大小为10；运行
结果：
图1.4
Exercise 1.5
要求：画出两个变量的散点图；
案例：输入变量X、Y，画散点图；
命令：X<-c(10,20,30,40,50,60,70,80,90,100)
Y<-c(100,201,300,410,500,1200,710,820,900,1000)
win.graph(width=4.875, height=4.5,pointsize=10)
plot(X,Y)
运行结果：
图1.5
Exercise 1.6
要求：在Exercise 1.5的图1.5中加入X轴、Y轴标注；
案例：将散点图1.5的X轴标注为“整十倍数”、Y轴标注为“奇异值”；命令：plot(X,Y,xlab='整十倍数',ylab='奇异值')
运行结果：
图1.6
Exercise 1.7
要求：从“*.TXT”文件中读取数据；
案例：将1978-2012年我国的GDP时间序列数据赋给变量GDP；
操作：先将数据源转换为“*.TXT”文件，比如：“GDP78-12.TXT”；保存于路径“E:\”下；在R 中修改当前目录于“E:\”；
命令：点击“文件”/“改变当前目录”
GDP<-read.table(“GDP78-12.txt”,head=TRUE)
GDP
运行结果：
GDP
1978年3,645.22
1979年4,062.58
1980年4,545.62
1981年4,891.56
1982年5,323.35
.
.
.
Exercise 1.8
要求：利用scan( )函数从“*.TXT”文件中读取数据；
案例：将100，200，300，400，500，600，700，1800，900，1000存于文件“1000.TXT”
中，利用scan( )函数从文件“1000.TXT”中读取数据，赋于变量X1；
操作：先将数据源保存于“1000.TXT”文件，比如：保存于路径“E:\”下（为方便起见，其它类同）；在R中修改当前目录于“E:\”；
命令：点击“文件”/“改变当前目录”
X1<-scan(“1000.txt”);X1
运行结果：
[1] 100 200 300 400 500 600 700 1800 900 1000
Exercise 1.9
要求：在R中作出P1中图E1.1；
操作：先打开library(TSA)---TSA是R的外部程序包，设定图宽、高等，调用数据--data(larain)，作图plot(larain,ylab='Inches',xlab='Year',type='o')
命令：library(TSA)
win.graph(width=4.875, height=2.5,pointsize=8)
data(larain); plot(larain,ylab='Inches',xlab='Year',type='o')
运行结果：
Exercise 1.10
要求：在R中作出P2中图E1.2；
操作：若第一次作图，先打开library(TSA)，第二次即不需要，定图宽、高等，调用数据--data(larain)，作图plot(larain,ylab='Inches',xlab='Year',type='o')
命令：library(TSA)
win.graph(width=3,height=3,pointsize=8)
plot(y=larain,x=zlag(larain),ylab='Inches', xlab='Previous Year Inches')
运行结果：
注：字符用汉字：
library(TSA)
win.graph(width=3,height=3,pointsize=8)
plot(y=larain,x=zlag(larain),ylab='英寸', xlab='上一年英寸')
Exercise 1.11
要求：在R中作出P2中图E1.3；
命令：library(TSA)
win.graph(width=4.875, height=2.5,pointsize=8)
data(color)
plot(color,ylab='Color Property',xlab='Batch',type='o') 运行结果：
Exercise 1.12
要求：在R中作出P3中图E1.4；
命令：library(TSA)
win.graph(width=3,height=3,pointsize=8)
plot(y=color,x=zlag(color),ylab='Color Property', xlab='Previous Batch Color Property') 运行结果：。