-
Notifications
You must be signed in to change notification settings - Fork 21
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Miscellaneous contributions from Omar #47
base: main
Are you sure you want to change the base?
Conversation
set deterministic ut0 shock
- remove any debugging code - unify the function signature between `irfs` and `irfs2` - call `irfs2` properly from `stoch_simul_core!`
display_autocorrelation2(simulation_results, endogenous_names, m, options) | ||
end | ||
|
||
function long_second_order_simulation(context; periods = 100_000, burning = 100) |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
long_second_order_simulation
needs a better name but I didn't know what to name it given the context
data = Matrix{Any}(undef, original_endo_nbr + 1, 4) | ||
data[1, :] = ["VARIABLE", "MEAN", "STD. DEV.", "VARIANCE"] | ||
data[2:end, 1] = endogenous_names[1:original_endo_nbr] | ||
data[2:end, 2] .= map(mean, eachcol(sim_results)) |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
The mean calculated in display_mean_sd_variance2
differs significantly that the theoretical moments calculated in display_mean_sd_variance
for first order solution.
Notice: the simulated 2nd order mean changes if the shock in the mod file is modified.
Michel argued that this should not be happening.
except of function names, most of this code is probably generalization to k-order solutions
this is just a draft PR to keep track any issues or WIP tasks that need to be done before @Omar-Elrefaei 's work is merged