power anaylysis для ановы со средними и сигмами |
Здравствуйте, гость ( Вход | Регистрация )
power anaylysis для ановы со средними и сигмами |
16.11.2017 - 16:49
Сообщение
#1
|
|
Группа: Пользователи Сообщений: 76 Регистрация: 27.04.2014 Пользователь №: 26375 |
Подскажите, пожалуйста, как мне рассчитать мощность для ановы с тремя группами, имея только средние и стандартные отклонения?
пример М=40, S=4 M2=35 S=10 M3=45 S=8 Как мне мощность рассчитать в R |
|
16.11.2017 - 16:58
Сообщение
#2
|
|
Группа: Пользователи Сообщений: 1091 Регистрация: 26.08.2010 Пользователь №: 22699 |
Подскажите, пожалуйста, как мне рассчитать мощность для ановы с тремя группами, имея только средние и стандартные отклонения? пример М=40, S=4 M2=35 S=10 M3=45 S=8 Как мне мощность рассчитать в R https://cran.r-project.org/web/packages/pwr2/pwr2.pdf |
|
16.11.2017 - 17:02
Сообщение
#3
|
|
Группа: Пользователи Сообщений: 76 Регистрация: 27.04.2014 Пользователь №: 26375 |
Код Description Calculate power for one-way ANOVA models. Usage pwr.1way(k=k, n=n, alpha=alpha, f=NULL, delta=delta, sigma=sigma) Arguments k Number of groups n Sample size per group f Effect size alpha Significant level (Type I error probability) delta The smallest difference among k groups sigma Standard deviation, i.e. square root of variance И где в этих аргументах вводить среднее и сигмы, у меня три группы и все имеют разные средние и сигмы)) а тут дана только одна сигма Сообщение отредактировал nastushka - 16.11.2017 - 17:04 |
|
16.11.2017 - 22:02
Сообщение
#4
|
|
Группа: Пользователи Сообщений: 1091 Регистрация: 26.08.2010 Пользователь №: 22699 |
Код Description Calculate power for one-way ANOVA models. Usage pwr.1way(k=k, n=n, alpha=alpha, f=NULL, delta=delta, sigma=sigma) Arguments k Number of groups n Sample size per group f Effect size alpha Significant level (Type I error probability) delta The smallest difference among k groups sigma Standard deviation, i.e. square root of variance И где в этих аргументах вводить среднее и сигмы, у меня три группы и все имеют разные средние и сигмы)) а тут дана только одна сигма Ну так посчитайте-оцените размер эффекта https://cran.r-project.org/web/packages/pwr...r-vignette.html |
|
17.11.2017 - 15:33
Сообщение
#5
|
|
Группа: Пользователи Сообщений: 76 Регистрация: 27.04.2014 Пользователь №: 26375 |
p2004r, там в любом случае надо вводить сигму, но только одну,а у меня их три)
Прям как в том анекдоте про блондинке в машине, ножек 2, а педалек 3)) |
|
17.11.2017 - 22:19
Сообщение
#6
|
|
Группа: Пользователи Сообщений: 1202 Регистрация: 13.01.2008 Из: Челябинск Пользователь №: 4704 |
Во-первых, классический дисперсионный анализ (как и классический t-критерий) требует равенства (однородности) дисперсий в группах. Это прописано во всех учебниках. Неправильно требовать от пакетов решений для неклассических задач. То, что на практике это требование может нарушаться и придуманы какие-то манёвры для его обхода - из другой области. Хотя можно попытаться использовать что-то из этих наработок для вашего случая. Например, модифицировать (уменьшить) степени свободы для F-критерия, т.е. попробовать посчитать мощность для Welch ANOVA. Но пока не очень себе представляю получится ли сделать это в рамках готовых пакетов типа рекомендованного pwr2.
Во-вторых, вы не указали количества наблюдений в группах. Без этого не понятно, отличаются ли дисперсии в группах (квадраты стандартных отклонений) статистически значимо или это может быть случайными выборочными вариациями одной и той же дисперсии. Если не отличаются - всё проще: нужно специальным образом усреднить ст. отклонения и считать мощность как для обычного случая. Всё-таки, приведите ещё и n - вдруг руки дойдут... Сообщение отредактировал nokh - 17.11.2017 - 22:30 |
|
18.11.2017 - 12:42
Сообщение
#7
|
|
Группа: Пользователи Сообщений: 76 Регистрация: 27.04.2014 Пользователь №: 26375 |
nokh, спасибо Вам, в каждой группе 25 чел. А как усреднить стандартные отклонения?
|
|
18.11.2017 - 19:27
Сообщение
#8
|
|
Группа: Пользователи Сообщений: 1202 Регистрация: 13.01.2008 Из: Челябинск Пользователь №: 4704 |
nokh, спасибо Вам, в каждой группе 25 чел. А как усреднить стандартные отклонения? Нужно гуглить, формулы видел. Также у Урбаха (Биометрические методы) есть на стр. 118 для дисперсий. Но теперь стало понятно, что вам это не подойдёт, т.к. ваши дисперисии сильно отличаются, например, для крайних вариантов 4 и 10, (т.е. дисперсий 16 и 100) F(24; 24)=100/16=6,25; P=0,000014. Т.о. усреднять ст. отклонения нельзя. Но столь сильные различия дисперсий наводят на мысль о ненормальном, возможно, сильно асимметричном распределении ошибки модели anova. Вы проверяли остатки на нормальность, использовали нормализующие преобразования в дисперсионном анализе? Сообщение отредактировал nokh - 19.11.2017 - 05:59 |
|
20.11.2017 - 15:18
Сообщение
#9
|
|
Группа: Пользователи Сообщений: 76 Регистрация: 27.04.2014 Пользователь №: 26375 |
nokh, все в порядке. Тут правда остатки распределены не нормально, а значит и не надо парится с анализом мощности, под эти данные, он не подходит.
|
|