Displaying 1 result from an estimated 1 matches for "o_j".
Did you mean:
o__
2007 Dec 31
0
Optimize jackknife code
...e indizes of the block boundaries
a = 1+k*(i-1);
b = k*i;
# sum of block to delete
delete_block_sums[i] = sum(O[a:b]);
}
v = (total_sum - delete_block_sums) / (N-k)
# The Jackknife error is given by
# eps^2 = (N_B-1)/N_B sum( (O_J - \bar{O})^2 )
# I don't understand the prefactor
#c(k, (N_B-1)**2/N_B * var(v));
print(c(k, total_sum/N, sqrt( (N_B-1)/N_B * sum( (v-mean(v))**2 )) ));
}
--