@@ -6,7 +6,7 @@ using ABCredit, Statistics, Plots, DelimitedFiles
6
6
7
7
# install the HPFilter and KernelDensity packages
8
8
using Pkg
9
- Pkg. add (url= " https://github.com/sdBrinkmann/HPFilter.jl" )
9
+ Pkg. add (url = " https://github.com/sdBrinkmann/HPFilter.jl" )
10
10
Pkg. add (" KernelDensity" )
11
11
using HPFilter, KernelDensity
12
12
@@ -23,10 +23,10 @@ T = 1000
23
23
d = ABCredit. run_one_sim! (model, T; seed = 100 , burn_in = 300 )
24
24
25
25
# select some variables
26
- Y = d. Y_real;
27
- P = d. gdp_deflator;
28
- I = d. Investment;
29
- C = d. consumption;
26
+ Y = d. Y_real;
27
+ P = d. gdp_deflator;
28
+ I = d. Investment;
29
+ C = d. consumption;
30
30
U = d. Un;
31
31
32
32
# load the data from FRED saved in the package as an example
@@ -38,16 +38,21 @@ c_real = data[:, 4];
38
38
u_real = data[:, 5 ];
39
39
40
40
# apply HP filter to the simulated data
41
- y = log .(Y) - HP (log .(Y),1600 );
42
- c = log .(C) - HP (log .(C),1600 );
43
- invest = log .(I) - HP (log .(I),1600 );
41
+ y = log .(Y) - HP (log .(Y), 1600 );
42
+ c = log .(C) - HP (log .(C), 1600 );
43
+ invest = log .(I) - HP (log .(I), 1600 );
44
44
u = U;
45
45
pi = diff (log .(P)) .- mean ((diff (log .(P))));
46
46
47
47
# plot the histograms of real and simulated data
48
48
pdf_sim = kde (u);
49
49
pdf_real = kde (u_real);
50
- p1 = plot ([pdf_sim. density, pdf_real. density], title = " unemployment rate" , titlefont = 10 , labels= [" simulated" " real" ])
50
+ p1 = plot (
51
+ [pdf_sim. density, pdf_real. density],
52
+ title = " unemployment rate" ,
53
+ titlefont = 10 ,
54
+ labels = [" simulated" " real" ],
55
+ )
51
56
52
57
pdf_sim = kde (y);
53
58
pdf_real = kde (y_real);
@@ -65,4 +70,4 @@ pdf_sim = kde(invest);
65
70
pdf_real = kde (invest_real);
66
71
p5 = plot ([pdf_sim. density, pdf_real. density], title = " investment gap" , titlefont = 10 , legend = :none )
67
72
68
- plot (p1, p2, p3, p4, p5, layout = (2 , 3 ))
73
+ plot (p1, p2, p3, p4, p5, layout = (2 , 3 ))
0 commit comments