クロ ータ分析演習
6.Time Series Data Shige Okajima
Excel の ータをStataのフォー ッ に
• TimeDatePractice.xlsxをイン ー
• 例え Excel フ ー ッ 日次 以下 う
• mdy(m,d,y) 関数利用
• gen date=mdy(month, day, year)
• 注意month,day, year 変数 名前
例題 (日 ータ
• 日次 mdy(m,d,y) 関数利用 ータMIKwaterを用い
gen date=mdy(month, day, year)をdoし 以下 表を く さい
例題 (日 ータ
format date %tdを用い
例題 (日 ータ
時系列 日時を横軸 フを書く
graph twoway line water date
2004006008001000
01jan1983 01feb198301mar198301apr198301may198301jun1983 01jul1983 01aug1983 date
例題 4 半期 ータ
• 四半期 yq(y,q) 関数利用
http://wps.aw.com/aw_stock_ie_3/178/45691/11696965.cw/ USMacro_Quarterly Data (Excel Dataset)を利用
gen date=yq(Year,Quarter)をdoし 以下 表を く さい
例題 4 半期 ータ
format date %tqを用い
例題 4 半期 ータ
• 時系列 四半期を横軸 フを書く
graph twoway line Oil date
0
50001000015000
1950q1 1960q1 1970q1 1980q1 1990q1 2000q1 2010q1 date
まとめ と
• mdy(month,day year) 日次 ータ format %td
• yw(year,week) 週 ータ format %td %tw
• ym(year,month) 月次 ータ format %td %tm
• yq(year,quarter) 四半期 ータ format %td %tq
• yh(year,half-year) 半年 ータ format %td %th
• 年次 ータ そ format %ty
問題
http://wps.aw.com/aw_stock_ie_3/178/45691/11696965.cw/ USMacro_Monthly Data (Excel Dataset)を利用し 例題 手 順 図を書い ください
Excel の ータをStataのフォー ッ に
• TimeDatePractice2.xlsxをイン ー
• 例え Excel フ ー ッ 日次 (文字列 以下 う
• codebook studyday 確認 文字列 う を
• date(変数名, ) 関数利用
• gen daily=date(studyday,"MDY")
• 注 studyday 変数 名前 D 日M 月Y 年を表し い
文字列を
• format date %tdを用い
文字列を
graph twoway line time daily
文字列を
• weekly(変数名１ っwy、) 週 ータを
• montly (変数名１ っmy、) 月 ータを
• quartely(変数名１ っqy、) 4半期 ータを
• halfyearly(変数名１ っhy、) 半年 ータを
• yearly(変数名１ っy、) 年 ータを
文字列を
• TimeDatePractice2xlsxをイン ー
• 例え Excel フ ー ッ 日次 (文字列 以下 う
• codebook studyday 確認 文字列 う を
文字列を
• 数字 そ 大丈夫
• graph twoway line grade classday
0510
20-May-99 23-May-99 26-May-99 29-May-99 1-Jun-99 classday
時系列 ータを描く
ータ 時系列だ Stata 認識
ータArcticを用い
tsset year
時間設定を ア
tsset, clear
再度tsset year
時系列 フを書く
tsline area 3
456
area
1980 1990 2000 2010
year
時系列 ータを描く 二本の線
• tsline area extent, ytitle("Million km{superscript:2}") /// title("Arctic sea ice, 1979-2011")
345678
2
1980 1990 2000 2010
year
area extent
Arctic sea ice, 1979-2011
Km 2乗を
時系列 ータを描く 二本の線 つ右と
左のｙ軸 違う値
tsline area, ytitle("Ice area") ///
|| tsline tempN, yaxis(2) ytitle("Arctic temperature", axis(2))
-.50.511.52
Arctic temperature
3456
1980 1990 2000 2010
year
area tempN
時系列分析
• 時系列分析 場合被説明変数 を説明変数 入 こ あ
くあ そ 場合被説明変数 作 方
時間をStata 認識さ tsset year
時点過去を
generate area_Lag_1 =L1.area
時系列の変数の作 方
* 時点過去を *
generate area_Lag_1 =L1.area
* 時点過去を *
generate area_Lag_2 = L2.area
* 時点先を *
generate area_Lead_1 = F1.area
* 時点先を *
generate area_Lead_2 = F2.area
*一次 差分ｘ ｔ －ｘ ｔ－ *
generate area_Dif_1 = D1.area
*二次 差分ｘ ｔ －ｘ ｔ－ － ｘ ｔ－ －ｘ ｔ－ *
generate area_Dif_2 = D2.area
問題
• 問題 フをtslineを使 図を書い ください
• イ 一期前 変数を作成し ください
総合
• MIKwaterをイン ー し ください
• date変数を作成し 下さい
• フ ー ッ を変え ください
• そ 変数をtime series 読 込 ください
• そし water 時系列 図を書く
総合
300400500600700800900
01jan1983 01mar1983 01may1983 01jul1983
tsline water, ylabel(300(100)900) ttitle("") ///
tlabel(01jan1983 01mar1983 01may1983 01jul1983, grid) /// ttick(01feb1983 01apr1983 01jun1983 01aug1983)
ylabel(300(100)900)
300 始め 900
100 上昇
横軸 タイ し
ここ
総合
• Water 過去三ヶ月 平均 変数water3aを作 そ 変数 water 変数両方を図 描く
• generate water_Lag_1 =L1.water
• generate water_Lag_2 =L2.water
• generate water_Lag_3 =L3.water
• gen water3a=(water_Lag_1+water_Lag_2+water_Lag_3)/3
総合
tsline water,clwidth(medium) ///
|| tsline water3a, clwidth(medthick) ///
||, ylabel(300(100)900) ytitle(っWaがer uかe in 1000 gallonか、) ttitle("") ///
tlabel(01jan1983 01mar1983 01may1983 01jul1983, grid) /// ttick(01feb1983 01apr1983 01jun1983 01aug1983) ///
legend(position(4) ring(0) rows(2) label(1 "daily water use") label(2 "3 day average"))
総合
• legend(position(4) ring(0) rows(2) label(1 "daily water use") label(2 "3 day average"))
legend ここ
Position 時計方向 ４ ここだ 4時
Ring 枠 内 外 だ 内 Rows 何行 ？ だ 2行 label 名前を け
300400500600700800900
01jan1983 01mar1983 01may1983 01jul1983
daily water use 3 day average
