def test1(self):
df = DataFrame()
fname = 'error~subjectXtimeofdayXcourseXmodel.csv'
df.read_tbl(fname)
aov = df.anova('ERROR', wfactors=['TIMEOFDAY', 'COURSE',
'MODEL']) #,transform='windsor05')
aov.output2html(fname[:-4] + 'RESULTS.htm')
def test02(self):
"""using loftus and masson error bars"""
# a simple plot
df=DataFrame()
df.read_tbl('data/words~ageXcondition.csv')
aov = df.anova('WORDS', wfactors=['AGE','CONDITION'])
aov.plot('WORDS','AGE', seplines='CONDITION',
errorbars='ci', output_dir='output')
def test02(self):
"""using loftus and masson error bars"""
# a simple plot
df = DataFrame()
df.read_tbl('data/words~ageXcondition.csv')
aov = df.anova('WORDS', wfactors=['AGE', 'CONDITION'])
aov.plot('WORDS',
'AGE',
seplines='CONDITION',
errorbars='ci',
output_dir='output')
def test0(self):
## Within test
R = """\
ERROR ~ TIMEOFDAY * COURSE * MODEL
TESTS OF WITHIN SUBJECTS EFFECTS
Measure: ERROR
Source Type III eps df MS F Sig. et2_G Obs. SE 95% CI lambda Obs.
SS Power
=====================================================================================================================================================
TIMEOFDAY Sphericity Assumed 140.167 - 1 140.167 120.143 0.008 3.391 27 0.456 0.894 1621.929 1
Greenhouse-Geisser 140.167 1 1 140.167 120.143 0.008 3.391 27 0.456 0.894 1621.929 1
Huynh-Feldt 140.167 1 1 140.167 120.143 0.008 3.391 27 0.456 0.894 1621.929 1
Box 140.167 1 1 140.167 120.143 0.008 3.391 27 0.456 0.894 1621.929 1
-----------------------------------------------------------------------------------------------------------------------------------------------------
Error(TIMEOFDAY) Sphericity Assumed 2.333 - 2 1.167
Greenhouse-Geisser 2.333 1 2 1.167
Huynh-Feldt 2.333 1 2 1.167
Box 2.333 1 2 1.167
-----------------------------------------------------------------------------------------------------------------------------------------------------
COURSE Sphericity Assumed 56.778 - 2 28.389 1022.000 3.815e-06 1.374 18 0.056 0.109 9198.000 1
Greenhouse-Geisser 56.778 0.501 1.002 56.667 1022.000 9.664e-04 1.374 18 0.056 0.109 9198.000 1
Huynh-Feldt 56.778 0.504 1.008 56.336 1022.000 9.349e-04 1.374 18 0.056 0.109 9198.000 1
Box 56.778 0.500 1 56.778 1022.000 9.770e-04 1.374 18 0.056 0.109 9198.000 1
-----------------------------------------------------------------------------------------------------------------------------------------------------
Error(COURSE) Sphericity Assumed 0.111 - 4 0.028
Greenhouse-Geisser 0.111 0.501 2.004 0.055
Huynh-Feldt 0.111 0.504 2.016 0.055
Box 0.111 0.500 2 0.056
-----------------------------------------------------------------------------------------------------------------------------------------------------
MODEL Sphericity Assumed 51.444 - 2 25.722 92.600 4.470e-04 1.245 18 0.176 0.345 833.400 1
Greenhouse-Geisser 51.444 0.507 1.013 50.770 92.600 0.010 1.245 18 0.176 0.345 833.400 1.000
Huynh-Feldt 51.444 0.527 1.054 48.817 92.600 0.009 1.245 18 0.176 0.345 833.400 1.000
Box 51.444 0.500 1 51.444 92.600 0.011 1.245 18 0.176 0.345 833.400 1.000
-----------------------------------------------------------------------------------------------------------------------------------------------------
Error(MODEL) Sphericity Assumed 1.111 - 4 0.278
Greenhouse-Geisser 1.111 0.507 2.027 0.548
Huynh-Feldt 1.111 0.527 2.108 0.527
Box 1.111 0.500 2 0.556
-----------------------------------------------------------------------------------------------------------------------------------------------------
TIMEOFDAY * Sphericity Assumed 5.444 - 2 2.722 2.085 0.240 0.132 9 0.540 1.057 9.383 0.446
COURSE Greenhouse-Geisser 5.444 0.814 1.628 3.345 2.085 0.255 0.132 9 0.540 1.057 9.383 0.373
Huynh-Feldt 5.444 1 2 2.722 2.085 0.240 0.132 9 0.540 1.057 9.383 0.446
Box 5.444 0.500 1 5.444 2.085 0.286 0.132 9 0.540 1.057 9.383 0.244
-----------------------------------------------------------------------------------------------------------------------------------------------------
Error(TIMEOFDAY * Sphericity Assumed 5.222 - 4 1.306
COURSE) Greenhouse-Geisser 5.222 0.814 3.255 1.604
Huynh-Feldt 5.222 1 4 1.306
Box 5.222 0.500 2 2.611
-----------------------------------------------------------------------------------------------------------------------------------------------------
TIMEOFDAY * Sphericity Assumed 16.778 - 2 8.389 37.750 0.003 0.406 9 0.223 0.436 169.875 1.000
MODEL Greenhouse-Geisser 16.778 0.540 1.079 15.545 37.750 0.021 0.406 9 0.223 0.436 169.875 0.993
Huynh-Feldt 16.778 0.571 1.142 14.697 37.750 0.018 0.406 9 0.223 0.436 169.875 0.996
Box 16.778 0.500 1 16.778 37.750 0.025 0.406 9 0.223 0.436 169.875 0.985
-----------------------------------------------------------------------------------------------------------------------------------------------------
Error(TIMEOFDAY * Sphericity Assumed 0.889 - 4 0.222
MODEL) Greenhouse-Geisser 0.889 0.540 2.159 0.412
Huynh-Feldt 0.889 0.571 2.283 0.389
Box 0.889 0.500 2 0.444
-----------------------------------------------------------------------------------------------------------------------------------------------------
COURSE * Sphericity Assumed 8.778 - 4 2.194 3.762 0.052 0.212 6 0.367 0.719 11.286 0.504
MODEL Greenhouse-Geisser 8.778 0.354 1.415 6.204 3.762 0.157 0.212 6 0.367 0.719 11.286 0.223
Huynh-Feldt 8.778 0.354 1.415 6.204 3.762 0.157 0.212 6 0.367 0.719 11.286 0.223
Box 8.778 0.500 2 4.389 3.762 0.120 0.212 6 0.367 0.719 11.286 0.292
-----------------------------------------------------------------------------------------------------------------------------------------------------
Error(COURSE * Sphericity Assumed 4.667 - 8 0.583
MODEL) Greenhouse-Geisser 4.667 0.354 2.830 1.649
Huynh-Feldt 4.667 0.354 2.830 1.649
Box 4.667 0.500 4 1.167
-----------------------------------------------------------------------------------------------------------------------------------------------------
TIMEOFDAY * Sphericity Assumed 2.778 - 4 0.694 1.923 0.200 0.067 3 0.408 0.800 2.885 0.152
COURSE * Greenhouse-Geisser 2.778 0.290 1.159 2.397 1.923 0.293 0.067 3 0.408 0.800 2.885 0.087
MODEL Huynh-Feldt 2.778 0.290 1.159 2.397 1.923 0.293 0.067 3 0.408 0.800 2.885 0.087
Box 2.778 0.500 2 1.389 1.923 0.260 0.067 3 0.408 0.800 2.885 0.109
-----------------------------------------------------------------------------------------------------------------------------------------------------
Error(TIMEOFDAY * Sphericity Assumed 2.889 - 8 0.361
COURSE * Greenhouse-Geisser 2.889 0.290 2.318 1.246
MODEL) Huynh-Feldt 2.889 0.290 2.318 1.246
Box 2.889 0.500 4 0.722
TABLES OF ESTIMATED MARGINAL MEANS
Estimated Marginal Means for TIMEOFDAY
TIMEOFDAY Mean Std. Error 95% Lower Bound 95% Upper Bound
==================================================================
T1 5.778 0.457 4.882 6.674
T2 2.556 0.229 2.108 3.003
Estimated Marginal Means for COURSE
COURSE Mean Std. Error 95% Lower Bound 95% Upper Bound
===============================================================
C1 5.222 0.608 4.031 6.414
C2 4.500 0.562 3.399 5.601
C3 2.778 0.432 1.931 3.625
Estimated Marginal Means for MODEL
MODEL Mean Std. Error 95% Lower Bound 95% Upper Bound
==============================================================
M1 5.333 0.686 3.989 6.678
M2 4.222 0.558 3.129 5.315
M3 2.944 0.328 2.301 3.588
Estimated Marginal Means for TIMEOFDAY * COURSE
TIMEOFDAY COURSE Mean Std. Error 95% Lower Bound 95% Upper Bound
===========================================================================
T1 C1 7.222 0.641 5.966 8.478
T1 C2 6.111 0.790 4.564 7.659
T1 C3 4 0.577 2.868 5.132
T2 C1 3.222 0.401 2.437 4.007
T2 C2 2.889 0.261 2.378 3.400
T2 C3 1.556 0.294 0.979 2.132
Estimated Marginal Means for TIMEOFDAY * MODEL
TIMEOFDAY MODEL Mean Std. Error 95% Lower Bound 95% Upper Bound
==========================================================================
T1 M1 7.444 0.835 5.807 9.081
T1 M2 6.111 0.512 5.107 7.115
T1 M3 3.778 0.465 2.867 4.689
T2 M1 3.222 0.434 2.372 4.073
T2 M2 2.333 0.408 1.533 3.133
T2 M3 2.111 0.261 1.600 2.622
Estimated Marginal Means for COURSE * MODEL
COURSE MODEL Mean Std. Error 95% Lower Bound 95% Upper Bound
=======================================================================
C1 M1 6.667 1.085 4.540 8.794
C1 M2 5.167 1.195 2.825 7.509
C1 M3 3.833 0.601 2.656 5.011
C2 M1 6.167 1.195 3.825 8.509
C2 M2 4.167 0.792 2.614 5.720
C2 M3 3.167 0.477 2.231 4.102
C3 M1 3.167 0.872 1.457 4.877
C3 M2 3.333 0.882 1.605 5.062
C3 M3 1.833 0.307 1.231 2.436
Estimated Marginal Means for TIMEOFDAY * COURSE * MODEL
TIMEOFDAY COURSE MODEL Mean Std. Error 95% Lower Bound 95% Upper Bound
===================================================================================
T1 C1 M1 9 0.577 7.868 10.132
T1 C1 M2 7.667 0.333 7.013 8.320
T1 C1 M3 5 0.577 3.868 6.132
T1 C2 M1 8.667 0.882 6.938 10.395
T1 C2 M2 5.667 0.882 3.938 7.395
T1 C2 M3 4 0.577 2.868 5.132
T1 C3 M1 4.667 1.202 2.311 7.022
T1 C3 M2 5 0.577 3.868 6.132
T1 C3 M3 2.333 0.333 1.680 2.987
T2 C1 M1 4.333 0.333 3.680 4.987
T2 C1 M2 2.667 0.882 0.938 4.395
T2 C1 M3 2.667 0.333 2.013 3.320
T2 C2 M1 3.667 0.333 3.013 4.320
T2 C2 M2 2.667 0.333 2.013 3.320
T2 C2 M3 2.333 0.333 1.680 2.987
T2 C3 M1 1.667 0.333 1.013 2.320
T2 C3 M2 1.667 0.882 -0.062 3.395
T2 C3 M3 1.333 0.333 0.680 1.987
"""
df = DataFrame()
fname = 'data/error~subjectXtimeofdayXcourseXmodel.csv'
df.read_tbl(fname)
aov = df.anova('ERROR', wfactors=['TIMEOFDAY', 'COURSE', 'MODEL'])
## print(aov)
self.assertEqual(str(aov), R)
def test1(self):
df=DataFrame()
fname='error~subjectXtimeofdayXcourseXmodel.csv'
df.read_tbl(fname)
aov=df.anova('ERROR',wfactors=['TIMEOFDAY','COURSE','MODEL'])#,transform='windsor05')
aov.output2html(fname[:-4]+'RESULTS.htm')
def test1(self):
## Within test with Windsor transformation
R = """\
WINDSOR05_ERROR ~ TIMEOFDAY * COURSE * MODEL
TESTS OF WITHIN SUBJECTS EFFECTS
Measure: WINDSOR05_ERROR
Source Type III eps df MS F Sig. et2_G Obs. SE 95% CI lambda Obs.
SS Power
======================================================================================================================================================
TIMEOFDAY Sphericity Assumed 130.667 - 1 130.667 92.235 0.019 3.698 27 0.622 1.218 1512.000 1
Greenhouse-Geisser 130.667 1 1 130.667 92.235 0.019 3.698 27 0.622 1.218 1512.000 1
Huynh-Feldt 130.667 1 1 130.667 92.235 0.019 3.698 27 0.622 1.218 1512.000 1
Box 130.667 1 1 130.667 92.235 0.019 3.698 27 0.622 1.218 1512.000 1
------------------------------------------------------------------------------------------------------------------------------------------------------
Error(TIMEOFDAY) Sphericity Assumed 2.333 - 1.647 1.417
Greenhouse-Geisser 2.333 1 1.647 1.417
Huynh-Feldt 2.333 1 1.647 1.417
Box 2.333 1 1.647 1.417
------------------------------------------------------------------------------------------------------------------------------------------------------
COURSE Sphericity Assumed 51.370 - 2 25.685 1142.235 2.087e-05 1.454 18 0.055 0.107 12483.000 1
Greenhouse-Geisser 51.370 0.501 1.003 51.229 1142.235 0.002 1.454 18 0.055 0.107 12483.000 1
Huynh-Feldt 51.370 0.506 1.011 50.810 1142.235 0.002 1.454 18 0.055 0.107 12483.000 1
Box 51.370 0.500 1 51.370 1142.235 0.002 1.454 18 0.055 0.107 12483.000 1
------------------------------------------------------------------------------------------------------------------------------------------------------
Error(COURSE) Sphericity Assumed 0.074 - 3.294 0.022
Greenhouse-Geisser 0.074 0.501 1.652 0.045
Huynh-Feldt 0.074 0.506 1.665 0.044
Box 0.074 0.500 1.647 0.045
------------------------------------------------------------------------------------------------------------------------------------------------------
MODEL Sphericity Assumed 47.148 - 2 23.574 123.336 8.004e-04 1.334 18 0.159 0.312 1347.882 1
Greenhouse-Geisser 47.148 0.504 1.008 46.773 123.336 0.015 1.334 18 0.159 0.312 1347.882 1.000
Huynh-Feldt 47.148 0.516 1.032 45.671 123.336 0.014 1.334 18 0.159 0.312 1347.882 1.000
Box 47.148 0.500 1 47.148 123.336 0.015 1.334 18 0.159 0.312 1347.882 1.000
------------------------------------------------------------------------------------------------------------------------------------------------------
Error(MODEL) Sphericity Assumed 0.630 - 3.294 0.191
Greenhouse-Geisser 0.630 0.504 1.660 0.379
Huynh-Feldt 0.630 0.516 1.700 0.370
Box 0.630 0.500 1.647 0.382
------------------------------------------------------------------------------------------------------------------------------------------------------
TIMEOFDAY * Sphericity Assumed 5.444 - 2 2.722 1.834 0.292 0.154 9 0.627 1.229 10.023 0.410
COURSE Greenhouse-Geisser 5.444 0.720 1.441 3.779 1.834 0.313 0.154 9 0.627 1.229 10.023 0.302
Huynh-Feldt 5.444 0.933 1.867 2.916 1.834 0.296 0.154 9 0.627 1.229 10.023 0.385
Box 5.444 0.500 1 5.444 1.834 0.332 0.154 9 0.627 1.229 10.023 0.217
------------------------------------------------------------------------------------------------------------------------------------------------------
Error(TIMEOFDAY * Sphericity Assumed 4.889 - 3.294 1.484
COURSE) Greenhouse-Geisser 4.889 0.720 2.373 2.060
Huynh-Feldt 4.889 0.933 3.075 1.590
Box 4.889 0.500 1.647 2.968
------------------------------------------------------------------------------------------------------------------------------------------------------
TIMEOFDAY * Sphericity Assumed 14.333 - 2 7.167 70.824 0.002 0.406 9 0.164 0.321 387.000 1.000
MODEL Greenhouse-Geisser 14.333 0.511 1.023 14.011 70.824 0.023 0.406 9 0.164 0.321 387.000 0.998
Huynh-Feldt 14.333 0.520 1.040 13.776 70.824 0.022 0.406 9 0.164 0.321 387.000 0.998
Box 14.333 0.500 1 14.333 70.824 0.024 0.406 9 0.164 0.321 387.000 0.997
------------------------------------------------------------------------------------------------------------------------------------------------------
Error(TIMEOFDAY * Sphericity Assumed 0.333 - 3.294 0.101
MODEL) Greenhouse-Geisser 0.333 0.511 1.685 0.198
Huynh-Feldt 0.333 0.520 1.714 0.195
Box 0.333 0.500 1.647 0.202
------------------------------------------------------------------------------------------------------------------------------------------------------
COURSE * Sphericity Assumed 8.296 - 4 2.074 3.804 0.064 0.235 6 0.368 0.722 13.856 0.545
MODEL Greenhouse-Geisser 8.296 0.327 1.307 6.347 3.804 0.188 0.235 6 0.368 0.722 13.856 0.211
Huynh-Feldt 8.296 0.327 1.307 6.347 3.804 0.188 0.235 6 0.368 0.722 13.856 0.211
Box 8.296 0.500 2 4.148 3.804 0.139 0.235 6 0.368 0.722 13.856 0.304
------------------------------------------------------------------------------------------------------------------------------------------------------
Error(COURSE * Sphericity Assumed 3.593 - 6.588 0.545
MODEL) Greenhouse-Geisser 3.593 0.327 2.153 1.669
Huynh-Feldt 3.593 0.327 2.153 1.669
Box 3.593 0.500 3.294 1.091
------------------------------------------------------------------------------------------------------------------------------------------------------
TIMEOFDAY * Sphericity Assumed 2.222 - 4 0.556 1.318 0.355 0.063 3 0.458 0.898 2.400 0.125
COURSE * Greenhouse-Geisser 2.222 0.336 1.343 1.654 1.318 0.387 0.063 3 0.458 0.898 2.400 0.080
MODEL Huynh-Feldt 2.222 0.336 1.343 1.654 1.318 0.387 0.063 3 0.458 0.898 2.400 0.080
Box 2.222 0.500 2 1.111 1.318 0.380 0.063 3 0.458 0.898 2.400 0.093
------------------------------------------------------------------------------------------------------------------------------------------------------
Error(TIMEOFDAY * Sphericity Assumed 2.778 - 6.588 0.422
COURSE * Greenhouse-Geisser 2.778 0.336 2.212 1.256
MODEL) Huynh-Feldt 2.778 0.336 2.212 1.256
Box 2.778 0.500 3.294 0.843
TABLES OF ESTIMATED MARGINAL MEANS
Estimated Marginal Means for TIMEOFDAY
TIMEOFDAY Mean Std. Error 95% Lower Bound 95% Upper Bound
==================================================================
T1 5.704 0.433 4.855 6.552
T2 2.593 0.215 2.171 3.014
Estimated Marginal Means for COURSE
COURSE Mean Std. Error 95% Lower Bound 95% Upper Bound
===============================================================
C1 5.167 0.584 4.021 6.312
C2 4.444 0.532 3.403 5.486
C3 2.833 0.414 2.021 3.645
Estimated Marginal Means for MODEL
MODEL Mean Std. Error 95% Lower Bound 95% Upper Bound
==============================================================
M1 5.222 0.645 3.959 6.485
M2 4.278 0.535 3.229 5.327
M3 2.944 0.328 2.301 3.588
Estimated Marginal Means for TIMEOFDAY * COURSE
TIMEOFDAY COURSE Mean Std. Error 95% Lower Bound 95% Upper Bound
===========================================================================
T1 C1 7.111 0.588 5.959 8.263
T1 C2 6 0.726 4.576 7.424
T1 C3 4 0.577 2.868 5.132
T2 C1 3.222 0.401 2.437 4.007
T2 C2 2.889 0.261 2.378 3.400
T2 C3 1.667 0.236 1.205 2.129
Estimated Marginal Means for TIMEOFDAY * MODEL
TIMEOFDAY MODEL Mean Std. Error 95% Lower Bound 95% Upper Bound
==========================================================================
T1 M1 7.222 0.760 5.733 8.711
T1 M2 6.111 0.512 5.107 7.115
T1 M3 3.778 0.465 2.867 4.689
T2 M1 3.222 0.434 2.372 4.073
T2 M2 2.444 0.338 1.782 3.107
T2 M3 2.111 0.261 1.600 2.622
Estimated Marginal Means for COURSE * MODEL
COURSE MODEL Mean Std. Error 95% Lower Bound 95% Upper Bound
=======================================================================
C1 M1 6.500 0.992 4.556 8.444
C1 M2 5.167 1.195 2.825 7.509
C1 M3 3.833 0.601 2.656 5.011
C2 M1 6 1.095 3.853 8.147
C2 M2 4.167 0.792 2.614 5.720
C2 M3 3.167 0.477 2.231 4.102
C3 M1 3.167 0.872 1.457 4.877
C3 M2 3.500 0.764 2.003 4.997
C3 M3 1.833 0.307 1.231 2.436
Estimated Marginal Means for TIMEOFDAY * COURSE * MODEL
TIMEOFDAY COURSE MODEL Mean Std. Error 95% Lower Bound 95% Upper Bound
===================================================================================
T1 C1 M1 8.667 0.333 8.013 9.320
T1 C1 M2 7.667 0.333 7.013 8.320
T1 C1 M3 5 0.577 3.868 6.132
T1 C2 M1 8.333 0.667 7.027 9.640
T1 C2 M2 5.667 0.882 3.938 7.395
T1 C2 M3 4 0.577 2.868 5.132
T1 C3 M1 4.667 1.202 2.311 7.022
T1 C3 M2 5 0.577 3.868 6.132
T1 C3 M3 2.333 0.333 1.680 2.987
T2 C1 M1 4.333 0.333 3.680 4.987
T2 C1 M2 2.667 0.882 0.938 4.395
T2 C1 M3 2.667 0.333 2.013 3.320
T2 C2 M1 3.667 0.333 3.013 4.320
T2 C2 M2 2.667 0.333 2.013 3.320
T2 C2 M3 2.333 0.333 1.680 2.987
T2 C3 M1 1.667 0.333 1.013 2.320
T2 C3 M2 2 0.577 0.868 3.132
T2 C3 M3 1.333 0.333 0.680 1.987
"""
df=DataFrame()
fname='data/error~subjectXtimeofdayXcourseXmodel.csv'
df.read_tbl(fname)
aov=df.anova('ERROR',wfactors=['TIMEOFDAY','COURSE','MODEL'],transform='windsor05')
## print(aov)
self.assertEqual(str(aov),R)
def test1(self):
## Within test with Windsor transformation
R = """\
WINDSOR05_ERROR ~ TIMEOFDAY * COURSE * MODEL
TESTS OF WITHIN SUBJECTS EFFECTS
Measure: WINDSOR05_ERROR
Source Type III eps df MS F Sig. et2_G Obs. SE 95% CI lambda Obs.
SS Power
======================================================================================================================================================
TIMEOFDAY Sphericity Assumed 130.667 - 1 130.667 92.235 0.019 3.698 27 0.622 1.218 1512.000 1
Greenhouse-Geisser 130.667 1 1 130.667 92.235 0.019 3.698 27 0.622 1.218 1512.000 1
Huynh-Feldt 130.667 1 1 130.667 92.235 0.019 3.698 27 0.622 1.218 1512.000 1
Box 130.667 1 1 130.667 92.235 0.019 3.698 27 0.622 1.218 1512.000 1
------------------------------------------------------------------------------------------------------------------------------------------------------
Error(TIMEOFDAY) Sphericity Assumed 2.333 - 1.647 1.417
Greenhouse-Geisser 2.333 1 1.647 1.417
Huynh-Feldt 2.333 1 1.647 1.417
Box 2.333 1 1.647 1.417
------------------------------------------------------------------------------------------------------------------------------------------------------
COURSE Sphericity Assumed 51.370 - 2 25.685 1142.235 2.087e-05 1.454 18 0.055 0.107 12483.000 1
Greenhouse-Geisser 51.370 0.501 1.003 51.229 1142.235 0.002 1.454 18 0.055 0.107 12483.000 1
Huynh-Feldt 51.370 0.506 1.011 50.810 1142.235 0.002 1.454 18 0.055 0.107 12483.000 1
Box 51.370 0.500 1 51.370 1142.235 0.002 1.454 18 0.055 0.107 12483.000 1
------------------------------------------------------------------------------------------------------------------------------------------------------
Error(COURSE) Sphericity Assumed 0.074 - 3.294 0.022
Greenhouse-Geisser 0.074 0.501 1.652 0.045
Huynh-Feldt 0.074 0.506 1.665 0.044
Box 0.074 0.500 1.647 0.045
------------------------------------------------------------------------------------------------------------------------------------------------------
MODEL Sphericity Assumed 47.148 - 2 23.574 123.336 8.004e-04 1.334 18 0.159 0.312 1347.882 1
Greenhouse-Geisser 47.148 0.504 1.008 46.773 123.336 0.015 1.334 18 0.159 0.312 1347.882 1.000
Huynh-Feldt 47.148 0.516 1.032 45.671 123.336 0.014 1.334 18 0.159 0.312 1347.882 1.000
Box 47.148 0.500 1 47.148 123.336 0.015 1.334 18 0.159 0.312 1347.882 1.000
------------------------------------------------------------------------------------------------------------------------------------------------------
Error(MODEL) Sphericity Assumed 0.630 - 3.294 0.191
Greenhouse-Geisser 0.630 0.504 1.660 0.379
Huynh-Feldt 0.630 0.516 1.700 0.370
Box 0.630 0.500 1.647 0.382
------------------------------------------------------------------------------------------------------------------------------------------------------
TIMEOFDAY * Sphericity Assumed 5.444 - 2 2.722 1.834 0.292 0.154 9 0.627 1.229 10.023 0.410
COURSE Greenhouse-Geisser 5.444 0.720 1.441 3.779 1.834 0.313 0.154 9 0.627 1.229 10.023 0.302
Huynh-Feldt 5.444 0.933 1.867 2.916 1.834 0.296 0.154 9 0.627 1.229 10.023 0.385
Box 5.444 0.500 1 5.444 1.834 0.332 0.154 9 0.627 1.229 10.023 0.217
------------------------------------------------------------------------------------------------------------------------------------------------------
Error(TIMEOFDAY * Sphericity Assumed 4.889 - 3.294 1.484
COURSE) Greenhouse-Geisser 4.889 0.720 2.373 2.060
Huynh-Feldt 4.889 0.933 3.075 1.590
Box 4.889 0.500 1.647 2.968
------------------------------------------------------------------------------------------------------------------------------------------------------
TIMEOFDAY * Sphericity Assumed 14.333 - 2 7.167 70.824 0.002 0.406 9 0.164 0.321 387.000 1.000
MODEL Greenhouse-Geisser 14.333 0.511 1.023 14.011 70.824 0.023 0.406 9 0.164 0.321 387.000 0.998
Huynh-Feldt 14.333 0.520 1.040 13.776 70.824 0.022 0.406 9 0.164 0.321 387.000 0.998
Box 14.333 0.500 1 14.333 70.824 0.024 0.406 9 0.164 0.321 387.000 0.997
------------------------------------------------------------------------------------------------------------------------------------------------------
Error(TIMEOFDAY * Sphericity Assumed 0.333 - 3.294 0.101
MODEL) Greenhouse-Geisser 0.333 0.511 1.685 0.198
Huynh-Feldt 0.333 0.520 1.714 0.195
Box 0.333 0.500 1.647 0.202
------------------------------------------------------------------------------------------------------------------------------------------------------
COURSE * Sphericity Assumed 8.296 - 4 2.074 3.804 0.064 0.235 6 0.368 0.722 13.856 0.545
MODEL Greenhouse-Geisser 8.296 0.327 1.307 6.347 3.804 0.188 0.235 6 0.368 0.722 13.856 0.211
Huynh-Feldt 8.296 0.327 1.307 6.347 3.804 0.188 0.235 6 0.368 0.722 13.856 0.211
Box 8.296 0.500 2 4.148 3.804 0.139 0.235 6 0.368 0.722 13.856 0.304
------------------------------------------------------------------------------------------------------------------------------------------------------
Error(COURSE * Sphericity Assumed 3.593 - 6.588 0.545
MODEL) Greenhouse-Geisser 3.593 0.327 2.153 1.669
Huynh-Feldt 3.593 0.327 2.153 1.669
Box 3.593 0.500 3.294 1.091
------------------------------------------------------------------------------------------------------------------------------------------------------
TIMEOFDAY * Sphericity Assumed 2.222 - 4 0.556 1.318 0.355 0.063 3 0.458 0.898 2.400 0.125
COURSE * Greenhouse-Geisser 2.222 0.336 1.343 1.654 1.318 0.387 0.063 3 0.458 0.898 2.400 0.080
MODEL Huynh-Feldt 2.222 0.336 1.343 1.654 1.318 0.387 0.063 3 0.458 0.898 2.400 0.080
Box 2.222 0.500 2 1.111 1.318 0.380 0.063 3 0.458 0.898 2.400 0.093
------------------------------------------------------------------------------------------------------------------------------------------------------
Error(TIMEOFDAY * Sphericity Assumed 2.778 - 6.588 0.422
COURSE * Greenhouse-Geisser 2.778 0.336 2.212 1.256
MODEL) Huynh-Feldt 2.778 0.336 2.212 1.256
Box 2.778 0.500 3.294 0.843
TABLES OF ESTIMATED MARGINAL MEANS
Estimated Marginal Means for TIMEOFDAY
TIMEOFDAY Mean Std. Error 95% Lower Bound 95% Upper Bound
==================================================================
T1 5.704 0.433 4.855 6.552
T2 2.593 0.215 2.171 3.014
Estimated Marginal Means for COURSE
COURSE Mean Std. Error 95% Lower Bound 95% Upper Bound
===============================================================
C1 5.167 0.584 4.021 6.312
C2 4.444 0.532 3.403 5.486
C3 2.833 0.414 2.021 3.645
Estimated Marginal Means for MODEL
MODEL Mean Std. Error 95% Lower Bound 95% Upper Bound
==============================================================
M1 5.222 0.645 3.959 6.485
M2 4.278 0.535 3.229 5.327
M3 2.944 0.328 2.301 3.588
Estimated Marginal Means for TIMEOFDAY * COURSE
TIMEOFDAY COURSE Mean Std. Error 95% Lower Bound 95% Upper Bound
===========================================================================
T1 C1 7.111 0.588 5.959 8.263
T1 C2 6 0.726 4.576 7.424
T1 C3 4 0.577 2.868 5.132
T2 C1 3.222 0.401 2.437 4.007
T2 C2 2.889 0.261 2.378 3.400
T2 C3 1.667 0.236 1.205 2.129
Estimated Marginal Means for TIMEOFDAY * MODEL
TIMEOFDAY MODEL Mean Std. Error 95% Lower Bound 95% Upper Bound
==========================================================================
T1 M1 7.222 0.760 5.733 8.711
T1 M2 6.111 0.512 5.107 7.115
T1 M3 3.778 0.465 2.867 4.689
T2 M1 3.222 0.434 2.372 4.073
T2 M2 2.444 0.338 1.782 3.107
T2 M3 2.111 0.261 1.600 2.622
Estimated Marginal Means for COURSE * MODEL
COURSE MODEL Mean Std. Error 95% Lower Bound 95% Upper Bound
=======================================================================
C1 M1 6.500 0.992 4.556 8.444
C1 M2 5.167 1.195 2.825 7.509
C1 M3 3.833 0.601 2.656 5.011
C2 M1 6 1.095 3.853 8.147
C2 M2 4.167 0.792 2.614 5.720
C2 M3 3.167 0.477 2.231 4.102
C3 M1 3.167 0.872 1.457 4.877
C3 M2 3.500 0.764 2.003 4.997
C3 M3 1.833 0.307 1.231 2.436
Estimated Marginal Means for TIMEOFDAY * COURSE * MODEL
TIMEOFDAY COURSE MODEL Mean Std. Error 95% Lower Bound 95% Upper Bound
===================================================================================
T1 C1 M1 8.667 0.333 8.013 9.320
T1 C1 M2 7.667 0.333 7.013 8.320
T1 C1 M3 5 0.577 3.868 6.132
T1 C2 M1 8.333 0.667 7.027 9.640
T1 C2 M2 5.667 0.882 3.938 7.395
T1 C2 M3 4 0.577 2.868 5.132
T1 C3 M1 4.667 1.202 2.311 7.022
T1 C3 M2 5 0.577 3.868 6.132
T1 C3 M3 2.333 0.333 1.680 2.987
T2 C1 M1 4.333 0.333 3.680 4.987
T2 C1 M2 2.667 0.882 0.938 4.395
T2 C1 M3 2.667 0.333 2.013 3.320
T2 C2 M1 3.667 0.333 3.013 4.320
T2 C2 M2 2.667 0.333 2.013 3.320
T2 C2 M3 2.333 0.333 1.680 2.987
T2 C3 M1 1.667 0.333 1.013 2.320
T2 C3 M2 2 0.577 0.868 3.132
T2 C3 M3 1.333 0.333 0.680 1.987
"""
df = DataFrame()
fname = 'data/error~subjectXtimeofdayXcourseXmodel.csv'
df.read_tbl(fname)
aov = df.anova('ERROR',
wfactors=['TIMEOFDAY', 'COURSE', 'MODEL'],
transform='windsor05')
## print(aov)
self.assertEqual(str(aov), R)
def test3(self):
## Mixed Between/Within test
R = """\
SUPPRESSION ~ CYCLE * PHASE * GROUP
TESTS OF BETWEEN-SUBJECTS EFFECTS
Measure: SUPPRESSION
Source Type III df MS F Sig. et2_G Obs. SE 95% CI lambda Obs.
SS Power
=========================================================================================================
Between Subjects 2.034 23
GROUP 0.462 2 0.231 3.083 0.067 0.150 8 0.103 0.201 2.349 0.229
=========================================================================================================
Error 1.572 21 0.075
TESTS OF WITHIN SUBJECTS EFFECTS
Measure: SUPPRESSION
Source Type III eps df MS F Sig. et2_G Obs. SE 95% CI lambda Obs.
SS Power
======================================================================================================================================================
CYCLE Sphericity Assumed 0.273 - 3 0.091 12.027 2.535e-06 0.094 48 0.013 0.025 27.491 0.995
Greenhouse-Geisser 0.273 0.668 2.005 0.136 12.027 7.281e-05 0.094 48 0.013 0.025 27.491 0.967
Huynh-Feldt 0.273 0.668 2.005 0.136 12.027 7.281e-05 0.094 48 0.013 0.025 27.491 0.967
Box 0.273 0.333 1 0.273 12.027 0.002 0.094 48 0.013 0.025 27.491 0.823
------------------------------------------------------------------------------------------------------------------------------------------------------
CYCLE * GROUP Sphericity Assumed 0.105 - 6 0.017 2.309 0.044 0.038 16 0.022 0.043 3.519 0.218
Greenhouse-Geisser 0.105 0.572 3.434 0.030 2.309 0.085 0.038 16 0.022 0.043 3.519 0.167
Huynh-Feldt 0.105 0.572 3.434 0.030 2.309 0.085 0.038 16 0.022 0.043 3.519 0.167
Box 0.105 0.167 1 0.105 2.309 0.158 0.038 16 0.022 0.043 3.519 0.107
------------------------------------------------------------------------------------------------------------------------------------------------------
Error(CYCLE) Sphericity Assumed 0.476 - 63 0.008
Greenhouse-Geisser 0.476 0.572 36.061 0.013
Huynh-Feldt 0.476 0.572 36.061 0.013
Box 0.476 0.167 10.500 0.045
------------------------------------------------------------------------------------------------------------------------------------------------------
PHASE Sphericity Assumed 1.170 - 1 1.170 129.855 1.878e-10 0.308 96 0.010 0.020 593.625 1
Greenhouse-Geisser 1.170 1 1 1.170 129.855 1.878e-10 0.308 96 0.010 0.020 593.625 1
Huynh-Feldt 1.170 1 1 1.170 129.855 1.878e-10 0.308 96 0.010 0.020 593.625 1
Box 1.170 1 1 1.170 129.855 1.878e-10 0.308 96 0.010 0.020 593.625 1
------------------------------------------------------------------------------------------------------------------------------------------------------
PHASE * GROUP Sphericity Assumed 0.405 - 2 0.203 22.493 6.012e-06 0.134 32 0.018 0.035 68.551 1.000
Greenhouse-Geisser 0.405 0.820 1.641 0.247 22.493 3.263e-05 0.134 32 0.018 0.035 68.551 1.000
Huynh-Feldt 0.405 0.820 1.641 0.247 22.493 3.263e-05 0.134 32 0.018 0.035 68.551 1.000
Box 0.405 0.500 1 0.405 22.493 6.896e-04 0.134 32 0.018 0.035 68.551 1.000
------------------------------------------------------------------------------------------------------------------------------------------------------
Error(PHASE) Sphericity Assumed 0.189 - 21 0.009
Greenhouse-Geisser 0.189 0.820 17.228 0.011
Huynh-Feldt 0.189 0.820 17.228 0.011
Box 0.189 0.500 10.500 0.018
------------------------------------------------------------------------------------------------------------------------------------------------------
CYCLE * Sphericity Assumed 0.074 - 3 0.025 4.035 0.011 0.027 24 0.016 0.032 4.612 0.386
PHASE Greenhouse-Geisser 0.074 0.686 2.057 0.036 4.035 0.024 0.027 24 0.016 0.032 4.612 0.313
Huynh-Feldt 0.074 0.686 2.057 0.036 4.035 0.024 0.027 24 0.016 0.032 4.612 0.313
Box 0.074 0.333 1 0.074 4.035 0.058 0.027 24 0.016 0.032 4.612 0.220
------------------------------------------------------------------------------------------------------------------------------------------------------
CYCLE * PHASE * GROUP Sphericity Assumed 0.127 - 6 0.021 3.466 0.005 0.046 8 0.028 0.055 2.641 0.169
Greenhouse-Geisser 0.127 0.605 3.629 0.035 3.466 0.019 0.046 8 0.028 0.055 2.641 0.137
Huynh-Feldt 0.127 1 6 0.021 3.466 0.005 0.046 8 0.028 0.055 2.641 0.169
Box 0.127 0.167 1 0.127 3.466 0.091 0.046 8 0.028 0.055 2.641 0.093
------------------------------------------------------------------------------------------------------------------------------------------------------
Error(CYCLE * Sphericity Assumed 0.386 - 63 0.006
PHASE) Greenhouse-Geisser 0.386 0.605 38.105 0.010
Huynh-Feldt 0.386 1 63 0.006
Box 0.386 0.167 10.500 0.037
TABLES OF ESTIMATED MARGINAL MEANS
Estimated Marginal Means for CYCLE
CYCLE Mean Std. Error 95% Lower Bound 95% Upper Bound
==============================================================
1 0.220 0.022 0.177 0.262
2 0.306 0.022 0.263 0.349
3 0.307 0.024 0.259 0.354
4 0.308 0.026 0.257 0.359
Estimated Marginal Means for PHASE
PHASE Mean Std. Error 95% Lower Bound 95% Upper Bound
==============================================================
I 0.207 0.014 0.180 0.234
II 0.363 0.016 0.332 0.394
Estimated Marginal Means for GROUP
GROUP Mean Std. Error 95% Lower Bound 95% Upper Bound
==============================================================
AA 0.222 0.018 0.187 0.256
AB 0.292 0.022 0.250 0.334
LAB 0.341 0.020 0.302 0.381
Estimated Marginal Means for CYCLE * PHASE
CYCLE PHASE Mean Std. Error 95% Lower Bound 95% Upper Bound
======================================================================
1 I 0.173 0.028 0.119 0.228
1 II 0.266 0.031 0.206 0.326
2 I 0.224 0.026 0.173 0.275
2 II 0.387 0.027 0.335 0.439
3 I 0.223 0.027 0.170 0.276
3 II 0.391 0.032 0.327 0.454
4 I 0.207 0.031 0.146 0.269
4 II 0.408 0.030 0.350 0.466
Estimated Marginal Means for CYCLE * GROUP
CYCLE GROUP Mean Std. Error 95% Lower Bound 95% Upper Bound
======================================================================
1 AA 0.193 0.046 0.104 0.282
1 AB 0.177 0.025 0.128 0.225
1 LAB 0.289 0.035 0.221 0.358
2 AA 0.253 0.033 0.187 0.318
2 AB 0.323 0.035 0.254 0.392
2 LAB 0.341 0.044 0.256 0.427
3 AA 0.219 0.036 0.149 0.289
3 AB 0.331 0.039 0.255 0.406
3 LAB 0.371 0.044 0.285 0.456
4 AA 0.223 0.026 0.172 0.273
4 AB 0.337 0.057 0.225 0.449
4 LAB 0.364 0.040 0.287 0.442
Estimated Marginal Means for PHASE * GROUP
PHASE GROUP Mean Std. Error 95% Lower Bound 95% Upper Bound
======================================================================
I AA 0.209 0.024 0.162 0.255
I AB 0.179 0.023 0.134 0.225
I LAB 0.233 0.026 0.183 0.283
II AA 0.235 0.026 0.184 0.286
II AB 0.404 0.023 0.359 0.450
II LAB 0.450 0.016 0.419 0.481
Estimated Marginal Means for CYCLE * PHASE * GROUP
CYCLE PHASE GROUP Mean Std. Error 95% Lower Bound 95% Upper Bound
==============================================================================
1 I AA 0.177 0.060 0.060 0.295
1 I AB 0.126 0.036 0.056 0.196
1 I LAB 0.216 0.047 0.124 0.309
1 II AA 0.209 0.072 0.067 0.351
1 II AB 0.228 0.025 0.179 0.276
1 II LAB 0.363 0.038 0.288 0.437
2 I AA 0.224 0.039 0.148 0.300
2 I AB 0.235 0.049 0.139 0.331
2 I LAB 0.214 0.053 0.110 0.317
2 II AA 0.281 0.055 0.173 0.389
2 II AB 0.411 0.026 0.360 0.463
2 II LAB 0.469 0.026 0.417 0.520
3 I AA 0.231 0.057 0.120 0.342
3 I AB 0.200 0.031 0.140 0.260
3 I LAB 0.238 0.054 0.133 0.342
3 II AA 0.208 0.047 0.115 0.300
3 II AB 0.461 0.024 0.414 0.508
3 II LAB 0.504 0.016 0.472 0.535
4 I AA 0.203 0.039 0.126 0.279
4 I AB 0.156 0.062 0.036 0.277
4 I LAB 0.264 0.058 0.149 0.378
4 II AA 0.242 0.034 0.176 0.309
4 II AB 0.517 0.031 0.457 0.578
4 II LAB 0.465 0.020 0.425 0.505
"""
df = DataFrame()
fname = 'data/suppression~subjectXgroupXcycleXphase.csv'
df.read_tbl(fname)
df['SUPPRESSION'] = [.01 * x for x in df['SUPPRESSION']]
aov = df.anova('SUPPRESSION',
wfactors=['CYCLE', 'PHASE'],
bfactors=['GROUP'])
## print(aov)
self.assertEqual(str(aov), R)
mse = aov[ss][-1] / aov['df'][-1]
aov['omega_sq'] = 'NaN'
aov['omega_sq'] = (aov[:-1][ss] -
(aov[:-1]['df'] * mse)) / (sum(aov[ss]) + mse)
return aov
eta_squared(aov_table, 'SS')
omega_squared(aov_table, 'SS')
print(aov_table)
# a short cut
formula = 'len ~ C(supp) + C(dose) + C(supp):C(dose)'
model = ols(formula, data).fit()
aov_table = anova_lm(model, typ=2)
eta_squared(aov_table, 'sum_sq')
omega_squared(aov_table, 'sum_sq')
print(aov_table)
# qqplot
res = model.resid
fig = sm.api.qqplot(res, line='s')
plt.show()
# use library pyvttbl, library may need to be installed
from pyvttbl import DataFrame
df = DataFrame()
df.read_tbl(datafile)
df['id'] = range(len(df['len']))
print(df.anova('len', sub='id', bfactors=['supp', 'dose']))
def test0(self):
## Within test
R = """\
ERROR ~ TIMEOFDAY * COURSE * MODEL
TESTS OF WITHIN SUBJECTS EFFECTS
Measure: ERROR
Source Type III eps df MS F Sig. et2_G Obs. SE 95% CI lambda Obs.
SS Power
=====================================================================================================================================================
TIMEOFDAY Sphericity Assumed 140.167 - 1 140.167 120.143 0.008 3.391 27 0.456 0.894 1621.929 1
Greenhouse-Geisser 140.167 1 1 140.167 120.143 0.008 3.391 27 0.456 0.894 1621.929 1
Huynh-Feldt 140.167 1 1 140.167 120.143 0.008 3.391 27 0.456 0.894 1621.929 1
Box 140.167 1 1 140.167 120.143 0.008 3.391 27 0.456 0.894 1621.929 1
-----------------------------------------------------------------------------------------------------------------------------------------------------
Error(TIMEOFDAY) Sphericity Assumed 2.333 - 2 1.167
Greenhouse-Geisser 2.333 1 2 1.167
Huynh-Feldt 2.333 1 2 1.167
Box 2.333 1 2 1.167
-----------------------------------------------------------------------------------------------------------------------------------------------------
COURSE Sphericity Assumed 56.778 - 2 28.389 1022.000 3.815e-06 1.374 18 0.056 0.109 9198.000 1
Greenhouse-Geisser 56.778 0.501 1.002 56.667 1022.000 9.664e-04 1.374 18 0.056 0.109 9198.000 1
Huynh-Feldt 56.778 0.504 1.008 56.336 1022.000 9.349e-04 1.374 18 0.056 0.109 9198.000 1
Box 56.778 0.500 1 56.778 1022.000 9.770e-04 1.374 18 0.056 0.109 9198.000 1
-----------------------------------------------------------------------------------------------------------------------------------------------------
Error(COURSE) Sphericity Assumed 0.111 - 4 0.028
Greenhouse-Geisser 0.111 0.501 2.004 0.055
Huynh-Feldt 0.111 0.504 2.016 0.055
Box 0.111 0.500 2 0.056
-----------------------------------------------------------------------------------------------------------------------------------------------------
MODEL Sphericity Assumed 51.444 - 2 25.722 92.600 4.470e-04 1.245 18 0.176 0.345 833.400 1
Greenhouse-Geisser 51.444 0.507 1.013 50.770 92.600 0.010 1.245 18 0.176 0.345 833.400 1.000
Huynh-Feldt 51.444 0.527 1.054 48.817 92.600 0.009 1.245 18 0.176 0.345 833.400 1.000
Box 51.444 0.500 1 51.444 92.600 0.011 1.245 18 0.176 0.345 833.400 1.000
-----------------------------------------------------------------------------------------------------------------------------------------------------
Error(MODEL) Sphericity Assumed 1.111 - 4 0.278
Greenhouse-Geisser 1.111 0.507 2.027 0.548
Huynh-Feldt 1.111 0.527 2.108 0.527
Box 1.111 0.500 2 0.556
-----------------------------------------------------------------------------------------------------------------------------------------------------
TIMEOFDAY * Sphericity Assumed 5.444 - 2 2.722 2.085 0.240 0.132 9 0.540 1.057 9.383 0.446
COURSE Greenhouse-Geisser 5.444 0.814 1.628 3.345 2.085 0.255 0.132 9 0.540 1.057 9.383 0.373
Huynh-Feldt 5.444 1 2 2.722 2.085 0.240 0.132 9 0.540 1.057 9.383 0.446
Box 5.444 0.500 1 5.444 2.085 0.286 0.132 9 0.540 1.057 9.383 0.244
-----------------------------------------------------------------------------------------------------------------------------------------------------
Error(TIMEOFDAY * Sphericity Assumed 5.222 - 4 1.306
COURSE) Greenhouse-Geisser 5.222 0.814 3.255 1.604
Huynh-Feldt 5.222 1 4 1.306
Box 5.222 0.500 2 2.611
-----------------------------------------------------------------------------------------------------------------------------------------------------
TIMEOFDAY * Sphericity Assumed 16.778 - 2 8.389 37.750 0.003 0.406 9 0.223 0.436 169.875 1.000
MODEL Greenhouse-Geisser 16.778 0.540 1.079 15.545 37.750 0.021 0.406 9 0.223 0.436 169.875 0.993
Huynh-Feldt 16.778 0.571 1.142 14.697 37.750 0.018 0.406 9 0.223 0.436 169.875 0.996
Box 16.778 0.500 1 16.778 37.750 0.025 0.406 9 0.223 0.436 169.875 0.985
-----------------------------------------------------------------------------------------------------------------------------------------------------
Error(TIMEOFDAY * Sphericity Assumed 0.889 - 4 0.222
MODEL) Greenhouse-Geisser 0.889 0.540 2.159 0.412
Huynh-Feldt 0.889 0.571 2.283 0.389
Box 0.889 0.500 2 0.444
-----------------------------------------------------------------------------------------------------------------------------------------------------
COURSE * Sphericity Assumed 8.778 - 4 2.194 3.762 0.052 0.212 6 0.367 0.719 11.286 0.504
MODEL Greenhouse-Geisser 8.778 0.354 1.415 6.204 3.762 0.157 0.212 6 0.367 0.719 11.286 0.223
Huynh-Feldt 8.778 0.354 1.415 6.204 3.762 0.157 0.212 6 0.367 0.719 11.286 0.223
Box 8.778 0.500 2 4.389 3.762 0.120 0.212 6 0.367 0.719 11.286 0.292
-----------------------------------------------------------------------------------------------------------------------------------------------------
Error(COURSE * Sphericity Assumed 4.667 - 8 0.583
MODEL) Greenhouse-Geisser 4.667 0.354 2.830 1.649
Huynh-Feldt 4.667 0.354 2.830 1.649
Box 4.667 0.500 4 1.167
-----------------------------------------------------------------------------------------------------------------------------------------------------
TIMEOFDAY * Sphericity Assumed 2.778 - 4 0.694 1.923 0.200 0.067 3 0.408 0.800 2.885 0.152
COURSE * Greenhouse-Geisser 2.778 0.290 1.159 2.397 1.923 0.293 0.067 3 0.408 0.800 2.885 0.087
MODEL Huynh-Feldt 2.778 0.290 1.159 2.397 1.923 0.293 0.067 3 0.408 0.800 2.885 0.087
Box 2.778 0.500 2 1.389 1.923 0.260 0.067 3 0.408 0.800 2.885 0.109
-----------------------------------------------------------------------------------------------------------------------------------------------------
Error(TIMEOFDAY * Sphericity Assumed 2.889 - 8 0.361
COURSE * Greenhouse-Geisser 2.889 0.290 2.318 1.246
MODEL) Huynh-Feldt 2.889 0.290 2.318 1.246
Box 2.889 0.500 4 0.722
TABLES OF ESTIMATED MARGINAL MEANS
Estimated Marginal Means for TIMEOFDAY
TIMEOFDAY Mean Std. Error 95% Lower Bound 95% Upper Bound
==================================================================
T1 5.778 0.457 4.882 6.674
T2 2.556 0.229 2.108 3.003
Estimated Marginal Means for COURSE
COURSE Mean Std. Error 95% Lower Bound 95% Upper Bound
===============================================================
C1 5.222 0.608 4.031 6.414
C2 4.500 0.562 3.399 5.601
C3 2.778 0.432 1.931 3.625
Estimated Marginal Means for MODEL
MODEL Mean Std. Error 95% Lower Bound 95% Upper Bound
==============================================================
M1 5.333 0.686 3.989 6.678
M2 4.222 0.558 3.129 5.315
M3 2.944 0.328 2.301 3.588
Estimated Marginal Means for TIMEOFDAY * COURSE
TIMEOFDAY COURSE Mean Std. Error 95% Lower Bound 95% Upper Bound
===========================================================================
T1 C1 7.222 0.641 5.966 8.478
T1 C2 6.111 0.790 4.564 7.659
T1 C3 4 0.577 2.868 5.132
T2 C1 3.222 0.401 2.437 4.007
T2 C2 2.889 0.261 2.378 3.400
T2 C3 1.556 0.294 0.979 2.132
Estimated Marginal Means for TIMEOFDAY * MODEL
TIMEOFDAY MODEL Mean Std. Error 95% Lower Bound 95% Upper Bound
==========================================================================
T1 M1 7.444 0.835 5.807 9.081
T1 M2 6.111 0.512 5.107 7.115
T1 M3 3.778 0.465 2.867 4.689
T2 M1 3.222 0.434 2.372 4.073
T2 M2 2.333 0.408 1.533 3.133
T2 M3 2.111 0.261 1.600 2.622
Estimated Marginal Means for COURSE * MODEL
COURSE MODEL Mean Std. Error 95% Lower Bound 95% Upper Bound
=======================================================================
C1 M1 6.667 1.085 4.540 8.794
C1 M2 5.167 1.195 2.825 7.509
C1 M3 3.833 0.601 2.656 5.011
C2 M1 6.167 1.195 3.825 8.509
C2 M2 4.167 0.792 2.614 5.720
C2 M3 3.167 0.477 2.231 4.102
C3 M1 3.167 0.872 1.457 4.877
C3 M2 3.333 0.882 1.605 5.062
C3 M3 1.833 0.307 1.231 2.436
Estimated Marginal Means for TIMEOFDAY * COURSE * MODEL
TIMEOFDAY COURSE MODEL Mean Std. Error 95% Lower Bound 95% Upper Bound
===================================================================================
T1 C1 M1 9 0.577 7.868 10.132
T1 C1 M2 7.667 0.333 7.013 8.320
T1 C1 M3 5 0.577 3.868 6.132
T1 C2 M1 8.667 0.882 6.938 10.395
T1 C2 M2 5.667 0.882 3.938 7.395
T1 C2 M3 4 0.577 2.868 5.132
T1 C3 M1 4.667 1.202 2.311 7.022
T1 C3 M2 5 0.577 3.868 6.132
T1 C3 M3 2.333 0.333 1.680 2.987
T2 C1 M1 4.333 0.333 3.680 4.987
T2 C1 M2 2.667 0.882 0.938 4.395
T2 C1 M3 2.667 0.333 2.013 3.320
T2 C2 M1 3.667 0.333 3.013 4.320
T2 C2 M2 2.667 0.333 2.013 3.320
T2 C2 M3 2.333 0.333 1.680 2.987
T2 C3 M1 1.667 0.333 1.013 2.320
T2 C3 M2 1.667 0.882 -0.062 3.395
T2 C3 M3 1.333 0.333 0.680 1.987
"""
df=DataFrame()
fname='data/error~subjectXtimeofdayXcourseXmodel.csv'
df.read_tbl(fname)
aov=df.anova('ERROR',wfactors=['TIMEOFDAY','COURSE','MODEL'])
## print(aov)
self.assertEqual(str(aov),R)
def test3(self):
## Mixed Between/Within test
R = """\
SUPPRESSION ~ CYCLE * PHASE * GROUP
TESTS OF BETWEEN-SUBJECTS EFFECTS
Measure: SUPPRESSION
Source Type III df MS F Sig. et2_G Obs. SE 95% CI lambda Obs.
SS Power
=========================================================================================================
Between Subjects 2.034 23
GROUP 0.462 2 0.231 3.083 0.067 0.150 8 0.103 0.201 2.349 0.229
=========================================================================================================
Error 1.572 21 0.075
TESTS OF WITHIN SUBJECTS EFFECTS
Measure: SUPPRESSION
Source Type III eps df MS F Sig. et2_G Obs. SE 95% CI lambda Obs.
SS Power
======================================================================================================================================================
CYCLE Sphericity Assumed 0.273 - 3 0.091 12.027 2.535e-06 0.094 48 0.013 0.025 27.491 0.995
Greenhouse-Geisser 0.273 0.668 2.005 0.136 12.027 7.281e-05 0.094 48 0.013 0.025 27.491 0.967
Huynh-Feldt 0.273 0.668 2.005 0.136 12.027 7.281e-05 0.094 48 0.013 0.025 27.491 0.967
Box 0.273 0.333 1 0.273 12.027 0.002 0.094 48 0.013 0.025 27.491 0.823
------------------------------------------------------------------------------------------------------------------------------------------------------
CYCLE * GROUP Sphericity Assumed 0.105 - 6 0.017 2.309 0.044 0.038 16 0.022 0.043 3.519 0.218
Greenhouse-Geisser 0.105 0.572 3.434 0.030 2.309 0.085 0.038 16 0.022 0.043 3.519 0.167
Huynh-Feldt 0.105 0.572 3.434 0.030 2.309 0.085 0.038 16 0.022 0.043 3.519 0.167
Box 0.105 0.167 1 0.105 2.309 0.158 0.038 16 0.022 0.043 3.519 0.107
------------------------------------------------------------------------------------------------------------------------------------------------------
Error(CYCLE) Sphericity Assumed 0.476 - 63 0.008
Greenhouse-Geisser 0.476 0.572 36.061 0.013
Huynh-Feldt 0.476 0.572 36.061 0.013
Box 0.476 0.167 10.500 0.045
------------------------------------------------------------------------------------------------------------------------------------------------------
PHASE Sphericity Assumed 1.170 - 1 1.170 129.855 1.878e-10 0.308 96 0.010 0.020 593.625 1
Greenhouse-Geisser 1.170 1 1 1.170 129.855 1.878e-10 0.308 96 0.010 0.020 593.625 1
Huynh-Feldt 1.170 1 1 1.170 129.855 1.878e-10 0.308 96 0.010 0.020 593.625 1
Box 1.170 1 1 1.170 129.855 1.878e-10 0.308 96 0.010 0.020 593.625 1
------------------------------------------------------------------------------------------------------------------------------------------------------
PHASE * GROUP Sphericity Assumed 0.405 - 2 0.203 22.493 6.012e-06 0.134 32 0.018 0.035 68.551 1.000
Greenhouse-Geisser 0.405 0.820 1.641 0.247 22.493 3.263e-05 0.134 32 0.018 0.035 68.551 1.000
Huynh-Feldt 0.405 0.820 1.641 0.247 22.493 3.263e-05 0.134 32 0.018 0.035 68.551 1.000
Box 0.405 0.500 1 0.405 22.493 6.896e-04 0.134 32 0.018 0.035 68.551 1.000
------------------------------------------------------------------------------------------------------------------------------------------------------
Error(PHASE) Sphericity Assumed 0.189 - 21 0.009
Greenhouse-Geisser 0.189 0.820 17.228 0.011
Huynh-Feldt 0.189 0.820 17.228 0.011
Box 0.189 0.500 10.500 0.018
------------------------------------------------------------------------------------------------------------------------------------------------------
CYCLE * Sphericity Assumed 0.074 - 3 0.025 4.035 0.011 0.027 24 0.016 0.032 4.612 0.386
PHASE Greenhouse-Geisser 0.074 0.686 2.057 0.036 4.035 0.024 0.027 24 0.016 0.032 4.612 0.313
Huynh-Feldt 0.074 0.686 2.057 0.036 4.035 0.024 0.027 24 0.016 0.032 4.612 0.313
Box 0.074 0.333 1 0.074 4.035 0.058 0.027 24 0.016 0.032 4.612 0.220
------------------------------------------------------------------------------------------------------------------------------------------------------
CYCLE * PHASE * GROUP Sphericity Assumed 0.127 - 6 0.021 3.466 0.005 0.046 8 0.028 0.055 2.641 0.169
Greenhouse-Geisser 0.127 0.605 3.629 0.035 3.466 0.019 0.046 8 0.028 0.055 2.641 0.137
Huynh-Feldt 0.127 1 6 0.021 3.466 0.005 0.046 8 0.028 0.055 2.641 0.169
Box 0.127 0.167 1 0.127 3.466 0.091 0.046 8 0.028 0.055 2.641 0.093
------------------------------------------------------------------------------------------------------------------------------------------------------
Error(CYCLE * Sphericity Assumed 0.386 - 63 0.006
PHASE) Greenhouse-Geisser 0.386 0.605 38.105 0.010
Huynh-Feldt 0.386 1 63 0.006
Box 0.386 0.167 10.500 0.037
TABLES OF ESTIMATED MARGINAL MEANS
Estimated Marginal Means for CYCLE
CYCLE Mean Std. Error 95% Lower Bound 95% Upper Bound
==============================================================
1 0.220 0.022 0.177 0.262
2 0.306 0.022 0.263 0.349
3 0.307 0.024 0.259 0.354
4 0.308 0.026 0.257 0.359
Estimated Marginal Means for PHASE
PHASE Mean Std. Error 95% Lower Bound 95% Upper Bound
==============================================================
I 0.207 0.014 0.180 0.234
II 0.363 0.016 0.332 0.394
Estimated Marginal Means for GROUP
GROUP Mean Std. Error 95% Lower Bound 95% Upper Bound
==============================================================
AA 0.222 0.018 0.187 0.256
AB 0.292 0.022 0.250 0.334
LAB 0.341 0.020 0.302 0.381
Estimated Marginal Means for CYCLE * PHASE
CYCLE PHASE Mean Std. Error 95% Lower Bound 95% Upper Bound
======================================================================
1 I 0.173 0.028 0.119 0.228
1 II 0.266 0.031 0.206 0.326
2 I 0.224 0.026 0.173 0.275
2 II 0.387 0.027 0.335 0.439
3 I 0.223 0.027 0.170 0.276
3 II 0.391 0.032 0.327 0.454
4 I 0.207 0.031 0.146 0.269
4 II 0.408 0.030 0.350 0.466
Estimated Marginal Means for CYCLE * GROUP
CYCLE GROUP Mean Std. Error 95% Lower Bound 95% Upper Bound
======================================================================
1 AA 0.193 0.046 0.104 0.282
1 AB 0.177 0.025 0.128 0.225
1 LAB 0.289 0.035 0.221 0.358
2 AA 0.253 0.033 0.187 0.318
2 AB 0.323 0.035 0.254 0.392
2 LAB 0.341 0.044 0.256 0.427
3 AA 0.219 0.036 0.149 0.289
3 AB 0.331 0.039 0.255 0.406
3 LAB 0.371 0.044 0.285 0.456
4 AA 0.223 0.026 0.172 0.273
4 AB 0.337 0.057 0.225 0.449
4 LAB 0.364 0.040 0.287 0.442
Estimated Marginal Means for PHASE * GROUP
PHASE GROUP Mean Std. Error 95% Lower Bound 95% Upper Bound
======================================================================
I AA 0.209 0.024 0.162 0.255
I AB 0.179 0.023 0.134 0.225
I LAB 0.233 0.026 0.183 0.283
II AA 0.235 0.026 0.184 0.286
II AB 0.404 0.023 0.359 0.450
II LAB 0.450 0.016 0.419 0.481
Estimated Marginal Means for CYCLE * PHASE * GROUP
CYCLE PHASE GROUP Mean Std. Error 95% Lower Bound 95% Upper Bound
==============================================================================
1 I AA 0.177 0.060 0.060 0.295
1 I AB 0.126 0.036 0.056 0.196
1 I LAB 0.216 0.047 0.124 0.309
1 II AA 0.209 0.072 0.067 0.351
1 II AB 0.228 0.025 0.179 0.276
1 II LAB 0.363 0.038 0.288 0.437
2 I AA 0.224 0.039 0.148 0.300
2 I AB 0.235 0.049 0.139 0.331
2 I LAB 0.214 0.053 0.110 0.317
2 II AA 0.281 0.055 0.173 0.389
2 II AB 0.411 0.026 0.360 0.463
2 II LAB 0.469 0.026 0.417 0.520
3 I AA 0.231 0.057 0.120 0.342
3 I AB 0.200 0.031 0.140 0.260
3 I LAB 0.238 0.054 0.133 0.342
3 II AA 0.208 0.047 0.115 0.300
3 II AB 0.461 0.024 0.414 0.508
3 II LAB 0.504 0.016 0.472 0.535
4 I AA 0.203 0.039 0.126 0.279
4 I AB 0.156 0.062 0.036 0.277
4 I LAB 0.264 0.058 0.149 0.378
4 II AA 0.242 0.034 0.176 0.309
4 II AB 0.517 0.031 0.457 0.578
4 II LAB 0.465 0.020 0.425 0.505
"""
df=DataFrame()
fname='data/suppression~subjectXgroupXcycleXphase.csv'
df.read_tbl(fname)
df['SUPPRESSION']=[.01*x for x in df['SUPPRESSION']]
aov=df.anova('SUPPRESSION',wfactors=['CYCLE','PHASE'],bfactors=['GROUP'])
## print(aov)
self.assertEqual(str(aov),R)
