《应用时间序列分析：R软件陪同》——1.5 习题

本节书摘来自华章计算机《应用时间序列分析：R软件陪同》一书中的第1章，第1.5节,作者：吴喜之，刘苗著， 更多章节内容可以访问云栖社区“华章计算机”公众号查看。

1.5 习题

1.从网上下载R软件.

2.在联网的情况下,打开R,点击“帮助-CRAN主页-Packages”看看今天有多少可供使用的程序包(package),再点击按照时间顺序排列的程序包列表,看看今天又增加了多少新程序包(包括更新的).

3.在联网的情况下,打开R,点击“帮助-CRAN主页-TaskViews”看看有多少你感兴趣的领域或方向,再点击你感兴趣的领域,看看有多少软件包可用.

4.在R中(不必要联网),通过语句?lm来看线性模型(linearmodel)函数lm()的用法,然后看说明后面的例子,把这些例子的代码逐行粘贴到R的运行界面,查看结果.

5.在R中,用语句mylm=fix(lm)或者mylm=edit(lm)就可以把函数lm()变成你自己的名为mylm()的函数,而且还可以修改.

6.在联网状态,在R菜单中点击“程序包–安装程序包”后,R会要求你选择镜像网站,选择完了就会出现一个很长的程序包列表,你可以选择一个程序包或多个程序包(按住Ctrl键)来安装,这种安装会自动安装这些程序包以及它们所需要的支持程序包(即使你没有选).

7.你也可以通过键入命令来安装程序包.你可以试着键入下面的命令安装一个和多个程序包:
install.packages("Ecdat"); install.packages(c("TSA","vars"))

这时,R会要求你选择镜像网站,选择完了就会自动安装这些程序包以及它们所需要的支持程序包(即使你没有选).

8.根据上题在R中安装了程序包Ecdat,用library(Ecdat);?Macrodat打开程序包Ecdat所提供的数据Macrodat的说明,然后用plot(Macrodat)语句画出有关的7个时间序列的图.当然,你也可以用语句plot(Macrodat,plot.type="single")把它们画到一张图上,看看有什么区别.

9.在上一题最后一个画图语句plot(Macrodat,plot.type="single")中增加一个选项:type="o",也就是换成plot(Macrodat,plot.type="single",type="o"),看看会产生何种不同的点图.由此,根据帮助?plot来了解相应于type的不同选项所产生的不同图形模式.

10.建立工作目录是很方便的,如果有了工作目录,存取文件都在这个目录文件夹之下,不用再指名路径.有两种方式建立工作路径:

(a)敲入诸如setwd("D:/mywork")设立你自己的工作目录(当然这个文件夹必须存在).

(b)在R的菜单中,点击“文件–改变工作目录”后,会出现路径窗口,你自己选择适合的工作目录文件夹.
如果你不知道目前的工作目录是什么,可以用getwd()来获得.

11.存取文件.在选择工作目录,并且用library(Ecdat)载入程序包verbEcdat+之后,先点击“文件–新建程序脚本”,就会出现一个编辑窗口.然后用语
句write.csv(Macrodat,"Macrodat.csv",row.names=F)把数据Macrodat存入硬盘的目录中.你再到该目录下寻找这个文件,并打开查看,你会发现没有时间存入(因为原数据就没有时间列).如果要提取这个文件的数据,加上数据说明中注明的时间段,可以用语句w=read.csv("Macrodat.csv");w=ts(w,start=c(1959,1),frequency=4)得到.这时点击w就可以看到有时间列的数据.在关闭程序脚本前,可以把这个脚本存为以“.R”结尾的文件.

12.你可以用命令history()来查看工作时所用命令的记录,可以把这个记录存到你的硬盘上,以后还可以重复原先的工作.