power anaylysis для ановы со средними и сигмами Опции

[nastushka]

Подскажите, пожалуйста, как мне рассчитать мощность для ановы с тремя группами, имея только средние и стандартные отклонения?
пример
М=40, S=4
M2=35 S=10
M3=45 S=8
Как мне мощность рассчитать в R

[nokh]

Во-первых, классический дисперсионный анализ (как и классический t-критерий) требует равенства (однородности) дисперсий в группах. Это прописано во всех учебниках. Неправильно требовать от пакетов решений для неклассических задач. То, что на практике это требование может нарушаться и придуманы какие-то манёвры для его обхода - из другой области. Хотя можно попытаться использовать что-то из этих наработок для вашего случая. Например, модифицировать (уменьшить) степени свободы для F-критерия, т.е. попробовать посчитать мощность для Welch ANOVA. Но пока не очень себе представляю получится ли сделать это в рамках готовых пакетов типа рекомендованного pwr2.

Во-вторых, вы не указали количества наблюдений в группах. Без этого не понятно, отличаются ли дисперсии в группах (квадраты стандартных отклонений) статистически значимо или это может быть случайными выборочными вариациями одной и той же дисперсии. Если не отличаются - всё проще: нужно специальным образом усреднить ст. отклонения и считать мощность как для обычного случая. Всё-таки, приведите ещё и n - вдруг руки дойдут...

Сообщение отредактировал nokh - 17.11.2017 - 22:30