test_that("Stuart_marginal_homogeneity", {
  n <- vision_data
  result <- Stuart_marginal_homogeneity(n)
  chisq <- result$chisq
  expect_true(abs(11.957 - chisq) <= 0.0005)
  df <- result$df
  expect_identical(3, df)
}
)