Power analysis of HLM LMM
1.0.0
نظرًا لقيمة استنساخ علم النفس بشكل متزايد ، فإن الطلب على حسابات الطاقة للتحليل الإحصائي يزداد أيضًا. يمكن إكمال حسابات الطاقة لأكثر الطرق الإحصائية الأساسية ( t test ، ANOVA ، إلخ) في بعض البرامج (مثل jamovi ). ومع ذلك ، كطريقة إحصائية متقدمة ، لا يمكن إجراء حساب الطاقة للنموذج المختلط الخطي (LMM) في برامج مثل jamovi مع واجهة مرئية ، ولا تزال حزم ووظائف البيانات ذات الصلة في R مطلوبة. تقدم هذه المقالة بشكل أساسي كيفية إجراء حساب الطاقة لـ LMM في R ، بما في ذلك: مقدمة للبيانات ووظائفها ، السيناريوهات المطبقة والمقيدة ، ومثلات مثال محددة.
حزمة بيانات simr هي حزمة أكثر ملاءمة لتحليل الطاقة LMM في R في التحليل ، يتم محاكاة بيانات النموذج عن طريق محاكاة Monte Carlo ويتم حساب النسبة الكبيرة في عدد المحاكاة المحددة. الوظيفة الرئيسية هي powerSim() .
احسب قوة التأثير الخلفي ، أي حساب الطاقة بناءً على التأثيرات في النموذج المبني على البيانات الموجودة ؛
احسب قوة التأثير المسبق ، أي حساب الطاقة بناءً على حجم التأثير في الأبحاث الحالية ؛
تقدير عدد الموضوعات/العناصر ، لأن حجم الطاقة يتأثر بعدد الموضوعات أو العناصر. النموذج المحدد هو أنه كلما زاد عدد العدد ، كلما زادت الطاقة ، والتي توفر لنا وظيفة تقدير عدد الموضوعات/العناصر المطلوبة للتجارب الرسمية بناءً على تأثيرات ما قبل التجريبية ؛
هنا نستخدم البيانات من تجربة علم النفس (عنوان تنزيل البيانات) لإظهار. جمعت هذه التجربة تفاعلات 10 مواضيع ( subj ) عند إجراء مهمة إدراكية ( DV ). كان للتجربة شرطين ( CondA و CondB ) ، والتي تم تصميمها لمواضيع 2*2 وداخل المشروع ( item ). وهذا هو ، يعتبر عاملان عشوائيان في نفس الوقت.
نظرًا لأننا نركز على التأثيرات الثابتة والتأثيرات غير العشوائية هنا ، فإننا ننظر فقط في اعتراضات عشوائية في جزء التأثيرات العشوائية عند النمذجة ، ويبحث جزء التأثيرات الثابتة في التأثيرين الرئيسيين وتفاعلاتها.
النمذجة كما يلي:
> Model = lmer(data = DemoData, DV~CondA*CondB+(1|subj)+(1|item))
جزء التأثير الثابت من النموذج هو:
> summary(Model)$coef
# Estimate Std. Error df t value Pr(>|t|)
# (Intercept) 279.43090 23.37537 11.71977 11.9540740 6.422210e-08
# CondAA2 24.29565 13.49694 498.43142 1.8000866 7.245150e-02
# CondBB2 12.18484 13.36445 509.31328 0.9117351 3.623395e-01
# CondAA2:CondBB2 -32.82881 19.29629 509.76684 -1.7013020 8.949610e-02
هنا نأخذ التأثير الثابت لفحص التفاعل بين الاثنين ( CondAA2:CondBB2 ) كمثال:
> PowerAB_ttest = simr::powerSim(fit = Model, # 要考察的模型
test = fixed('CondAA2:CondBB2', # 要考察的固定效应的名称
method = 't'), # 选取检验方法，因为固定效应为t检验，因此method设置为t
nsim=50) # 设置模拟次数，建议设置为500 (此时可以获取到较稳定的power)
> PowerAB_ttest
# Power for predictor 'CondAA2:CondBB2', (95% confidence interval):
# 46.00% (31.81, 60.68)
#
# Test: t-test with Satterthwaite degrees of freedom (package lmerTest)
# Effect size for CondAA2:CondBB2 is -33.
#
# Based on 50 simulations, (2 warnings, 0 errors)
# alpha = 0.05, nrow = 553
#
# Time elapsed: 0 h 0 m 8 s
#
# nb: result might be an observed power calculation
كانت قوة التأثير الثابت للتفاعل 0.46 (فاصل الثقة 0.318 – 0.607 ).
وبالمثل ، يمكن تنفيذ الطاقة التي تحسب التأثير الثابت CondA من خلال الكود التالي:
> PowerA_ttest = simr::powerSim(fit = Model, test = fixed('CondAA2', method = 't'), nsim=50)
جزء التأثير/التفاعل الرئيسي من النموذج هو:
> anova(Model)
# Type III Analysis of Variance Table with Satterthwaite's method
# Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
# CondA 8546 8546 1 489.18 0.6694 0.4137
# CondB 2453 2453 1 509.78 0.1921 0.6613
# CondA:CondB 36951 36951 1 509.77 2.8944 0.0895 .
# ---
# Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
هنا نأخذ فحص التأثير الرئيسي CondA كمثال:
> PowerA_Ftest = simr::powerSim(fit = Model, # 要考察的模型
test = fixed('CondA', # 要考察的主效应/交互作用的名称
method = 'f'), # 选取检验方法，因为主效应为f检验，因此method设置为f
nsim=50) # 设置模拟次数，建议设置为500 (此时可以获取到较稳定的power)
> PowerA_Ftest
# Power for predictor 'CondA', (95% confidence interval):
# 12.00% ( 4.53, 24.31)
#
# Test: Type-II F-test (package car)
#
# Based on 50 simulations, (50 warnings, 0 errors)
# alpha = 0.05, nrow = 553
#
# Time elapsed: 0 h 0 m 58 s
#
# nb: result might be an observed power calculation
قوة التأثير الرئيسي CondA هي 0.12 (فاصل الثقة هو 0.045 – 0.243 ).
إذا تم الحصول على حجم تأثير مستقر نسبيًا في الدراسات السابقة ، على سبيل المثال ، فقد وجد أن التأثير الثابت CondA حوالي 50 ، يمكن أيضًا حساب الطاقة بناءً على التأثيرات السابقة. الخطوات كما يلي:
> fixef(Model) # 查看模型固定效应
# (Intercept) CondAA2 CondBB2 CondAA2:CondBB2
# 279.43090 24.29565 12.18484 -32.82881
> fixef(Model)[2] = 50 # 修改CondA的效应
> fixef(Model) # 查看修改后的固定效应
# (Intercept) CondAA2 CondBB2 CondAA2:CondBB2
# 279.43090 50.00000 12.18484 -32.82881
> PowerA_ttest = simr::powerSim(fit = Model, test = fixed('CondAA2', method = 't'), nsim=50)
> PowerA_ttest
# Power for predictor 'CondAA2', (95% confidence interval):
# 88.00% (75.69, 95.47)
#
# Test: t-test with Satterthwaite degrees of freedom (package lmerTest)
# Effect size for CondAA2 is 50.
#
# Based on 50 simulations, (3 warnings, 0 errors)
# alpha = 0.05, nrow = 553
#
# Time elapsed: 0 h 0 m 9 s
في هذا الوقت ، فإن قوة التأثير الثابتة CondA هي 0.88 .
> PowerA_Ftest = simr::powerSim(fit = Model, test = fixed('CondA', method = 'f'), nsim=50)
> PowerA_Ftest
# Power for predictor 'CondA', (95% confidence interval):
# 94.00% (83.45, 98.75)
#
# Test: Type-II F-test (package car)
#
# Based on 50 simulations, (50 warnings, 0 errors)
# alpha = 0.05, nrow = 553
#
# Time elapsed: 0 h 1 m 4 s
قوة التأثير الرئيسية CondA هي 0.94 .
> Model = lmer(data = DemoData, DV~CondA*CondB+(1|subj)+(1|item)) # 这里仍以后验power为例
قوة التأثير الثابت للتفاعل هي 0.46 ، ونريد اختبار عدد الموضوعات ، يمكن أن تصل الطاقة إلى 0.8 . الخطوات كما يلي:
> Model2 = extend(object = Model, along = 'subj', n = 30)
> PowerA_ttest2 = powerSim(Model2, fixed('CondAA2', method = 't'), nsim=50)
> PowerA_ttest2
# Power for predictor 'CondAA2', (95% confidence interval):
# 84.00% (70.89, 92.83)
#
# Test: t-test with Satterthwaite degrees of freedom (package lmerTest)
# Effect size for CondAA2 is 24.
#
# Based on 50 simulations, (6 warnings, 0 errors)
# alpha = 0.05, nrow = 1659
#
# Time elapsed: 0 h 0 m 11 s
#
# nb: result might be an observed power calculation
في هذا الوقت ، ارتفعت الطاقة إلى 0.84 .
powerCurve() > Pcurve = powerCurve(fit = Model2, test = fixed('CondAA2', method = 't'), along = 'subj', nsim=50)
> Pcurve
# Power for predictor 'CondAA2', (95% confidence interval), by largest value of subj:
# 11: 22.00% (11.53, 35.96) - 180 rows
# 14: 36.00% (22.92, 50.81) - 325 rows
# 17: 50.00% (35.53, 64.47) - 507 rows
# 2: 58.00% (43.21, 71.81) - 667 rows
# 22: 64.00% (49.19, 77.08) - 837 rows
# 25: 58.00% (43.21, 71.81) - 989 rows
# 28: 74.00% (59.66, 85.37) - 1166 rows
# 30: 82.00% (68.56, 91.42) - 1318 rows
# 6: 86.00% (73.26, 94.18) - 1489 rows
# 9: 90.00% (78.19, 96.67) - 1659 rows
#
# Time elapsed: 0 h 1 m 39 s
لاحظ أن القيمة أمام كل قوة غير دقيقة في هذا الوقت ، والقيمة الحقيقية هي:
> Pcurve$nlevels
# [1] 3 6 9 12 15 18 21 24 27 30
يمكن ملاحظة أنه عندما يصل حوالي 24 موضوعًا إلى الطاقة ، يمكن أن تصل الطاقة إلى أكثر من 0.8 .
يمكن ضبط نطاق حجم العينة يدويًا:
> Pcurve = powerCurve(fit = Model2, test = fixed('CondAA2', method = 't'), along = 'subj', nsim=50,
breaks=c(21, 22, 23, 24, 25, 30)) # 考察样本量为21、22、23、24、25、30时的power
> Pcurve
# Power for predictor 'CondAA2', (95% confidence interval), by largest value of subj:
# 28: 72.00% (57.51, 83.77) - 1166 rows
# 29: 74.00% (59.66, 85.37) - 1220 rows
# 3: 70.00% (55.39, 82.14) - 1262 rows
# 30: 74.00% (59.66, 85.37) - 1318 rows
# 4: 76.00% (61.83, 86.94) - 1369 rows
# 9: 80.00% (66.28, 89.97) - 1659 rows
#
# Time elapsed: 0 h 1 m 12 s
> Pcurve$nlevels
# [1] 21 22 23 24 25 30
إن حساب الطاقة للتأثير في النموذج المعمم يشبه ما سبق ، والاختلافات الرئيسية هي:
method في fixed() على 'z' لأن التأثير الثابت في النموذج المعمم هو اختبار z ؛method في fixed() على 'chisq' ، لأن التأثير الرئيسي في النموذج المعمم هو اختبار تشي مربع
