As with the one-way case, testing the omnibus hypothesis via two-way ANOVA is simple process in R. This tutorial will explore how R can be used to perform a two-way ANOVA to test the difference between two (or more) group means.
Tutorial FilesBefore we begin, you may want to download the sample data (.csv) used in this tutorial. Be sure to right-click and save the file to your R working directory. This dataset contains a hypothetical sample of 60 participants who are divided by gender (male and female) and treatment group (control and treatment). The values represent a scale that ranges from 1 to 5. For instance, this dataset could be conceptualized as a comparison between two professional training programs, where the control group participated the company’s longstanding program and the treatment group participated in an experimental program. The values could represent the attitudes of employees towards the training programs on a scale from 1 (poor) to 5 (excellent).
Beginning StepsTo begin, we need to read our dataset into R and store its contents in a variable.
- > #read the two-way ANOVA dataset into an R variable using the read.csv(file) function
- > dataTwoWay <- read.csv("dataset_ANOVA_TwoWay.csv")
- > #display the data
- > dataTwoWay
The first ten rows of our two-way ANOVA dataset.
Two-Way ANOVANow that our data are ready, we can conduct a two-way omnibus ANOVA test using the anova(object) function. Note that the only step necessary to add a second independent variable into our ANOVA model is to incorporate it into our lm(model, dataset) function using the * operator. Whereas our one-way model was lm(Values ~ Group), our two-way model becomes lm(Values ~ Group * Gender). As you can see from the results below, adding a second independent variable in this manner also gives us information about the interaction between our variables.
- > #use anova(object) to test the omnibus hypothesis in two-way ANOVA
- > #Are the differences between the group means for treatment and gender statistically significant?
- > #Is there a statistically significant interaction between treatment and gender?
- > anova(lm(Values ~ Group * Gender, dataTwoWay))
Our two-way ANOVA table.
The output of our ANOVA test indicates that the difference between our treatment group means is statistically significant (p < .001) and that the difference between genders is not (p = .585). However, in light of the statistically significant interaction between treatment group and gender (p = .032), we would generally elect to forgo these main effects. Subsequently, a series of follow-up procedures could be carried out to examine the simple main effects for treatment group and gender.