thanks, jerry. there're typos and corrected for those who may also need
this. on the other hand, i don't know what "single degree of freedom to test
for non-additivity" is and i wonder if a ranking saying which pairs are
significantly different, thanks again!

Dim rAveD(), cAveD() ' arrays of deviations of row/column averages
from grand mean

Dim rAveD(), cAveD() ' arrays of deviations of row/column averages from
grand mean

SSr As Double, SSc As Double, SSTukey As Double, SStot As Double,

Tukey1df = WorksheetFunction.FDist(FTukey, 1, dfe) ' p-value for
non-additivity

Tukey1df = WorksheetFunction.FDist(FTukey, 1, dfe) ' p-value for
non-additivity