Saya harus menemukan 95% CI pada median dan persentil lainnya. Saya tidak tahu bagaimana mendekati ini. Saya terutama menggunakan R sebagai alat pemrograman.
r
confidence-interval
median
Dominic Comtois
sumber
sumber
library(boot)
muncul untuk mengonfirmasi ini:> boot.ci (boot (x, function (x, i) median (x [i]), R = 1000)) Interval: Level Normal Dasar 95% (74.42, 78.22) (75.00 , 78.49) Level Persentil BCa 95% (73.51, 77.00) (73.00, 77.00)Pendekatan lain didasarkan pada kuantil dari distribusi binomial.
misalnya:
sumber
Lihatlah bootstrap resampling. Cari bantuan R untuk fungsi boot. Bergantung pada data Anda dengan resampling, Anda dapat memperkirakan interval kepercayaan untuk apa saja.
sumber
wilcox.test(..., conf.int=TRUE)
fungsi R.Dan ada pendekatan lain: Salah satunya didasarkan pada uji Jumlah Wilcoxon Rank diterapkan untuk satu sampel dengan koreksi kontinuitas. Dalam R ini dapat diberikan sebagai:
Dan ada CI David Olive untuk median dibahas di sini:
CI untuk Median
sumber
Hasil berdasarkan pendekatan qbinom tidak benar untuk sampel kecil. Misalkan x memiliki 10 komponen. Kemudian qbinom (c (.025, .975), 10, .5) memberikan 2 dan 8. Interval yang dihasilkan tidak memperlakukan statistik pesanan di ekor bawah secara simetris dengan yang dari ekor atas; Anda harus mendapatkan 2 dan 9, atau 3 dan 8. Jawaban yang tepat adalah 2 dan 9. Anda dapat memeriksa terhadap proc univariate di SAS. Tangkap di sini adalah Anda tidak perlu lebih dari 0,025 probabilitas di bawah dan di atas; kuantil yang lebih rendah tidak melakukan ini, karena memberikan setidaknya 0,025 pada atau di bawah. Anda bisa diselamatkan di bagian bawah karena hitungan yang seharusnya 1 harus dipetakan ke statistik urutan kedua, menghitung 0, dan begitu "off per satu" dibatalkan. Pembatalan yang tidak disengaja ini tidak terjadi di atas, sehingga Anda mendapatkan jawaban yang salah di sini. Pengurutan kode (x) [qbinom (c (.025, .975), panjang (x) ,. 5) + c (0,1)] hampir berfungsi, dan 0,5 dapat digantikan oleh nilai-nilai kuantil lain untuk mendapatkan interval kepercayaan untuk kuantil lain, tetapi itu tidak akan benar ketika ada sedemikian sehingga P [X <= a ] =. 025. Lihat, misalnya, Higgins, Statisitcs Nonparametric.
sumber