Misalkan saya punya variabel acak normal independen
dan . Bagaimana saya menandai kerapatan jika distribusi setiap masing -masing terpotong ke dalam ? Dengan kata lain, saya mengambil sampel dari distribusi normal independen, membuang sampel yang tidak berada dalam dari setiap rata-rata, dan menjumlahkannya.
Saat ini, saya melakukan ini dengan kode R di bawah ini:
x_mu <- c(12, 18, 7)
x_sd <- c(1.5, 2, 0.8)
a <- x_mu - 2 * x_sd
b <- x_mu + 2 * x_sd
samples <- sapply(1:3, function(i) {
return(rtruncnorm(100000, a[i], b[i], x_mu[i], x_sd[i]))
})
y <- rowSums(samples)
Apakah ada metode untuk menghasilkan kerapatan secara langsung?
Jawaban:
Anda bisa menggunakan pendekatan dengan metode saddlepoint, untuk jumlah normals terpotong. Saya tidak akan memberikan detailnya sekarang, Anda dapat melihat jawaban saya untuk jumlah umum distribusi Gamma untuk petunjuk. Yang kita butuhkan adalah menemukan fungsi penghasil momen untuk normal terpotong, yang mudah. Saya akan melakukannya di sini untuk standar normal terpotong± 2 , yang memiliki kepadatan
Fungsi pembangkit momen dapat dihitung sebagai
sumber
Saya ingin tahu mengapa, tapi ya, ada cara sederhana untuk menghasilkan pdf dari jumlah distribusi ini:
sumber