国土技術政策総合研究所資料

Size: px
Start display at page:

Download "国土技術政策総合研究所資料"

Transcription

1 ISSN 国総研資料第 652 号平成 23 年 9 月 国土技術政策総合研究所資料 TECHNICAL NOTE of Naional Insiue for Land and Infrasrucure Managemen No.652 Sepember 2011 航空需要予測における計量時系列分析手法の適用性に関する基礎的研究 ~ 季節変動自己回帰移動平均モデル及びベクトル誤差修正モデルの適用性 ~ 井上岳 丹生清輝 Sudy on Applicaion of Time-series Analysis o Air Transpor Demand Esimaion Gaku INOUE, Kiyoeru TANSEI 国土交通省国土技術政策総合研究所 Naional Insiue for Land and Infrasrucure Managemen Minisry of Land, Infrasrucure, Transpor and Tourism, Japan

2 No (YSK-N-238) * ** (ARIMA) 5-10 (VECM) GDP GDP GDP GDP : (ARIMA) (VECM) * ** : Fax: inoue-g23i@ysk.nilim.go.jp i

3 Technical Noe of NILIM No.652 Sepember 2011 (YSK-N-238) Sudy on Applicaion of Time-series Analysis o Air Transpor Demand Esimaion Gaku INOUE* Kiyoeru TANSEI** Synopsis Time-series analysis echniques have been applied o he various fields, including economic and financial analyses. This sudy examined he general applicabiliy of Auoregressive Inegraed Moving Average model(arima), one of ime-series analysis echniques, o air ranspor demand esimaion. Furhermore, his sudy also analyzed long-erm equilibrium beween air ranspor demand and GDP by using Vecor Error Correcion Model (VECM). These analysis comes o be basic daa for fuure research on air ranspor demand esimaion. Key Words: Air Transpor Demand Esimaion, Time-series Analysis, Auoregressive Inegraed Moving Average model(arima), Vecor Error Correcion Model (VECM) * Senior Researcher, Airpor Deparmen ** Direcor of Airpor Planning Division, Airpor Deparmen Nagase, Yokosuka Japan Phone: Fax: inoue-g23i@ysk.nilim.go.jp ii

4 (ARIMA) (VECM) GDP (VECM) A iii

5 iv

6 No ) 1 2) GDP GDP (ARIMA) (VECM) GDP 2 (ARIMA) 3 (VECM) GDP 4 2. (ARIMA) 2.1 (1) (ARIMA) (ARIMA) 2) (2010) 3) Hamilon(1994) 4) y ϕ(l)φ(l s ) D s d y = θ(l)θ(l s )ε (1) ϕ(l) = 1 ϕ 1L ϕ 2L 2 ϕ pl p Φ(L s ) = 1 Φ 1 L s Φ 2 L 2s Φ P L sp θ(l) = 1 + θ 1 L + θ 2 L θ q L q Θ(L s ) = 1 + Θ 1 L s + Θ 2 L 2s + + Θ Q L sq ϕ i, Φ i, θ i, Θ i (1) L : Ly = y 1 : y = y y 1. s : s y = y y s ε E(ε ) = 0 { σ 2 (k = 0) E(ε ε k ) = 0 (k 0) (1) ARIMA(p, d, q) (P, D, Q) s (2) p, d, q p : d : y q : P : D : Q : s : ARIMA - 1 -

7 / k E(y ) = µ Var(y ) = E[(y µ) 2 ] = γ 0 (3) Cov(y, y k ) = E[(y µ)(y k µ)] = γ k (4) (5) (1) ϵ ρ k = Corr(y, y k ) = Cov(y, y k ) Var(y )Var(y k ). (5) ρ k = Corr(y, y k ) = γ k γ 0 (6) k y y k k y y k y 1, y 2,, y k+1 5) k 0 95% ϵ ε (Pormaneau es) H 0 : ρ 1 = ρ 2 = = ρ m = 0 H 1 : 1 k [1, m] ρ k 0 Ljung and Box Q Q Ljung and Box Q(m) = T (T + 2) m k=1 ˆρ k 2 T k χ2 (m) (7) χ 2 (m) m T Q(m) χ 2 (m) 95% 5% ε p, d, q (AIC) (BIC) AIC AIC = 2 L(ˆθ) + 2n (8) L(ˆθ) n 2 BIC BIC = 2 L(ˆθ) + log(t )n (9) (2) (3) - 2 -

8 No.652 domesic d_domesic Auocorrelaions of d_domesic Lag Barle s formula for MA(q) 95% confidence bands Parial auocorrelaions of d_domesic Lag 95% Confidence bands [se = 1/sqr(n)] - 3 m m E(y ) = µ % 5% ρ k ARIMA(p, d, q) (P, D, Q) s d = D = 1 s = 12 0 p, q, P, Q 2 81 (p, q, P, Q) ϕ, Φ, θ, Θ 2.1 AIC ϕ, Φ, θ, Θ - 3 -

9 / y predicion, one sep domesic - 1 ARIMA(1, 1, 1) (0, 1, 1) 12 p AR(1) MA(1) MA(12) AIC Q p (0.4185) - 5 Auocorrelaions of e Lag Barle s formula for MA(q) 95% confidence bands Parial auocorrelaions of e Lag 95% Confidence bands [se = 1/sqr(n)] ARIMA(1, 1, 1) (0, 1, 1) 12-1 ϕ 1 = AR(1) θ 1 = MA(1) Θ 1 = MA(12) -5 ARIMA % -1 (Pormaneau es) Q χ 2 (m) 95% P 0.42 (1) ( ) (2) ( ) -2-3 (1) (2) ARIMA(0, 1, 1) (0, 1, 1) ARIMA(1, 1, 1) (0, 1, 1) Q - 4 -

10 No m1 2000m1 2005m1 2010m1 prediced(1994m12) domesic prediced(1999m12) m1 2000m1 2005m1 2010m1 observaion upper(95%) lower(68%) predicion upper(68%) lower(95%) % p MA(1) MA(12) AIC Q p (0.4583) p MA(1) MA(12) AIC Q p (0.2589) χ 2 (m) 95% P (1) (2) (3) % % 40% (MSE) MSE MSE (2010) 3) MSE h ŷ +1,, ŷ +h 1 y ˆε 2 1 y ˆε h ε ε 0 0 ε ˆε ε +h 0 ε (k) +1,, ε(k) +h N(0, σ2 ) y, y 1,, ε, ε 1, ε y (k) +h - 5 -

11 / N y (k), k = 1, 2,, N +h MSE y (k), k = 1, 2,, N +h N = 1, 000 MSE h 95% (±1.96σ) h 68% (±σ) % 68% 68% 95% 95% % 5 100% 68% E(y ) = µ % 5% ARIMA(p, d, q) (P, D, Q) s d = D = 1 s = 12 0 p, q, P, Q 2 81 (p, q, P, Q) ϕ, Φ, θ, Θ 2.1 AIC ϕ, Φ, θ, Θ ARIMA(1, 1, 2) (0, 1, 1) 12-4 ϕ 1 = AR(1) θ 1 = MA(1) θ 2 = MA(2) Θ 1 = MA(12) -15 ARIMA

12 No.652 inl d_inl d_ln_inl Auocorrelaions of d_inl Lag Barle s formula for MA(q) 95% confidence bands Parial auocorrelaions of d_inl Lag 95% Confidence bands [se = 1/sqr(n)] % (Pormaneau es) Q χ 2 (m) 95% P ARIMA(1, 1, 1) (2, 1, 2) 12 ARIMA(1, 1, 2) (0, 1, 1) 12 ϕ 1 = AR(1) θ 1 = MA(1) Φ 1 = AR(12) Φ 2 = AR(13) Θ 1 = MA(12) Θ 2 = MA(13) Φ 1 Φ 2 Θ 1 40 Q χ 2 (m) 95% - 14 P SARS 10% % 40-50% % 68% 68% 95% (MSE) 95% 68% 95% 95% 5 50% % 95% SARS - 7 -

13 / observaion y predicion, one sep ARIMA(1, 1, 2) (0, 1, 1) 12 p AR(1) MA(1) MA(2) MA(12) AIC Q p ( ) Auocorrelaions of e Lag Barle s formula for MA(q) 95% confidence bands Parial auocorrelaions of e Lag 95% Confidence bands [se = 1/sqr(n)] SARS SARS SARS SARS SARS SARS SARS SARS 3 (150) SARS 3 0.3% SARS : SARS :

14 No observaion predicion p AR(1) MA(1) AR(12) AR(13) MA(12) MA(13) AIC Q p ( ) m1 2002m1 2004m1 2006m1 2008m1 2010m1 observaion upper(95%) lower(68%) predicion upper(68%) lower(95%) residual, one sep m1 2002m1 2004m1 2006m1 2008m1 2010m % SARS SARS SARS SARS 2000m1 2002m1 2004m1 2006m1 2008m1 2010m1 observaion errorism predicion errorism&sars - 21 SARS - 9 -

15 / inlcargo d_inlcargo d_ln_inlcargo Auocorrelaions of d_ln_inlcargo Lag Barle s formula for MA(q) 95% confidence bands Parial auocorrelaions of d_ln_inlcargo Lag 95% Confidence bands [se = 1/sqr(n)] - 25 SARS E(y ) = µ % 5% ARIMA(p, d, q)

16 No observaion predicion - 6 ARIMA(0, 1, 1) (0, 1, 1) 12 p MA(1) MA(12) AIC Q p (0.9964) - 27 Auocorrelaions of e Lag Barle s formula for MA(q) 95% confidence bands Parial auocorrelaions of e Lag 95% Confidence bands [se = 1/sqr(n)] - 28 (P, D, Q) s d = D = 1 s = 12 0 p, q, P, Q 2 81 (p, q, P, Q) ϕ, Φ, θ, Θ 2.1 AIC ϕ, Φ, θ, Θ ARIMA(0, 1, 1) (0, 1, 1) 12-6 θ 1 = MA(1) Θ 1 = MA(12) -27 ARIMA % (Pormaneau es) Q χ 2 (m) 95% P ARIMA(0, 1, 1) (0, 1, 1) θ 1 = MA(1) Θ 1 = MA(12) 40 Q

17 / observaion predicion p MA(1) MA(12) AIC Q p ( ) χ 2 (m) 95% P (MSE) 95% 68% 68% 95% % % m1 2002m1 2004m1 2006m1 2008m1 2010m1 observaion upper(95%) lower(68%) predicion upper(68%) lower(95%) % (VECM) GDP 3.1 1) (1) (2) (3) (4) ( ) (5) GDP

18 No random1 random2 random3 random4 random5 α (152.67) (240.50) (91.99) (228.14) (226.37) β (15.79) (16.62) (8.74) (13.92) (18.87) R random1 random2 random3 random4 random5 α (86.64) (141.89) (59.13) (127.53) (137.00) β (15.11) (16.23) (9.86) (12.63) (19.15) R random1 random2 random3 random4 random5 α (89.00) (162.32) (61.99) (138.05) (176.93) β (16.78) (20.16) (11.47) (14.91) (26.95) R year random1 random2 random3 random4 random5-32 GDP (OLS) y y = y y 1 y x x y GDP (y ) 5 (x ) y = α + βx x = x ε, ε N(0, 2 2 ) (10) -A.1-32 R GDP GDP (VECM) 3.2 (VECM) (VECM) (2010) 3) Hamilon(1994) 4) GDP Chang 6)

19 / - 11 Phillips-Perron uni roo ess Variables Wih a ime rend DP IP IC Y DP ** IP ** IC ** Y ** ** : 1% - 12 Maximum likelihood coinegraion ess Saiics Saiics * 14.28* * 15.31* * 95% ** 99% 2 x y ψ y ψx 2 2 x y 2 2 x y y = α 1 + λ 11 y λ 1p y p + β 11 x β 1p x p + η 1 (y ψx ) + ε 1 x = α 1 + λ 21 y λ 2p y p + β 21 x β 2p x p + η 2(y ψx ) + ε 2 (11) (Vecor Error Correcion Model) y ψx (x, y ) 2 x x x, y x (MSE) y x (11) λ 2i = 0 i p VECM F ( (11) η i ) GDP 2000( 12) GDP 1979 GDP 68SNA 93SNA 3.4 (1) GDP Phillips and Perron PP -11 DP IP IC Y GDP 1 PP GDP (2) GDP H 0 : 5% 95% VECM GDP

20 No VECM ( - 14 VECM ( - 15 VECM ( ) DP Y ( 0.00) (-0.04) DP(-1) ( 2.42)* ( 1.55) Y(-1) (-0.67) ( 2.08)* ECT(-1) (-3.18)** (-1.84) R F DP-0.992Y(-6.07) ) IP Y ( 0.00) (-0.18) IP(-1) ( 0.46) ( 1.41) Y(-1) (-0.96) ( 2.87)** ECT(-1) (-3.78)** (-1.12) R F IP-1.591Y(-9.16) ) IC Y ( 0.01) (-0.68) IC(-1) ( 0.19) ( 2.31)* Y(-1) (-0.99) ( 3.61)** ECT(-1) (-3.74)** (-0.46) R F IC-2.052Y(-8.89) F T ( DP Y ECT 1 DP ** Y ** 1% F T ( IP Y ECT 1 IP ** Y ** 1% F T ( IC Y ECT 1 IC ** Y 5.33 * ** 1% *5% (AIC) (BIC) AIC AIC BIC 2 BIC 2 AIC BIC 2 2 (3) VECM VECM ECT (11) η i ( ) ( ) GDP GDP GDP GDP GDP Y OLS GDP 5% GDP GDP GDP DP Y GDP GDP GDP GDP GDP Y OLS GDP 5% GDP

21 / GDP GDP IP Y GDP GDP GDP GDP GDP Y OLS GDP 5% GDP GDP GDP GDP IP Y GDP GDP GDP GDP GDP GDP GDP GDP GDP VECM GDP GDP Y OLS GDP 4. ARIMA VECM GDP GDP GDP GDP GDP GDP GDP

22 No.652 ( ) 1) (2007): hp:// juyou1.hml (2011/8/31 ) 2) (2003): No ) (2010): 4) Hamilon,J.(1994), Time Series Analysis, Princeon Universiy Press. 5) A.C.(1985): 6) Chang,Y. and Chang,Y.(2009), Air cargo expansion and economic growh: Finding he empirical link, Journal of Air Transpor Managemen Vol.15, pp

23 / A A year random1 random2 random3 random4 random x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x = x x

24 国土技術政策総合研究所資料 TECHNICAL NOTE of N I L I M No. 652 Sepember 2011 編集 発行 国土技術政策総合研究所 本資料の転載 複写のお問い合わせは 神奈川県横須賀市長瀬 管理調整部企画調整課電話 :

国土技術政策総合研究所 研究資料

国土技術政策総合研究所 研究資料 ISSN TECHNICAL NOTE of National Institute for Land and Infrastructure Management No256 September 2005 Experimental Study on Seismic Behavior of Seawalls for Controlled Waste Disposal Shingo KANO, Katsuya

More information

²�ËÜËܤǻþ·ÏÎó²òÀÏÊÙ¶¯²ñ - Â裱¾Ï¤ÈÂ裲¾ÏÁ°È¾

²�ËÜËܤǻþ·ÏÎó²òÀÏÊÙ¶¯²ñ - Â裱¾Ï¤ÈÂ裲¾ÏÁ°È¾ Kano Lab. Yuchi MATSUOKA December 22, 2016 1 / 32 1 1.1 1.2 1.3 1.4 2 ARMA 2.1 ARMA 2 / 32 1 1.1 1.2 1.3 1.4 2 ARMA 2.1 ARMA 3 / 32 1.1.1 - - - 4 / 32 1.1.2 - - - - - 5 / 32 1.1.3 y t µ t = E(y t ), V

More information

国土技術政策総合研究所 研究資料

国土技術政策総合研究所 研究資料 ISSN TECHNICAL NOTE of National Institute for Land and Infrastructure Management No255 September 2005 A Result of Investigation into the Actual Conditions of Japanese Ports Pavement and Cargo Handling

More information

N cos s s cos ψ e e e e 3 3 e e 3 e 3 e

N cos s s cos ψ e e e e 3 3 e e 3 e 3 e 3 3 5 5 5 3 3 7 5 33 5 33 9 5 8 > e > f U f U u u > u ue u e u ue u ue u e u e u u e u u e u N cos s s cos ψ e e e e 3 3 e e 3 e 3 e 3 > A A > A E A f A A f A [ ] f A A e > > A e[ ] > f A E A < < f ; >

More information

Microsoft Word - 表紙.docx

Microsoft Word - 表紙.docx 黒住英司 [ 著 ] サピエンティア 計量経済学 訂正および練習問題解答 (206/2/2 版 ) 訂正 練習問題解答 3 .69, 3.8 4 (X i X)U i i i (X i μ x )U i ( X μx ) U i. i E [ ] (X i μ x )U i i E[(X i μ x )]E[U i ]0. i V [ ] (X i μ x )U i i 2 i j E [(X i

More information

橡表紙参照.PDF

橡表紙参照.PDF CIRJE-J-58 X-12-ARIMA 2000 : 2001 6 How to use X-12-ARIMA2000 when you must: A Case Study of Hojinkigyo-Tokei Naoto Kunitomo Faculty of Economics, The University of Tokyo Abstract: We illustrate how to

More information

医系の統計入門第 2 版 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 第 2 版 1 刷発行時のものです.

医系の統計入門第 2 版 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます.   このサンプルページの内容は, 第 2 版 1 刷発行時のものです. 医系の統計入門第 2 版 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/009192 このサンプルページの内容は, 第 2 版 1 刷発行時のものです. i 2 t 1. 2. 3 2 3. 6 4. 7 5. n 2 ν 6. 2 7. 2003 ii 2 2013 10 iii 1987

More information

McCain & McCleary (1979) The Statistical Analysis of the Simple Interrupted Time-Series Quasi-Experiment

McCain & McCleary (1979) The Statistical Analysis of the Simple Interrupted Time-Series Quasi-Experiment Quasi-Experimenaion Ch.6 005/8/7 ypo rep: The Saisical Analysis of he Simple Inerruped Time-Series Quasi-Experimen INTRODUCTION () THE PROBLEM WITH ORDINAR LEAST SQUARE REGRESSION OLS (Ordinary Leas Square)

More information

chap9.dvi

chap9.dvi 9 AR (i) (ii) MA (iii) (iv) (v) 9.1 2 1 AR 1 9.1.1 S S y j = (α i + β i j) D ij + η j, η j = ρ S η j S + ε j (j =1,,T) (1) i=1 {ε j } i.i.d(,σ 2 ) η j (j ) D ij j i S 1 S =1 D ij =1 S>1 S =4 (1) y j =

More information

時系列データ解析による予測と最適化 ~エネルギー需要、発電、価格のモデリング~

時系列データ解析による予測と最適化 ~エネルギー需要、発電、価格のモデリング~ MATLAB による時系列データ解析と予測 MahWorks Japan アプリケーションエンジニアリング部テクニカルコンピューティング 中川慶子 2015 The MahWorks, Inc. 1 アジェンダ 需要予測 : 時系列データモデリング 1. データの準備 データの取得 生データの前処理 2. 機械学習 非線形重回帰 ニューラルネットワーク RNN 3. 自己回帰系モデル ARIMA/GARCHモデル

More information

chap10.dvi

chap10.dvi . q {y j } I( ( L y j =Δy j = u j = C l ε j l = C(L ε j, {ε j } i.i.d.(,i q ( l= y O p ( {u j } q {C l } A l C l

More information

y = x x R = 0. 9, R = σ $ = y x w = x y x x w = x y α ε = + β + x x x y α ε = + β + γ x + x x x x' = / x y' = y/ x y' =

y = x x R = 0. 9, R = σ $ = y x w = x y x x w = x y α ε = + β + x x x y α ε = + β + γ x + x x x x' = / x y' = y/ x y' = y x = α + β + ε =,, ε V( ε) = E( ε ) = σ α $ $ β w ( 0) σ = w σ σ y α x ε = + β + w w w w ε / w ( w y x α β ) = α$ $ W = yw βwxw $β = W ( W) ( W)( W) w x x w x x y y = = x W y W x y x y xw = y W = w w

More information

fiúŁÄ”s‘ê‡ÌŁª”U…−…X…N…v…„…~…A…•‡Ì ”s‘ê™´›ß…−…^†[…fiŠ‚ª›Âfl’«

fiúŁÄ”s‘ê‡ÌŁª”U…−…X…N…v…„…~…A…•‡Ì ”s‘ê™´›ß…−…^†[…fiŠ‚ª›Âfl’« 2016/3/11 Realized Volatility RV 1 RV 1 Implied Volatility IV Volatility Risk Premium VRP 1 (Fama and French(1988) Campbell and Shiller(1988)) (Hodrick(1992)) (Lettau and Ludvigson (2001)) VRP (Bollerslev

More information

LA-VAR Toda- Yamamoto(1995) VAR (Lag Augmented vector autoregressive model LA-VAR ) 2 2 Nordhaus(1975) 3 1 (D2)

LA-VAR Toda- Yamamoto(1995) VAR (Lag Augmented vector autoregressive model LA-VAR ) 2 2 Nordhaus(1975) 3 1 (D2) LA-VAR 1 1 1973 4 2000 4 Toda- Yamamoto(1995) VAR (Lag Augmented vector autoregressive model LA-VAR ) 2 2 Nordhaus(1975) 3 1 (D2) E-mail b1215@yamaguchi-u.ac.jp 2 Toda, Hiro Y. and Yamamoto,T.(1995) 3

More information

研究シリーズ第40号

研究シリーズ第40号 165 PEN WPI CPI WAGE IIP Feige and Pearce 166 167 168 169 Vector Autoregression n (z) z z p p p zt = φ1zt 1 + φ2zt 2 + + φ pzt p + t Cov( 0 ε t, ε t j )= Σ for for j 0 j = 0 Cov( ε t, zt j ) = 0 j = >

More information

微分積分 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 初版 1 刷発行時のものです.

微分積分 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます.   このサンプルページの内容は, 初版 1 刷発行時のものです. 微分積分 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. ttp://www.morikita.co.jp/books/mid/00571 このサンプルページの内容は, 初版 1 刷発行時のものです. i ii 014 10 iii [note] 1 3 iv 4 5 3 6 4 x 0 sin x x 1 5 6 z = f(x, y) 1 y = f(x)

More information

金融政策の波及経路と政策手段

金融政策の波及経路と政策手段 Krugman(988) Woodford(999) (2000) (2000) 4 rae-of-reurn dominance 405 4 406 (i) a 2 cash good credi good b King and Wolman(999) (ii) 407 3 4 90 (iii) (iv) 408 λ κ (2.8) π x π λ = x κ Svensson 999 sric

More information

untitled

untitled 3 4 4 2.1 4 2.2 5 2.3 6 6 7 4.1 RC 7 4.2 RC 8 4.3 9 10 5.1 10 5.2 10 11 12 13-1 - Bond Behavior Between Corroded Rebar and Concrete Ema KATO* Mitsuyasu IWANAMI** Hiroshi YOKOTA*** Hajime ITO**** Fuminori

More information

.2 ρ dv dt = ρk grad p + 3 η grad (divv) + η 2 v.3 divh = 0, rote + c H t = 0 dive = ρ, H = 0, E = ρ, roth c E t = c ρv E + H c t = 0 H c E t = c ρv T

.2 ρ dv dt = ρk grad p + 3 η grad (divv) + η 2 v.3 divh = 0, rote + c H t = 0 dive = ρ, H = 0, E = ρ, roth c E t = c ρv E + H c t = 0 H c E t = c ρv T NHK 204 2 0 203 2 24 ( ) 7 00 7 50 203 2 25 ( ) 7 00 7 50 203 2 26 ( ) 7 00 7 50 203 2 27 ( ) 7 00 7 50 I. ( ν R n 2 ) m 2 n m, R = e 2 8πε 0 hca B =.09737 0 7 m ( ν = ) λ a B = 4πε 0ħ 2 m e e 2 = 5.2977

More information

ii 3.,. 4. F. (), ,,. 8.,. 1. (75%) (25%) =7 20, =7 21 (. ). 1.,, (). 3.,. 1. ().,.,.,.,.,. () (12 )., (), 0. 2., 1., 0,.

ii 3.,. 4. F. (), ,,. 8.,. 1. (75%) (25%) =7 20, =7 21 (. ). 1.,, (). 3.,. 1. ().,.,.,.,.,. () (12 )., (), 0. 2., 1., 0,. 24(2012) (1 C106) 4 11 (2 C206) 4 12 http://www.math.is.tohoku.ac.jp/~obata,.,,,.. 1. 2. 3. 4. 5. 6. 7.,,. 1., 2007 (). 2. P. G. Hoel, 1995. 3... 1... 2.,,. ii 3.,. 4. F. (),.. 5... 6.. 7.,,. 8.,. 1. (75%)

More information

Stata 11 Stata ts (ARMA) ARCH/GARCH whitepaper mwp 3 mwp-083 arch ARCH 11 mwp-051 arch postestimation 27 mwp-056 arima ARMA 35 mwp-003 arima postestim

Stata 11 Stata ts (ARMA) ARCH/GARCH whitepaper mwp 3 mwp-083 arch ARCH 11 mwp-051 arch postestimation 27 mwp-056 arima ARMA 35 mwp-003 arima postestim TS001 Stata 11 Stata ts (ARMA) ARCH/GARCH whitepaper mwp 3 mwp-083 arch ARCH 11 mwp-051 arch postestimation 27 mwp-056 arima ARMA 35 mwp-003 arima postestimation 49 mwp-055 corrgram/ac/pac 56 mwp-009 dfgls

More information

基礎数学I

基礎数学I I & II ii ii........... 22................. 25 12............... 28.................. 28.................... 31............. 32.................. 34 3 1 9.................... 1....................... 1............

More information

Lecture 12. Properties of Expanders

Lecture 12. Properties of Expanders Lecture 12. Properties of Expanders M2 Mitsuru Kusumoto Kyoto University 2013/10/29 Preliminalies G = (V, E) L G : A G : 0 = λ 1 λ 2 λ n : L G ψ 1,..., ψ n : L G µ 1 µ 2 µ n : A G ϕ 1,..., ϕ n : A G (Lecture

More information

201711grade1ouyou.pdf

201711grade1ouyou.pdf 2017 11 26 1 2 52 3 12 13 22 23 32 33 42 3 5 3 4 90 5 6 A 1 2 Web Web 3 4 1 2... 5 6 7 7 44 8 9 1 2 3 1 p p >2 2 A 1 2 0.6 0.4 0.52... (a) 0.6 0.4...... B 1 2 0.8-0.2 0.52..... (b) 0.6 0.52.... 1 A B 2

More information

1 Tokyo Daily Rainfall (mm) Days (mm)

1 Tokyo Daily Rainfall (mm) Days (mm) ( ) r-taka@maritime.kobe-u.ac.jp 1 Tokyo Daily Rainfall (mm) 0 100 200 300 0 10000 20000 30000 40000 50000 Days (mm) 1876 1 1 2013 12 31 Tokyo, 1876 Daily Rainfall (mm) 0 50 100 150 0 100 200 300 Tokyo,

More information

10

10 z c j = N 1 N t= j1 [ ( z t z ) ( )] z t j z q 2 1 2 r j /N j=1 1/ N J Q = N(N 2) 1 N j j=1 r j 2 2 χ J B d z t = z t d (1 B) 2 z t = (z t z t 1 ) (z t 1 z t 2 ) (1 B s )z t = z t z t s _ARIMA CONSUME

More information

LCC LCC INOUE, Gaku TANSEI, Kiyoteru KIDO, Motohiro IMAMURA, Takahiro LCC 7 LCC Ryanair 1 Ryanair Number of Passengers 2,000,000 1,800,000 1,

LCC LCC INOUE, Gaku TANSEI, Kiyoteru KIDO, Motohiro IMAMURA, Takahiro LCC 7 LCC Ryanair 1 Ryanair Number of Passengers 2,000,000 1,800,000 1, LCC LCC INOUE, Gaku TANSEI, Kiyoteru KIDO, Motohiro IMAMURA, Takahiro 1 23 LCC 7 LCC Ryanair 1 Ryanair Number of Passengers 2,, 1,8, 1,6, 1,4, 1,2, 1,, 8, 6, 4, 2, 1985 1987 199 1993 1996 1998 Barrett

More information

D = [a, b] [c, d] D ij P ij (ξ ij, η ij ) f S(f,, {P ij }) S(f,, {P ij }) = = k m i=1 j=1 m n f(ξ ij, η ij )(x i x i 1 )(y j y j 1 ) = i=1 j

D = [a, b] [c, d] D ij P ij (ξ ij, η ij ) f S(f,, {P ij }) S(f,, {P ij }) = = k m i=1 j=1 m n f(ξ ij, η ij )(x i x i 1 )(y j y j 1 ) = i=1 j 6 6.. [, b] [, d] ij P ij ξ ij, η ij f Sf,, {P ij } Sf,, {P ij } k m i j m fξ ij, η ij i i j j i j i m i j k i i j j m i i j j k i i j j kb d {P ij } lim Sf,, {P ij} kb d f, k [, b] [, d] f, d kb d 6..

More information

時系列データ解析による予測と最適化 ~エネルギー需要、発電、価格のモデリング~

時系列データ解析による予測と最適化 ~エネルギー需要、発電、価格のモデリング~ 時系列データ解析による予測と最適化 ~ エネルギー需要 発電 価格のモデリング ~ MahWorks Japan アプリケーションエンジニアリング部テクニカルコンピューティング 中川慶子 2015 The MahWorks, Inc. 1 アジェンダ 需要予測 : 時系列データモデリング 1. 機械学習ニューラルネットワーク データの取得 生データの前処理 ニューラルネットワークによるモデリング 2.

More information

研修コーナー

研修コーナー l l l l l l l l l l l α α β l µ l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l

More information

Myers, Montgomery & Anderson-Cook (2009) Response Surface Methodology

Myers, Montgomery & Anderson-Cook (2009) Response Surface Methodology Myers, R.H., Montgomery, D.C. & Anderson-Cook, C.M. (2009) Response Surface Methodology, Third Edition. Chapter 7. Experimantal Designs for Fitting Response Surfaces - I. (response surface methodology)

More information

seminar0220a.dvi

seminar0220a.dvi 1 Hi-Stat 2 16 2 20 16:30-18:00 2 2 217 1 COE 4 COE RA E-MAIL: ged0104@srv.cc.hit-u.ac.jp 2004 2 25 S-PLUS S-PLUS S-PLUS S-code 2 [8] [8] [8] 1 2 ARFIMA(p, d, q) FI(d) φ(l)(1 L) d x t = θ(l)ε t ({ε t }

More information

COE-RES Discussion Paper Series Center of Excellence Project The Normative Evaluation and Social Choice of Contemporary Economic Systems Graduate Scho

COE-RES Discussion Paper Series Center of Excellence Project The Normative Evaluation and Social Choice of Contemporary Economic Systems Graduate Scho COE-RES Discussion Paper Series Center of Excellence Project The Normative Evaluation and Social Choice of Contemporary Economic Systems Graduate School of Economics and Institute of Economic Research

More information

renshumondai-kaito.dvi

renshumondai-kaito.dvi 3 1 13 14 1.1 1 44.5 39.5 49.5 2 0.10 2 0.10 54.5 49.5 59.5 5 0.25 7 0.35 64.5 59.5 69.5 8 0.40 15 0.75 74.5 69.5 79.5 3 0.15 18 0.90 84.5 79.5 89.5 2 0.10 20 1.00 20 1.00 2 1.2 1 16.5 20.5 12.5 2 0.10

More information

On model selection problems in terms of prediction mean squared error and interpretaion of AIC (slides)

On model selection problems in terms of prediction mean squared error and interpretaion of AIC (slides) Applications in Econometrics and Finance by Long Memory Processes 2007 11 6 13:30-16:10 Table of Contents PART1 PART2 PART3 PART1 1 {y t } ρ(h) = ( h ) = Cov[y t y t+h ]/ Var[y t ] (yt y)(y t+h y) ρ(h)

More information

II III II 1 III ( ) [2] [3] [1] 1 1:

II III II 1 III ( ) [2] [3] [1] 1 1: 2015 4 16 1. II III II 1 III () [2] [3] 2013 11 18 [1] 1 1: [5] [6] () [7] [1] [1] 1998 4 2008 8 2014 8 6 [1] [1] 2 3 4 5 2. 2.1. t Dt L DF t A t (2.1) A t = Dt L + Dt F (2.1) 3 2 1 2008 9 2008 8 2008

More information

,,,,., = (),, (1) (4) :,,,, (1),. (2),, =. (3),,. (4),,,,.. (1) (3), (4).,,., () : = , ( ) : = F 1 + F 2 + F 3 + ( ) : = i Fj j=1 2

,,,,., = (),, (1) (4) :,,,, (1),. (2),, =. (3),,. (4),,,,.. (1) (3), (4).,,., () : = , ( ) : = F 1 + F 2 + F 3 + ( ) : = i Fj j=1 2 6 2 6.1 2 2, 2 5.2 R 2, 2 (R 2, B, µ)., R 2,,., 1, 2, 3,., 1, 2, 3,,. () : = 1 + 2 + 3 + (6.1.1).,,, 1 ,,,,., = (),, (1) (4) :,,,, (1),. (2),, =. (3),,. (4),,,,.. (1) (3), (4).,,., () : = 1 + 2 + 3 +,

More information

ii 3.,. 4. F. ( ), ,,. 8.,. 1. (75% ) (25% ) =7 24, =7 25, =7 26 (. ). 1.,, ( ). 3.,...,.,.,.,.,. ( ) (1 2 )., ( ), 0., 1., 0,.

ii 3.,. 4. F. ( ), ,,. 8.,. 1. (75% ) (25% ) =7 24, =7 25, =7 26 (. ). 1.,, ( ). 3.,...,.,.,.,.,. ( ) (1 2 )., ( ), 0., 1., 0,. (1 C205) 4 10 (2 C206) 4 11 (2 B202) 4 12 25(2013) http://www.math.is.tohoku.ac.jp/~obata,.,,,..,,. 1. 2. 3. 4. 5. 6. 7. 8. 1., 2007 ( ).,. 2. P. G., 1995. 3. J. C., 1988. 1... 2.,,. ii 3.,. 4. F. ( ),..

More information

Dirac 38 5 Dirac 4 4 γ µ p µ p µ + m 2 = ( p µ γ µ + m)(p ν γ ν + m) (5.1) γ = p µ p ν γ µ γ ν p µ γ µ m + mp ν γ ν + m 2 = 1 2 p µp ν {γ µ, γ ν } + m

Dirac 38 5 Dirac 4 4 γ µ p µ p µ + m 2 = ( p µ γ µ + m)(p ν γ ν + m) (5.1) γ = p µ p ν γ µ γ ν p µ γ µ m + mp ν γ ν + m 2 = 1 2 p µp ν {γ µ, γ ν } + m Dirac 38 5 Dirac 4 4 γ µ p µ p µ + m 2 p µ γ µ + mp ν γ ν + m 5.1 γ p µ p ν γ µ γ ν p µ γ µ m + mp ν γ ν + m 2 1 2 p µp ν {γ µ, γ ν } + m 2 5.2 p m p p µ γ µ {, } 10 γ {γ µ, γ ν } 2η µν 5.3 p µ γ µ + mp

More information

% 10%, 35%( 1029 ) p (a) 1 p 95% (b) 1 Std. Err. (c) p 40% 5% (d) p 1: STATA (1). prtesti One-sample test of pr

% 10%, 35%( 1029 ) p (a) 1 p 95% (b) 1 Std. Err. (c) p 40% 5% (d) p 1: STATA (1). prtesti One-sample test of pr 1 1. 2014 6 2014 6 10 10% 10%, 35%( 1029 ) p (a) 1 p 95% (b) 1 Std. Err. (c) p 40% 5% (d) p 1: STATA (1). prtesti 1029 0.35 0.40 One-sample test of proportion x: Number of obs = 1029 Variable Mean Std.

More information

untitled

untitled 17 5 13 1 2 1.1... 2 1.2... 2 1.3... 3 2 3 2.1... 3 2.2... 5 3 6 3.1... 6 3.2... 7 3.3 t... 7 3.4 BC a... 9 3.5... 10 4 11 1 1 θ n ˆθ. ˆθ, ˆθ, ˆθ.,, ˆθ.,.,,,. 1.1 ˆθ σ 2 = E(ˆθ E ˆθ) 2 b = E(ˆθ θ). Y 1,,Y

More information

SO(3) 7 = = 1 ( r ) + 1 r r r r ( l ) (5.17) l = 1 ( sin θ ) + sin θ θ θ ϕ (5.18) χ(r)ψ(θ, ϕ) l ψ = αψ (5.19) l 1 = i(sin ϕ θ l = i( cos ϕ θ l 3 = i ϕ

SO(3) 7 = = 1 ( r ) + 1 r r r r ( l ) (5.17) l = 1 ( sin θ ) + sin θ θ θ ϕ (5.18) χ(r)ψ(θ, ϕ) l ψ = αψ (5.19) l 1 = i(sin ϕ θ l = i( cos ϕ θ l 3 = i ϕ SO(3) 71 5.7 5.7.1 1 ħ L k l k l k = iϵ kij x i j (5.117) l k SO(3) l z l ± = l 1 ± il = i(y z z y ) ± (z x x z ) = ( x iy) z ± z( x ± i y ) = X ± z ± z (5.118) l z = i(x y y x ) = 1 [(x + iy)( x i y )

More information

6.1 (P (P (P (P (P (P (, P (, P.

6.1 (P (P (P (P (P (P (, P (, P. (011 30 7 0 ( ( 3 ( 010 1 (P.3 1 1.1 (P.4.................. 1 1. (P.4............... 1 (P.15.1 (P.16................. (P.0............3 (P.18 3.4 (P.3............... 4 3 (P.9 4 3.1 (P.30........... 4 3.

More information

..3. Ω, Ω F, P Ω, F, P ). ) F a) A, A,..., A i,... F A i F. b) A F A c F c) Ω F. ) A F A P A),. a) 0 P A) b) P Ω) c) [ ] A, A,..., A i,... F i j A i A

..3. Ω, Ω F, P Ω, F, P ). ) F a) A, A,..., A i,... F A i F. b) A F A c F c) Ω F. ) A F A P A),. a) 0 P A) b) P Ω) c) [ ] A, A,..., A i,... F i j A i A .. Laplace ). A... i),. ω i i ). {ω,..., ω } Ω,. ii) Ω. Ω. A ) r, A P A) P A) r... ).. Ω {,, 3, 4, 5, 6}. i i 6). A {, 4, 6} P A) P A) 3 6. ).. i, j i, j) ) Ω {i, j) i 6, j 6}., 36. A. A {i, j) i j }.

More information

量子力学 問題

量子力学 問題 3 : 203 : 0. H = 0 0 2 6 0 () = 6, 2 = 2, 3 = 3 3 H 6 2 3 ϵ,2,3 (2) ψ = (, 2, 3 ) ψ Hψ H (3) P i = i i P P 2 = P 2 P 3 = P 3 P = O, P 2 i = P i (4) P + P 2 + P 3 = E 3 (5) i ϵ ip i H 0 0 (6) R = 0 0 [H,

More information

ii 3.,. 4. F. (), ,,. 8.,. 1. (75% ) (25% ) =9 7, =9 8 (. ). 1.,, (). 3.,. 1. ( ).,.,.,.,.,. ( ) (1 2 )., ( ), 0. 2., 1., 0,.

ii 3.,. 4. F. (), ,,. 8.,. 1. (75% ) (25% ) =9 7, =9 8 (. ). 1.,, (). 3.,. 1. ( ).,.,.,.,.,. ( ) (1 2 )., ( ), 0. 2., 1., 0,. 23(2011) (1 C104) 5 11 (2 C206) 5 12 http://www.math.is.tohoku.ac.jp/~obata,.,,,.. 1. 2. 3. 4. 5. 6. 7.,,. 1., 2007 ( ). 2. P. G. Hoel, 1995. 3... 1... 2.,,. ii 3.,. 4. F. (),.. 5.. 6.. 7.,,. 8.,. 1. (75%

More information

評論・社会科学 123号(P)☆/1.福田

評論・社会科学 123号(P)☆/1.福田 VECM 2002 1 2007 12 VECM 55 VECM 1 2 3 3-1ECM 3-2 3-3VECM 3-4 4 1 2017 9 28 2017 10 16 2 1 1 2015, 79-85 BLUEBest Linear Unbiased Estimator 1 Vector Error-Correction Model VECM II III VECM VECM 55 15-34

More information

Hi-Stat Discussion Paper Series No.248 東京圏における 1990 年代以降の住み替え行動 住宅需要実態調査 を用いた Mixed Logit 分析 小林庸平行武憲史 March 2008 Hitotsubashi University Research Unit

Hi-Stat Discussion Paper Series No.248 東京圏における 1990 年代以降の住み替え行動 住宅需要実態調査 を用いた Mixed Logit 分析 小林庸平行武憲史 March 2008 Hitotsubashi University Research Unit Hi-Stat Discussion Paper Series No.248 東京圏における 1990 年代以降の住み替え行動 住宅需要実態調査 を用いた Logit 分析 小林庸平行武憲史 March 2008 Hitotsubashi University Research Unit for Statistical Analysis in Social Sciences A 21st-Century

More information

JFE.dvi

JFE.dvi ,, Department of Civil Engineering, Chuo University Kasuga 1-13-27, Bunkyo-ku, Tokyo 112 8551, JAPAN E-mail : atsu1005@kc.chuo-u.ac.jp E-mail : kawa@civil.chuo-u.ac.jp SATO KOGYO CO., LTD. 12-20, Nihonbashi-Honcho

More information

, 1), 2) (Markov-Switching Vector Autoregression, MSVAR), 3) 3, ,, , TOPIX, , explosive. 2,.,,,.,, 1

, 1), 2) (Markov-Switching Vector Autoregression, MSVAR), 3) 3, ,, , TOPIX, , explosive. 2,.,,,.,, 1 2016 1 12 4 1 2016 1 12, 1), 2) (Markov-Switching Vector Autoregression, MSVAR), 3) 3, 1980 1990.,, 225 1986 4 1990 6, TOPIX,1986 5 1990 2, explosive. 2,.,,,.,, 1986 Q2 1990 Q2,,. :, explosive, recursiveadf,

More information

1. 1 A : l l : (1) l m (m 3) (2) m (3) n (n 3) (4) A α, β γ α β + γ = 2 m l lm n nα nα = lm. α = lm n. m lm 2β 2β = lm β = lm 2. γ l 2. 3

1. 1 A : l l : (1) l m (m 3) (2) m (3) n (n 3) (4) A α, β γ α β + γ = 2 m l lm n nα nα = lm. α = lm n. m lm 2β 2β = lm β = lm 2. γ l 2. 3 1. 1 A : l l : (1) l m (m 3) (2) m (3) n (n 3) (4) A 2 1 2 1 2 3 α, β γ α β + γ = 2 m l lm n nα nα = lm. α = lm n. m lm 2β 2β = lm β = lm 2. γ l 2. 3 4 P, Q R n = {(x 1, x 2,, x n ) ; x 1, x 2,, x n R}

More information

Auerbach and Kotlikoff(1987) (1987) (1988) 4 (2004) 5 Diamond(1965) Auerbach and Kotlikoff(1987) 1 ( ) ,

Auerbach and Kotlikoff(1987) (1987) (1988) 4 (2004) 5 Diamond(1965) Auerbach and Kotlikoff(1987) 1 ( ) , ,, 2010 8 24 2010 9 14 A B C A (B Negishi(1960) (C) ( 22 3 27 ) E-mail:fujii@econ.kobe-u.ac.jp E-mail:082e527e@stu.kobe-u.ac.jp E-mail:iritani@econ.kobe-u.ac.jp 1 1 1 2 3 Auerbach and Kotlikoff(1987) (1987)

More information

○松本委員

○松本委員 CIRJE-J-100 2003 11 CIRJE hp://www.e.u-okyo.ac.jp/cirje/research/03research02dp_j.hml Credi Risk Modeling Approaches 2003 11 17 Absrac This aricle originaes from a speech given by he auhor in he seminar

More information

80 X 1, X 2,, X n ( λ ) λ P(X = x) = f (x; λ) = λx e λ, x = 0, 1, 2, x! l(λ) = n f (x i ; λ) = i=1 i=1 n λ x i e λ i=1 x i! = λ n i=1 x i e nλ n i=1 x

80 X 1, X 2,, X n ( λ ) λ P(X = x) = f (x; λ) = λx e λ, x = 0, 1, 2, x! l(λ) = n f (x i ; λ) = i=1 i=1 n λ x i e λ i=1 x i! = λ n i=1 x i e nλ n i=1 x 80 X 1, X 2,, X n ( λ ) λ P(X = x) = f (x; λ) = λx e λ, x = 0, 1, 2, x! l(λ) = n f (x i ; λ) = n λ x i e λ x i! = λ n x i e nλ n x i! n n log l(λ) = log(λ) x i nλ log( x i!) log l(λ) λ = 1 λ n x i n =

More information

waseda2010a-jukaiki1-main.dvi

waseda2010a-jukaiki1-main.dvi November, 2 Contents 6 2 8 3 3 3 32 32 33 5 34 34 6 35 35 7 4 R 2 7 4 4 9 42 42 2 43 44 2 5 : 2 5 5 23 52 52 23 53 53 23 54 24 6 24 6 6 26 62 62 26 63 t 27 7 27 7 7 28 72 72 28 73 36) 29 8 29 8 29 82 3

More information

( ) Note (e ) (µ ) (τ ) ( (ν e,e ) e- (ν µ, µ ) µ- (ν τ,τ ) τ- ) ( ) ( ) (SU(2) ) (W +,Z 0,W ) * 1) 3 * 2) [ ] [ ] [ ] ν e ν µ ν τ e

( ) Note (e ) (µ ) (τ ) ( (ν e,e ) e- (ν µ, µ ) µ- (ν τ,τ ) τ- ) ( ) ( ) (SU(2) ) (W +,Z 0,W ) * 1) 3 * 2) [ ] [ ] [ ] ν e ν µ ν τ e ( ) Note 3 19 12 13 8 8.1 (e ) (µ ) (τ ) ( (ν e,e ) e- (ν µ, µ ) µ- (ν τ,τ ) τ- ) ( ) ( ) (SU(2) ) (W +,Z 0,W ) * 1) 3 * 2) [ ] [ ] [ ] ν e ν µ ν τ e µ τ, e R, µ R, τ R (1a) L ( ) ) * 3) W Z 1/2 ( - )

More information

AR(1) y t = φy t 1 + ɛ t, ɛ t N(0, σ 2 ) 1. Mean of y t given y t 1, y t 2, E(y t y t 1, y t 2, ) = φy t 1 2. Variance of y t given y t 1, y t

AR(1) y t = φy t 1 + ɛ t, ɛ t N(0, σ 2 ) 1. Mean of y t given y t 1, y t 2, E(y t y t 1, y t 2, ) = φy t 1 2. Variance of y t given y t 1, y t 87 6.1 AR(1) y t = φy t 1 + ɛ t, ɛ t N(0, σ 2 ) 1. Mean of y t given y t 1, y t 2, E(y t y t 1, y t 2, ) = φy t 1 2. Variance of y t given y t 1, y t 2, V(y t y t 1, y t 2, ) = σ 2 3. Thus, y t y t 1,

More information

151021slide.dvi

151021slide.dvi : Mac I 1 ( 5 Windows (Mac Excel : Excel 2007 9 10 1 4 http://asakura.co.jp/ books/isbn/978-4-254-12172-8/ (1 1 9 1/29 (,,... (,,,... (,,, (3 3/29 (, (F7, Ctrl + i, (Shift +, Shift + Ctrl (, a i (, Enter,

More information

総研大恒星進化概要.dvi

総研大恒星進化概要.dvi The Structure and Evolution of Stars I. Basic Equations. M r r =4πr2 ρ () P r = GM rρ. r 2 (2) r: M r : P and ρ: G: M r Lagrange r = M r 4πr 2 rho ( ) P = GM r M r 4πr. 4 (2 ) s(ρ, P ) s(ρ, P ) r L r T

More information

autocorrelataion cross-autocorrelataion Lo/MacKinlay [1988, 1990] (A)

autocorrelataion cross-autocorrelataion Lo/MacKinlay [1988, 1990] (A) Discussion Paper Series A No.425 2002 2 186-8603 iwaisako@ier.hit-u.ac.jp 14 1 24 autocorrelataion cross-autocorrelataion Lo/MacKinlay [1988, 1990] 1990 12 13 (A) 12370027 13 1 1980 Lo/MacKinlay [1988]

More information

財政赤字の経済分析:中長期的視点からの考察

財政赤字の経済分析:中長期的視点からの考察 000 364 645 (ax smoohing (Barro(979 Barro(979 (consumpion smoohing Hall(978 Campbell(987 Barro(979 Hall(978 Campbell(987 Campbell(987 B ( r B G T ( ( r E B G T Y T/Y r ( No Ponzi Game ( [ Ω ] Y ( lim (

More information

6.1 (P (P (P (P (P (P (, P (, P.101

6.1 (P (P (P (P (P (P (, P (, P.101 (008 0 3 7 ( ( ( 00 1 (P.3 1 1.1 (P.3.................. 1 1. (P.4............... 1 (P.15.1 (P.15................. (P.18............3 (P.17......... 3.4 (P................ 4 3 (P.7 4 3.1 ( P.7...........

More information

( ) s n (n = 0, 1,...) n n = δ nn n n = I n=0 ψ = n C n n (1) C n = n ψ α = e 1 2 α 2 n=0 α, β α n n! n (2) β α = e 1 2 α 2 1

( ) s n (n = 0, 1,...) n n = δ nn n n = I n=0 ψ = n C n n (1) C n = n ψ α = e 1 2 α 2 n=0 α, β α n n! n (2) β α = e 1 2 α 2 1 (3.5 3.8) 03032s 2006.7.0 n (n = 0,,...) n n = δ nn n n = I n=0 ψ = n C n n () C n = n ψ α = e 2 α 2 n=0 α, β α n n (2) β α = e 2 α 2 2 β 2 n=0 =0 = e 2 α 2 β n α 2 β 2 n=0 = e 2 α 2 2 β 2 +β α β n α!

More information

カルマンフィルターによるベータ推定( )

カルマンフィルターによるベータ推定( ) β TOPIX 1 22 β β smoothness priors (the Capital Asset Pricing Model, CAPM) CAPM 1 β β β β smoothness priors :,,. E-mail: koiti@ism.ac.jp., 104 1 TOPIX β Z i = β i Z m + α i (1) Z i Z m α i α i β i (the

More information

( )/2 hara/lectures/lectures-j.html 2, {H} {T } S = {H, T } {(H, H), (H, T )} {(H, T ), (T, T )} {(H, H), (T, T )} {1

( )/2   hara/lectures/lectures-j.html 2, {H} {T } S = {H, T } {(H, H), (H, T )} {(H, T ), (T, T )} {(H, H), (T, T )} {1 ( )/2 http://www2.math.kyushu-u.ac.jp/ hara/lectures/lectures-j.html 1 2011 ( )/2 2 2011 4 1 2 1.1 1 2 1 2 3 4 5 1.1.1 sample space S S = {H, T } H T T H S = {(H, H), (H, T ), (T, H), (T, T )} (T, H) S

More information

1 n 1 1 2 2 3 3 3.1............................ 3 3.2............................. 6 3.2.1.............. 6 3.2.2................. 7 3.2.3........................... 10 4 11 4.1..........................

More information

44 4 I (1) ( ) (10 15 ) ( 17 ) ( 3 1 ) (2)

44 4 I (1) ( ) (10 15 ) ( 17 ) ( 3 1 ) (2) (1) I 44 II 45 III 47 IV 52 44 4 I (1) ( ) 1945 8 9 (10 15 ) ( 17 ) ( 3 1 ) (2) 45 II 1 (3) 511 ( 451 1 ) ( ) 365 1 2 512 1 2 365 1 2 363 2 ( ) 3 ( ) ( 451 2 ( 314 1 ) ( 339 1 4 ) 337 2 3 ) 363 (4) 46

More information

i ii i iii iv 1 3 3 10 14 17 17 18 22 23 28 29 31 36 37 39 40 43 48 59 70 75 75 77 90 95 102 107 109 110 118 125 128 130 132 134 48 43 43 51 52 61 61 64 62 124 70 58 3 10 17 29 78 82 85 102 95 109 iii

More information

1 (Contents) (1) Beginning of the Universe, Dark Energy and Dark Matter Noboru NAKANISHI 2 2. Problem of Heat Exchanger (1) Kenji

1 (Contents) (1) Beginning of the Universe, Dark Energy and Dark Matter Noboru NAKANISHI 2 2. Problem of Heat Exchanger (1) Kenji 8 4 2018 6 2018 6 7 1 (Contents) 1. 2 2. (1) 22 3. 31 1. Beginning of the Universe, Dark Energy and Dark Matter Noboru NAKANISHI 2 2. Problem of Heat Exchanger (1) Kenji SETO 22 3. Editorial Comments Tadashi

More information

オーストラリア研究紀要 36号(P)☆/3.橋本

オーストラリア研究紀要 36号(P)☆/3.橋本 36 p.9 202010 Tourism Demand and the per capita GDP : Evidence from Australia Keiji Hashimoto Otemon Gakuin University Abstract Using Australian quarterly data1981: 2 2009: 4some time-series econometrics

More information

untitled

untitled 2 : n =1, 2,, 10000 0.5125 0.51 0.5075 0.505 0.5025 0.5 0.4975 0.495 0 2000 4000 6000 8000 10000 2 weak law of large numbers 1. X 1,X 2,,X n 2. µ = E(X i ),i=1, 2,,n 3. σi 2 = V (X i ) σ 2,i=1, 2,,n ɛ>0

More information

Autumn II III Zon and Muysken 2005 Zon and Muysken 2005 IV II 障害者への所得移転の経済効果 分析に用いるデータ

Autumn II III Zon and Muysken 2005 Zon and Muysken 2005 IV II 障害者への所得移転の経済効果 分析に用いるデータ 212 Vol. 44 No. 2 I はじめに 2008 1 2 Autumn 08 213 II III Zon and Muysken 2005 Zon and Muysken 2005 IV II 障害者への所得移転の経済効果 17 18 1 分析に用いるデータ 1 2005 10 12 200 2 2006 9 12 1 1 2 129 35 113 3 1 2 6 1 2 3 4 4 1

More information

I II III IV V

I II III IV V I II III IV V N/m 2 640 980 50 200 290 440 2m 50 4m 100 100 150 200 290 390 590 150 340 4m 6m 8m 100 170 250 µ = E FRVβ β N/mm 2 N/mm 2 1.1 F c t.1 3 1 1.1 1.1 2 2 2 2 F F b F s F c F t F b F s 3 3 3

More information

( ) 2002 1 1 1 1.1....................................... 1 1.1.1................................. 1 1.1.2................................. 1 1.1.3................... 3 1.1.4......................................

More information

02.„o“φiflì„㙃fic†j

02.„o“φiflì„㙃fic†j X-12-ARIMA Band-PassDECOMP HP X-12-ARIMADECOMP HPBeveridge and Nelson DECOMP X-12-ARIMA Band-PassHodrick and PrescottHP DECOMPBeveridge and Nelson M CD X ARIMA DECOMP HP Band-PassDECOMP Kiyotaki and Moore

More information

デフォルト相関係数のインプライド推計( )

デフォルト相関係数のインプライド推計( ) ,,,.,.,.,,.,,.,,.,,. :,, OU, SUR,.,.,,.,.., 62 1, BIS,.,.,,,,.,,.,.,, [33],,,.,,.,,,,.,,.,,,.. 2,. 3,. 4,,. 5. 2,,.,. 2.1,.,,.,. Giesecke [10] dependen defaul.,,,, GDP,.,,.,...,,,.,.,.,,.,,. Giesecke [10],

More information

1 (Berry,1975) 2-6 p (S πr 2 )p πr 2 p 2πRγ p p = 2γ R (2.5).1-1 : : : : ( ).2 α, β α, β () X S = X X α X β (.1) 1 2

1 (Berry,1975) 2-6 p (S πr 2 )p πr 2 p 2πRγ p p = 2γ R (2.5).1-1 : : : : ( ).2 α, β α, β () X S = X X α X β (.1) 1 2 2005 9/8-11 2 2.2 ( 2-5) γ ( ) γ cos θ 2πr πρhr 2 g h = 2γ cos θ ρgr (2.1) γ = ρgrh (2.2) 2 cos θ θ cos θ = 1 (2.2) γ = 1 ρgrh (2.) 2 2. p p ρgh p ( ) p p = p ρgh (2.) h p p = 2γ r 1 1 (Berry,1975) 2-6

More information

ファイナンスのための数学基礎 第1回 オリエンテーション、ベクトル

ファイナンスのための数学基礎 第1回 オリエンテーション、ベクトル 時系列分析 変量時系列モデルとその性質 担当 : 長倉大輔 ( ながくらだいすけ 時系列モデル 時系列モデルとは時系列データを生み出すメカニズムとなるものである これは実際には未知である 私たちにできるのは観測された時系列データからその背後にある時系列モデルを推測 推定するだけである 以下ではいくつかの代表的な時系列モデルを考察する 自己回帰モデル (Auoregressive Model もっとも頻繁に使われる時系列モデルは自己回帰モデル

More information

# GDP 1 GDP 1992 GDP Hamilton and Flavin (1986) Doi and Ihori(2003) (2007) 2 Bohn(1998) /GDP Bohn # 1 2 Broda

# GDP 1 GDP 1992 GDP Hamilton and Flavin (1986) Doi and Ihori(2003) (2007) 2 Bohn(1998) /GDP Bohn # 1 2 Broda # masaya@econ.keio.ac.jp GDP 1 GDP 1992 GDP 1 2006 Hamilon and Flavin (1986) Doi and Ihori(2003) 1965 2003 (2007) 2 Bohn(1998) /GDP Bohn # 1 2 Broda and Weinsein (2005) (2007) 1 Hamilon and Flavin(1986)

More information

t χ 2 F Q t χ 2 F 1 2 µ, σ 2 N(µ, σ 2 ) f(x µ, σ 2 ) = 1 ( exp (x ) µ)2 2πσ 2 2σ 2 0, N(0, 1) (100 α) z(α) t χ 2 *1 2.1 t (i)x N(µ, σ 2 ) x µ σ N(0, 1

t χ 2 F Q t χ 2 F 1 2 µ, σ 2 N(µ, σ 2 ) f(x µ, σ 2 ) = 1 ( exp (x ) µ)2 2πσ 2 2σ 2 0, N(0, 1) (100 α) z(α) t χ 2 *1 2.1 t (i)x N(µ, σ 2 ) x µ σ N(0, 1 t χ F Q t χ F µ, σ N(µ, σ ) f(x µ, σ ) = ( exp (x ) µ) πσ σ 0, N(0, ) (00 α) z(α) t χ *. t (i)x N(µ, σ ) x µ σ N(0, ) (ii)x,, x N(µ, σ ) x = x+ +x N(µ, σ ) (iii) (i),(ii) z = x µ N(0, ) σ N(0, ) ( 9 97.

More information

22 / ( ) OD (Origin-Destination)

22 / ( ) OD (Origin-Destination) 23 2 15 22 / ( ) OD (Origin-Destination) 1 1 2 3 2.1....................................... 3 2.2......................................... 3 2.3.......................................... 5 2.4............................

More information

Milnor 1 ( ), IX,. [KN].,. 2 : (1),. (2). 1 ; 1950, Milnor[M1, M2]. Milnor,,. ([Hil, HM, IO, St] ).,.,,, ( 2 5 )., Milnor ( 4.1)..,,., [CEGS],. Ω m, P

Milnor 1 ( ), IX,. [KN].,. 2 : (1),. (2). 1 ; 1950, Milnor[M1, M2]. Milnor,,. ([Hil, HM, IO, St] ).,.,,, ( 2 5 )., Milnor ( 4.1)..,,., [CEGS],. Ω m, P Milnor 1 ( ), IX,. [KN].,. 2 : (1),. (2). 1 ; 1950, Milnor[M1, M2]. Milnor,,. ([Hil, HM, IO, St] ).,.,,, ( 2 5 )., Milnor ( 4.1)..,,., [CEGS],. Ω m, PC ( 4 5 )., 5, Milnor Milnor., ( 6 )., (I) Z modulo

More information

1990年代以降の日本の経済変動

1990年代以降の日本の経済変動 1990 * kenichi.sakura@boj.or.jp ** hitoshi.sasaki@boj.or.jp *** masahiro.higo@boj.or.jp No.05-J-10 2005 12 103-8660 30 * ** *** 1990 2005 12 1990 1990 1990 2005 11 2425 BIS E-mail: kenichi.sakura@boj.or.jp

More information

K E N Z U 01 7 16 HP M. 1 1 4 1.1 3.......................... 4 1.................................... 4 1..1..................................... 4 1...................................... 5................................

More information

untitled

untitled 2 book conference 1990 2003 14 Repeated Cross-Section Data 1 M1,M2 M1 Sekine(1998) Repeated Cross-Section Data 1 1. (1989), Yoshida and Rasche(1990), Rasche(1990), 19921997, Fujiki and Mulligan(1996),

More information

* n x 11,, x 1n N(µ 1, σ 2 ) x 21,, x 2n N(µ 2, σ 2 ) H 0 µ 1 = µ 2 (= µ ) H 1 µ 1 µ 2 H 0, H 1 *2 σ 2 σ 2 0, σ 2 1 *1 *2 H 0 H

* n x 11,, x 1n N(µ 1, σ 2 ) x 21,, x 2n N(µ 2, σ 2 ) H 0 µ 1 = µ 2 (= µ ) H 1 µ 1 µ 2 H 0, H 1 *2 σ 2 σ 2 0, σ 2 1 *1 *2 H 0 H 1 1 1.1 *1 1. 1.3.1 n x 11,, x 1n Nµ 1, σ x 1,, x n Nµ, σ H 0 µ 1 = µ = µ H 1 µ 1 µ H 0, H 1 * σ σ 0, σ 1 *1 * H 0 H 0, H 1 H 1 1 H 0 µ, σ 0 H 1 µ 1, µ, σ 1 L 0 µ, σ x L 1 µ 1, µ, σ x x H 0 L 0 µ, σ 0

More information

q π =0 Ez,t =ε σ {e ikz ωt e ikz ωt } i/ = ε σ sinkz ωt 5.6 x σ σ *105 q π =1 Ez,t = 1 ε σ + ε π {e ikz ωt e ikz ωt } i/ = 1 ε σ + ε π sinkz ωt 5.7 σ

q π =0 Ez,t =ε σ {e ikz ωt e ikz ωt } i/ = ε σ sinkz ωt 5.6 x σ σ *105 q π =1 Ez,t = 1 ε σ + ε π {e ikz ωt e ikz ωt } i/ = 1 ε σ + ε π sinkz ωt 5.7 σ H k r,t= η 5 Stokes X k, k, ε, ε σ π X Stokes 5.1 5.1.1 Maxwell H = A A *10 A = 1 c A t 5.1 A kη r,t=ε η e ik r ωt 5. k ω ε η k η = σ, π ε σ, ε π σ π A k r,t= q η A kη r,t+qηa kηr,t 5.3 η q η E = 1 c A

More information