test_that("test_Clayton", {
  n <- tonsils

  n <- n[,3:1]
  log_odds_est <- Clayton_marginal_location(n[1,], n[2,])
  expect_true(abs(exp(log_odds_est$log_theta_hat) - 1.77) <= 0.005, info="pooled")
  expect_true(abs(1.77 - log_odds_est$theta_hat) <= 0.005)
  expect_true(abs(6.46 - log_odds_est$chisq_theta_hat) < 0.005)
  expect_true(abs(exp(log_odds_est$log_mh_theta_hat) - 1.79) <= 0.01,
              info=paste("mh", 1.78, exp(log_odds_est$log_mh_theta_hat)))
  expect_true(abs(1.79 - log_odds_est$mh_theta_hat) <= 0.005)
  expect_true(abs(6.65 - log_odds_est$chisq_mh_theta_hat) <= 0.005)
  expect_true(abs(0.051 - log_odds_est$var_log_theta_hat) <= 0.0005)
  expect_equal(1, log_odds_est$df)
}
)