时间序列实验报告

重庆交通大学学生实验报告实验课程名称时间序列剖析开课实验中心数统学院实验教课中心开课学院数学与统计学院专业年级应用统计学2015级1班姓名XXXX学号6315XXXXXXXX任课老师XXXXX开课时间2017—2018学年第1学期此页空页！实验一R语言简介:基本操作一实验目的1、认识软件R：安装、启动、退出、帮助等。2、熟习R的操作界面。二、实验内容及要求：1、实验内容：1）R的安装；2）启动与退出；3）包的安装及R的更新；4）帮助及移除多个对象等；5）常有命令2、实验要求：1）熟习R的操作环境；2）熟习包的安装与帮助；3）学习常有命令，熟习R的操作界面。三、实验过程及结果1、（1）R的安装（2）启动与退出；（3）包的安装及R的更新；A、包的安装chooseCRANmirror( )( )B、R的更新("installr")library(installr)updateR( )（4）帮助及移除多个对象等；重点字重点字( )#帮助rm( )rm(list=ls( ))#移除多个对象（5）常有命令四、实验心得认识了R的一些基本使用及其常有的命令，为自己深入学习r的使用打下了基础。实验二R语言简介:数据集创立与办理一实验目的1、掌握R数据集的不一样创立形式。2、熟习并掌握利用R对时间序列数据集进行变换与办理。二、实验内容及要求1、实验内容：（1）利用函数创立数据集；（2）读取型数据框；（3）读取excel数据及对某变量数据进行某些办理（4）导出R中数据集（5）时间序列数据输入（6）对已有数据集中数据的办理2、实验要求：熟习R数据集的不一样创立方法，掌握利用R对时序数据集进行变换与办理三、实验过程及结果1、实验内容：（1）利用函数创立数据集；>mydata<(col1,col2,col3,...)#col1，col2，col3，...能够是任何种类的列向量（2）读取型数据框；m>w<("")（3）读取excel数据及对某变量数据进行某些办理library(readxl)dataset<-read_excel(NULL)>View(dataset)（4）导出R中数据集(x,file="",sep="",=TRUE,=TRUE,quote=TRUE)#x：需要导出的数据#file：导出的文件路径#sep：分开符，默以为空格（""），也就是以空格为切割列#：能否导出行序号，默以为TRUE，也就是导出行序号#：能否导出列名，默以为TRUE，也就是导出列名#quote：字符串能否使用引号表示，默以为TRUE，也就是使用引号表示（5）时间序列数据输入（6）对已有数据集中数据的办理（5）（6）合library(readxl)X2_7<-read_excel("C:/Users/Administrator/Desktop/")summary(X2_7)Min.:Min.:Min.:Min.:Min.:1stQu.:1stQu.:1stQu.:1stQu.:1stQu.:Median:Median:Median:Median:Median:Mean:Mean:Mean:Mean:Mean:3rdQu.:3rdQu.:3rdQu.:3rdQu.:3rdQu.:Max.:Max.:Max.:Max.:Max.:Min.:1stQu.:Median:Mean:3rdQu.:Max.:四、实验心得经过本次实验，第一，我知道了文件其余格式的文件怎样导入R，了解乐数据集的创立，使用及一些简单的办理。实验三安稳时间序列模型预办理一实验目的1、掌握安稳时间序列的特色2、学会经过序列图形对时间序列进行安稳性查验和纯随机性查验预办理二、实验内容及要求：1、实验内容：（1）导入课后作业（41页第2题）数据并绘制时序图，并进行安稳性查验和纯随机性查验。（2）导入课后作业（41页第3题）数据并绘制时序图，并进行安稳性查验和纯随机性查验。（3）导入课后作业（41页第5题）数据并绘制时序图，并进行安稳性查验和纯随机性查验。（4）导入课后作业（41页第6题）数据并绘制时序图，进行安稳性查验和纯随机性查验，并对序列进行差分，再做安稳性查验与纯随机性查验。2、实验要求：熟习经过序列图形对时间序列进行安稳性查验和纯随机性查验预办理三、实验过程及结果、实验内容：1）导入课后作业（41页第1题）数据并绘制时序图，并进行安稳性查验和纯随机性查验。a<-scan('3(1).txt')###在工作目录下调用txt文件3（1）###ppm<-ts(a,start=c(1975,1),frequency=12)####转变时序数据###plot(ppm)####绘制时序图######acf(ppm)###安稳性查验####for(iin1:2)print(ppm,lag=6*i))####纯随机性查验###时序图安稳性查验Seriesppm00148380336603F4mp4A3322033000332802319751976197719781979198019810.00.51.01.5TimeLag纯随机性查验结果Box-Piercetestdata:ppmX-squared=,df=6,p-value<Box-Piercetestdata:ppmX-squared=,df=12,p-value<（2）导入课后作业（41页第3题）数据并绘制时序图，并进行安稳性查验和纯随机性查验。z<-scan('3(2).txt')###在工作目录下调用txt文件3（2）###water<-ts(z,start=c(1945,1),frequency=12)####转变时序数据###plot(water)####绘制时序图######acf(water)###安稳性查验####for(iin1:2)print(water,lag=6*i))####纯随机性查验###时序图安稳性查验Serieswater010082006501rF4C00wA0210005020-019451946194719481949195019510.00.51.01.5TimeLag纯随机性查验结果Box-Piercetestdata:waterX-squared=,df=6,p-value=Box-Piercetestdata:waterX-squared=,df=12,p-value=（3）导入课后作业（41页第5题）数据并绘制时序图，并进行安稳性查验和纯随机性查验。d<-scan('3(3).txt')###在工作目录下调用txt文件3（3）###sales<-ts(d,start=c(2000,1),frequency=12)####转变时序数据###plot(sales)####绘制时序图######acf(sales)###安稳性查验####for(iin1:2)print(sales,lag=6*i))####纯随机性查验###时序图安稳性查验Seriessales00031800526.0004s2Fe.CasA025.100000120-200020012002200320040.00.20.40.60.81.01.2TimeLag纯随机性查验结果Box-Piercetestdata:salesX-squared=,df=6,p-value<Box-Piercetestdata:salesX-squared=,df=12,p-value<（4）导入课后作业（41页第6题）数据并绘制时序图，进行安稳性查验和纯随机性查验，并对序列进行差分，再做安稳性查验与纯随机性查验。b<-scan('3(4).txt')

###在工作目录下调用

txt

文件

3（4）###qb<-ts(b,start=c(2000,1),frequency=15)

####转变时序数据

###plot(qb)

####绘制时序图

######acf(qb)###安稳性查验

####for(iin1:2)print(qb,lag=6*i))####纯随机性查验###时序图安稳性查验Seriesqb05.13800356.20b0F4C.q20A5210000520-196919701971197219730.00.20.40.60.81.01.2TimeLag纯随机性查验结果Box-Piercetestdata:qbX-squared=,df=6,p-value=Box-Piercetestdata:qbX-squared=,df=12,p-value=四、实验心得经过本次实验，我学到了用R做时序图，以实时序图安稳性查验，纯随机性查验，了解了更多的R语言使用的知识，为此后的《时间序列剖析》