Форум врачей-аспирантов

Здравствуйте, гость ( Вход | Регистрация )

 
Добавить ответ в эту темуОткрыть тему
> power anaylysis для ановы со средними и сигмами
nastushka
сообщение 16.11.2017 - 16:49
Сообщение #1


Дух форума
*

Группа: Пользователи
Сообщений: 68
Регистрация: 27.04.2014
Пользователь №: 26375



Подскажите, пожалуйста, как мне рассчитать мощность для ановы с тремя группами, имея только средние и стандартные отклонения?
пример
М=40, S=4
M2=35 S=10
M3=45 S=8
Как мне мощность рассчитать в R
Вернуться в начало страницы
 
+Ответить с цитированием данного сообщения
 
p2004r
сообщение 16.11.2017 - 16:58
Сообщение #2


Дух форума
*

Группа: Пользователи
Сообщений: 1015
Регистрация: 26.08.2010
Пользователь №: 22699



Цитата(nastushka @ 16.11.2017 - 16:49) *
Подскажите, пожалуйста, как мне рассчитать мощность для ановы с тремя группами, имея только средние и стандартные отклонения?
пример
М=40, S=4
M2=35 S=10
M3=45 S=8
Как мне мощность рассчитать в R


https://cran.r-project.org/web/packages/pwr2/pwr2.pdf


Signature
Вернуться в начало страницы
 
+Ответить с цитированием данного сообщения
 
nastushka
сообщение 16.11.2017 - 17:02
Сообщение #3


Дух форума
*

Группа: Пользователи
Сообщений: 68
Регистрация: 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
Вернуться в начало страницы
 
+Ответить с цитированием данного сообщения
 
p2004r
сообщение 16.11.2017 - 22:02
Сообщение #4


Дух форума
*

Группа: Пользователи
Сообщений: 1015
Регистрация: 26.08.2010
Пользователь №: 22699



Цитата(nastushka @ 16.11.2017 - 17:02) *
Код
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


Signature
Вернуться в начало страницы
 
+Ответить с цитированием данного сообщения
 
nastushka
сообщение 17.11.2017 - 15:33
Сообщение #5


Дух форума
*

Группа: Пользователи
Сообщений: 68
Регистрация: 27.04.2014
Пользователь №: 26375



p2004r, там в любом случае надо вводить сигму, но только одну,а у меня их три)
Прям как в том анекдоте про блондинке в машине, ножек 2, а педалек 3))
Вернуться в начало страницы
 
+Ответить с цитированием данного сообщения
 
nokh
сообщение 17.11.2017 - 22:19
Сообщение #6


Дух форума
*

Группа: Пользователи
Сообщений: 1021
Регистрация: 13.01.2008
Из: Челябинск
Пользователь №: 4704



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

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

Сообщение отредактировал nokh - 17.11.2017 - 22:30
Вернуться в начало страницы
 
+Ответить с цитированием данного сообщения
 
nastushka
сообщение 18.11.2017 - 12:42
Сообщение #7


Дух форума
*

Группа: Пользователи
Сообщений: 68
Регистрация: 27.04.2014
Пользователь №: 26375



nokh, спасибо Вам, в каждой группе 25 чел. А как усреднить стандартные отклонения?
Вернуться в начало страницы
 
+Ответить с цитированием данного сообщения
 
nokh
сообщение 18.11.2017 - 19:27
Сообщение #8


Дух форума
*

Группа: Пользователи
Сообщений: 1021
Регистрация: 13.01.2008
Из: Челябинск
Пользователь №: 4704



Цитата(nastushka @ 18.11.2017 - 14:42) *
nokh, спасибо Вам, в каждой группе 25 чел. А как усреднить стандартные отклонения?

Нужно гуглить, формулы видел. Также у Урбаха (Биометрические методы) есть на стр. 118 для дисперсий. Но теперь стало понятно, что вам это не подойдёт, т.к. ваши дисперисии сильно отличаются, например, для крайних вариантов 4 и 10, (т.е. дисперсий 16 и 100) F(24; 24)=100/16=6,25; P=0,000014. Т.о. усреднять ст. отклонения нельзя. Но столь сильные различия дисперсий наводят на мысль о ненормальном, возможно, сильно асимметричном распределении ошибки модели anova. Вы проверяли остатки на нормальность, использовали нормализующие преобразования в дисперсионном анализе?

Сообщение отредактировал nokh - 19.11.2017 - 05:59
Вернуться в начало страницы
 
+Ответить с цитированием данного сообщения
 
nastushka
сообщение 20.11.2017 - 15:18
Сообщение #9


Дух форума
*

Группа: Пользователи
Сообщений: 68
Регистрация: 27.04.2014
Пользователь №: 26375



nokh, все в порядке. Тут правда остатки распределены не нормально, а значит и не надо парится с анализом мощности, под эти данные, он не подходит.
Вернуться в начало страницы
 
+Ответить с цитированием данного сообщения
 

Добавить ответ в эту темуОткрыть тему