tests/testthat/test-apa_corr.R

In stp4/stp25stat: Functions for statistical computations

context("test-apa_corr")

test_that("corr formula works", {
  set.seed(0815)
  n <- 2 * 20
  e <- rnorm(n)
  dat <- stp25aggregate::Label(
    data.frame(
      a = rnorm(n) + e / 2,
      b = rnorm(n) + e,
      c = rnorm(n),
      d = rnorm(n) + e * 10,
      g = gl(2, 20, labels = c("Control", "Treat"))
    ),
    a = "Alpha",
    b = "Beta",
    c = "Gamma"
  )


  expect_equal(APA_Correlation(~ a + b + c,
                               dat, output = FALSE)[1, 3],
               ".45**")
  expect_equal(APA_Correlation(a ~ c,
                               dat, output = FALSE)[1, 4],
               ".706")
  expect_equal(APA_Correlation(a + b + c ~ d,
                               dat, output = FALSE)[1, 3],
               ".49")
  expect_equal(APA_Correlation(a + b ~ c +  d,
                               dat, output = FALSE)[1, 6],
               ".49")
  expect_equal(APA_Correlation(a + b + c ~ d,
                               dat,
                               groups = ~ g, output =FALSE)[2, 3],
               ".75")
  
  
  
  
  x<-APA_Correlation( ~ a + b + c, dat, include.mean=TRUE, output=FALSE)
  
  expect_equal(x[,2],
 c( "-0.03 (1.19)", "-0.17 (1.34)", "0.23 (1.01)" )
  )
  

})