Starting preprocessing of the model file ... Substitution of endo lags >= 2: added 2 auxiliary variables and equations. Found 53 equation(s). Evaluating expressions...done Computing static model derivatives (order 1). Computing dynamic model derivatives (order 2). Processing outputs ... done Preprocessing completed. You did not declare endogenous variables after the estimation/calib_smoother command. Smoothed variables will be computed for the 51 endogenous variables of your model, this can take a long time .... Choose one of the following options: [1] Consider all the endogenous variables. [2] Consider all the observed endogenous variables. [3] Stop Dynare and change the mod file. options [default is 1] = 1 Initial value of the log posterior (or likelihood): -2874.7786 ========================================================== Change in the posterior covariance matrix = 10. Change in the posterior mean = 0.55739. Mode improvement = 1882.5151 New value of jscale = 0.00038099 ========================================================== ========================================================== Change in the posterior covariance matrix = 0.068778. Change in the posterior mean = 2.8057. Mode improvement = 656.3488 New value of jscale = 0.41249 ========================================================== ========================================================== Change in the posterior covariance matrix = 0.062431. Change in the posterior mean = 1.7122. Mode improvement = 480.4617 New value of jscale = 0.34248 ========================================================== Optimal value of the scale parameter = 0.34248 Final value of minus the log posterior (or likelihood):1168.150585 MODE CHECK Fval obtained by the minimization routine (minus the posterior/likelihood)): 1168.150585 [警告: 行列が特異なため、正確に処理できません。] [> evaluate_steady_state (行 261) 内 resol (行 61) 内 compute_decision_rules (行 35) 内 dynare_resolve (行 69) 内 dsge_var_likelihood (行 156) 内 mode_check (行 161) 内 dynare_estimation_1 (行 322) 内 dynare_estimation (行 105) 内 jpmodel5.driver (行 842) 内 dynare (行 293) 内 ] [警告: 行列は、特異行列に近いか、正しくスケーリングされていませ ん。結果は不正確になる可能性があります。RCOND = 5.239063e-21。] [> evaluate_steady_state (行 261) 内 resol (行 61) 内 compute_decision_rules (行 35) 内 dynare_resolve (行 69) 内 dsge_var_likelihood (行 156) 内 mode_check (行 161) 内 dynare_estimation_1 (行 322) 内 dynare_estimation (行 105) 内 jpmodel5.driver (行 842) 内 dynare (行 293) 内 ] [警告: 行列は、特異行列に近いか、正しくスケーリングされていませ ん。結果は不正確になる可能性があります。RCOND = 7.727338e-22。] [> evaluate_steady_state (行 261) 内 resol (行 61) 内 compute_decision_rules (行 35) 内 dynare_resolve (行 69) 内 dsge_var_likelihood (行 156) 内 mode_check (行 161) 内 dynare_estimation_1 (行 322) 内 dynare_estimation (行 105) 内 jpmodel5.driver (行 842) 内 dynare (行 293) 内 ] [警告: 行列は、特異行列に近いか、正しくスケーリングされていませ ん。結果は不正確になる可能性があります。RCOND = 4.696307e-22。] [> evaluate_steady_state (行 261) 内 resol (行 61) 内 compute_decision_rules (行 35) 内 dynare_resolve (行 69) 内 dsge_var_likelihood (行 156) 内 mode_check (行 161) 内 dynare_estimation_1 (行 322) 内 dynare_estimation (行 105) 内 jpmodel5.driver (行 842) 内 dynare (行 293) 内 ] [警告: 行列が特異なため、正確に処理できません。] [> evaluate_steady_state (行 261) 内 resol (行 61) 内 compute_decision_rules (行 35) 内 dynare_resolve (行 69) 内 dsge_var_likelihood (行 156) 内 mode_check (行 161) 内 dynare_estimation_1 (行 322) 内 dynare_estimation (行 105) 内 jpmodel5.driver (行 842) 内 dynare (行 293) 内 ] [警告: 行列が特異なため、正確に処理できません。] [> evaluate_steady_state (行 261) 内 resol (行 61) 内 compute_decision_rules (行 35) 内 dynare_resolve (行 69) 内 dsge_var_likelihood (行 156) 内 mode_check (行 161) 内 dynare_estimation_1 (行 322) 内 dynare_estimation (行 105) 内 jpmodel5.driver (行 842) 内 dynare (行 293) 内 ] [警告: 行列が特異なため、正確に処理できません。] [> evaluate_steady_state (行 261) 内 resol (行 61) 内 compute_decision_rules (行 35) 内 dynare_resolve (行 69) 内 dsge_var_likelihood (行 156) 内 mode_check (行 161) 内 dynare_estimation_1 (行 322) 内 dynare_estimation (行 105) 内 jpmodel5.driver (行 842) 内 dynare (行 293) 内 ] [警告: 行列が特異なため、正確に処理できません。] [> evaluate_steady_state (行 261) 内 resol (行 61) 内 compute_decision_rules (行 35) 内 dynare_resolve (行 69) 内 dsge_var_likelihood (行 156) 内 mode_check (行 161) 内 dynare_estimation_1 (行 322) 内 dynare_estimation (行 105) 内 jpmodel5.driver (行 842) 内 dynare (行 293) 内 ] RESULTS FROM POSTERIOR ESTIMATION parameters prior mean mode s.d. prior pstdev dsge_prior_weight 1.000 0.9867 0.0083 unif 0.5774 v_gm 0.000 1.1063 0.0322 norm 1.5000 v_gp 0.000 0.4141 0.0238 norm 1.5000 v 0.250 0.1678 0.0159 norm 1.0000 omega 0.250 0.1737 0.0042 beta 0.1000 sigma 1.000 1.2520 0.0363 gamm 0.3750 theta 0.700 0.4970 0.0089 beta 0.1500 chi 2.000 2.6843 0.0446 gamm 0.7500 zetainv 4.000 4.3963 0.0452 gamm 1.5000 mu 1.000 1.0480 0.0908 gamm 1.0000 phi 0.075 0.0779 0.0005 gamm 0.0125 gammaw 0.500 0.4326 0.0104 beta 0.2500 xiw 0.375 0.5426 0.0042 beta 0.1000 gammap 0.500 0.9705 0.0155 beta 0.2500 xip 0.375 0.5602 0.0046 beta 0.1000 lambdap 0.150 0.2427 0.0025 gamm 0.0500 zstar 0.190 0.1909 0.0048 gamm 0.0500 lstar 0.000 -0.0245 0.0023 norm 0.0500 pistar 0.175 0.1267 0.0010 gamm 0.0500 rstar 0.498 0.5201 0.0021 gamm 0.0500 phi_r 0.800 0.7449 0.0031 beta 0.1000 phi_r_pi 0.500 0.4592 0.0126 norm 0.2000 phi_r_y 0.125 0.0823 0.0008 gamm 0.0500 phi_gm 0.125 0.7445 0.0017 beta 0.0500 phi_gm_y 0.000 1.4087 0.0210 norm 0.5000 phi_gm_b 0.000 1.0212 0.0259 norm 0.5000 phi_gp 0.800 0.8391 0.0033 beta 0.1000 phi_gp_y 0.000 0.2792 0.0168 norm 0.5000 phi_gp_b 0.000 -0.0598 0.0088 norm 0.5000 phi_gi 0.800 0.7668 0.0024 beta 0.1000 phi_gi_y 0.000 0.2803 0.0100 norm 0.5000 phi_gi_b 0.000 -0.1518 0.0149 norm 0.5000 phi_T 0.800 0.9178 0.0026 beta 0.1000 phi_T_y 0.000 -0.2708 0.0127 norm 0.5000 phi_T_b 0.000 -0.0680 0.0187 norm 0.5000 phi_tc 0.800 0.8605 0.0023 beta 0.1000 phi_tc_y 0.000 -0.0034 0.0380 norm 0.5000 phi_tc_b 0.000 -0.0025 0.0247 norm 0.5000 phi_tw 0.800 0.8002 0.0019 beta 0.1000 phi_tw_y 0.000 -0.1736 0.0218 norm 0.5000 phi_tw_b 0.000 0.0597 0.0053 norm 0.5000 phi_tk 0.800 0.8212 0.0028 beta 0.1000 phi_tk_y 0.000 0.0345 0.0099 norm 0.5000 phi_tk_b 0.000 0.0383 0.0136 norm 0.5000 rho_z 0.500 0.2234 0.0061 beta 0.2000 rho_b 0.500 0.2939 0.0096 beta 0.2000 rho_i 0.500 0.4579 0.0039 beta 0.2000 rho_w 0.500 0.6570 0.0040 beta 0.2000 rho_p 0.500 0.5564 0.0020 beta 0.2000 rho_x 0.500 0.8254 0.0056 beta 0.2000 rho_r 0.500 0.4978 0.0064 beta 0.2000 rho_gm 0.500 0.5579 0.0090 beta 0.2000 rho_gp 0.500 0.0576 0.0060 beta 0.2000 rho_gi 0.500 0.4397 0.0029 beta 0.2000 rho_T 0.500 0.1034 0.0069 beta 0.2000 standard deviation of shocks prior mean mode s.d. prior pstdev e_z 0.100 0.5437 0.0669 invg Inf e_b 0.100 1.4655 0.0973 invg Inf e_i 0.100 1.9719 0.1785 invg Inf e_w 0.100 0.5048 0.0244 invg Inf e_p 0.100 0.3061 0.0418 invg Inf e_x 0.100 4.0029 0.2757 invg Inf e_r 0.100 0.2062 0.0192 invg Inf e_gm 0.100 0.4992 0.0310 invg Inf e_gp 0.100 0.7118 0.0419 invg Inf e_gi 0.100 1.6939 0.0547 invg Inf e_T 0.100 0.0475 0.0363 invg Inf e_tc 0.100 0.0078 0.0007 invg Inf e_tw 0.100 0.4474 0.0011 invg Inf e_tk 0.100 0.0091 0.0008 invg Inf Log data density [Laplace approximation] is -1601.676070. Estimation::mcmc: Multiple chains mode. Estimation::mcmc: Old metropolis.log file successfully erased! Estimation::mcmc: Creation of a new metropolis.log file. Estimation::mcmc: Searching for initial values... Estimation::mcmc: Initial values found! Estimation::mcmc: Write details about the MCMC... Ok! Estimation::mcmc: Details about the MCMC are available in jpmodel5/metropolis/jpmodel5_mh_history_0.mat Estimation::mcmc: Number of mh files: 2 per block. Estimation::mcmc: Total number of generated files: 4. Estimation::mcmc: Total number of iterations: 200000. Estimation::mcmc: Current acceptance ratio per chain: Chain 1: 41.0335% Chain 2: 37.1695% Estimation::mcmc: Total number of MH draws per chain: 200000. Estimation::mcmc: Total number of generated MH files: 2. Estimation::mcmc: I'll use mh-files 1 to 2. Estimation::mcmc: In MH-file number 1 I'll start at line 60001. Estimation::mcmc: Finally I keep 140000 draws per chain. MCMC Inefficiency factors per block Parameter Block 1 Block 2 SE_e_z 658.169 639.700 SE_e_b 505.658 565.177 SE_e_i 470.724 620.827 SE_e_w 564.861 423.889 SE_e_p 195.697 259.421 SE_e_x 599.790 610.127 SE_e_r 135.408 180.089 SE_e_gm 368.410 383.779 SE_e_gp 371.486 408.434 SE_e_gi 647.870 631.643 SE_e_T 666.373 535.502 SE_e_tc 69.623 69.465 SE_e_tw 736.467 740.024 SE_e_tk 67.272 70.246 dsge_prior_weight 716.526 733.405 v_gm 722.972 725.805 v_gp 733.582 738.244 v 539.570 550.173 omega 734.608 708.307 sigma 709.750 712.156 theta 709.900 690.669 chi 738.189 724.542 zetainv 705.253 712.366 mu 638.562 698.666 phi 711.625 732.711 gammaw 715.592 694.481 xiw 734.417 717.501 gammap 714.430 548.124 xip 717.841 719.004 lambdap 738.366 697.009 zstar 704.842 680.990 lstar 690.328 696.593 pistar 730.653 723.263 rstar 693.708 684.423 phi_r 735.175 696.372 phi_r_pi 717.386 731.037 phi_r_y 729.118 721.860 phi_gm 731.807 739.592 phi_gm_y 717.904 730.761 phi_gm_b 734.028 669.719 phi_gp 728.608 719.974 phi_gp_y 725.866 687.360 phi_gp_b 743.091 739.894 phi_gi 709.957 719.839 phi_gi_y 706.419 722.722 phi_gi_b 719.132 726.417 phi_T 734.136 717.433 phi_T_y 726.922 721.440 phi_T_b 701.688 648.803 phi_tc 715.955 698.505 phi_tc_y 611.100 244.922 phi_tc_b 682.841 286.077 phi_tw 720.765 705.876 phi_tw_y 680.100 721.863 phi_tw_b 722.494 694.189 phi_tk 740.121 718.644 phi_tk_y 518.177 439.538 phi_tk_b 576.529 343.513 rho_z 715.631 703.641 rho_b 658.470 720.789 rho_i 699.857 685.580 rho_w 731.975 690.318 rho_p 742.539 736.552 rho_x 722.000 689.502 rho_r 718.687 660.529 rho_gm 722.655 739.500 rho_gp 727.308 690.830 rho_gi 738.043 731.700 rho_T 676.747 721.867 Estimation::mcmc::diagnostics: Univariate convergence diagnostic, Brooks and Gelman (1998): Parameter 1... Done! Parameter 2... Done! Parameter 3... Done! Parameter 4... Done! Parameter 5... Done! Parameter 6... Done! Parameter 7... Done! Parameter 8... Done! Parameter 9... Done! Parameter 10... Done! Parameter 11... Done! Parameter 12... Done! Parameter 13... Done! Parameter 14... Done! Parameter 15... Done! Parameter 16... Done! Parameter 17... Done! Parameter 18... Done! Parameter 19... Done! Parameter 20... Done! Parameter 21... Done! Parameter 22... Done! Parameter 23... Done! Parameter 24... Done! Parameter 25... Done! Parameter 26... Done! Parameter 27... Done! Parameter 28... Done! Parameter 29... Done! Parameter 30... Done! Parameter 31... Done! Parameter 32... Done! Parameter 33... Done! Parameter 34... Done! Parameter 35... Done! Parameter 36... Done! Parameter 37... Done! Parameter 38... Done! Parameter 39... Done! Parameter 40... Done! Parameter 41... Done! Parameter 42... Done! Parameter 43... Done! Parameter 44... Done! Parameter 45... Done! Parameter 46... Done! Parameter 47... Done! Parameter 48... Done! Parameter 49... Done! Parameter 50... Done! Parameter 51... Done! Parameter 52... Done! Parameter 53... Done! Parameter 54... Done! Parameter 55... Done! Parameter 56... Done! Parameter 57... Done! Parameter 58... Done! Parameter 59... Done! Parameter 60... Done! Parameter 61... Done! Parameter 62... Done! Parameter 63... Done! Parameter 64... Done! Parameter 65... Done! Parameter 66... Done! Parameter 67... Done! Parameter 68... Done! Parameter 69... Done! Estimation::marginal density: I'm computing the posterior mean and covariance... Done! Estimation::marginal density: I'm computing the posterior log marginal density (modified harmonic mean)... Done! ESTIMATION RESULTS Log data density is -1421.256629. parameters prior mean post. mean 90% HPD interval prior pstdev dsge_prior_weight 1.000 1.0809 0.9486 1.1986 unif 0.5774 v_gm 0.000 1.0836 0.5501 1.6952 norm 1.5000 v_gp 0.000 0.5640 0.1424 0.9720 norm 1.5000 v 0.250 0.2302 0.1837 0.2755 norm 1.0000 omega 0.250 0.1473 0.1173 0.1756 beta 0.1000 sigma 1.000 1.4648 1.3227 1.5958 gamm 0.3750 theta 0.700 0.4218 0.3404 0.5009 beta 0.1500 chi 2.000 2.2444 1.8512 2.5427 gamm 0.7500 zetainv 4.000 4.0956 3.6185 4.7292 gamm 1.5000 mu 1.000 1.1135 0.9092 1.3253 gamm 1.0000 phi 0.075 0.0725 0.0694 0.0763 gamm 0.0125 gammaw 0.500 0.3551 0.2585 0.4317 beta 0.2500 xiw 0.375 0.6080 0.5641 0.6600 beta 0.1000 gammap 0.500 0.9561 0.8928 0.9993 beta 0.2500 xip 0.375 0.6098 0.5767 0.6438 beta 0.1000 lambdap 0.150 0.2791 0.2349 0.3285 gamm 0.0500 zstar 0.190 0.1780 0.1637 0.1922 gamm 0.0500 lstar 0.000 -0.0228 -0.0351 -0.0106 norm 0.0500 pistar 0.175 0.1432 0.1329 0.1543 gamm 0.0500 rstar 0.498 0.5391 0.5229 0.5550 gamm 0.0500 phi_r 0.800 0.7650 0.7462 0.7968 beta 0.1000 phi_r_pi 0.500 0.4988 0.4294 0.5752 norm 0.2000 phi_r_y 0.125 0.0858 0.0686 0.1024 gamm 0.0500 phi_gm 0.125 0.7693 0.7485 0.7859 beta 0.0500 phi_gm_y 0.000 1.2123 0.9964 1.4297 norm 0.5000 phi_gm_b 0.000 0.7509 0.5076 0.9669 norm 0.5000 phi_gp 0.800 0.8509 0.8243 0.8747 beta 0.1000 phi_gp_y 0.000 -0.0370 -0.2569 0.1600 norm 0.5000 phi_gp_b 0.000 0.2389 0.0509 0.4622 norm 0.5000 phi_gi 0.800 0.7917 0.7700 0.8122 beta 0.1000 phi_gi_y 0.000 0.3006 0.1972 0.4173 norm 0.5000 phi_gi_b 0.000 -0.2766 -0.4559 -0.0893 norm 0.5000 phi_T 0.800 0.8873 0.8658 0.9065 beta 0.1000 phi_T_y 0.000 -0.0859 -0.1954 0.0252 norm 0.5000 phi_T_b 0.000 0.1097 -0.0804 0.2825 norm 0.5000 phi_tc 0.800 0.9413 0.8881 0.9954 beta 0.1000 phi_tc_y 0.000 0.0106 -0.0299 0.0650 norm 0.5000 phi_tc_b 0.000 -0.0442 -0.1193 0.0120 norm 0.5000 phi_tw 0.800 0.8162 0.7949 0.8369 beta 0.1000 phi_tw_y 0.000 -0.1826 -0.2665 -0.0823 norm 0.5000 phi_tw_b 0.000 0.0334 -0.0504 0.1293 norm 0.5000 phi_tk 0.800 0.8747 0.8298 0.9327 beta 0.1000 phi_tk_y 0.000 0.0247 -0.0017 0.0490 norm 0.5000 phi_tk_b 0.000 0.0412 0.0202 0.0620 norm 0.5000 rho_z 0.500 0.1747 0.1243 0.2205 beta 0.2000 rho_b 0.500 0.3251 0.2718 0.3666 beta 0.2000 rho_i 0.500 0.4034 0.3524 0.4463 beta 0.2000 rho_w 0.500 0.5654 0.4653 0.6524 beta 0.2000 rho_p 0.500 0.4810 0.3858 0.5611 beta 0.2000 rho_x 0.500 0.8357 0.7883 0.8822 beta 0.2000 rho_r 0.500 0.5607 0.5245 0.6039 beta 0.2000 rho_gm 0.500 0.4724 0.4215 0.5712 beta 0.2000 rho_gp 0.500 0.1292 0.0605 0.2143 beta 0.2000 rho_gi 0.500 0.4605 0.3932 0.5130 beta 0.2000 rho_T 0.500 0.0763 0.0197 0.1379 beta 0.2000 standard deviation of shocks prior mean post. mean 90% HPD interval prior pstdev e_z 0.100 0.8713 0.6148 1.1362 invg Inf e_b 0.100 1.4480 1.0546 1.8461 invg Inf e_i 0.100 2.0133 1.6385 2.4147 invg Inf e_w 0.100 0.5201 0.4243 0.6141 invg Inf e_p 0.100 0.2818 0.2391 0.3219 invg Inf e_x 0.100 4.1148 3.5547 4.6631 invg Inf e_r 0.100 0.2149 0.1844 0.2449 invg Inf e_gm 0.100 0.5302 0.4506 0.6063 invg Inf e_gp 0.100 0.7289 0.6288 0.8307 invg Inf e_gi 0.100 1.7593 1.4829 2.0377 invg Inf e_T 0.100 0.1027 0.0227 0.2279 invg Inf e_tc 0.100 0.0078 0.0070 0.0086 invg Inf e_tw 0.100 0.4146 0.3992 0.4391 invg Inf e_tk 0.100 0.0093 0.0083 0.0104 invg Inf Estimation::mcmc: Smoothed variables Estimation::mcmc: Smoothed variables, done! Estimation::mcmc: Smoothed shocks Estimation::mcmc: Smoothed shocks, done! Estimation::mcmc: Trend_coefficients Estimation::mcmc: Trend_coefficients, done! Estimation::mcmc: Smoothed constant Estimation::mcmc: Smoothed constant, done! Estimation::mcmc: Smoothed trend Estimation::mcmc: Smoothed trend, done! Estimation::mcmc: Updated Variables Estimation::mcmc: Updated Variables, done! MODEL SUMMARY Number of variables: 53 Number of stochastic shocks: 14 Number of state variables: 31 Number of jumpers: 11 Number of static variables: 19 MATRIX OF COVARIANCE OF EXOGENOUS SHOCKS Variables e_b e_w e_z e_p e_i e_x e_r e_gm e_gp e_gi e_T e_tc e_tw e_tk e_b 2.096836 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 e_w 0.000000 0.270499 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 e_z 0.000000 0.000000 0.759137 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 e_p 0.000000 0.000000 0.000000 0.079426 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 e_i 0.000000 0.000000 0.000000 0.000000 4.053321 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 e_x 0.000000 0.000000 0.000000 0.000000 0.000000 16.931574 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 e_r 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.046161 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 e_gm 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.281137 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 e_gp 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.531339 0.000000 0.000000 0.000000 0.000000 0.000000 e_gi 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 3.095251 0.000000 0.000000 0.000000 0.000000 e_T 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.010542 0.000000 0.000000 0.000000 e_tc 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000061 0.000000 0.000000 e_tw 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.171925 0.000000 e_tk 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000087 POLICY AND TRANSITION FUNCTIONS dy_obs dhc_obs dhi_obs dgpc_obs dgic_obs dgi_obs dw_obs l_obs pi_obs rn_obs tax_c_obs tax_w_obs tax_k_obs Constant 0.178004 0.178004 0.178004 0.178004 0.178004 0.178004 0.178004 -0.022762 0.143158 0.682254 0.061000 0.273000 0.446000 c(-1) 0 -1.000000 0 0 0 0 0 0 0 0 0 0 0 rn(-1) 0.163722 0.225694 0.248163 0 0 0 0.330991 0.183729 0.361487 0.778867 0 0 0 k(-1) -0.067635 -0.050702 -0.176610 0 0 0 -0.110924 -0.395212 -0.183351 -0.006738 0 0 0 y(-1) -1.013395 -0.044085 -0.023038 -0.041129 0.106430 0.120276 -0.010880 -0.024236 -0.012126 -0.000626 -0.003217 0.039690 0.002063 b(-1) 0.213961 0.300991 0.259418 0.035607 0.173213 -0.057638 0.289920 0.353646 0.217942 0.010703 0.002592 -0.006133 -0.005159 k_g(-1) -0.130231 -0.203070 -0.156390 0 0 0 -0.224551 -0.560142 -0.280467 -0.010846 0 0 0 gm(-1) 0.129982 0.006834 0.092071 0 -0.230664 0.791654 0.131276 0.187160 0.092955 0.005346 0 0 0 gp(-1) 0.091027 0.037460 0.067979 -0.149052 0 0 0.094396 0.134299 0.068766 0.003851 0 0 0 gi(-1) 0 0 0 0 0 -1.000000 0 0 0 0 0 0 0 t(-1) 0.261445 0.425276 0.283294 0 0 0 0.329266 0.420991 0.246370 0.012494 0 0 0 t_c(-1) -0.784317 -1.400952 -0.631855 0 0 0 -0.874758 -1.200800 -0.655072 -0.035017 -0.058692 0 0 t_w(-1) -0.467617 -0.744097 -0.535513 0 0 0 -0.375769 -0.763958 -0.353354 -0.019787 0 -0.183785 0 z_w(-1) 0.027416 0.083175 -0.037493 0 0 0 0.445126 -0.027310 0.144501 0.004788 0 0 0 z_p(-1) -0.114565 -0.144053 -0.197824 0 0 0 -0.382981 -0.065314 0.320150 0.007072 0 0 0 z_x(-1) 0.057355 -0.029543 -0.051213 0 0 0 -0.003623 0.046695 -0.009524 0.000878 0 0 0 z_r(-1) -0.055330 -0.104680 -0.034385 0 0 0 0.054715 -0.157371 0.190059 0.565202 0 0 0 z_gm(-1) 0.084105 0.048997 0.090510 0 0.472400 0 0.117602 0.137951 0.088059 0.004277 0 0 0 z_gp(-1) 0.014991 0.007558 0.011657 0.129172 0 0 0.016074 0.022377 0.011772 0.000647 0 0 0 z_gi(-1) 0.030480 0.013647 0.005423 0 0 0.460528 0.015753 0.035615 0.008813 0.000873 0 0 0 z_T(-1) 0.023938 0.038835 0.026120 0 0 0 0.030322 0.038637 0.022706 0.001148 0 0 0 ypot(-1) -0.013000 0.034401 -0.009112 0.005522 -0.279643 -0.062638 -0.025401 -0.017684 -0.015117 -0.000705 0.000624 -0.033557 0.003097 pi(-2) -0.003007 -0.006142 -0.001079 0 0 0 0.003168 -0.009076 0.010470 0.029553 0 0 0 pi(-3) -0.000913 -0.002514 0.000801 0 0 0 0.002963 -0.004904 0.006543 0.029480 0 0 0 c_e(-1) 0.171396 0.425128 -0.069168 0 0 0 0.004038 0.157001 -0.018093 0.002927 0 0 0 pi(-1) -0.161115 -0.228805 -0.232532 0 0 0 -0.386970 -0.173168 0.523378 0.041396 0 0 0 i(-1) 0.139031 -0.105161 -0.272240 0 0 0 -0.030125 0.089745 -0.044729 0.001493 0 0 0 w(-1) 0.016988 0.055172 -0.029562 0 0 0 -0.475929 -0.089373 0.105702 0.003440 0 0 0 t_k(-1) -0.456872 -0.625453 -0.700093 0 0 0 -0.533928 -0.371186 -0.303935 -0.018122 0 0 -0.125260 z_b(-1) 0.048606 0.128008 -0.032584 0 0 0 0.039447 0.033857 0.003041 0.001069 0 0 0 z_i(-1) -0.033994 0.062048 -0.241233 0 0 0 0.024105 -0.004874 0.025364 0.000058 0 0 0 z_z(-1) 0.038472 0.028736 0.005133 0.026043 0.040303 0.036403 0.009328 -0.044812 -0.032976 -0.001446 0 0 0 e_b 0.149521 0.393773 -0.100235 0 0 0 0.121345 0.104151 0.009355 0.003290 0 0 0 e_w 0.048492 0.147116 -0.066317 0 0 0 0.787317 -0.048305 0.255587 0.008469 0 0 0 e_z 0.220186 0.164463 0.029379 0.149052 0.230664 0.208346 0.053387 -0.256473 -0.188734 -0.008277 0 0 0 e_p -0.238198 -0.299509 -0.411306 0 0 0 -0.796275 -0.135798 0.665640 0.014703 0 0 0 e_i -0.084261 0.153802 -0.597955 0 0 0 0.059750 -0.012081 0.062872 0.000143 0 0 0 e_x 0.068633 -0.035353 -0.061284 0 0 0 -0.004335 0.055876 -0.011397 0.001050 0 0 0 e_r -0.098671 -0.186679 -0.061319 0 0 0 0.097575 -0.280645 0.338939 1.007943 0 0 0 e_gm 0.178038 0.103719 0.191596 0 1.000000 0 0.248946 0.292022 0.186407 0.009054 0 0 0 e_gp 0.116058 0.058510 0.090241 1.000000 0 0 0.124436 0.173231 0.091132 0.005012 0 0 0 e_gi 0.066184 0.029633 0.011775 0 0 1.000000 0.034206 0.077336 0.019136 0.001896 0 0 0 e_T 0.313623 0.508789 0.342204 0 0 0 0.397268 0.506204 0.297486 0.015044 0 0 0 e_tc -0.833220 -1.488304 -0.671252 0 0 0 -0.929300 -1.275672 -0.695917 -0.037201 1.000000 0 0 e_tw -0.572908 -0.911643 -0.656093 0 0 0 -0.460380 -0.935976 -0.432917 -0.024243 0 1.000000 0 e_tk -0.522295 -0.715015 -0.800343 0 0 0 -0.610385 -0.424339 -0.347457 -0.020717 0 0 1.000000 THEORETICAL MOMENTS VARIABLE MEAN STD. DEV. VARIANCE dy_obs 0.1780 0.6576 0.4325 dhc_obs 0.1780 0.9545 0.9110 dhi_obs 0.1780 1.4459 2.0905 dgpc_obs 0.1780 0.8046 0.6474 dgic_obs 0.1780 0.6730 0.4530 dgi_obs 0.1780 2.1344 4.5557 dw_obs 0.1780 0.7539 0.5684 l_obs -0.0228 2.2539 5.0803 pi_obs 0.1432 1.9504 3.8041 rn_obs 0.6823 1.7331 3.0035 tax_c_obs 0.0610 0.0104 0.0001 tax_w_obs 0.2730 0.4399 0.1935 tax_k_obs 0.4460 0.0150 0.0002 VARIANCE DECOMPOSITION (in percent) e_b e_w e_z e_p e_i e_x e_r e_gm e_gp e_gi e_T e_tc e_tw e_tk dy_obs 17.46 2.31 19.32 3.11 7.68 24.00 0.46 3.13 2.16 4.10 0.32 0.01 15.93 0.01 dhc_obs 46.93 2.12 10.19 2.48 12.58 3.26 0.51 1.07 0.49 0.41 0.39 0.02 19.57 0.01 dhi_obs 1.89 1.30 2.95 1.84 80.64 5.19 0.18 0.79 0.34 0.03 0.09 0.00 4.75 0.00 dgpc_obs 0.06 0.02 11.11 0.01 0.02 0.06 0.00 0.04 88.60 0.00 0.00 0.00 0.07 0.00 dgic_obs 0.55 1.71 16.30 2.03 0.79 0.10 0.36 76.74 0.24 0.14 0.06 0.00 0.99 0.00 dgi_obs 0.10 0.13 1.86 0.15 0.05 0.06 0.03 5.16 0.02 92.34 0.01 0.00 0.11 0.00 dw_obs 8.28 40.29 12.12 16.17 4.24 0.20 0.65 4.55 2.14 0.83 0.43 0.01 10.06 0.01 l_obs 0.84 2.46 82.51 1.13 1.40 1.50 0.30 1.16 0.62 0.50 0.11 0.00 7.46 0.00 pi_obs 0.44 1.91 60.61 1.79 4.66 1.21 18.75 2.00 0.86 0.20 0.19 0.01 7.39 0.00 rn_obs 0.10 0.14 70.71 0.05 0.72 0.30 25.52 0.43 0.17 0.09 0.04 0.00 1.74 0.00 tax_c_obs 2.09 0.70 30.25 0.31 0.89 2.26 0.10 1.70 0.73 0.10 0.16 57.60 3.10 0.00 tax_w_obs 0.07 0.02 0.32 0.03 0.08 0.10 0.00 0.03 0.02 0.01 0.00 0.00 99.31 0.00 tax_k_obs 1.88 1.96 46.46 1.58 0.38 1.19 0.37 1.68 0.70 0.04 0.17 0.01 1.97 41.60 MATRIX OF CORRELATIONS Variables dy_obs dhc_obs dhi_obs dgpc_obs dgic_obs dgi_obs dw_obs l_obs pi_obs rn_obs tax_c_obs tax_w_obs tax_k_obs dy_obs 1.0000 0.5248 0.3919 0.2840 0.3051 0.2513 0.5290 0.1337 -0.1167 -0.0395 0.1007 -0.4028 -0.0289 dhc_obs 0.5248 1.0000 -0.1487 0.1604 0.1479 0.0772 0.6813 0.0876 -0.0970 -0.0458 0.0687 -0.4465 -0.0232 dhi_obs 0.3919 -0.1487 1.0000 0.1069 0.0985 0.0250 0.0503 0.0288 -0.0616 -0.0124 0.0416 -0.1986 0.0145 dgpc_obs 0.2840 0.1604 0.1069 1.0000 0.0598 0.0296 0.2369 0.0611 -0.0676 -0.0482 -0.0221 -0.0227 0.1077 dgic_obs 0.3051 0.1479 0.0985 0.0598 1.0000 0.1101 0.2822 0.0782 -0.0322 -0.0276 0.1205 0.0164 -0.1875 dgi_obs 0.2513 0.0772 0.0250 0.0296 0.1101 1.0000 0.1136 0.0570 -0.0061 -0.0109 0.0186 0.0024 -0.0120 dw_obs 0.5290 0.6813 0.0503 0.2369 0.2822 0.1136 1.0000 0.1235 -0.1426 -0.0558 0.0736 -0.2954 -0.0306 l_obs 0.1337 0.0876 0.0288 0.0611 0.0782 0.0570 0.1235 1.0000 -0.6071 -0.7403 -0.1571 -0.0878 0.0911 pi_obs -0.1167 -0.0970 -0.0616 -0.0676 -0.0322 -0.0061 -0.1426 -0.6071 1.0000 0.9072 0.0356 0.0007 -0.0494 rn_obs -0.0395 -0.0458 -0.0124 -0.0482 -0.0276 -0.0109 -0.0558 -0.7403 0.9072 1.0000 0.0682 0.0262 -0.0626 tax_c_obs 0.1007 0.0687 0.0416 -0.0221 0.1205 0.0186 0.0736 -0.1571 0.0356 0.0682 1.0000 -0.0533 -0.4217 tax_w_obs -0.4028 -0.4465 -0.1986 -0.0227 0.0164 0.0024 -0.2954 -0.0878 0.0007 0.0262 -0.0533 1.0000 0.0089 tax_k_obs -0.0289 -0.0232 0.0145 0.1077 -0.1875 -0.0120 -0.0306 0.0911 -0.0494 -0.0626 -0.4217 0.0089 1.0000 COEFFICIENTS OF AUTOCORRELATION Order 1 2 3 4 5 dy_obs 0.1727 0.0161 -0.0464 -0.0639 -0.0583 dhc_obs 0.1267 -0.0887 -0.1248 -0.1018 -0.0632 dhi_obs 0.3122 0.0387 -0.0673 -0.1028 -0.1068 dgpc_obs 0.1374 0.0202 0.0001 -0.0056 -0.0078 dgic_obs 0.3785 0.0411 -0.0819 -0.1261 -0.1315 dgi_obs -0.2176 -0.1068 -0.0578 -0.0336 -0.0216 dw_obs 0.3173 -0.0402 -0.1824 -0.1937 -0.1420 l_obs 0.9497 0.8962 0.8538 0.8259 0.8106 pi_obs 0.9672 0.9079 0.8484 0.8017 0.7702 rn_obs 0.9897 0.9708 0.9492 0.9276 0.9070 tax_c_obs 0.3797 0.3457 0.3117 0.2825 0.2593 tax_w_obs -0.0936 -0.0833 -0.0688 -0.0540 -0.0411 tax_k_obs 0.5286 0.4761 0.4214 0.3732 0.3338 {エラー: dseries/subsref (行 212) dseries::subsref: Indices are out of bounds! Subsample cannot end after 152Y. エラー: makedataset (行 253) DynareDataset = DynareDataset(FIRSTOBS:lastobs); エラー: dynare_estimation_init (行 553) [dataset_, dataset_info, newdatainterfaceflag] = makedataset(options_, options_.dsge_var*options_.dsge_varlag, gsa_flag); エラー: evaluate_smoother (行 64) [dataset_,dataset_info,xparam1, hh, M_, options_, oo_, estim_params_,bayestopt_] = dynare_estimation_init(var_list, M_.fname, [], M_, options_, oo_, estim_params_, bayestopt_); エラー: shock_decomposition (行 87) [oo_, M_, ~, ~, Smoothed_Variables_deviation_from_mean] = evaluate_smoother(parameter_set, varlist, M_, oo_, options_, bayestopt_, estim_params_); エラー: jpmodel5.driver (行 847) oo_ = shock_decomposition(M_,oo_,options_,var_list_,bayestopt_,estim_params_); エラー: dynare (行 293) evalin('base',[fname '.driver']) ; }