Masalah saya: uji coba acak kelompok paralel memiliki distribusi hasil primer yang sangat miring. Saya tidak ingin menganggap normalitas dan menggunakan 95% CI berbasis normal (yaitu menggunakan 1,96 X SE).
Saya merasa nyaman mengekspresikan ukuran kecenderungan sentral sebagai median, tetapi pertanyaan saya adalah bagaimana membangun 95% CI dari perbedaan median antara kedua kelompok.
Hal pertama yang terlintas dalam pikiran adalah bootstrap (resample dengan penggantian, tentukan median di masing-masing dua kelompok dan kurangi satu dari yang lain, ulangi 1000 kali, dan gunakan Bias-koreksi 95% CI). Apakah ini pendekatan yang benar? Ada saran lain?
wilcox.test()
(bawahDetails
), ini terkait erat dengan perbedaan median, tetapi tidak persis sama.Jawaban:
Prosedur bootstrap yang Anda uraikan harus valid. Namun, penting untuk diingat bahwa, seperti CI 95% berbasis normal, interval kepercayaan bootstrap hanya dijamin memiliki cakupan yang benar tanpa gejala. Satu hal yang menyenangkan tentang bekerja dengan median atau kuantil lainnya adalah bahwa Anda dapat membuat interval kepercayaan sampel hingga yang tepat di bawah asumsi yang sangat lemah. Ide dasarnya adalah bahwa di bawah nol bahwa median adalah , indikator untuk adalah variabel acak Bernoulli 0,5. Anda dapat menggunakan pengamatan ini untuk membuat statistik uji dengan distribusi sampel terbatas yang diketahui. Lihat Chernozhukov, Hansen, Jansson (2009) untuk lebih jelasnya.y m y< m
sumber
Anda juga dapat mencoba metode yang disarankan di http://www.ncbi.nlm.nih.gov/pubmed/12243307 (Bonett, Price; 2002) sebagai alternatif yang lebih sederhana (setidaknya secara komputasi, menurut saya). Pertanyaan yang bagus.
sumber