| RootStock | R Documentation |
In a classic experiment carried out from 1918 to 1934, growth of apple trees of six different rootstocks were compared on four measures of size. How do the measures of size vary with the type of rootstock?
A data frame with 48 observations on the following 5 variables.
rootstocka factor with levels 1 2 3 4 5 6
girth4a numeric vector: trunk girth at 4 years (mm x 100)
ext4a numeric vector: extension growth at 4 years (m)
girth15a numeric vector: trunk girth at 15 years (mm x 100)
weight15a numeric vector: weight of tree above ground at 15 years (lb x 1000)
This is a balanced, one-way MANOVA design, with n=8 trees for each rootstock.
Andrews, D. and Herzberg, A. (1985). Data: A Collection of Problems from Many Fields for the Student and Research Worker Springer-Verlag, pp. 357–360.
Rencher, A. C. (1995). Methods of Multivariate Analysis. New York: Wiley, Table 6.2
library(car)
data(RootStock)
str(RootStock)
root.mod <- lm(cbind(girth4, ext4, girth15, weight15) ~ rootstock, data=RootStock)
car::Anova(root.mod)
pairs(root.mod)
# test two orthogonal contrasts among the rootstocks
hyp <- matrix(c(2,-1,-1,-1,-1,2,
1, 0,0,0,0,-1), 2, 6, byrow=TRUE)
car::linearHypothesis(root.mod, hyp)
heplot(root.mod, hypotheses=list(Contrasts=hyp, C1=hyp[1,], C2=hyp[2,]))
heplot1d(root.mod, hypotheses=list(Contrasts=hyp, C1=hyp[1,], C2=hyp[2,]))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.