If I wanted a 95% credible interval instead of a 90% credible interval for the QTL peak, would it suffice to change this code in mp_inference.py as follows?
cumul, mean = 0.0, 0.0
left, right = None, None
for i,val in enumerate(temp):
cumul += val
if cumul >= 0.05 and left is None:
left = i-1
if cumul >= 0.95 and right is None:
right = i
mean += val*bins[i]
cumul, mean = 0.0, 0.0
left, right = None, None
for i,val in enumerate(temp):
cumul += val
if cumul >= 0.025 and left is None:
left = i-1
if cumul >= 0.975 and right is None:
right = i
mean += val*bins[i]
Thank you for the great software.
If I wanted a 95% credible interval instead of a 90% credible interval for the QTL peak, would it suffice to change this code in mp_inference.py as follows?
cumul, mean = 0.0, 0.0
left, right = None, None
for i,val in enumerate(temp):
cumul += val
if cumul >= 0.05 and left is None:
left = i-1
if cumul >= 0.95 and right is None:
right = i
mean += val*bins[i]
cumul, mean = 0.0, 0.0
left, right = None, None
for i,val in enumerate(temp):
cumul += val
if cumul >= 0.025 and left is None:
left = i-1
if cumul >= 0.975 and right is None:
right = i
mean += val*bins[i]
Thank you for the great software.