Commit f8aa1c1a authored by LANORE Vincent's avatar LANORE Vincent
Browse files

Modified analysis for diffsel

parent e67f59cb
......@@ -89,8 +89,8 @@ if rddata.shape[1] != 20:
MESSAGE(failure() + "Expected 20 columns!")
exit(1)
MESSAGE("Computing maximum probability per site")
rddata2 = 0.5 + abs(0.5 - rddata)
maxes = rddata2.max(axis=1)
rddata2 = 2 * abs(0.5 - rddata)
maxes = rddata2.mean(axis=1)
#===================================================================================================
print(step("Writing result to file"))
......
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