Distribusi polinomial derajat kedua dari variabel acak Gaussian

8

Saya ingin menghitung

P(Y=aX2+bX+c<0)

di mana . Saya bisa melakukannya dengan mudah menggunakan Monte Carlo. Namun, saya telah diminta untuk menemukan analitik dari dan kemudian menghitungXN(0,σ)fY(y)Y

I=0fY(y)dy

Saya kira akan sedemikian rupa sehingga hanya dapat dihitung secara numerik. Namun, karena ini merupakan integral univariat, metode numerik tersedia untuk menghitungnya dengan akurasi yang sangat tinggi. Apakah ada ekspresi (relatif sederhana) untuk , sehingga saya dapat melakukan integrasi numerik? Atau apakah ada kemungkinan lain untuk menghitung , selain dari Monte Carlo (yang menurut saya pendekatan yang paling masuk akal)? fY(y)IfY(y)I

DeltaIV
sumber
1
Apakah Anda harus menemukan pdf terlebih dahulu dan kemudian mengintegrasikannya ke garis real negatif, atau dapatkah Anda menggunakan metode yang ditunjukkan oleh mpiktas yang menghindari menemukan pdf ? YY
Dilip Sarwate
@DilipSarwate, terima kasih atas pertanyaannya. Saya secara khusus diminta untuk 1. menemukan dan 2. mengintegrasikan lebih dari . Jadi, jawaban yang benar-benar bagus. Di sisi lain, saya dapat menunjukkan bahwa permintaan tersebut tidak masuk akal dan saya sudah memiliki dua metode yang sangat bagus (MC dan @mpiktas) yang berfungsi dengan baik. Dengan demikian, jawaban atas pertanyaan Anda adalah: Saya tidak harus (saya tidak dipecat jika tidak melakukannya), tetapi saya sangat menghargai dapat melakukannya (sehingga menghindari diskusi lain dengan pemohon) . fY(y)[,0]
DeltaIV
3
Baiklah begini. Perhatikan bahwa yang dapat diekspresikan dalam istilah standar Gaussian CDF menggunakan metode yang dijelaskan dalam jawaban @ mpkitas. Mengambil wrt derivatif, maka akan memberikan pdf . Juga, beri tahu pemohon Anda bahwa Anda sebenarnya tidak perlu untuk secara eksplisit mengintegrasikan pdf untuk menemukan sejak yang nilainya telah Anda tentukan. FY(y)=P{Yy}=P{aX2+bX+cy0}Φ()yfY(y)sayasaya=FY(0)
Dilip Sarwate
@DilipSarwate fantastis! Dengan kata lain, dan dalam jawaban mpkitas menjadi fungsi , dan kemudian saya hanya menerapkan aturan rantai untuk derivasi. Terima kasih banyak! x1x2y
DeltaIV

Jawaban:

16

Perhatikan bahwa , di mana dan adalah akar dari polinomial . Kita harus mengasumsikan bahwa dan adalah nyata dan tidak sama, jika tidak, probabilitas yang dimaksud adalah nol atau satu.Sebuahx2+bx+c=Sebuah(x-x1)(x-x2)x1x2Sebuahx2+bx+cx1x2

Kami punya dua kasus.

  1. Sebuah>0 , lalu .P(SebuahX2+bX+c<0)=P(x1<X<x2)

  2. Sebuah<0 , laluP(SebuahX2+bX+c<0)=P(X<x1X>x2)=1-P(x1<X<x2).

Karena adalah normal, probabilitas dapat dihitung menggunakan fungsi distribusi kumulatif dari variabel normal.X

mpiktas
sumber
luar biasa! Dan tentu saja kami memiliki ekspresi analitik untuk dan (akar persamaan kuadrat), jadi semuanya sangat sederhana. Terima kasih! PS tentu saja dan selalu nyata dan berbeda dalam kasus saya. Saya lupa menentukannya. x1x2x1x2
DeltaIV
Solusi yang sangat pintar. Saya akan menyarankan untuk menurunkan distribusi tetapi sama sekali tidak perlu. Sudah selesai dilakukan dengan baik! Y
ramhiser
@ JohnA.Ramey Lihat komentar di pertanyaan utama.
Dilip Sarwate