Saya bertanya-tanya apakah ada metode untuk menghitung ukuran sampel dalam model campuran? Saya menggunakan lmer
di R agar sesuai dengan model (saya punya lereng acak dan penyadapan).
r
mixed-model
lme4-nlme
power-analysis
Nikita Kuznetsov
sumber
sumber
Jawaban:
The
longpower
paket alat perhitungan ukuran sampel di Liu dan Liang (1997) dan Diggle et al (2002). Dokumentasi memiliki kode contoh. Ini salah satunya, menggunakanlmmpower()
fungsi:Juga periksa
liu.liang.linear.power()
yang " melakukan perhitungan ukuran sampel untuk model campuran linier"Liu, G., & Liang, KY (1997). Perhitungan ukuran sampel untuk studi dengan observasi berkorelasi. Biometrik, 53 (3), 937-47.
Diggle PJ, Heagerty PJ, Liang K, Zeger SL. Analisis data longitudinal. Edisi kedua. Oxford. Ilmu Statistik Melayani. 2002
Sunting: Cara lain adalah dengan "memperbaiki" efek pengelompokan. Dalam model linier biasa, setiap pengamatan adalah independen, tetapi di hadapan pengamatan clustering tidak independen yang dapat dianggap memiliki lebih sedikit pengamatan independen - ukuran sampel efektif lebih kecil. Hilangnya efektivitas ini dikenal sebagai efek desain :
sumber
Untuk apa pun di luar tes 2 sampel sederhana, saya lebih suka menggunakan simulasi untuk ukuran sampel atau studi kekuatan. Dengan rutinitas yang telah dikemas sebelumnya, Anda terkadang dapat melihat perbedaan besar antara hasil dari program berdasarkan pada asumsi yang mereka buat (dan Anda mungkin tidak dapat mengetahui apa asumsi itu, apalagi jika itu masuk akal untuk studi Anda). Dengan simulasi Anda mengontrol semua asumsi.
Berikut ini tautan ke sebuah contoh:
https://stat.ethz.ch/pipermail/r-sig-mixed-models/2009q1/001790.html
sumber