library(ANOFA)
options(superb.feedback = 'none') # shut down 'warnings' and 'design' interpretation messages
# Lets run the analysis
L <- anofa( obsfreq ~ vocation * gender, LightMargolin1971)
summary(L)
# a quick plot
anofaPlot(L)
# Some simple effects.
e <- emFrequencies(L, ~ gender | vocation )
summary(e)
# some contrasts:
e <- emFrequencies(L, ~ vocation | gender )
f <- contrastFrequencies(e, list(
"teacher college vs. gymnasium"=c( 0, 0, 1,-1, 0),
"vocational vs. university" = c( 0, 1, 0, 0,-1),
"another" = c( 0, 1,-1,-1,+1)/2,
"to exhaust the df" = c( 4,-1,-1,-1,-1)/4
)
)
Run the code above in your browser using DataLab