r语言时间序列分析

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

1、R语言时间序列初探！说明：作者：AvrilCoghlanEmail：alc@sanger.ac.uk原文：http://a-little-book-of-r-for-time-series.readthedocs.org/en/latest/src/timeseries.html本文档由梁德明、赵华蕾合译；联系方式：hunt4game@gmail.com；内容遵循许可协议CC3.0BY本文档PDF版本地址为：http://doc.datapanda.net/a-Little-Book-of-R-f

2、or-Time-Series.pdf本文档中所提及数据集，如出现无法访问，只需将robjhyndman.com替换为doc.datapanda.net即可这是一份关于使用R统计软件进行时间序列分析的入门文档。内容大纲：使用R进行实现序列分析o时间序列分析o读取时间序列数据o绘制时间曲线图o分解时间序列分解非季节性数据分解季节性数据季节性的修正o使用指数平滑法进行预测简单指数平滑法Holt指数平滑法Holt-Winters指数平滑法oARIMA模型时间序列的差分选择一个合适的ARIMA模

3、型使用ARIMA模型进行预测o链接与拓展阅读时间序列分析这本小册子将告诉你如何使用R软件实现对常见时间序列数据进行简单分析。这本小册子假定读者已经有一定时间序列分析的基础，并且本册的侧重点并不在于解释时间序列分析，而是在于如何使用R软件实现这些分析。如果你是一位学习时间序列的新手，并且很想学习更多关于本册中提及的概念，我强烈推荐英国公开大学的《Timeseries》（产品编号M249/02）,你可以从英国公开大学商店购买到。在本册中，我使用的时间序列数据是由RobHyndman的时间序列数据库（http

4、://robjhyndman.com/TSDL/）提供的。本册的PDF版本（英文）在：https://media.readthedocs.org/pdf/a-little-book-of-r-for-time-series/latest/a-little-book-of-r-for-time-series.pdf如果你喜欢本册，也许你应该瞧瞧这本《R语言与生物统计》（英文）（http://a-little-book-of-r-for-biomedical-statistics.readthedocs.or

5、g/）是否感兴趣，还有我的另一本《R语言与多元统计分析》（英文）（http://little-book-of-r-for-multivariate-analysis.readthedocs.org/）读取时间序列数据当你想要分析时间序列数据时，第一件事就是将数据读入R软件，并且绘制时间序列图。假设你的数据是以连续的时间点的形式存储在简单的文本文件的一列里，你可以使用scan()函数将数据读入到R。例如，这个文件（http://robjhyndman.com/tsdldata/misc/kings.dat）

6、包含着从威廉一世开始的英国国王的去世年龄数据。（原始出处：HipelandMcleod,1994）数据集如下：AgeofDeathofSuccessiveKingsofEngland#startingwithWilliamtheConqueror#Source:McNeill,"InteractiveDataAnalysis"604367505642506568436534...这里仅显示了文件的一些行。第一个三行包含了一些数据的注解内容，在我们读入数据到R的时候，直接忽略这些。我们可以使用sacn()函

7、数的”skip”参数指定文件中从顶部开始有多少行需要忽略。为了将数据读入到R，并且忽略掉文件中的前三行，我们输入以下代码：>kings<-scan("http://robjhyndman.com/tsdldata/misc/kings.dat",skip=3)Read42items>kings[1]60436750564250656843653447344941133553561643695948[26]5986556851334967778167718168707756在这个案例中有42位英国国王的去世

8、年龄数据被读入到变量”kings”中。一旦你将时间序列数据读入到R，下一步便是将数据存储到R中的一个时间序列对象里，以便你能使用R的很多函数分析时间序列数据。我们在R中使用ts()函数将数据存储到一个时间序列对象中去。例如存储“kings”这个变量中的数据到时间序列对象中，我们输入以下代码：>kingstimeseries<-ts(kings)>kingstimeseriesTimeSeries:Start=1End=42Fr