《时间序列实验报告.doc》由会员分享,可在线阅读,更多相关《时间序列实验报告.doc(15页珍藏版)》请在装配图网上搜索。
1、重 庆 交 通 大 学学 生 实 验 报 告实验课程名称 时间序列分析 开课实验中心 数统学院实验教学中心 开课学院 数学与统计学院 专业年级 应用统计学2015级1班 姓 名 XXXX 学 号 6315XXXXXXXX 任课老师 XXXXX 开课时间 20172018学年第1学期 1此页空页实验一 R语言简介 基本操作一 实验目的1、了解软件R安装、启动、退出、帮助等。2、熟悉R的操作界面。二、实验内容及要求1、 实验内容(1)R的安装;(2)启动与退出;(3)包的安装及R的更新;(4)帮助及移除多个对象等;(5)常见命令2、实验要求(1)熟悉R的操作环境;(2)熟悉包的安装与帮助;(3)学。
2、习常见命令,熟悉 R 的操作界面。三、实验过程及结果1、(1)R的安装(2)启动与退出;(3)包的安装及R的更新;A、包的安装 chooseCRANmirror install.packages B、R的更新 install.packagesinstallr libraryinstallr updateR (4)帮助及移除多个对象等;关键字关键字 help.start帮助rmrmlistls移除多个对象 (5)常见命令四、实验心得 了解了R的一些基本使用及其常见的命令,为自己深入学习r的使用打下了基础。实验二 R语言简介 数据集创建与处理一 实验目的1、掌握R数据集的不同创建形式。2、熟悉并掌。
3、握利用R对时间序列数据集进行变换与处理。二、实验内容及要求1、实验内容(1)利用data.frame函数创建数据集;(2)读取 d.txt 型数据框;(3)读取 excel 数据及对某变量数据进行某些处理(4)导出 R 中数据集(5)时间序列数据输入(6)对已有数据集中数据的处理2、实验要求熟悉R数据集的不同创建方法,掌握利用R对时序数据集进行变换与处理3、 实验过程及结果1、实验内容(1)利用data.frame函数创建数据集; mydata-data.framecol1,col2,col3,.col1,col2,col3,.可以是任何类型的列向量(2)读取 d.txt 型数据框;m w-r。
4、ead.tabletext.data.txt(3)读取 excel 数据及对某变量数据进行某些处理 libraryreadxl dataset - read_excelNULL Viewdataset(4)导出 R 中数据集write.table x, file , sep , row.names TRUE, col.names TRUE, quote TRUEx需要导出的数据file导出的文件路径sep分隔符,默认为空格( ),也就是以空格为分割列row.names是否导出行序号,默认为TRUE,也就是导出行序号col.names是否导出列名,默认为TRUE,也就是导出列名quote字符串是。
5、否使用引号表示,默认为TRUE,也就是使用引号表示(5)时间序列数据输入(6)对已有数据集中数据的处理(5)(6)合 libraryreadxl X2_7- read_excelC/Users/Administrator/Desktop/2.7.xlsx summaryX2_7 330.45 330.97 331.64 332.87 333.61 Min. 331.6 Min. 330.1 Min. 328.6 Min. 328.3 Min. 329.4 1st Qu.332.9 1st Qu.332.4 1st Qu.331.9 1st Qu.331.5 1st Qu.332.8 Media。
6、n 334.7 Median 334.4 Median 333.7 Median 334.4 Median 335.1 Mean 335.0 Mean 334.2 Mean 333.9 Mean 334.3 Mean 335.2 3rd Qu.336.8 3rd Qu.336.1 3rd Qu.335.9 3rd Qu.337.0 3rd Qu.337.7 Max. 339.2 Max. 338.2 Max. 339.9 Max. 340.6 Max. 341.2 333.55 Min. 330.6 1st Qu.333.9 Median 336.0 Mean 335.7 3rd Qu.338。
7、.0 Max. 340.9 四、实验心得 通过本次实验,首先,我知道了文件其他格式的文件如何导入R,知晓乐数据集的创建,使用及一些简单的处理。实验三 平稳时间序列模型预处理一 实验目的1、掌握平稳时间序列的特点2、学会通过序列图形对时间序列进行平稳性检验和纯随机性检验预处理二、实验内容及要求1、实验内容(1)导入课后作业(41页第2题)数据并绘制时序图,并进行平稳性检验和纯随机性检验。(2)导入课后作业(41页第3题)数据并绘制时序图,并进行平稳性检验和纯随机性检验。(3)导入课后作业(41页第5题)数据并绘制时序图,并进行平稳性检验和纯随机性检验。(4)导入课后作业(41页第6题)数据并绘制。
8、时序图,进行平稳性检验和纯随机性检验,并对序列进行差分,再做平稳性检验与纯随机性检验。2、实验要求熟悉通过序列图形对时间序列进行平稳性检验和纯随机性检验预处理 三、实验过程及结果 1、实验内容(1)导入课后作业(41页第1题)数据并绘制时序图,并进行平稳性检验和纯随机性检验。a-scan31.txt在工作目录下调用txt文件3(1)ppm-tsa,startc1975,1,frequency12转化时序数据plotppm绘制时序图acfppm平稳性检验fori in 12 printBox.testppm,lag6*i纯随机性检验时序图 平稳性检验纯随机性检验结果 Box-Pierce tes。
9、tdata ppmX-squared 131.88, df 6, p-value 2.2e-16 Box-Pierce testdata ppmX-squared 217.26, df 12, p-value 2.2e-16(2)导入课后作业(41页第3题)数据并绘制时序图,并进行平稳性检验和纯随机性检验。z-scan32.txt在工作目录下调用txt文件3(2)water-tsz,startc1945,1,frequency12转化时序数据plotwater绘制时序图acfwater平稳性检验fori in 12 printBox.testwater,lag6*i纯随机性检验时序。
10、图 平稳性检验纯随机性检验结果 Box-Pierce testdata waterX-squared 7.7505, df 6, p-value 0.257 Box-Pierce testdata waterX-squared 19.914, df 12, p-value 0.06873(3)导入课后作业(41页第5题)数据并绘制时序图,并进行平稳性检验和纯随机性检验。d-scan33.txt在工作目录下调用txt文件3(3)sales-tsd,startc2000,1,frequency12转化时序数据plotsales绘制时序图acfsales平稳性检验fori in 12 p。
11、rintBox.testsales,lag6*i纯随机性检验时序图 平稳性检验纯随机性检验结果 Box-Pierce testdata salesX-squared 134.49, df 6, p-value 2.2e-16 Box-Pierce testdata salesX-squared 163.85, df 12, p-value 2.2e-16(4)导入课后作业(41页第6题)数据并绘制时序图,进行平稳性检验和纯随机性检验,并对序列进行差分,再做平稳性检验与纯随机性检验。 b-scan34.txt在工作目录下调用txt文件3(4)qb-tsb,startc2000,1,f。
12、requency15转化时序数据plotqb绘制时序图acfqb平稳性检验fori in 12 printBox.testqb,lag6*i纯随机性检验时序图 平稳性检验纯随机性检验结果 Box-Pierce testdata qbX-squared 60.084, df 6, p-value 4.327e-11 Box-Pierce testdata qbX-squared 81.066, df 12, p-value 2.583e-12四、实验心得通过本次实验,我学到了用R做时序图,以及时序图平稳性检验,纯随机性检验,知晓了更多的R语言使用的知识,为以后的时间序列分析的学习打下了一定的r语言实践基础。实验四 平稳线性ARMA模型初步识别与定阶一 实验目的1、 掌握AR模型平稳性、MA模型的可逆性及ARMA模型平稳可逆性的判别条件;2、 掌握利用序列的自相关图和偏自相关图识别模型并进行初步定阶。二、实验内容及要求对上次实验数据进行模型初步定阶三、实验过程及结果1. 图从上图可看出,该序列样本自相关系数五、实验心得通过本次实验,学会。15。