1 2 Octave/MATLAB Dynare Dynare Octave/MATLAB 1.1 Dynare Dynare Dynare DSGE 3 4 Dynare Octave MAT- LAB Dynare stable release
|
|
- ひろと まるこ
- 7 years ago
- Views:
Transcription
1 Dynare Ver Dynare Dynare Octave Dynare RBC.mod tmodel2.mod NK Linear.mod, NK Linear stoch.mod Dynare Dynare NK Linear EST2.mod mod
2 1 2 Octave/MATLAB Dynare Dynare Octave/MATLAB 1.1 Dynare Dynare Dynare DSGE 3 4 Dynare Octave MAT- LAB Dynare stable release Octave Octave GNU Dynare Octave MATLAB *1 Octave Windows *2 MinGW Microsoft Visual Studio 2 MinGW Octave Octave Dynare *1 MATLAB MathWorks *2 MATLAB 2
3 1 Octave Octave exe - Windows Octave Windows Octave exit 3
4 1.3 Dynare Octave/MATLAB.txt Dynare.mod Octave/MATLAB.m *3 DSGE Dynare mod mod Dynare Octave/MATLAB Dynare Octave/MATLAB *4 2 Octave Dynare mod RBC mod Dynare RBC mod RBC.mod C:\work\dsge *5 main.m addpath C:\dynare\4.3.1\matlab cd C:\work\dsge dynare RBC m main.m RBC.mod Octave 3 Octave Octave/Dynare RBC.mod *3 Windows 7 *4 Octave/MATLAB *5 \ Yen Y= 4
5 2 Octave main dynare RBC 3 Dynare 3.1 RBC.mod DSGE a a Dynare 4 3 RBC w t = (γ + 1) µl γ t C t (1) C t+1 = β (r t+1 δ + 1) C t (2) Y t = A t Kt α L 1 α t (3) w t = (1 α) A t Kt α L α t (4) r t = αa t Kt α 1 L 1 α t (5) K t+1 = Y t + (1 δ) K t C t (6) ln (A t+1 ) = ρ ln (A t ) + e t+1 (7) *6 *6 4 5
6 C t C L t L K t K Y t Y w t w r t r A t A e t e 2 RBC 1 RBC 2. 3 α alpha 0.3 β beta 0.99 δ delta µ mu 1.0 γ gamma 1.0 ρ AR(1) rho RBC 3. (1) (7) 7 6
7 4. *7 A = 1 (8) r = β 1 + δ 1 (9) K ( ) 1 r L = α 1 A (10) α Y ( ) K α L = A L (11) C L = Y L δ K L (12) ( ) K w = (1 α)a α (13) { w L = (γ + 1)µ L } 1 γ+1 ( C L ) 1 γ+1 Dynare *8 (14) 5. t = 0 t = 1 e 1 = % * mod RBC.mod 1. var varexo 2. parameters 3. model *7 4 *8 Octave fsolve * % 7
8 4. initval steady check 5. simul m % mod // * var varexo // 1. var C L K Y w r A; varexo e; 2. 2 parameters // 2. parameters alpha beta delta mu gamma rho; // alpha = 0.3; beta = 0.99; delta = 0.025; mu = 1.0; gamma = 1.0; rho = 0.9; *10 mod Octave/MATLAB mod % 8
9 3. model model model; end; t + 1 C(+1) t 1 A( 1) t + 1 K t A t 2 (6) (7) K t = Y t 1 + (1 δ) K t 1 C t 1 (15) ln (A t ) = ρ ln (A t 1 ) + e t (16) K t+1 A t+1 model * 11 // 3. model; w/c = (gamma+1)*mu*l^gamma; C(+1)/C = beta*(r(+1)-delta+1); Y = A*K^alpha*L^(1-alpha); w = (1-alpha)*A*K^alpha*L^(-alpha); r = alpha*a*k^(alpha-1)*l^(1-alpha); K = Y(-1)+(1-delta)*K(-1)-C(-1); log(a) = rho*log(a(-1)) + e; end; 4. initval Dynare star *11 log exp 9
10 // 4. Astar = 1; rstar = 1/beta + delta - 1; K_L = (rstar/alpha/astar)^(1/(alpha-1)); Y_L = Astar*K_L^alpha; C_L = Y_L-delta*K_L; wstar = (1-alpha)*Astar*K_L^alpha; Lstar = (wstar/(gamma+1)/mu)^(1/(gamma+1))*c_l^(-1/(gamma+1)); Kstar = K_L*Lstar; Ystar = Y_L*Lstar; Cstar = C_L*Lstar; // Dynare initval; C = Cstar; L = Lstar; K = Kstar; Y = Ystar; w = wstar; r = rstar; A = Astar; end; steady Dynare * 12 Dynare *12 initval steady Dynare steady histval 10
11 // Dynare steady; // [Cstar; Lstar; Kstar; Ystar; wstar; rstar; Astar] ans = Dynare ans = STEADY-STATE RESULTS: C L K Y w r A 1 ans = check // check; EIGENVALUES: check 1 11
12 C(+1) r(+1) 2 rank * 13 The rank condition is verified. EIGENVALUES: Modulus Real Imaginary 3.204e e e e Inf Inf 0 There are 2 eigenvalue(s) larger than 1 in modulus for 2 forward-looking variable(s) The rank condition is verified. 5. t = 0 t = 1 e 1 = % simul perfect foresight solution periods=150 t = 151 periods var C t C t = *13 Blanchard and Kahn [1980] 12
13 // 5. // shocks; var e; periods 1; values -0.05; end; // simul(periods=150); Dynare Octave/MATLAB simul 13
14 // C1 = (C./Cstar-1)*100; L1 = (L./Lstar-1)*100; K1 = (K./Kstar-1)*100; Y1 = (Y./Ystar-1)*100; w1 = (w./wstar-1)*100; r1 = (r-rstar)*100; A1 = (A./Astar-1)*100; I1 = ((Y-C)./(Ystar-Cstar)-1)*100; // figure(1) subplot(2,2,1) plot(0:50, A1(1:51)); title( A ) subplot(2,2,2) plot(0:50, Y1(1:51)); title( Y ) subplot(2,2,3) plot(0:50, C1(1:51)); title( C ) subplot(2,2,4) plot(0:50, K1(1:51)); title( K ) figure(2) subplot(2,2,1) plot(0:50, L1(1:51)); title( L ) subplot(2,2,2) plot(0:50, I1(1:51)); title( I ) subplot(2,2,3) plot(0:50, w1(1:51)); title( w ) subplot(2,2,4) plot(0:50, r1(1:51)); title( r ) 14
15 2 e 1 = RBC.mod 2 * 14 csv *14 print( filename.eps, -depsc2 ) 15
16 csvwrite( rbc_rslt.csv,[c, L, K, Y, w, r, A]); RBC.mod rbc rslt.csv csv csv EXCEL 3.2 tmodel2.mod 3 mod C t C K t K τ c,t tauc 0 0 τ k,t tauk g t g α alpha 0.3 β beta 0.99 δ delta 0.25 A t At
17 3. (1 + τ c,t+1 )C t+1 = β [ (1 τ k,t+1 )αa t+1 Kt+1 α 1 (1 + τ c,t )C t (17) K t+1 = A t Kt α + (1 δ)k t C t g t (18) 4. Dynare 5. t = 0 t = 10 τ k,t 0.1 t = 1 τ c,t mod tmodel2.mod 1. var varexo 2. parameters 3. model 4. initval 5. endval steady 6. check 7. simul RBC initval t = 0 initval steady Dynare initval Dynare * 15 Dynare * t = 0 *12 17
18 // 4. initval; C = 1; K = 1; tauc = 0; tauk = 0; g = 0.1; end; steady; endval Dynare steady Dynare // 5. endval; C = 1; K = 1; tauc = 0; tauk = 0.1; g = 0.1; end; steady; check
19 t = 0 t = 10 τ k,t 0.1 t = 1 // 7. shocks; var tauk; periods 1:9; values 0; end; simul(periods=31); 3.3 NK Linear.mod, NK Linear stoch.mod 5 Dynare ˆx t GDP x π t ppi â t a î t ii v t vv ν t â t+1 â t nu 7 e t e ε t eps 8 19
20 2. 9 β beta 0.99 γ gamma 1 ϱ varrho 0.8 κ κ = (1 ϱ)(1 ϱβ)(γ+1) ϱ kappa ϕ π phi pi 1.5 ϕ y GDP phi y 0.5 ρ A AR(1) rho A 0.9 ρ v AR(1) rho v π t = βπ t+1 + κˆx t (19) ˆx t = ˆx t+1 (î t π t+1 ) + ν t (20) î t = ϕ π π t + ϕ y ˆx t + v t (21) v t+1 = ρ v v t + ε t+1 (22) â t+1 = ρ A â t + e t+1 (23) ν t = â t+1 â t (24) 4. * t = 0 t = 1 +5% *16 20
21 initval steady simul stoch simul simul shocks; var e; periods 1; values 0.05; end; simul(periods=150); stoch simul shocks; var e = 5^2; end; stoch_simul(order=1, irf = 100) x ii ppi a v; 6 * 17 stoch simul RBC stoch.mod * 18 *17 1 mod simul NK Linear.mod stoch simul NK Linear stoch.mod stoch simul 1 *18 Dynare 21
22 3.4 Dynare simul stoch simu 2 Deterministic Simulation Stochastic Simulation 22
23 4 Dynare NK Linear EST2.mod DSGE Dynare * 19 Dynare MCMC M-H 4.1 DSGE π t = βπ t+1 + κˆx t (25) ˆx t = ˆx t+1 (î t π t+1 ) + (ρ A 1)â t (26) î t = (1 + ϕ π )π t + ϕ y ˆx t + v t (27) v t+1 = ρ v v t + u t+1 (28) â t+1 = ρ A â t + ε t+1 (29) x obs t x obs = ˆx t (30) π obs t π obs = π t + ϵ π,t (31) (i obs t ī obs )/4 = î t (32) *19 Dynare 23
24 * ˆx t GDP x π t ppi â t a î t ii v t v 10 e t e u t u ϵ π,t errppi β beta γ ϱ ϕ π ϕ y GDP ρ A AR(1) ρ v AR(1) σ ε ε t σ u u t σ ϵπ ϵ π,t 12 *
25 5. GDP x obs t i obs t 3 GDP ; πt obs 6. MCMC 20, * x act t π act t i act t = x obs t x obs (33) = π obs t π obs (34) = (i obs t ī obs )/4 (35) x act t, πt act, i act t mod script dataset.m *21 DSGE 25
26 script dataset.m EO90 = csvread( EO90.csv, 1, 1); GAP0 = EO90(:,1); PGDP0 = EO90(:,2); IRS0 = EO90(:,3); PC_PGDP0 = [NaN; (PGDP0(2:end)./PGDP0(1:end-1)-1)*100]; tt = 2:77; MGAP = mean(gap0(tt)); MPC_PGDP = mean(pc_pgdp0(tt)); MIRS = mean(irs0(tt)); GAP = GAP0(tt)-MGAP; PC_PGDP = PC_PGDP0(tt)-MPC_PGDP; IRS4 = (IRS0(tt)-MIRS)*0.25; x = GAP; ppiact = PC_PGDP; ii = IRS4; save dset.mat x ppiact ii; csv EO90.csv GDP (33) (35) x act t, πt act, i act t x, ppiact, ii dset.mat mat MATLAB mod 4.3 mod mod NK Linear EST2.mod 1. var varexo 2. parameters 26
27 3. model 4. estimated params 5. varobs 6. estimation 1. 3 var x ppi a ii v ppiact; varexo e u errppi; ppiact // 2. parameters beta gamma varrho phi_pi phi_y rho_a rho_v; // beta = 0.99; 3. (31) π act t = π t + ϵ π,t (36) π act t ppiact κ kappa # model 27
28 model(linear); # kappa = (1-varrho)*(1-varrho*beta)*(gamma+1)/varrho; ppi = beta*ppi(+1)+kappa*x; x = x(+1)-(ii-ppi(+1))+(rho_a-1)*a; ii = (phi_pi+1)*ppi + phi_y*x + v; a = rho_a*a(-1) + e; v = rho_v*v(-1) + u; ppiact = ppi + errppi; end; ˆx t, πt act, î t mod x, ppiact, ii dset.mat // 4. estimated_params; gamma, gamma_pdf, 1, 0.5; varrho, beta_pdf, 0.8, 0.1; phi_pi, gamma_pdf, 0.5, 0.25; phi_y, gamma_pdf, 0.5, 0.25; rho_a, beta_pdf, 0.8, 0.05; rho_v, beta_pdf, 0.8, 0.1; stderr e, inv_gamma_pdf, 0.5, 0.5; stderr u, inv_gamma_pdf, 0.5, 0.5; stderr errppi, inv_gamma_pdf, 0.5, 0.5; end; ˆx t, πt act, î t 28
29 // 5. varobs x ppiact ii; 6. estimation datafile dset mh replic MCMC 20,000 mh nblocks 2 mh drop 0.25 mh jscale M-H acceptation rate // 6. estimation(datafile=dset, mh_replic=20000, mh_nblocks=2, mh_drop = 0.25, mh_jscale=0.7); 4.4 mat mod 2 29
30 ESTIMATION RESULTS Log data density is parameters prior mean post. mean conf. interval prior pstdev gamma gamma varrho beta phi_pi gamma phi_y gamma rho_a beta rho_v beta standard deviation of shocks prior mean post. mean conf. interval prior pstdev e invg u invg errppi invg Total computing time : 0h02m00s 30
DSGE Dynamic Stochastic General Equilibrium Model DSGE 5 2 DSGE DSGE ω 0 < ω < 1 1 DSGE Blanchard and Kahn VAR 3 MCMC 2 5 4 1 1 1.1 1. 2. 118
7 DSGE 2013 3 7 1 118 1.1............................ 118 1.2................................... 123 1.3.............................. 125 1.4..................... 127 1.5...................... 128 1.6..............
More information受賞講演要旨2012cs3
アハ ート アハ ート アハ ート アハ ート アハ ート アハ ート アハ ート アハ ート アハ ート アハ ート アハ ート アハ ート アハ ート アハ ート アハ ート アハ ート アハ ート アハ ート アハ ート アハ ート アハ ート アハ ート アハ ート アハ ート アハ ート アハ ート アハ ート アハ ート アハ ート アハ ート アハ ート α β α α α α α
More information第86回日本感染症学会総会学術集会後抄録(II)
χ μ μ μ μ β β μ μ μ μ β μ μ μ β β β α β β β λ Ι β μ μ β Δ Δ Δ Δ Δ μ μ α φ φ φ α γ φ φ γ φ φ γ γδ φ γδ γ φ φ φ φ φ φ φ φ φ φ φ φ φ α γ γ γ α α α α α γ γ γ γ γ γ γ α γ α γ γ μ μ κ κ α α α β α
More information一般演題(ポスター)
6 5 13 : 00 14 : 00 A μ 13 : 00 14 : 00 A β β β 13 : 00 14 : 00 A 13 : 00 14 : 00 A 13 : 00 14 : 00 A β 13 : 00 14 : 00 A β 13 : 00 14 : 00 A 13 : 00 14 : 00 A β 13 : 00 14 : 00 A 13 : 00 14 : 00 A
More information基礎数学I
I & II ii ii........... 22................. 25 12............... 28.................. 28.................... 31............. 32.................. 34 3 1 9.................... 1....................... 1............
More information467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 B =(1+R ) B +G τ C C G τ R B C = a R +a W W ρ W =(1+R ) B +(1+R +δ ) (1 ρ) L B L δ B = λ B + μ (W C λ B )
More informationAR(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第85 回日本感染症学会総会学術集会後抄録(III)
β β α α α µ µ µ µ α α α α γ αβ α γ α α γ α γ µ µ β β β β β β β β β µ β α µ µ µ β β µ µ µ µ µ µ γ γ γ γ γ γ µ α β γ β β µ µ µ µ µ β β µ β β µ α β β µ µµ β µ µ µ µ µ µ λ µ µ β µ µ µ µ µ µ µ µ
More information診療ガイドライン外来編2014(A4)/FUJGG2014‐01(大扉)
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
More information日本糖尿病学会誌第58巻第1号
α β β β β β β α α β α β α l l α l μ l β l α β β Wfs1 β β l l l l μ l l μ μ l μ l Δ l μ μ l μ l l ll 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 information136 pp p µl µl µl
135 2006 PCB C 12 H 10-n Cl n n 1 10 CAS No. 42 PCB: 53469-21-9, 54 PCB: 11097-69-1 0.01 mg/m 3 PCB PCB 25 µg/l 136 pp p µl µl µl 137 1 γ 138 1 γ γ γ µl µl µl µl µl µl µl l µl µl µl µl µl l 139 µl µl µl
More information日本糖尿病学会誌第58巻第3号
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 information330
330 331 332 333 334 t t P 335 t R t t i R +(P P ) P =i t P = R + P 1+i t 336 uc R=uc P 337 338 339 340 341 342 343 π π β τ τ (1+π ) (1 βτ )(1 τ ) (1+π ) (1 βτ ) (1 τ ) (1+π ) (1 τ ) (1 τ ) 344 (1 βτ )(1
More information80 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 informationuntitled
Y = Y () x i c C = i + c = ( x ) x π (x) π ( x ) = Y ( ){1 + ( x )}( 1 x ) Y ( )(1 + C ) ( 1 x) x π ( x) = 0 = ( x ) R R R R Y = (Y ) CS () CS ( ) = Y ( ) 0 ( Y ) dy Y ( ) A() * S( π ), S( CS) S( π ) =
More information日本糖尿病学会誌第58巻第2号
β γ Δ Δ β β β 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第89回日本感染症学会学術講演会後抄録(I)
! ! ! β !!!!!!!!!!! !!! !!! μ! μ! !!! β! β !! β! β β μ! μ! μ! μ! β β β β β β μ! μ! μ!! β ! β ! ! β β ! !! ! !!! ! ! ! β! !!!!! !! !!!!!!!!! μ! β !!!! β β! !!!!!!!!! !! β β β β β β β β !!
More information(interval estimation) 3 (confidence coefficient) µ σ/sqrt(n) 4 P ( (X - µ) / (σ sqrt N < a) = α a α X α µ a σ sqrt N X µ a σ sqrt N 2
7 2 1 (interval estimation) 3 (confidence coefficient) µ σ/sqrt(n) 4 P ( (X - µ) / (σ sqrt N < a) = α a α X α µ a σ sqrt N X µ a σ sqrt N 2 (confidence interval) 5 X a σ sqrt N µ X a σ sqrt N - 6 P ( X
More informationStata 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 information24.15章.微分方程式
m d y dt = F m d y = mg dt V y = dy dt d y dt = d dy dt dt = dv y dt dv y dt = g dv y dt = g dt dt dv y = g dt V y ( t) = gt + C V y ( ) = V y ( ) = C = V y t ( ) = gt V y ( t) = dy dt = gt dy = g t dt
More information第88回日本感染症学会学術講演会後抄録(III)
!!!! β! !!μ μ!!μ μ!!μ! !!!! α!!! γδ Φ Φ Φ Φ! Φ Φ Φ Φ Φ! α!! ! α β α α β α α α α α α α α β α α β! β β μ!!!! !!μ !μ!μ!!μ!!!!! !!!!!!!!!! !!!!!!μ! !!μ!!!μ!!!!!! γ γ γ γ γ γ! !!!!!! β!!!! β !!!!!! β! !!!!μ!!!!!!
More informationweb04.dvi
4 MATLAB 1 visualization MATLAB 2 Octave gnuplot Octave copyright c 2004 Tatsuya Kitamura / All rights reserved. 35 4 4.1 1 1 y =2x x 5 5 x y plot 4.1 Figure No. 1 figure window >> x=-5:5;ψ >> y=2*x;ψ
More informationMicrosoft Word - ■3中表紙(2006版).doc
18 Annual Report on Research Activity by Faculty of Medicine, University of the Ryukyus 2006 FACULTY OF MEDICINE UNIVERSITY OF THE RYUKYUS α αγ α β α βγ β α β α β β β γ κα κ κ βγ ε α γδ β
More information1 913 10301200 A B C D E F G H J K L M 1A1030 10 : 45 1A1045 11 : 00 1A1100 11 : 15 1A1115 11 : 30 1A1130 11 : 45 1A1145 12 : 00 1B1030 1B1045 1C1030
1 913 9001030 A B C D E F G H J K L M 9:00 1A0900 9:15 1A0915 9:30 1A0930 9:45 1A0945 10 : 00 1A1000 10 : 15 1B0900 1B0915 1B0930 1B0945 1B1000 1C0900 1C0915 1D0915 1C0930 1C0945 1C1000 1D0930 1D0945 1D1000
More informationron04-02/ky768450316800035946
β α β α β β β α α α Bugula neritina α β β β γ γ γ γ β β γ β β β β γ β β β β β β β β! ! β β β β μ β μ β β β! β β β β β μ! μ! μ! β β α!! β γ β β β β!! β β β β β β! β! β β β!! β β β β β β β β β β β β!
More informationQ.5-1 Ans.
5 Q.5-1 Q.5-2 Q.5-3 Q.5-4 Q.5-5 Q.5-6 Q.5-7 Q.5-8 Q.5-9 Q.5-10 Q.5-11 Q.5-12 Q.5-13 Q.5-14 Q.5-15 Q.5-1 Ans. Q.4-5 1 Q.5-3 Check Q.5-2 200 203 197 199 Q.5-2 Ans. Check Q.5-3 Ans. Q.5-2 12 Q.6-4 69 56
More information確率論と統計学の資料
5 June 015 ii........................ 1 1 1.1...................... 1 1........................... 3 1.3... 4 6.1........................... 6................... 7 ii ii.3.................. 8.4..........................
More information4
4 5 6 7 + 8 = ++ 9 + + + + ++ 10 + + 11 12 WS LC VA L WS = LC VA = LC L L VA = LC L VA L 13 i LC VA WS WS = LC = VA LC VA VA = VA α WS α = VA VA i WS = LC VA i t t+1 14 WS = α WS + WS α WS = WS WS WS =
More information1 1 2 GDP 3 1 GDP 2 GDP 3 GDP GDP GDP 4 GDP GDP GDP 1 GDP 2 CPI 2
Macroeconomics/ Olivier Blanchard, 1996 1 2 1........... 2 2............. 2 2 3 3............... 3 4...... 4 5.............. 4 6 IS-LM................. 5 3 6 7. 6 8....... 7 9........... 8 10..... 8 8
More information平成19年度
1 2 3 4 H 3 H CC N + 3 O H 3 C O CO CH 3 CH O CO O CH2 CH 3 P O O 5 H H H CHOH H H H N + CHOH CHOH N + CH CH COO- CHOH CH CHOH 6 1) 7 2 ) 8 3 ) 4 ) 9 10 11 12 13 14 15 16 17 18 19 20 A A 0 21 ) exp( )
More informationdicutil1_5_2.book
Kabayaki for Windows Version 1.5.2 ...1...1 1...3...3 2...5...5...5...7...7 3...9...9...9...10...10...11...12 1 2 Kabayaki ( ) Kabayaki Kabayaki ( ) Kabayaki Kabayaki Kabayaki 1 2 1 Kabayaki ( ) ( ) CSV
More information1 GDP Q GDP (a) (b) (c) (d) (e) (f) A (b) (e) (f) Q GDP A GDP GDP = Q 1990 GNP GDP 4095 3004 1091 GNP A Q 1995 7 A 2 2
/, 2001 1 GDP................................... 2 2.......................... 2 3.................................... 4 4........................................ 5 5.....................................
More information1 1 ( ) ( 1.1 1.1.1 60% mm 100 100 60 60% 1.1.2 A B A B A 1
1 21 10 5 1 E-mail: qliu@res.otaru-uc.ac.jp 1 1 ( ) ( 1.1 1.1.1 60% mm 100 100 60 60% 1.1.2 A B A B A 1 B 1.1.3 boy W ID 1 2 3 DI DII DIII OL OL 1.1.4 2 1.1.5 1.1.6 1.1.7 1.1.8 1.2 1.2.1 1. 2. 3 1.2.2
More information46 Y 5.1.1 Y Y Y 3.1 R Y Figures 5-1 5-3 3.2mm Nylon Glass Y (X > X ) X Y X Figure 5-1 X min Y Y d Figure 5-3 X =X min Y X =10 Y Y Y 5.1.2 Y Figure 5-
45 5 5.1 Y 3.2 Eq. (3) 1 R [s -1 ] ideal [s -1 ] Y [-] Y [-] ideal * [-] S [-] 3 R * ( ω S ) = ω Y = ω 3-1a ideal ideal X X R X R (X > X ) ideal * X S Eq. (3-1a) ( X X ) = Y ( X ) R > > θ ω ideal X θ =
More information10
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 information204 / CHEMISTRY & CHEMICAL INDUSTRY Vol.69-1 January 2016 047
9 π 046 Vol.69-1 January 2016 204 / CHEMISTRY & CHEMICAL INDUSTRY Vol.69-1 January 2016 047 β γ α / α / 048 Vol.69-1 January 2016 π π π / CHEMISTRY & CHEMICAL INDUSTRY Vol.69-1 January 2016 049 β 050 Vol.69-1
More informationJune 2016 i (statistics) F Excel Numbers, OpenOffice/LibreOffice Calc ii *1 VAR STDEV 1 SPSS SAS R *2 R R R R *1 Excel, Numbers, Microsoft Office, Apple iwork, *2 R GNU GNU R iii URL http://ruby.kyoto-wu.ac.jp/statistics/training/
More information需 要 予 測 のための 統 計 モテ ルの 研 究 異 常 値 検 知 のための 基 本 的 モテ ルの 考 察 平 成 26 年 3 月 東 京 大 学 大 学 院 情 報 理 工 学 系 研 究 科 教 授 博 士 課 程 修 士 課 程 竹 村 彰 通 小 川 光 紀 笹 井 健 行 特 定 非 営 利 活 動 法 人 ヒ ュー コミュニケーションス 副 理 事 長 小 松 秀 樹 主 任
More information36
36 37 38 P r R P 39 (1+r ) P =R+P g P r g P = R r g r g == == 40 41 42 τ R P = r g+τ 43 τ (1+r ) P τ ( P P ) = R+P τ ( P P ) n P P r P P g P 44 R τ P P = (1 τ )(r g) (1 τ )P R τ 45 R R σ u R= R +u u~ (0,σ
More information: (EQS) /EQUATIONS V1 = 30*V F1 + E1; V2 = 25*V *F1 + E2; V3 = 16*V *F1 + E3; V4 = 10*V F2 + E4; V5 = 19*V99
218 6 219 6.11: (EQS) /EQUATIONS V1 = 30*V999 + 1F1 + E1; V2 = 25*V999 +.54*F1 + E2; V3 = 16*V999 + 1.46*F1 + E3; V4 = 10*V999 + 1F2 + E4; V5 = 19*V999 + 1.29*F2 + E5; V6 = 17*V999 + 2.22*F2 + E6; CALIS.
More informationCopyrght 7 Mzuho-DL Fnancal Technology Co., Ltd. All rghts reserved.
766 Copyrght 7 Mzuho-DL Fnancal Technology Co., Ltd. All rghts reserved. Copyrght 7 Mzuho-DL Fnancal Technology Co., Ltd. All rghts reserved. 3 Copyrght 7 Mzuho-DL Fnancal Technology Co., Ltd. All rghts
More informationL A TEX ver L A TEX LATEX 1.1 L A TEX L A TEX tex 1.1 1) latex mkdir latex 2) latex sample1 sample2 mkdir latex/sample1 mkdir latex/sampl
L A TEX ver.2004.11.18 1 L A TEX LATEX 1.1 L A TEX L A TEX tex 1.1 1) latex mkdir latex 2) latex sample1 sample2 mkdir latex/sample1 mkdir latex/sample2 3) /staff/kaede work/www/math/takase sample1.tex
More informationチュートリアル:ノンパラメトリックベイズ
{ x,x, L, xn} 2 p( θ, θ, θ, θ, θ, } { 2 3 4 5 θ6 p( p( { x,x, L, N} 2 x { θ, θ2, θ3, θ4, θ5, θ6} K n p( θ θ n N n θ x N + { x,x, L, N} 2 x { θ, θ2, θ3, θ4, θ5, θ6} log p( 6 n logθ F 6 log p( + λ θ F θ
More information日本糖尿病学会誌第58巻第7号
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基礎から学ぶトラヒック理論 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 初版 1 刷発行時のものです.
基礎から学ぶトラヒック理論 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/085221 このサンプルページの内容は, 初版 1 刷発行時のものです. i +α 3 1 2 4 5 1 2 ii 3 4 5 6 7 8 9 9.3 2014 6 iii 1 1 2 5 2.1 5 2.2 7
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日本分子第4巻2号_10ポスター発表.indd
JSMI Report 62 63 JSMI Report γ JSMI Report 64 β α 65 JSMI Report JSMI Report 66 67 JSMI Report JSMI Report 68 69 JSMI Report JSMI Report 70 71 JSMI Report JSMI Report 72 73 JSMI Report JSMI Report 74
More information0.,,., m Euclid m m. 2.., M., M R 2 ψ. ψ,, R 2 M.,, (x 1 (),, x m ()) R m. 2 M, R f. M (x 1,, x m ), f (x 1,, x m ) f(x 1,, x m ). f ( ). x i : M R.,,
2012 10 13 1,,,.,,.,.,,. 2?.,,. 1,, 1. (θ, φ), θ, φ (0, π),, (0, 2π). 1 0.,,., m Euclid m m. 2.., M., M R 2 ψ. ψ,, R 2 M.,, (x 1 (),, x m ()) R m. 2 M, R f. M (x 1,, x m ), f (x 1,, x m ) f(x 1,, x m ).
More information<4D6963726F736F667420576F7264202D204850835483938376838B8379815B83578B6594BB2D834A836F815B82D082C88C60202E646F63>
例 題 で 学 ぶ Excel 統 計 入 門 第 2 版 サンプルページ この 本 の 定 価 判 型 などは, 以 下 の URL からご 覧 いただけます. http://www.morikita.co.jp/books/mid/084302 このサンプルページの 内 容 は, 第 2 版 発 行 当 時 のものです. i 2 9 2 Web 2 Excel Excel Excel 11 Excel
More informationJMP V4 による生存時間分析
V4 1 SAS 2000.11.18 4 ( ) (Survival Time) 1 (Event) Start of Study Start of Observation Died Died Died Lost End Time Censor Died Died Censor Died Time Start of Study End Start of Observation Censor
More information5989_4840JAJP.qxd
Agilent Application Note 1287-11 2 3 4 5 Zc Z T 1+ G 1 e - γ 1+ G 2 G i G 1 G 2 0 0 G2 G 1 G T 1+ G 2 e - γ 1+ G 1 a b [ T XI ] [ T L ] [ T XO ] [ G L ] Zc Zr ZT Zr Γ1 = Γ2 = Γ1ΓT = (1.1) Zc+ Zr ZT + Zr
More information( ) 1.1 Polychoric Correlation Polyserial Correlation Graded Response Model Partial Credit Model Tetrachoric Correlation ( ) 2 x y x y s r 1 x 2
1 (,2007) SPSSver8 1997 (2002) 1. 2. polychoric correlation coefficient (polyserial correlation coefficient) 3. (1999) M-plus R 1 ( ) 1.1 Polychoric Correlation Polyserial Correlation Graded Response Model
More information2 H23 BioS (i) data d1; input group patno t sex censor; cards;
H BioS (i) data d1; input group patno t sex censor; cards; 0 1 0 0 0 0 1 0 1 1 0 4 4 0 1 0 5 5 1 1 0 6 5 1 1 0 7 10 1 0 0 8 15 0 1 0 9 15 0 1 0 10 4 1 0 0 11 4 1 0 1 1 5 1 0 1 1 7 0 1 1 14 8 1 0 1 15 8
More informationuntitled
Global Quantitative Research / -2- -3- -4- -5- 35 35 SPC SPC REIT REIT -6- -7- -8- -9- -10- -11- -12- -13- -14- -15- -16- -17- 100m$110-18- Global Quantitative Research -19- -20- -21- -22- -23- -24- -25-
More information203 0 5% 8% 204 4 204 205 0 0% 8% Peroi202 203 Iwaa20 205 0 0% 0%. 4 2. 2
三 田 祭 論 文 消 費 税 増 税 時 の 政 策 オプション 3 本 の 矢 を 折 らないために 慶 應 義 塾 大 学 廣 瀬 康 生 研 究 会 荻 野 秀 明 川 邉 美 帆 財 津 薫 平 関 谷 裕 鏡 高 野 隼 一 橋 本 壮 広 橋 本 龍 一 郎 藤 村 和 輝 保 里 俊 介 203 年 月 203 0 5% 8% 204 4 204 205 0 0% 8% Peroi202
More informationECCS. ECCS,. ( 2. Mac Do-file Editor. Mac Do-file Editor Windows Do-file Editor Top Do-file e
1 1 2015 4 6 1. ECCS. ECCS,. (https://ras.ecc.u-tokyo.ac.jp/guacamole/) 2. Mac Do-file Editor. Mac Do-file Editor Windows Do-file Editor Top Do-file editor, Do View Do-file Editor Execute(do). 3. Mac System
More information<4D6963726F736F667420576F7264202D204850835483938376838B8379815B83578B6594BB2D834A836F815B82D082C88C60202E646F63>
誤 り 訂 正 技 術 の 基 礎 サンプルページ この 本 の 定 価 判 型 などは, 以 下 の URL からご 覧 いただけます http://wwwmorikitacojp/books/mid/081731 このサンプルページの 内 容 は, 第 1 版 発 行 時 のものです http://wwwmorikitacojp/support/ e mail editor@morikitacojp
More information: , 2.0, 3.0, 2.0, (%) ( 2.
2017 1 2 1.1...................................... 2 1.2......................................... 4 1.3........................................... 10 1.4................................. 14 1.5..........................................
More information木オートマトン•トランスデューサによる 自然言語処理
木オートマトン トランスデューサによる 自然言語処理 林 克彦 NTTコミュニケーション科学基礎研究所 hayashi.katsuhiko@lab.ntt.co.jp n I T 1 T 2 I T 1 Pro j(i T 1 T 2 ) (Σ,rk) Σ rk : Σ N {0} nσ (n) rk(σ) = n σ Σ n Σ (n) Σ (n)(σ,rk)σ Σ T Σ (A) A
More informationE B m e ( ) γma = F = e E + v B a m = 0.5MeV γ = E e m =957 E e GeV v β = v SPring-8 γ β γ E e [GeV] [ ] NewSUBARU.0 957 0.999999869 SPring-8 8.0 5656
SPring-8 PF( ) ( ) UVSOR( HiSOR( SPring-8.. 3. 4. 5. 6. 7. E B m e ( ) γma = F = e E + v B a m = 0.5MeV γ = E e m =957 E e GeV v β = v SPring-8 γ β γ E e [GeV] [ ] NewSUBARU.0 957 0.999999869 SPring-8
More information5 36 5................................................... 36 5................................................... 36 5.3..............................
9 8 3............................................. 3.......................................... 4.3............................................ 4 5 3 6 3..................................................
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
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卒業論文
Y = ax 1 b1 X 2 b2...x k bk e u InY = Ina + b 1 InX 1 + b 2 InX 2 +...+ b k InX k + u X 1 Y b = ab 1 X 1 1 b 1 X 2 2...X bk k e u = b 1 (ax b1 1 X b2 2...X bk k e u ) / X 1 = b 1 Y / X 1 X 1 X 1 q YX1
More informationuntitled
MCMC 2004 23 1 I. MCMC 1. 2. 3. 4. MH 5. 6. MCMC 2 II. 1. 2. 3. 4. 5. 3 I. MCMC 1. 2. 3. 4. MH 5. 4 1. MCMC 5 2. A P (A) : P (A)=0.02 A B A B Pr B A) Pr B A c Pr B A)=0.8, Pr B A c =0.1 6 B A 7 8 A, :
More information