tests/testthat/test-CvM.pvalue.R

In LiYao-sfu/EDFtest: Goodness of Fit Based on Empirical Distribution Function

context("test-CvM.pvalue.R")

test_that("CvM P-value for uniform sample", {
  x = c(0.3103205, -0.1386720, -0.5988813,  0.9483934, -0.5213117,
        0.7555062, -0.6821612,  0.7957394, -0.7387161,  0.1910647)
  wsq = CvM.uniform(x)
  p_value = CvM.uniform.pvalue(w=wsq)$P

  expect_equal(p_value,0.43739209)
})

test_that("CvM P-value for normal sample", {
  x = c(0.25024690, -0.33712454, -0.11335370, -0.09888291, 0.26408682,
        0.13898369, -0.24226950, 0.05903138, -0.17727187, 0.79468027)
  wsq = CvM.normal(x)
  p_value = CvM.normal.pvalue(w=wsq)$P

  expect_equal(p_value,0.49021528)
})

test_that("CvM for gamma sample", {
  x = c(0.5047757, 0.1538300, 0.5704100, 0.3013008, 1.2775724,
        1.0468233, 0.6525627, 0.4376768, 2.4700737, 1.0944885)
  wsq = CvM.gamma(x)
  shape = estimate.gamma(x)[2]
  p_value = CvM.gamma.pvalue(w=wsq,shape=shape)$P

  expect_equal(p_value,0.88483667)
})

test_that("CvM for logistic sample", {
  x = c(-1.1263960,  0.9562103, -3.3860294,  0.1980448,  0.7667096,
        -0.8461510, -0.4524666,  1.0070690,  3.2450939,  1.1559508)
  wsq = CvM.logistic(x)
  p_value = CvM.logistic.pvalue(w=wsq)$P

  expect_equal(p_value,0.54949935)
})

test_that("CvM for laplace sample", {
  x = c(-0.23539279,  0.16009027,  2.84634962,  0.35710312, -0.40466195,
        -0.41113889,  2.16169132, -0.27151351,  0.13770907,  0.02330074)
  wsq = CvM.laplace(x)
  p_value = CvM.laplace.pvalue(w=wsq)$P

  expect_equal(p_value,0.046049247)
})

test_that("CvM for weibull sample", {
  x = c(0.36218715, 0.16506700, 0.16757965, 0.93681048, 1.87396510,
        0.44718470, 1.24767735, 0.07435952, 1.86023456, 0.03682825)
  wsq = CvM.weibull(x)
  p_value = CvM.weibull.pvalue(w=wsq)$P

  expect_equal(p_value,0.56603163)
})

test_that("CvM for exponential sample", {
  x = c(13.0581121,  0.8301048,  0.5207504,  1.0923122,  0.7086793,
        0.1271974,  3.9326089,  0.0510448,  4.3839846,  3.4396530)
  wsq = CvM.exp(x)
  p_value = CvM.exp.pvalue(w=wsq)$P

  expect_equal(p_value,0.1653965)
})