# convert the data to a data.frame
dta <- data.frame(Detergent)
# run the anofa analysis
if (FALSE) {
w <- anofa( Freq ~ Temperature * M_User * Preference * Water_softness, dta)
# make a plot with all the factors
anofaPlot(w)
# ... or with just a few factors
anofaPlot(w, ~ Preference * M_User )
anofaPlot(w, ~ Temperature )
# extract simple effects
e <- emFrequencies(w, ~ M_User | Preference )
}
Run the code above in your browser using DataLab