From 14ef010b9e483cf4a396c3f5007a41e2547a3f72 Mon Sep 17 00:00:00 2001 From: Johannes Pfeifer Date: Sun, 16 Feb 2020 12:14:04 +0100 Subject: [PATCH] geweke_moments.m: correctly use parameter instead of hard-coded value --- matlab/convergence_diagnostics/geweke_moments.m | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/matlab/convergence_diagnostics/geweke_moments.m b/matlab/convergence_diagnostics/geweke_moments.m index 61c22979f..51f894edd 100644 --- a/matlab/convergence_diagnostics/geweke_moments.m +++ b/matlab/convergence_diagnostics/geweke_moments.m @@ -89,7 +89,7 @@ results_struct.rne_iid = results_vec(1,4); centered_window_means=window_means-total_mean; autocov_grouped_means=zeros(n_groups,1); for lag=0:n_groups-1 - autocov_grouped_means(lag+1)=centered_window_means(lag+1:n_groups,1)'*centered_window_means(1:n_groups-lag,1)/100; + autocov_grouped_means(lag+1)=centered_window_means(lag+1:n_groups,1)'*centered_window_means(1:n_groups-lag,1)/n_groups; end % numerical standard error with tapered autocovariance functions