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

应用时间序列分析上机实验教案 第一次实验：（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.2Exercise 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，点大小为； 命令：win.graph(width=4.875, height=2.5,pointsize=) plot(X,xlab=序号,ylab=10-100的整十倍数,type=’o’) 注解：“width=4.875”--宽度为4.875，“height=2.5”--高度为2.5，“pointsize=”--点大小为；运行结果： Exercise 1.5 要求：图 案例：图 命令：win.graph(width=4.875, height=4.5,pointsize=10) plot(X,Y) 运行结果： Exercise 1.6 要求：Exercise 1.5的图 案例：图 命令： 运行结果： 1.6 Exercise 1.7 要求： 案例： 操作：先将数据源转换为“*.TXT”文件，比如：“GDP78-12.TXT”；保存于路径“E:\”下；在R中修改当前目录于“E:\”； 命令：“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 要求： 案例： 操作：先将数据源保存于“1000.TXT”文件，比如：保存于路径“E:\”下（为方便起见，其它类同）；在R中修改当前目录于“E:\”； 命令：“1000.txt”);X1 运行结果：[1] 100 200 300 400 500 600 700 1800 900 1000 Exercise 1.9 要求： 操作：先打开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 要求： 操作：若第一次作图，先打开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) 运行结果