knitr::opts_chunk$set( collapse = TRUE, comment = "#>" )
library(SolowVariants) theme_set(theme_bw()) # Setting up the parameter grid --------------------------------- aux_np <- 200 aux_parameter_grid <- create_parameter_grid( c("alpha", "phi", "lambda", "rho", "s", "sR", "delta", "n"), c(1/3, 1, 0.2, 0.2, 0.2, 0.15, 0.15, 0), c(NA, NA, NA, NA, NA, NA, NA, NA), c(NA, NA, NA, NA, NA, NA, NA, NA), aux_np ) # Setting starting values for K and L --------------------------------- aux_startvalues <- list(A = 1, K = 1, L = 1) # Simulate the basic Solow growth model --------------------------------- aux_simulation <- SimulateExtendedSolowModelEndogenousGrowthRomer(aux_parameter_grid, aux_np, aux_startvalues) # Visualise some variables --------------------------------- VisualiseSimulation( aux_simulation, c("YpW", "KpEW", "gYpW", "gKpW", "YpEW"), "free" ) # Economy in SS in the End --------------------------------- steadystate_checker(aux_simulation, aux_parameter_grid, "ESEGRomer")
Parameters defining the ESEG
Remark: Those parameters can be found via the command: getRequiredParams("BS")
.
Starting value of endogenous variables defining the BS
Remark: Those parameters can be found via the command: getRequiredStartingValues("BS")
.
$$\begin{aligned} Y_t &= K_t^\alpha (A_tL_{Y,t})^{1-\alpha} \ A_{t+1} &= \rho * A_t^\phi L_{A,t}^\lambda + A_t \ K_{t+1}&= sY_t + (1-\delta)K_t \ L_{t+1}&=(1+n)L_t \ L_t &= L_{A,t} + L_{Y,t} \ L_{A, t} &= s_R * L_t (\rightarrow L_{Y,t} = (1-s_R * L_t)) \end{aligned}$$
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.