tests/testthat/test-regressionloglinear.R

In jasp-stats/jaspFrequencies: Frequencies Module for JASP

context("Log-Linear Regression")

test_that("Main table results match", {
  options <- jaspTools::analysisOptions("RegressionLogLinear")
  options$count <- "facFifty"
  options$factors <- c("contBinom", "facGender")
  options$modelTerms <- list(
    list(components="contBinom"),
    list(components="facGender"),
    list(components=c("contBinom", "facGender"))
  )
  results <- jaspTools::runAnalysis("RegressionLogLinear", "test.csv", options)
  table <- results[["results"]][["AnovaTable"]][["data"]]
  jaspTools::expect_equal_tables(table,
    list(" ", " ", "NULL", " ", 936.356249443911, 99, 9.73545292814663,
         1, "contBinom", 0.00180747470901472, 926.620796515764, 98, 7.02546792150429,
         1, "facGender", 0.00803584702758609, 919.59532859426, 97, 0.769509416901883,
         1, "contBinom<unicode><unicode><unicode>facGender",
         0.380368860922424, 918.825819177358, 96)
  )
})

test_that("Coefficients table matches", {
  options <- jaspTools::analysisOptions("RegressionLogLinear")
  options$count <- "facFifty"
  options$factors <- c("contBinom", "facFive")
  options$modelTerms <- list(
    list(components="contBinom"),
    list(components="facFive"),
    list(components=c("contBinom", "facFive"))
  )
  options$regressionCoefficientsEstimates <- TRUE
  options$regressionCoefficientsCi <- TRUE
  options$regressionCoefficientsCiLevel <- 0.95
  results <- jaspTools::runAnalysis("RegressionLogLinear", "test.csv", options)
  table <- results[["results"]][["CoefficientsTable"]][["data"]]
  jaspTools::expect_equal_tables(table,
    list("(Intercept)", 3.36441813015886, 0.0536828127084252, 3.25920175066154,
         3.46963450965618, 62.6721656414034, 0, "contBinom = 1", -0.63167531645938,
         0.104937579028612, -0.837349191980285, -0.426001440938476, -6.01953392013314,
         1.74920000092307e-09, "facFive = 2", -0.379859157202991, 0.0921944777734458,
         -0.560557013212423, -0.199161301193559, -4.1201942499901, 3.78553110594005e-05,
         "facFive = 3", -0.343993244014492, 0.0833479734033365, -0.507352270069434,
         -0.18063421795955, -4.1271938592897, 3.67216813648207e-05, "facFive = 4",
         -0.250902820948484, 0.085593742903271, -0.418663474340876, -0.0831421675560923,
         -2.93132199198285, 0.00337522739489266, "facFive = 5", -0.0685812641545293,
         0.0731503942633908, -0.211953402365681, 0.0747908740566221,
         -0.937537860802095, 0.348481958057894, "contBinom = 1*facFive = 2",
         0.992224258514625, 0.140836788661736, 0.716189225039344, 1.26825929198991,
         7.04520649713029, 1.85186666239764e-12, "contBinom = 1*facFive = 3",
         1.10775799178149, 0.137349243064596, 0.838558422071047, 1.37695756149194,
         8.06526462807306, 7.30773999647178e-16, "contBinom = 1*facFive = 4",
         0.935886690862371, 0.136875960028964, 0.667614738856258, 1.20415864286848,
         6.83748037759394, 8.05981239553115e-12, "contBinom = 1*facFive = 5",
         0.703134280441524, 0.142750071638859, 0.423349281238848, 0.982919279644201,
         4.92563171681182, 8.40882439228307e-07)
  )
})

test_that("Analysis handles errors - infinity", {
  options <- jaspTools::analysisOptions("RegressionLogLinear")
  options$factors <- "facGender"
  options$count <- "debInf"
  results <- jaspTools::runAnalysis("RegressionLogLinear", "test.csv", options)
  status <- results[["status"]]
  expect_identical(status, "validationError")
})

test_that("Analysis handles errors - missing values (factors)", {
  options <- jaspTools::analysisOptions("RegressionLogLinear")
  options$factors <- "debBinMiss20"
  options$modelTerms <- list(
    list(components="debBinMiss20")
  )
  results <- jaspTools::runAnalysis("RegressionLogLinear", "test.csv", options)
  status <- results[["status"]]
  expect_identical(status, "validationError")
})

test_that("Analysis handles errors - missing values (counts)", {
  options <- jaspTools::analysisOptions("RegressionLogLinear")
  options$factors <- "contBinom"
  options$count <- "debMiss30"
  results <- jaspTools::runAnalysis("RegressionLogLinear", "test.csv", options)
  status <- results[["status"]]
  expect_identical(status, "validationError")
})

test_that("Analysis handles errors - negatives", {
  options <- jaspTools::analysisOptions("RegressionLogLinear")
  options$factors <- "facGender"
  options$count <- "contNormal"
  results <- jaspTools::runAnalysis("RegressionLogLinear", "test.csv", options)
  status <- results[["status"]]
  expect_identical(status, "validationError")
})