
import spm1d




#(0) Load dataset:
dataset      = spm1d.data.uv0d.anova3onerm.NYUCaffeine()
dataset      = spm1d.data.uv0d.anova3onerm.Southampton3onerm()
y,A,B,C,SUBJ = dataset.get_data()
print( dataset )



#(1) Run ANOVA:
FF        = spm1d.stats.anova3onerm(y, A, B, C, SUBJ, equal_var=True)
FFi       = FF.inference(0.05)
print( FFi )

