tests, version, docs, news

DESCRIPTION

@@ -1,7 +1,7 @@
 Type: Package
 Package: datawizard
 Title: Easy Data Wrangling and Statistical Transformations
-Version: 0.8.0.4
+Version: 0.8.0.5
 Authors@R: c(
     person("Indrajeet", "Patil", , "[email protected]", role = "aut",
            comment = c(ORCID = "0000-0003-1995-6531", Twitter = "@patilindrajeets")),

NEWS.md

@@ -1,5 +1,9 @@
 # datawizard (devel)
 
+NEW FUNCTIONS
+
+* `contr.deviation()` for sum-deviation contrast coding of factors.
+
 CHANGES
 
 * `recode_into()` gains an `overwrite` argument to skip overwriting already

R/contrs.R

@@ -23,7 +23,7 @@
 #' \dontrun{
 #' data("mtcars")
 #'
-#' mtcars$cyl <- factor(mtcars$cyl)
+#' mtcars <- data_modify(mtcars, cyl = factor(cyl))
 #'
 #' c.treatment <- cbind(Intercept = 1, contrasts(mtcars$cyl))
 #' solve(c.treatment)
@@ -50,10 +50,9 @@
 #' #> 8         -1.000 0.000 1.000   # 3rd level - 1st level
 #'
 #' ## With Interactions -----------------------------------------
-#' mtcars$am <- factor(mtcars$am)
-#' contrasts(mtcars$am) <- contr.deviation
+#' mtcars <- data_modify(mtcars, am = factor(am))
+#' mtcars <- data_arrange(mtcars, select = c("cyl", "am"))
 #'
-#' mtcars <- mtcars[order(mtcars$cyl, mtcars$am),]
 #' mm <- unique(model.matrix(~ cyl * am, data = mtcars))
 #' rownames(mm) <- c("cyl4.am0", "cyl4.am1", "cyl6.am0",
 #'                   "cyl6.am1", "cyl8.am0", "cyl8.am1")

tests/testthat/test-contr.deviation.R

@@ -0,0 +1,8 @@
+test_that("contr.deviation", {
+  c.treatment <- solve(cbind(Intercept = 1, contr.treatment(3)))
+  c.sum <- solve(cbind(Intercept = 1, contr.sum(3)))
+  c.deviation <- solve(cbind(Intercept = 1, contr.deviation(3)))
+
+  expect_equal(c.deviation[1,], c.sum[1,])
+  expect_equal(c.deviation[-1,], c.treatment[-1,])
+})