bug in MultiComparison.tukeyhsd() where the sum of group nobs was being passed into stats.tukeyhsd() instead of the list of independent group nobs. This meant that df_total and thus the crit_int interval widths were being calculated incorrectly. Matters for pull request #809 so that now the interval widths of plot_intervals will match the results_table.