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

`contrasts`

, `C`

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)
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

All documentation is copyright its authors; we didn't write any of that.