# running the anofa
L <- anofa( obsfreq ~ provider * program, LandisBarrettGalvin2013)
# getting a plot
anofaPlot(L)
# the G table shows a significant interaction
summary(L)
# getting the simple effect
e <- emFrequencies(L, ~ program | provider )
## Getting some contrast by provider (i.e., on e)
f <- contrastFrequencies(e, list(
"(PBH & CBH) vs. BM"=c(1,1,-2)/2,
"PBH vs. CBH"=c(1,-1,0))
)
Run the code above in your browser using DataLab