RF ND-ANOVA: separate classes for different model types
ENH glm:  unbalanced ND-ANOVA
STY/DOC
FIX Dataset.from_r()
ENH Dataset:  add .to_r()
TEST glm:  unbalanced 3 by 2 ANOVA with Adler dataset