[R] Tukey HSD

Kurt Debono kurtdebono at hotmail.com
Thu Aug 9 23:41:23 CEST 2007


Hi,
I was wondering if you could help me:
The following are the first few lines of my data set:

subject group condition depvar
s1	c	ver	114.87
s1	c	feet	114.87
s1	c	body	114.87
s2	c	ver	73.54
s2	c	feet	64.32
s2	c	body	61.39
s3	a	ver	114.87
s3	a	feet	97.21
s3	a	body	103.31 etc.

I entered the following ANOVA command:

>dat <- read.table("mydata.txt", header=T)
>summary(aov(depvar ~ group * condition + Error(subject), data=dat))

Error: subject
          Df  Sum Sq Mean Sq F value Pr(>F)
group      1   443.3   443.3  1.0314 0.3185
Residuals 28 12035.3   429.8

Error: Within
                Df  Sum Sq Mean Sq F value   Pr(>F)
condition        2  615.82  307.91  6.6802 0.002501 **
group:condition  2   61.51   30.75  0.6672 0.517168
Residuals       56 2581.18   46.09
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
>


I cannot find a way to perform a Tukey HSD on the main effect of condition. 
since the ANOVA formula contains the 'error' command.

Could you help me please?
Kurt

_________________________________________________________________
[[trailing spam removed]]



More information about the R-help mailing list