Katie Surrence
2009-Oct-27 20:32 UTC
[R] anova interaction contrasts: crossing helmert and linear contrasts
I am new to statistics, R, and this list, so apologies in advance for
the errors etiquette I am certain to make (in spite of reading the
posting guide, help on
various commands, etc.). ?Any help is greatly appreciated.
Here is my data:
fatigue =
c(3,2,2,3,2,3,4,3,2,4,5,3,3,2,4,5,4,5,5,6,4,6,9,8,4,3,5,5,6,6,6,7,9,10,12,9)
n <- 3
train <- gl(3, 4*n, labels=c("6wks", "4wks",
"2wks"))
dist <-
rep(gl(4,n,labels=c("1mi","2mi","3mi","4mi")),3)
What I would like to do is test two interaction contrasts that cross linear
coefficients of the variable "dist" with the two sets of helmert
contrasts of
the variable "train".
I have tried so many things I think are wrong I won't reproduce them all
here.
Here is a representative example:
My code:
contrasts(dist) = c(-3, -1, 1, 3)
contrasts(train) = contr.helmert(3)
aov1 = aov(fatigue~dist*train)
summary(aov1, intercept = T, split=list(train=1:2), expand.split = T)
My output:
? ? ? ? ? ? ? Df Sum Sq Mean Sq ?F value ? ?Pr(>F)
(Intercept) ? ? ? 1 890.03 ?890.03 801.0250 < 2.2e-16 ***
dist ? ? ? ? ? ? ?3 ?85.19 ? 28.40 ?25.5583 1.199e-07 ***
train ? ? ? ? ? ? 2 ?88.39 ? 44.19 ?39.7750 2.403e-08 ***
?train: C1 ? ? ? 1 ?18.38 ? 18.38 ?16.5375 0.0004452 ***
?train: C2 ? ? ? 1 ?70.01 ? 70.01 ?63.0125 3.621e-08 ***
dist:train ? ? ? ?6 ?18.72 ? ?3.12 ? 2.8083 0.0325092 *
?dist:train: C1 ?3 ? 3.46 ? ?1.15 ? 1.0375 0.3938042
?dist:train: C2 ?3 ?15.26 ? ?5.09 ? 4.5792 0.0113201 *
Residuals ? ? ? ?24 ?26.67 ? ?1.11
---
Signif. codes: ?0 ?***? 0.001 ?**? 0.01 ?*? 0.05 ?.? 0.1 ? ? 1
Among the reasons this seems wrong to me is that I want 1 df
interaction contrasts. "dist:train: C1" and "dist:train: C2"
seem (to
my naive eye) to be testing respectively whether my first and second
helmert contrast of train differ significantly across levels of dist.
I want to test the significance of two contrasts of linear trend over
dist--the slope of linear trend at 6wks and 4wks against each other,
and the average of the slopes at 6wks and 4wks against the slope at
2wks. I don't know how to do this.
Thanks for any help you can offer.
