Hi, I found what I think is a bug in miser.c. In the function estimate_corrmc, line 689 in the file sigma_l[i] = sqrt (fsum2_l[i] - fsum_l[i] * fsum_l[i] / hits_l[i]); is supposed to be the variance in the left hand region with the division in the i-axis , Shouldn't this be: sqrt( ( fsum2_l[i] / hits_l[i] - fsum_l[i] * fsum_l[i] ) / (hits_l[i] - 1) ) and correspondingly for the right hand variance? Or is there some subtlety I am missing?
Rudy.
