Commit 72601757 by Eric Coissac

Change the variance formula to the variance estimator formula

parent 9e789f40
......@@ -755,7 +755,7 @@ def pairEndedConnected(self,assgraph,edge1,edge2,back=250):
if (len(delta)):
dmean = sum(delta) / len(delta)
dmean2= sum(x ** 2 for x in delta) / len(delta)
dmean2= sum(x ** 2 for x in delta) / (len(delta) - 1)
dsd = math.sqrt(dmean2 - dmean ** 2)
dmean+=ri.getReadSize()
else:
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment