This function performs analysis of variance of data from experiments using an incomplete split-plot design for the situation when blocks are complete with respect to main plot treatments and mainplots are incomplete with respect to subplot treatments
aov.ispd.cmis(obs, block, mp, sp, y)
observation numbers
block
main plot treatment
subplot treatment
response variable
Returns ANOVA table of incomplete split-plot design