tests/testthat/test-contingency_supp4_1.R

In cjvanlissa/gorica: Evaluation of Inequality Constrained Hypotheses Using GORICA

dat <- structure(c(7L, 5L, 8L, 0L), .Dim = c(2L, 2L))

class(dat) <- "table"
#Parameters used in evaluation (eta)
hyp <- "a:=x[1,1]/(x[1,1]+x[1,2]);
b:=x[2,1]/(x[2,1]+x[2,2]);
a>b"
#Number of equality and inequality constraints per model
# gorica:::parse_hypothesis(varnames = gorica:::rename_table_est(c("x[1,1]", "x[2,1]", "x[1,2]", "x[2,2]")), gorica:::rename_table_est(hyp))

test_that("contingency table supp 4_1 throws error", {
  expect_error(gorica(dat, hyp))
})
# LogLikelihood Penalty GORICA_value GORICA_weight
# H1     0.1956781  0.5004    0.6094439     0.6223653
# H2     0.1956781  1.0000    1.6086439     0.3776347
#
# some estimated as 0