diff --git a/meson.build b/meson.build index 126982f0441f0f009cf471e871ba72297d7f67e4..f5acef473978bda1f50a8e3a80a592e471af532a 100644 --- a/meson.build +++ b/meson.build @@ -866,7 +866,6 @@ mod_and_m_tests = [ { 'test' : [ 'fs2000/fs2000.mod' ], 'extra' : [ 'fs2000/fsdat_simul.m' ] }, { 'test' : [ 'ls2003/ls2003_hessian_zero.mod' ] }, - { 'test' : [ 'ep/rbc.mod' ] }, { 'test' : [ 'exogenous-observed-variables/preprocessor.mod' ] }, { 'test' : [ 'estimation/fs2000_with_weibull_prior.mod' ], 'extra' : [ 'estimation/fsdat_simul.m' ] }, diff --git a/tests/ep/rbc.mod b/tests/ep/rbc.mod deleted file mode 100644 index c1b29b77b28b91043de5ca09656fab74f3923c1d..0000000000000000000000000000000000000000 --- a/tests/ep/rbc.mod +++ /dev/null @@ -1,94 +0,0 @@ -var Capital, Output, Labour, Consumption, Efficiency, efficiency, ExpectedTerm; - -varexo EfficiencyInnovation; - -parameters beta, theta, tau, alpha, psi, delta, rho, effstar, sigma; - -/* -** Calibration -*/ - - -beta = 0.990; -theta = 0.357; -tau = 30.000; -alpha = 0.450; -psi = -5.000; -delta = 0.020; -rho = 0.950; -effstar = 1.000; -sigma = 0.010; - -model(use_dll); - - // Eq. n°1: - efficiency = rho*efficiency(-1) + sigma*EfficiencyInnovation; - - // Eq. n°2: - Efficiency = effstar*exp(efficiency-.5*sigma*sigma/(1-rho*rho)); - - // Eq. n°3: - Output = Efficiency*(alpha*(Capital(-1)^psi)+(1-alpha)*(Labour^psi))^(1/psi); - - // Eq. n°4: - Consumption + Capital - Output - (1-delta)*Capital(-1); - - // Eq. n°5: - ((1-theta)/theta)*(Consumption/(1-Labour)) - (1-alpha)*(Output/Labour)^(1-psi); - - // Eq. n°6: - (((Consumption^theta)*((1-Labour)^(1-theta)))^(1-tau))/Consumption - ExpectedTerm(1); - - // Eq. n°7: - ExpectedTerm = beta*((((Consumption^theta)*((1-Labour)^(1-theta)))^(1-tau))/Consumption)*(alpha*((Output/Capital(-1))^(1-psi))+1-delta); - -end; - -steady_state_model; -efficiency = 0; -Efficiency = effstar; -// Compute steady state ratios. -Output_per_unit_of_Capital=((1/beta-1+delta)/alpha)^(1/(1-psi)); -Consumption_per_unit_of_Capital=Output_per_unit_of_Capital-delta; -Labour_per_unit_of_Capital=(((Output_per_unit_of_Capital/Efficiency)^psi-alpha)/(1-alpha))^(1/psi); -Output_per_unit_of_Labour=Output_per_unit_of_Capital/Labour_per_unit_of_Capital; -Consumption_per_unit_of_Labour=Consumption_per_unit_of_Capital/Labour_per_unit_of_Capital; - -// Compute steady state share of capital. -ShareOfCapital=alpha/(alpha+(1-alpha)*Labour_per_unit_of_Capital^psi); - -/// Compute steady state of the endogenous variables. -Labour=1/(1+Consumption_per_unit_of_Labour/((1-alpha)*theta/(1-theta)*Output_per_unit_of_Labour^(1-psi))); -Consumption = Consumption_per_unit_of_Labour*Labour; -Capital = Labour/Labour_per_unit_of_Capital; -Output = Output_per_unit_of_Capital*Capital; -ExpectedTerm = beta*((((Consumption^theta)*((1-Labour)^(1-theta)))^(1-tau))/Consumption)*(alpha*((Output/Capital)^(1-psi))+1-delta); -end; - - -shocks; -var EfficiencyInnovation = 1; -end; - -steady(nocheck); - - -options_.ep.stochastic.order = 0; -options_.ep.stack_solve_algo=0; -ts0 = extended_path([], 10, [], options_, M_, oo_); - -options_.ep.stochastic.order = 1; - -options_.ep.stochastic.nodes = 3; -options_.ep.stochastic.IntegrationAlgorithm='Tensor-Gaussian-Quadrature'; -ts1_3 = extended_path([], 10, [], options_, M_, oo_); - -options_.ep.stochastic.nodes = 5; -ts1_5 = extended_path([], 10, [], options_, M_, oo_); - -options_.ep.stochastic.order = 2; -options_.ep.stochastic.nodes = 3; -ts2_3 = extended_path([], 10, [], options_, M_, oo_); - -options_.ep.stochastic.nodes = 5; -ts2_5 = extended_path([], 10, [], options_, M_, oo_);