test_that("Bowker chi-square is correct", {
  n <- vision_data
  result <- Bowker_symmetry(n)
  expect_true(abs(19.107 - result$chisq) <= 0.0005)
  expect_identical(6, result$df)
}
)