diff --git a/tests/testthat/tests-ll.r b/tests/testthat/tests-ll.r
index 73085038..28bfe3b5 100644
--- a/tests/testthat/tests-ll.r
+++ b/tests/testthat/tests-ll.r
@@ -4,29 +4,44 @@ test_that("Likelihoods can be calculated", {
   expect_lt(chain_ll(chains, "pois", "size", lambda = 0.5), 0)
   expect_lt(chain_ll(chains, "pois", "size", lambda = 0.5, exclude = 1), 0)
   expect_lt(chain_ll(chains, "pois", "size", lambda = 0.5, infinite = 5), 0)
-  expect_lt(chain_ll(chains, "pois", "size",
-    lambda = 0.5, obs_prob = 0.5,
-    nsim_obs = 1
+  expect_lt(chain_ll(
+    chains, "pois", "size", lambda = 0.5, obs_prob = 0.5, nsim_obs = 1
   ), 0)
-  expect_lt(chain_ll(chains, "pois", "length",
-    lambda = 0.5, obs_prob = 0.5,
-    nsim_obs = 1
+  expect_lt(chain_ll(
+    chains, "pois", "length", lambda = 0.5, obs_prob = 0.5, nsim_obs = 1
   ), 0)
-  expect_lt(chain_ll(chains, "pois", "size",
-    lambda = 0.5, infinite = 5,
-    obs_prob = 0.5, nsim_obs = 1
+  expect_lt(chain_ll(
+    chains, "pois", "size", lambda = 0.5, infinite = 5, obs_prob = 0.5,
+    nsim_obs = 1
   ), 0)
   expect_lt(chain_ll(chains, "binom", "size", size = 1, prob = 0.5), 0)
 })
 
-test_that("Analytical size or length distributions are implemented", {
-  expect_true(all(pois_size_ll(chains, lambda = 0.5) < 0))
-  expect_true(all(nbinom_size_ll(chains, mu = 0.5, size = 0.2) < 0))
-  expect_true(all(nbinom_size_ll(chains, prob = 0.5, size = 0.2) < 0))
-  expect_true(all(gborel_size_ll(chains, prob = 0.5, size = 0.2) < 0))
-  expect_true(all(gborel_size_ll(chains, prob = 0.5, size = 0.2) < 0))
-  expect_true(all(pois_length_ll(chains, lambda = 0.5) < 0))
-  expect_true(all(geom_length_ll(chains, prob = 0.5) < 0))
+test_that("Analytical size or length distributions are correctly calculated", {
+  expect_identical(
+    round(chain_ll(chains, "pois", "size", lambda = 0.5), 5),
+    -8.6072
+  )
+  expect_identical(
+    round(chain_ll(chains, "nbinom", "size", mu = 0.5, size = 0.2), 5),
+    -9.13437
+  )
+  expect_identical(
+    round(chain_ll(chains, "nbinom", "size", prob = 0.5, size = 0.2), 5),
+    -10.88944
+  )
+  expect_identical(
+    round(chain_ll(chains, "gborel", "size", prob = 0.5, size = 0.2), 5),
+    -11.21929
+  )
+  expect_identical(
+    round(chain_ll(chains, "pois", "length", lambda = 0.5), 5),
+    -9.39945
+  )
+  expect_identical(
+    round(chain_ll(chains, "geom", "length", prob = 0.5), 5),
+    -12.48639
+  )
 })
 
 test_that("Errors are thrown", {