diff options
| author | Justin Bedo <cu@cua0.org> | 2023-01-17 21:28:02 +1100 |
|---|---|---|
| committer | Justin Bedo <cu@cua0.org> | 2023-01-20 15:10:19 +1100 |
| commit | 5c4ffb1cd6c206884f6d0b1347f86964306a78c8 (patch) | |
| tree | a82d418b3ca9580a1c75e8fb0d1feea0340115af | |
| parent | 307d1afe9b82d24e9740cbb6b1bb5f9fa1cf8f2c (diff) | |
remove heterozygous state from likelihood
| -rw-r--r-- | bin/cluster.hs | 2 |
1 files changed, 1 insertions, 1 deletions
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] |