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]
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.