tests/testthat/test_contr.wec.weighted.R

In beckerbenj/eatWtd: Weighted Descriptive Statistics

### example from wec paper
dat <- data.frame ( group = as.factor(c(rep(1,3), rep(2,2))), wgt = c(2/3, 4/3, 2, 3/8, 5/8))
con <- contr.wec.weighted(x= dat[,"group"], omitted=1,weights=dat[,"wgt"])
con1<- contr.wec.weighted(x= dat[,"group"], omitted=1,weights=rep(1,nrow(dat)))


###
test_that("weighted contrasts", {
  expect_equal(dim(con), 2:1)
  expect_equal(as.vector(con), c(-0.25,1))
  expect_equal(dim(con1), 2:1)
  expect_equal(as.vector(con1), c(-2/3,1))
})


test_that("weighted contrasts errors", {
  dat2 <- data.frame ( group = c(rep(1,3), rep(2,2)), wgt = c(2/3, 4/3, 2, 3/8, 5/8))
  expect_error(contr.wec.weighted(x= dat2[,"group"], omitted=1,weights=rep(1,nrow(dat))), "Variable 'x' must be of class 'factor'.")
})

test_that("weighted contrasts errors", {
  expect_error(contr.wec.weighted(x= dat[,"group"], omitted=4,weights=rep(1,nrow(dat))), "Level declared to be omitted ('4') does not exist in 'x' or has zero observations.", fixed = TRUE)
})