用于R中重塑数据的3向ANOVA

Question

I have just discovered reshaping in R and am unsure of how to proceed with an ANOVA once the data is reshaped. I found this site which has the data organized in a way very similar to my own data. If I were using this hypothetical data, how would I conduct a 3-way ANOVA say between race, program and subject? Now that the subjects have been reshaped into a single column I'm having trouble seeing how to include this variable using the typical ANOVA code. Any help would be much appreciated!

Answer 1

Assuming the data are in 'long format' and 'score' is your dependent variable you could do something like:

mymodel = aov(score ~ prog + race + subj, data=l)
summary(my model)

Which in this case yields:

Df Sum Sq Mean Sq F value   Pr(>F)    
prog          1   2864    2864   31.32 2.82e-08 ***
race          1   5064    5064   55.39 2.14e-13 ***
subj          4    106      27    0.29    0.885    
Residuals   993  90780      91                     
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

nb this model contains only the main effects