tests/testthat/test-sgn.test.R

In ANSM5: Functions and Data for the Book "Applied Nonparametric Statistical Methods", 5th Edition

test_that("Example 3.1", {
  expect_equal(sgn.test(ch3$sampleI, 110)$pval.exact, 0.109375)
  expect_equal(sgn.test(ch3$sampleII, 110)$pval.exact, 0.041389465)
})

test_that("Example 3.2", {
  expect_equal(sgn.test(ch3$sampleI)$actualCIwidth.exact, 0.97851562)
  expect_equal(sgn.test(ch3$sampleI)$CI.exact.lower, 11)
  expect_equal(sgn.test(ch3$sampleI)$CI.exact.upper, 151)
  expect_equal(sgn.test(ch3$sampleII)$actualCIwidth.exact, 0.95861053)
  expect_equal(sgn.test(ch3$sampleII)$CI.exact.lower, 40)
  expect_equal(sgn.test(ch3$sampleII)$CI.exact.upper, 108)
})

test_that("Example 3.3", {
  expect_equal(sgn.test(ch3$sampleI, 110, do.asymp = TRUE)$pval.asymp, 0.113846298)
  expect_equal(sgn.test(ch3$sampleI, 110, do.asymp = TRUE, cont.corr = FALSE)$pval.asymp, 0.057779571)
  expect_equal(sgn.test(ch3$sampleII, 110, do.asymp = TRUE)$pval.asymp, 0.044171345)
  expect_equal(sgn.test(ch3$sampleII, 110, do.asymp = TRUE, cont.corr = FALSE)$pval.asymp, 0.025347319)
})

test_that("Example 3.12", {
  expect_equal(sgn.test(ch3$sampleA, 9)$pval.exact, 0.109375)
  expect_equal(sgn.test(ch3$sampleA, 9, do.asymp = TRUE, cont.corr = FALSE)$pval.asymp, 0.057779571)
  expect_equal(sgn.test(ch3$sampleA, 9, do.asymp = TRUE)$pval.asymp, 0.113846298)
  expect_equal(sgn.test(ch3$sampleB, 9)$pval.exact, 0.002576828)
  expect_equal(sgn.test(ch3$sampleB, 9, do.asymp = TRUE, cont.corr = FALSE)$pval.asymp, 0.0017451187)
  expect_equal(sgn.test(ch3$sampleB, 9, do.asymp = TRUE)$pval.asymp, 0.0036504344)
})

test_that("Example 3.13", {
  expect_equal(sgn.test(ch3$sampleB, do.asymp = TRUE)$CI.asymp.lower, 7)
  expect_equal(sgn.test(ch3$sampleB, do.asymp = TRUE)$CI.asymp.upper, 8.2)
  expect_equal(sgn.test(ch3$sampleB, 6.95, "greater")$pval.exact, 0.05765915)
  expect_equal(sgn.test(ch3$sampleB, 8.25, "less")$pval.exact, 0.05765915)
  expect_equal(sgn.test(ch3$sampleB, 7.05, "greater")$pval.exact, 0.131587982)
  expect_equal(sgn.test(ch3$sampleB, 8.15, "less")$pval.exact, 0.131587982)
})

test_that("Example 3.14", {
  expect_equal(sgn.test(ch3$heartrates2a)$CI.exact.lower, 1)
  expect_equal(sgn.test(ch3$heartrates2a)$CI.exact.upper, 17)
  expect_equal(sgn.test(ch3$heartrates2a, CI.width = 0.99)$CI.exact.lower, -2)
  expect_equal(sgn.test(ch3$heartrates2a, CI.width = 0.99)$CI.exact.upper, 20)
  expect_equal(sgn.test(ch3$heartrates2b)$CI.exact.lower, 1)
  expect_equal(sgn.test(ch3$heartrates2b)$CI.exact.upper, 17)
  expect_equal(sgn.test(ch3$heartrates2b, CI.width = 0.99)$CI.exact.lower, -2)
  expect_equal(sgn.test(ch3$heartrates2b, CI.width = 0.99)$CI.exact.upper, 20)
})

test_that("Exercise 3.1", {
  expect_equal(sgn.test(ch3$sampleIa, 110)$pval.exact, 0.34375)
})

test_that("Exercise 3.3", {
  expect_equal(sgn.test(ch3$parkingtime, 20)$pval.exact, 0.098737147)
  expect_equal(sgn.test(ch3$parkingtime, 20, do.asymp = TRUE)$pval.asymp, 0.100348246)
  expect_equal(sgn.test(ch3$parkingtime, 20, do.asymp = TRUE, cont.corr = FALSE)$pval.asymp, 0.067889155)
})

test_that("Exercise 3.10", {
  expect_equal(sgn.test(ch3$sleeptime, 2, alternative = "greater")$pval.exact, 0.032714844)
})

test_that("Exercise 3.11", {
  expect_equal(sgn.test(ch3$fishlengths)$actualCIwidth, 0.97118328)
  expect_equal(sgn.test(ch3$fishlengths)$CI.exact.lower, 69)
  expect_equal(sgn.test(ch3$fishlengths)$CI.exact.upper, 73)
  expect_equal(sgn.test(ch3$fishlengths, CI.width = 0.90)$actualCIwidth, 0.93475466)
  expect_equal(sgn.test(ch3$fishlengths, CI.width = 0.90)$CI.exact.lower, 70)
  expect_equal(sgn.test(ch3$fishlengths, CI.width = 0.90)$CI.exact.upper, 72)
})

test_that("Exercise 3.12", {
  expect_equal(sgn.test(ch3$weightloss, 5)$pval.exact, 0.17956543)
})

test_that("Exercise 3.13", {
  expect_equal(sgn.test(ch3$plants, 50)$pval.exact, 0.33678364)
})

test_that("Exercise 3.14", {
  expect_equal(sgn.test(ch3$birthprops, CI.width = 0.99)$CI.exact.lower, 0.27700)
  expect_equal(sgn.test(ch3$birthprops, CI.width = 0.99)$CI.exact.upper, 0.29800)
})

test_that("Exercise 3.15", {
  expect_equal(sgn.test(ch3$assembly)$CI.exact.lower, 11.8)
  expect_equal(sgn.test(ch3$assembly)$CI.exact.upper, 14)
})

test_that("Exercise 3.16", {
  expect_equal(sgn.test(ch3$weightchange)$CI.exact.lower, -2)
  expect_equal(sgn.test(ch3$weightchange)$CI.exact.upper, -0.1)
})

test_that("Exercise 3.17", {
  expect_equal(sgn.test(ch3$sampleI, 110)$pval.exact, 0.109375)
  expect_equal(sgn.test(ch3$sampleI, 110, do.asymp = TRUE,
                         cont.corr = FALSE)$pval.asymp, 0.057779571)
  expect_equal(sgn.test(ch3$sampleI, 110, do.asymp = TRUE)$pval.asymp,
               0.113846298)
  expect_equal(sgn.test(ch3$sampleII, 110)$pval.exact, 0.041389465)
  expect_equal(sgn.test(ch3$sampleII, 110, do.asymp = TRUE,
                         cont.corr = FALSE)$pval.asymp, 0.025347319)
  expect_equal(sgn.test(ch3$sampleII, 110, do.asymp = TRUE)$pval.asymp,
               0.044171345)
})

test_that("Exercise 4.15", {
  expect_equal(sgn.test(ch4$arrow.angles, 145)$pval.exact, 0.109375)
})

test_that("Example 5.1", {
  expect_equal(sgn.test(ch5$LVF - ch5$RVF, 0)$pval.exact, 0.0063476563)
  expect_equal(sgn.test(ch5$LVF - ch5$RVF)$CI.exact.lower, 7)
  expect_equal(sgn.test(ch5$LVF - ch5$RVF)$CI.exact.upper, 22)
})

test_that("Example 5.3", {
  expect_equal(sgn.test(ch5$bp)$CI.exact.lower, 15)
  expect_equal(sgn.test(ch5$bp)$CI.exact.upper, 40)
  expect_equal(sgn.test(ch5$bp.incorrect)$CI.exact.lower, 15)
  expect_equal(sgn.test(ch5$bp.incorrect)$CI.exact.upper, 40)
})

test_that("Example 5.4", {
  expect_equal(sgn.test(c(rep(0, 14), rep(1, 9)), 0.5)$pval.exact, 0.4048729)
  expect_equal(sgn.test(c(rep(0, 14), rep(1, 9)), 0.5, do.asymp = TRUE)$pval.asymp,
               0.4042485)
  expect_equal(sgn.test(c(rep(0, 14), rep(1, 9)), 0.5, do.asymp = TRUE,
                         cont.corr = FALSE)$pval.asymp, 0.29714653)
})

test_that("Exercise 5.1", {
  expect_equal(sgn.test(ch5$bp.diff, 0)$pval.exact, 0.2265625)
})

test_that("Exercise 5.4", {
  expect_equal(sgn.test(ch5$parent)$pval.exact, 0.049041748)
})

test_that("Exercise 5.6", {
  expect_equal(sgn.test(c(rep(0, 27), rep(1, 16)), 0.5)$pval.exact,
               0.126289474)
})

test_that("Exercise 6.2", {
  expect_equal(sgn.test(ch5$LVF - ch5$RVF, 0)$pval.exact, 0.0063476563)
})