http://tkustat.stat.tku.edu.tw
機率  |  統計  |  Data Mining  |  資料處理  |  Demo  |  Sample Data  |  Links  |  相關文件  ( ENGLISH )

時間數列資料初步檢視
Data: lynx (範例)
   

===========
輸入資料檔
===========

R 指令> library(tkustat)          
R 指令> data(lynx)          
R 指令> mydata <- lynx 

======================
未差分的時間數列 Lynx
======================
R 指令> z0.Lynx <- mydata$Lynx

==================================
先對原始時間數列 Lynx 作 Log 轉換
==================================
R 指令> n <- length(z0.Lynx)
R 指令> z0.Lynx <- log(abs(z0.Lynx)+1.0/(2*n))

==============================
 引用 R-package 的 ts library
==============================
R 指令> if(!exists("ts.plot", mode="function"))

R 指令> library(ts)

===========================
 Log Lynx 變數的敘述統計量
===========================
R 指令> z0.desc <- tkustat.desc(z0.Lynx)
R 指令> z0.desc

    樣本數     最小值         Q1     中位數     平均數         Q3     最大值 
 64.000000   3.664000   5.871000   6.615000   6.647000   7.674000   8.852000 
    變異數       總和     標準差        IQR        MAD 
  1.594680 425.422651   1.262806   1.802503   1.355358 

===========================
 Log Lynx 變數的時間數列圖
===========================

R 指令> ts.plot(z0.Lynx)
R 指令> abline(h=mean(z0.Lynx),col="red")
=======================================
 Log Lynx 變數的樣本 ACF 圖 及 PACF 圖
=======================================
R 指令> acf(z0.Lynx)
R 指令> pacf(z0.Lynx)
===========================
 Log Lynx 變數的 樣本 EACF
===========================
R 指令> eacf(z0.Lynx,ar.max=10,ma.max=10)

AR/MA
   0   1   2   3   4   5   6   7   8   9   10 
0  "X" "X" "O" "X" "X" "X" "O" "O" "X" "X" "X"
1  "X" "X" "O" "X" "X" "X" "O" "O" "X" "X" "X"
2  "O" "O" "O" "O" "O" "O" "O" "O" "O" "O" "O"
3  "X" "O" "O" "O" "O" "O" "O" "O" "O" "O" "O"
4  "X" "X" "O" "O" "O" "O" "O" "O" "O" "O" "O"
5  "X" "O" "O" "O" "O" "O" "O" "O" "O" "O" "O"
6  "X" "O" "O" "O" "O" "O" "O" "O" "O" "O" "O"
7  "X" "O" "O" "O" "O" "O" "O" "O" "O" "O" "O"
8  "X" "O" "O" "O" "O" "O" "O" "O" "O" "O" "O"
9  "X" "X" "X" "X" "O" "O" "O" "O" "O" "O" "O"
10 "X" "O" "X" "O" "O" "O" "O" "O" "O" "O" "O"

==============================
 對 Log Lynx 變數 作 1 次差分
==============================

R 指令> z1.Lynx <- diff(z0.Lynx, differences=1)

=======================================
 Log Lynx 變數 (差分 1 次): 敘述統計量
=======================================
R 指令> z1.desc <- tkustat.desc(z1.Lynx)
R 指令> z1.desc

    樣本數     最小值         Q1     中位數     平均數         Q3     最大值 
63.0000000 -1.8650000 -0.5420000  0.1304000 -0.0095670  0.7049000  1.2320000 
    變異數       總和     標準差        IQR        MAD 
 0.7124588 -0.6027216  0.8440728  1.2469845  0.8694722 

=======================================
 Log Lynx 變數 (差分 1 次): 時間數列圖
=======================================

R 指令> ts.plot(z1.Lynx)
R 指令> abline(h=mean(z1.Lynx),col="red")
===================================================
 Log Lynx 變數 (差分 1 次): 樣本 ACF 圖 及 PACF 圖
===================================================
R 指令> acf(z1.Lynx)
R 指令> pacf(z1.Lynx)
======================================
 Log Lynx 變數 (差分 1 次): 樣本 EACF
======================================
R 指令> eacf(z1.Lynx,ar.max=10,ma.max=10)

AR/MA
   0   1   2   3   4   5   6   7   8   9   10 
0  "X" "O" "X" "X" "X" "X" "O" "O" "X" "X" "X"
1  "X" "O" "O" "O" "X" "X" "O" "O" "X" "X" "X"
2  "X" "O" "O" "O" "O" "O" "O" "O" "O" "O" "O"
3  "X" "X" "O" "O" "O" "O" "O" "O" "O" "O" "O"
4  "X" "O" "O" "O" "O" "O" "O" "O" "O" "O" "O"
5  "X" "O" "O" "O" "O" "O" "O" "O" "O" "O" "O"
6  "X" "O" "O" "O" "O" "O" "O" "O" "O" "O" "O"
7  "X" "O" "O" "O" "O" "O" "O" "O" "O" "O" "O"
8  "X" "X" "X" "O" "O" "O" "O" "O" "O" "O" "O"
9  "X" "O" "X" "O" "O" "O" "O" "O" "O" "O" "O"
10 "X" "O" "X" "O" "O" "O" "O" "O" "O" "O" "O"