Bisakah standar deviasi dihitung untuk rata-rata harmonik? Saya mengerti bahwa standar deviasi dapat dihitung untuk rata-rata aritmatika, tetapi jika Anda memiliki rata-rata harmonis, bagaimana Anda menghitung standar deviasi atau CV?
sumber
Bisakah standar deviasi dihitung untuk rata-rata harmonik? Saya mengerti bahwa standar deviasi dapat dihitung untuk rata-rata aritmatika, tetapi jika Anda memiliki rata-rata harmonis, bagaimana Anda menghitung standar deviasi atau CV?
Rata-rata harmonik dari variabel acak didefinisikan sebagai
Mengambil momen pecahan adalah bisnis berantakan, jadi bukan saya lebih suka bekerja dengan . Sekarang
Dalam teorema limit pusat, kita segera mendapatkan itu
jika tentu saja dan adalah iid, karena kami bekerja sederhana dengan rata-rata aritmatika variabel .
Sekarang menggunakan metode delta untuk fungsi kita dapatkan itu
Hasil ini asimptotik, tetapi untuk aplikasi sederhana mungkin cukup.
Pembaruan Seperti yang ditunjukkan @whuber dengan tepat, aplikasi sederhana adalah nama yang salah. Teorema batas pusat hanya berlaku jika ada, yang merupakan asumsi yang sangat membatasi.
Perbarui 2 Jika Anda memiliki sampel, maka untuk menghitung standar deviasi, cukup masukkan momen sampel ke dalam rumus. Jadi untuk sampel , estimasi rata-rata harmonik adalah
contoh momen dan masing-masing adalah:
di sini berarti timbal balik.
Akhirnya rumus perkiraan untuk standar deviasi adalah
Saya menjalankan beberapa simulasi Monte-Carlo untuk variabel acak yang terdistribusi secara merata dalam interval . Ini kodenya:
hm <- function(x)1/mean(1/x)
sdhm <- function(x)sqrt((mean(1/x))^(-4)*var(1/x)/length(x))
n<-1000
nn <- c(10,30,50,100,500,1000,5000,10000)
N<-1000
mc<-foreach(n=nn,.combine=rbind) %do% {
rr <- matrix(runif(n*N,min=2,max=3),nrow=N)
c(n,mean(apply(rr,1,sdhm)),sd(apply(rr,1,sdhm)),sd(apply(rr,1,hm)))
}
colnames(mc) <- c("n","DeltaSD","sdDeltaSD","trueSD")
> mc
n DeltaSD sdDeltaSD trueSD
result.1 10 0.089879211 1.528423e-02 0.091677622
result.2 30 0.052870477 4.629262e-03 0.051738941
result.3 50 0.040915607 2.705137e-03 0.040257673
result.4 100 0.029017031 1.407511e-03 0.028284458
result.5 500 0.012959582 2.750145e-04 0.013200580
result.6 1000 0.009139193 1.357630e-04 0.009115592
result.7 5000 0.004094048 2.685633e-05 0.004070593
result.8 10000 0.002894254 1.339128e-05 0.002964259
Saya mensimulasikan N
sampel sampel n
berukuran. Untuk setiap n
sampel berukuran saya menghitung estimasi estimasi standar (fungsi sdhm
). Kemudian saya membandingkan rata-rata dan standar deviasi dari estimasi ini dengan standar deviasi sampel dari rata-rata harmonik yang diperkirakan untuk setiap sampel, yang seharusnya merupakan standar deviasi rata-rata rata-rata harmonik.
Seperti yang Anda lihat hasilnya cukup baik bahkan untuk ukuran sampel sedang. Tentu saja distribusi seragam adalah yang berperilaku sangat baik, sehingga tidak mengherankan bahwa hasilnya baik. Saya akan pergi ke orang lain untuk menyelidiki perilaku untuk distribusi lain, kode ini sangat mudah untuk beradaptasi.
Catatan: Di versi sebelumnya dari jawaban ini ada kesalahan dalam hasil metode delta, varians salah.
Kerugian utama adalah bahwa perhitungan tidak menghasilkan interval kepercayaan yang baik untuk distribusi dasar yang sangat miring. Itu kemungkinan akan menjadi masalah dengan metode tujuan umum: rata-rata harmonik peka terhadap keberadaan bahkan nilai kecil tunggal dalam dataset.
Berikut adalah
R
kode untuk simulasi dan gambar.sumber
Berikut adalah contoh untuk Exponential r.v's.
Varian (dan standar deviasi) rv ini sudah dikenal, lihat, misalnya di sini .
sumber
sumber
sumber