Description Usage Arguments Details Value Author(s) See Also Examples
View source: R/contrast_functions.R
named_contr_sum
creates sum contrasts for a factor which are named
with the levels of the factor rather than with numbers (e.g. if a factor
f1
has levels A
, B
, and C
, then rather than
creating contrast columns f11
and f12
, it creates columns
f1A
and f1B
). The absolute value of the non-zero elements
of the matrix can also be specified.
1 | named_contr_sum(x, scale = 1, return_contr = TRUE)
|
x |
An object coercible to factor or a numeric or character vector of levels. |
scale |
A positive number by which the entire contrast
matrix returned by |
return_contr |
A logical. If |
First, x
is coerced to factor, and its levels (excluding NA
)
are sorted alphabetically. If there are two unique values, and they are
equal to (ignoring case) "F" and "T", "FALSE" and "TRUE", "N" and "Y",
"NO" and "YES", or "0" and "1", then their order is reversed (this makes it
so the positive level gets the dummy coefficient rather than the negative
level, yielding a more intuitive interpretation for coefficients). Then
contr.sum
is called, and the column names of the
resulting contrast matrix are set using the character vector of unique values
(excluding the final element that gets coded as -1
for all dummy
variables). This entire matrix is then multiplied by scale
; with
the default value of 1
, this does not change the matrix; if, for
example, scale = 0.5
, then rather than each column containing values
in -1, 0, 1
, each column would contain values in -0.5, 0, 0.5
.
If return_contr = TRUE
, then this contrast matrix is
returned. If return_contr = FALSE
, then x
is converted to an
unordered factor with the named sum contrasts and returned. NA
is never
assigned as a level in the contrast matrix or in the factor returned by the
function, but NA
values in x
are not removed in the factor
returned when return_contr = FALSE
. See the examples.
If return_contr = TRUE
, a contrast matrix obtained from
contr.sum
with named columns rather than numbered
columns and deviations with magnitude scale
.
If return_contr = FALSE
, then x
is returned
as an unordered factor with the named sum contrasts applied.
Christopher D. Eager <eager.stats@gmail.com>
scaled_contr_poly
for ordered factors.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 | f <- factor(rep(c("a", "b", "c", NA), 2), levels = c("b", "c", "a"))
f <- addNA(f)
levels(f) # NA listed as factor level
contrasts(f) # NA included in contrast matrix
named_contr_sum(f) # named sum contrasts (NA dropped; levels alphabetized)
named_contr_sum(levels(f)) # same output
named_contr_sum(f, return_contr = FALSE) # factor with named sum contrasts
named_contr_sum(f, 0.5) # deviations of magniude 0.5
f <- c(TRUE, FALSE, FALSE, TRUE)
class(f) # logical
named_contr_sum(f) # TRUE gets the dummy variable
f <- named_contr_sum(f, return_contr = FALSE)
class(f) # factor
named_contr_sum(letters[1:5]) # character argument
named_contr_sum(rep(letters[1:5], 2), return_contr = FALSE) # creates factor
# ordered factors are converted to unordered factors, so use with caution
f <- factor(rep(1:3, 2), ordered = TRUE)
is.ordered(f) # TRUE
f
f <- named_contr_sum(f, return_contr = FALSE)
is.ordered(f) # FALSE
f
## Not run:
# error from stats::contr.sum because only one unique non-NA value
named_contr_sum(5)
named_contr_sum(rep(c("a", NA), 3))
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.