Как выполнить тест Tukey HSD с помощью команды Anova (автомобильный пакет)

Я имею дело с несбалансированным планом / образцом и изначально узнал aov () . Теперь я знаю, что для моих тестов ANOVA мне нужно использовать сумму квадратов типа III, которая включает использование подгонки с использованием lm () вместо использования aov () .

Проблема заключается в получении апостериорных тестов (в частности, HSD Тьюки) с использованием lm () . Все проведенные мной исследования показали, что использование simint в пакете multcomp будет работать, но теперь, когда он обновлен, эта команда кажется недоступной. Также кажется, что для вычисления требуется пройти через aov () .

По сути, все тесты Tukey HSD, которые я нашел для R, предполагают, что вы используете aov () для сравнения, а не lm () . Чтобы получить сумму квадратов типа III, необходимую для несбалансированного дизайна, я должен использовать:

mod<-lm(Snavg~StudentEthnicity*StudentGender)

Anova(mod, type="III")

Как мне использовать тест Tukey HSD с моим модом, используя lm () ? Или, наоборот, рассчитать мой дисперсионный анализ с использованием типа III и все еще иметь возможность запустить тест Tukey HSD?

Спасибо!