r语言时间序列中文教程

《r语言时间序列中文教程》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、39R语言时间序列中文教程2012特别声明：R语言是免费语言，其代码不带任何质量保证，使用R语言所产生的后果由使用者负全责。前言R语言是一种数据分析语言，它是科学的免费的数据分析语言，是凝聚了众多研究人员心血的成熟的使用范围广泛全面的语言，也是学习者能较快受益的语言。在R语言出现之前，数据分析的编程语言是SAS。当时SAS的功能比较有限。在贝尔实验室里，有一群科学家讨论提到，他们研究过程中需要用到数据分析软件。SAS的局限也限制了他们的研究。于是他们想，我们贝尔实验室的研究历史要比SAS长好几倍，技术力量也比SAS强好几

2、倍，且贝尔实验室里并不缺乏训练有素的专业编程人员，那么，我们贝尔实验室为什么不自己编写数据分析语言，来满足我们应用中所需要的特殊要求呢？于是，贝尔实验室研究出了S-PLUS语言。后来，新西兰奥克兰大学的两位教授非常青睐S-PLUS的广泛性能。他们决定重新编写与S-PLUS相似的语言，并且使之免费，提供给全世界所有相关研究人员使用。于是，在这两位教授努力下，一种叫做R的语言在奥克兰大学诞生了。R基本上是S-PLUS的翻版，但R是免费的语言，所有编程研究人员都可以对R语言做出贡献，且他们已经将大量研究成果写成了R命令或脚本，

3、因而R语言的功能比较强大，比较全面。研究人员可免费使用R语言，可通过阅读R语言脚本源代码，学习其他人的研究成果。笔者曾有幸在奥克兰大学受过几年熏陶，曾经向一位统计系的老师提请教过一个数据模拟方面的问题。那位老师只用一行R语句就解答了。R语言的强大功能非常令人惊讶。为了进一步推广R语言，为了方便更多研究人员学习使用R语言，我们收集了R语言时间序列分析实例，以供大家了解和学习使用。当然，这是非常简单的模仿练习，具体操作是，用复制粘贴把本材料中R代码放入R的编程环境；材料中蓝色背景的内容是相关代码和相应输出结果。经过反复模仿，

4、学习者便能熟悉和学会。需要提醒学习者的是：建议学习者安装了R语言编程，再继续阅读本材料；执行R命令时，请删除命令的中文注解，没使用过在命令中加入中文；如果学习者是初次接触R或者Splus，建议先阅读<>，如果学习者比较熟悉R语言，还可以阅读优秀时间序列读物EcomometricsinR，也可以上QuickR网站。39目录R语言时间序列分析1前言1目录21.运用R语言研究JJ数据32.运用R语言研究空气污染193.运用R语言研究自动回归移动平均集成模型224.运用R语言研究全球变暖理论315.运用

5、R语言研究非独立误差与线性回归326.运用R语言研究估计波动的频率367.运用R语言研究厄尔尼诺频率398.运用R语言研究太阳黑子周期频率40391.运用R语言研究JJ数据学习R言语时间序列分析程序操作，需要从最基础、最简单的学起，例如在命令窗口中，输入并执行2+2等于4的R语言命令。2+2[1]4执行完2+2等于4的R语言命令后，我们可以开始时间序列，即学着把玩johnson&Johnson数据。下载jj.dat或执行下面语句。这个数据已被人上传到因特网中。R所需要做的只是将网址进行扫描就可以将数据读取进入R的编程环境

6、中。下面有3种不同读取数据的方法：jj=scan("http://www.stat.pitt.edu/stoffer/tsa2/data/jj.dat")#readthedata读取数据jj<-scan("http://www.stat.pitt.edu/stoffer/tsa2/data/jj.dat")#readthedataanotherway第二种方法读取数据scan("http://www.stat.pitt.edu/stoffer/tsa2/data/jj.dat")->jj#andanother第三种方法读

7、取数据使用R语言的人，有的喜欢使用[<-]，有的喜欢使用[->]，大多数医疗系统的工作者喜欢用[=]，正因为如此才用了上面种不同读取数据的方法。读取数据后，键入并执行jj，数据在窗口便会有如下显示：jj[1]0.710.630.850.44[5]0.610.690.920.55..........[77]14.0412.9614.859.99[81]16.2014.6716.0211.61jj中有84个数据被称作对象。下面命令可以显示所有对象。objects()39如果使用matlab，你会认为jj是一个84行1列的向量

8、，但实际上不是这样。jj有次序，有长度，但没维度，R称这些对象为向量，要小心区别。在R里，矩阵有维度，但向量没维度。这都是程序语言的一些概念。jj[1]#thefirstelement列中第一个数据[1]0.71jj[84]#thelastelement列中最后一个数据[1]11.61jj[1:4]#thefirst