diff --git a/README.Rmd b/README.Rmd index fa1c78c..e48746c 100644 --- a/README.Rmd +++ b/README.Rmd @@ -93,26 +93,3 @@ irf = compute_impulse_responses(posterior, horizon = 60) plot(irf, probability = 0.68) ``` -A repliaction of Arias, Rubio-Ramírez and Waggoner (2018). - -```r -data(optimism) - -# optimism shock -# no contemporaneous effect on productivity -zero_irf = matrix(0, nrow = 5, ncol = 5) -zero_irf[1, 1] = 1 -# positive contemporaneous effect on stock prices -sign_irf = array(0, dim = c(5, 5, 1)) -sign_irf[2, 1, 1] = 1 - -specification = specify_bsvarSIGN$new(optimism * 100, - p = 4, - sign_irf = sign_irf, - zero_irf = zero_irf) - -posterior = estimate(specification, S = 100) -irf = compute_impulse_responses(posterior, horizon = 40) -plot(irf, probability = 0.68) -``` - diff --git a/README.md b/README.md index 798846a..75ed3f2 100644 --- a/README.md +++ b/README.md @@ -112,26 +112,3 @@ posterior = estimate(specification, S = 100) irf = compute_impulse_responses(posterior, horizon = 60) plot(irf, probability = 0.68) ``` - -A repliaction of Arias, Rubio-Ramírez and Waggoner (2018). - -``` r -data(optimism) - -# optimism shock -# no contemporaneous effect on productivity -zero_irf = matrix(0, nrow = 5, ncol = 5) -zero_irf[1, 1] = 1 -# positive contemporaneous effect on stock prices -sign_irf = array(0, dim = c(5, 5, 1)) -sign_irf[2, 1, 1] = 1 - -specification = specify_bsvarSIGN$new(optimism * 100, - p = 4, - sign_irf = sign_irf, - zero_irf = zero_irf) - -posterior = estimate(specification, S = 100) -irf = compute_impulse_responses(posterior, horizon = 40) -plot(irf, probability = 0.68) -```