Description Usage Arguments Details Value Author(s) See Also Examples
Helps you to construct a contrast matrix and naming the contrasts
1 |
contr |
a list of contrasts (given as vectors) or a matrix with contrasts given in (k-1) rows or (k-1) columns (k is the number of levels of a factor) |
contr.names |
a list or vector of contrast names |
For factors with more than two levels, comparisons among different levels are usually of interest (also comparisons among means of different levels). Therefore you need to be able to set your contrasts manually. Because it is also important to give your contrasts meaningful names. This is especially easy with this function. You can either provide them as row.names directly in the contrast matrix or in separate vector or list.
A (k times k-1) matrix
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 | data(pea)
trt.contr <- rbind("control-sugar"=c(1, -1/4, -1/4, -1/4, -1/4),
"pure-mixed"=c(0, 1/3, 1/3, -1, 1/3),
"monosaccharides-disaccharides"=c(0,1/2,1/2,0,-1),
"gluc-fruc"=c(0,1,-1,0,0))
model <- aov(length ~ trt,
contrasts=list(trt=mancontr(contr=trt.contr)), data=pea)
### ALTERNATIVE formulations:
## contrasts in a list:
# trt.contr <- list(c( 1, -1/4, -1/4, -1/4, -1/4),
# c( 0, 1/3, 1/3, -1, 1/3),
# c( 0, 1/2, 1/2, 0, -1),
# c( 0, 1, -1, 0, 0))
## contrasts in a matrix
# trt.contr <- rbind(c( 1, -1/4, -1/4,-1/4, -1/4),
# c( 0, 1/3, 1/3, -1, 1/3),
# c( 0, 1/2, 1/2, 0, -1),
# c( 0, 1, -1, 0, 0))
## names of contrasts in a list:
# trt.contr.names=list("control-sugar", "pure-mixed",
# "monosaccharides-disaccharides", "gluc-fruc")
## names of contrasts in a vector:
# trt.contr.names=c("control-sugar", "pure-mixed",
# "monosaccharides-disaccharides", "gluc-fruc")
#model <- aov(length ~ trt,
# contrasts=list(trt=mancontr(contr=trt.contr, contr.names=trt.contr.names)),
# data=pea)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.