From 5c4ffb1cd6c206884f6d0b1347f86964306a78c8 Mon Sep 17 00:00:00 2001 From: Justin Bedo Date: Tue, 17 Jan 2023 21:28:02 +1100 Subject: remove heterozygous state from likelihood --- bin/cluster.hs | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) (limited to 'bin') diff --git a/bin/cluster.hs b/bin/cluster.hs index dd8f2c6..0a7a438 100644 --- a/bin/cluster.hs +++ b/bin/cluster.hs @@ -76,7 +76,7 @@ model xs = do where likelihood ps cnts = product $ zipWith go ps (pairs cnts) where - go p (c, d) = max (pois (fromIntegral d * p) c) (pois (fromIntegral d * p / 2) c) + go p (c, d) = pois (fromIntegral d * p) c -- Tabulate list tabulate xs = M.elems $ M.fromListWith (+) [(c, 1) | c <- xs] -- cgit v1.2.3