data(splitplot)
splitplot.lm = lm(RT ~ SOA + Item + Subject +
SOA:Subject + Item:Subject, data = splitplot)
splitplot.aov = anova(splitplot.lm)
quasiF.fnc(splitplot.aov["SOA","Mean Sq"],
splitplot.aov["Item:Subject", "Mean Sq"],
splitplot.aov["SOA:Subject", "Mean Sq"],
splitplot.aov["Item", "Mean Sq"],
splitplot.aov["SOA","Df"],
splitplot.aov["Item:Subject", "Df"],
splitplot.aov["SOA:Subject", "Df"],
splitplot.aov["Item", "Df"])
# much simpler is
quasiF.split.plot.fnc(splitplot)$quasiF
Run the code above in your browser using DataLab