時(shí)間序列實(shí)驗(yàn)報(bào)告

1、重 慶 交 通 大 學(xué)學(xué) 生 實(shí) 驗(yàn) 報(bào) 告實(shí)驗(yàn)課程名稱(chēng) 時(shí)間序列分析 開(kāi)課實(shí)驗(yàn)中心 數(shù)統(tǒng)學(xué)院實(shí)驗(yàn)教學(xué)中心 開(kāi)課學(xué)院 數(shù)學(xué)與統(tǒng)計(jì)學(xué)院 專(zhuān)業(yè)年級(jí) 應(yīng)用統(tǒng)計(jì)學(xué)2015級(jí)1班 姓 名 XXXX 學(xué) 號(hào) 6315XXXXXXXX 任課老師 XXXXX 開(kāi)課時(shí)間 20172018學(xué)年第1學(xué)期 歡迎下載此頁(yè)空頁(yè)！實(shí)驗(yàn)一 R語(yǔ)言簡(jiǎn)介: 基本操作一 實(shí)驗(yàn)?zāi)康?、了解軟件R：安裝、啟動(dòng)、退出、幫助等。2、熟悉R的操作界面。二、實(shí)驗(yàn)內(nèi)容及要求：1、 實(shí)驗(yàn)內(nèi)容：（1）R的安裝；（2）啟動(dòng)與退出；（3）包的安裝及R的更新；（4）幫助及移除多個(gè)對(duì)象等；（5）常見(jiàn)命令2、實(shí)驗(yàn)要求：（1）熟悉R的操作環(huán)境；（2）熟悉包的安

2、裝與幫助；（3）學(xué)習(xí)常見(jiàn)命令，熟悉 R 的操作界面。三、實(shí)驗(yàn)過(guò)程及結(jié)果1、（1）R的安裝（2）啟動(dòng)與退出；（3）包的安裝及R的更新；A、包的安裝> chooseCRANmirror()> install.packages() B、R的更新> install.packages("installr")> library(installr)> updateR() （4）幫助及移除多個(gè)對(duì)象等；> ？關(guān)鍵字> ？關(guān)鍵字> help.start()#幫助> rm()> rm(list=l

3、s()#移除多個(gè)對(duì)象 （5）常見(jiàn)命令四、實(shí)驗(yàn)心得 了解了R的一些基本使用及其常見(jiàn)的命令，為自己深入學(xué)習(xí)r的使用打下了基礎(chǔ)。實(shí)驗(yàn)二 R語(yǔ)言簡(jiǎn)介: 數(shù)據(jù)集創(chuàng)建與處理一 實(shí)驗(yàn)?zāi)康?、掌握R數(shù)據(jù)集的不同創(chuàng)建形式。2、熟悉并掌握利用R對(duì)時(shí)間序列數(shù)據(jù)集進(jìn)行變換與處理。二、實(shí)驗(yàn)內(nèi)容及要求1、實(shí)驗(yàn)內(nèi)容：（1）利用data.frame函數(shù)創(chuàng)建數(shù)據(jù)集；（2）讀取 d.txt 型數(shù)據(jù)框；（3）讀取 excel 數(shù)據(jù)及對(duì)某變量數(shù)據(jù)進(jìn)行某些處理（4）導(dǎo)出 R 中數(shù)據(jù)集（5）時(shí)間序列數(shù)據(jù)輸入（6）對(duì)已有數(shù)據(jù)集中數(shù)據(jù)的處理2、實(shí)驗(yàn)要求：熟悉R數(shù)據(jù)集的不同創(chuàng)建方法，掌握利用R對(duì)時(shí)序數(shù)據(jù)集進(jìn)行變換與處理3、 實(shí)驗(yàn)過(guò)程及結(jié)果1、

4、實(shí)驗(yàn)內(nèi)容：（1）利用data.frame函數(shù)創(chuàng)建數(shù)據(jù)集；> mydata<-data.frame(col1,col2,col3,.)#col1，col2，col3，.可以是任何類(lèi)型的列向量（2）讀取 d.txt 型數(shù)據(jù)框；m> w<-read.table("text.data.txt")（3）讀取 excel 數(shù)據(jù)及對(duì)某變量數(shù)據(jù)進(jìn)行某些處理> library(readxl)> dataset <- read_excel(NULL)> View(dataset)（4）導(dǎo)出 R 中數(shù)據(jù)集>write.table (x, fi

5、le ="", sep ="", s =TRUE, s =TRUE, quote =TRUE)#x：需要導(dǎo)出的數(shù)據(jù)#file：導(dǎo)出的文件路徑#sep：分隔符，默認(rèn)為空格（" "），也就是以空格為分割列#s：是否導(dǎo)出行序號(hào)，默認(rèn)為T(mén)RUE，也就是導(dǎo)出行序號(hào)#s：是否導(dǎo)出列名，默認(rèn)為T(mén)RUE，也就是導(dǎo)出列名#quote：字符串是否使用引號(hào)表示，默認(rèn)為T(mén)RUE，也就是使用引號(hào)表示（5）時(shí)間序列數(shù)據(jù)輸入（6）對(duì)已有數(shù)據(jù)集中數(shù)據(jù)的處理（5）（6）合> library(read

6、xl)> X2_7<- read_excel("C:/Users/Administrator/Desktop/2.7.xlsx")> summary(X2_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 Median :334.7 Median :334.4 Median

7、: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.0 Ma

8、x. :340.9 四、實(shí)驗(yàn)心得 通過(guò)本次實(shí)驗(yàn)，首先，我知道了文件其他格式的文件如何導(dǎo)入R，知曉樂(lè)數(shù)據(jù)集的創(chuàng)建，使用及一些簡(jiǎn)單的處理。實(shí)驗(yàn)三 平穩(wěn)時(shí)間序列模型預(yù)處理一 實(shí)驗(yàn)?zāi)康?、掌握平穩(wěn)時(shí)間序列的特點(diǎn)2、學(xué)會(huì)通過(guò)序列圖形對(duì)時(shí)間序列進(jìn)行平穩(wěn)性檢驗(yàn)和純隨機(jī)性檢驗(yàn)預(yù)處理二、實(shí)驗(yàn)內(nèi)容及要求：1、實(shí)驗(yàn)內(nèi)容：（1）導(dǎo)入課后作業(yè)（41頁(yè)第2題）數(shù)據(jù)并繪制時(shí)序圖，并進(jìn)行平穩(wěn)性檢驗(yàn)和純隨機(jī)性檢驗(yàn)。（2）導(dǎo)入課后作業(yè)（41頁(yè)第3題）數(shù)據(jù)并繪制時(shí)序圖，并進(jìn)行平穩(wěn)性檢驗(yàn)和純隨機(jī)性檢驗(yàn)。（3）導(dǎo)入課后作業(yè)（41頁(yè)第5題）數(shù)據(jù)并繪制時(shí)序圖，并進(jìn)行平穩(wěn)性檢驗(yàn)和純隨機(jī)性檢驗(yàn)。（4）導(dǎo)入課后作業(yè)（41頁(yè)第6題）數(shù)據(jù)并繪制時(shí)序

9、圖，進(jìn)行平穩(wěn)性檢驗(yàn)和純隨機(jī)性檢驗(yàn)，并對(duì)序列進(jìn)行差分，再做平穩(wěn)性檢驗(yàn)與純隨機(jī)性檢驗(yàn)。2、實(shí)驗(yàn)要求：熟悉通過(guò)序列圖形對(duì)時(shí)間序列進(jìn)行平穩(wěn)性檢驗(yàn)和純隨機(jī)性檢驗(yàn)預(yù)處理 三、實(shí)驗(yàn)過(guò)程及結(jié)果 1、實(shí)驗(yàn)內(nèi)容：（1）導(dǎo)入課后作業(yè)（41頁(yè)第1題）數(shù)據(jù)并繪制時(shí)序圖，并進(jìn)行平穩(wěn)性檢驗(yàn)和純隨機(jī)性檢驗(yàn)。a<-scan('3(1).txt')#在工作目錄下調(diào)用txt文件3（1）#ppm<-ts(a,start=c(1975,1),frequency=12)#轉(zhuǎn)化時(shí)序數(shù)據(jù)#plot(ppm)#繪制時(shí)序圖#acf(ppm)#平穩(wěn)性檢驗(yàn)#for(i in 1:2) print(Box.test(ppm

10、,lag=6*i)#純隨機(jī)性檢驗(yàn)#時(shí)序圖 平穩(wěn)性檢驗(yàn)純隨機(jī)性檢驗(yàn)結(jié)果 Box-Pierce testdata: 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）導(dǎo)入課后作業(yè)（41頁(yè)第3題）數(shù)據(jù)并繪制時(shí)序圖，并進(jìn)行平穩(wěn)性檢驗(yàn)和純隨機(jī)性檢驗(yàn)。z<-scan('3(2).txt')#在工作目錄下調(diào)用txt文件3（2）#water<-ts(z,start=c(1945,

11、1),frequency=12)#轉(zhuǎn)化時(shí)序數(shù)據(jù)#plot(water)#繪制時(shí)序圖#acf(water)#平穩(wěn)性檢驗(yàn)#for(i in 1:2) print(Box.test(water,lag=6*i)#純隨機(jī)性檢驗(yàn)#時(shí)序圖 平穩(wěn)性檢驗(yàn)純隨機(jī)性檢驗(yàn)結(jié)果 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）導(dǎo)入課后作業(yè)（41頁(yè)第5題）數(shù)據(jù)并繪制時(shí)序圖，并進(jìn)行平

12、穩(wěn)性檢驗(yàn)和純隨機(jī)性檢驗(yàn)。d<-scan('3(3).txt')#在工作目錄下調(diào)用txt文件3（3）#sales<-ts(d,start=c(2000,1),frequency=12)#轉(zhuǎn)化時(shí)序數(shù)據(jù)#plot(sales)#繪制時(shí)序圖#acf(sales)#平穩(wěn)性檢驗(yàn)#for(i in 1:2) print(Box.test(sales,lag=6*i)#純隨機(jī)性檢驗(yàn)#時(shí)序圖 平穩(wěn)性檢驗(yàn)純隨機(jī)性檢驗(yàn)結(jié)果 Box-Pierce testdata: salesX-squared = 134.49, df = 6, p-value < 2.2e-16 Box-Pier

13、ce testdata: salesX-squared = 163.85, df = 12, p-value < 2.2e-16（4）導(dǎo)入課后作業(yè)（41頁(yè)第6題）數(shù)據(jù)并繪制時(shí)序圖，進(jìn)行平穩(wěn)性檢驗(yàn)和純隨機(jī)性檢驗(yàn)，并對(duì)序列進(jìn)行差分，再做平穩(wěn)性檢驗(yàn)與純隨機(jī)性檢驗(yàn)。 b<-scan('3(4).txt')#在工作目錄下調(diào)用txt文件3（4）#qb<-ts(b,start=c(2000,1),frequency=15)#轉(zhuǎn)化時(shí)序數(shù)據(jù)#plot(qb)#繪制時(shí)序圖#acf(qb)#平穩(wěn)性檢驗(yàn)#for(i in 1:2) print(Box.test(qb,lag=6*i)#純隨機(jī)性檢驗(yàn)#時(shí)序圖 平穩(wěn)性檢驗(yàn)純隨機(jī)性檢驗(yàn)結(jié)果 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四、實(shí)驗(yàn)心得通過(guò)本次實(shí)驗(yàn)，我學(xué)到了用R做時(shí)序圖，以及時(shí)序圖平穩(wěn)性檢驗(yàn)，純隨機(jī)性檢驗(yàn)