Saya memiliki dua variabel acak, mana adalah distribusi 0-1 yang seragam.
Kemudian, ini menghasilkan suatu proses, katakan:
Sekarang, saya bertanya-tanya apakah ada ekspresi bentuk tertutup untuk teoritis 75 persen kuantil dari untuk diberikan --Saya kira saya bisa melakukannya dengan komputer dan banyak realisasi , tetapi saya lebih suka bentuk tertutup--.
distributions
probability
stochastic-processes
pengguna603
sumber
sumber
quant = function(n,p,x) return( quantile(runif(n)*sin(x)+runif(n)*cos(x),p) )
danquant(100000,0.75,1)
.Jawaban:
Masalah ini dapat dengan cepat direduksi menjadi salah satu dari menemukan kuantil dari distribusi trapesium .
Mari kita menulis ulang proses sebagai mana dan adalah variabel acak iid ; dan, secara simetri, ini memiliki distribusi marjinal yang sama dengan proses Dua istilah pertama menentukan kepadatan trapesium simetris
Kuantitas distribusi trapesium
Misalkan mana dan independen dan distribusi. Asumsikan tanpa kehilangan keumuman bahwa . Kemudian, kepadatan terbentuk dengan menggabungkan kepadatan dan . Ini mudah dilihat sebagai trapesium dengan simpul , , dan .Y=X1+X2 X1 X2 U(−a,a) U(−b,b) a≥b Y X1 X2 (−a−b,0) (−a+b,1/2a) (a−b,1/2a) (a+b,0)
Quantum dari distribusi , untuk setiap adalah, dengan demikian, Dengan simetri, untuk , kita memiliki .Y p<1/2
Kembali ke kasing yang ada
Di atas sudah cukup memberikan ekspresi bentuk tertutup. Yang kita butuhkan adalah membobol dua kasusdanuntuk menentukan yang memainkan peran dan yang memainkan peran atas. (Faktor 2 di sini hanya untuk mengkompensasi divisi dengan dua dalam definisi .)|sinx|≥|cosx| |sinx|<|cosx| 2a 2b P¯¯¯¯(x)
Untuk , di, kita atur dan dan mendapatkan dan padaperannya terbalik. Demikian pula, untukp<1/2 |sinx|≥|cosx| a=|sinx|/2 b=|cosx|/2
Kuantil
Di bawah ini adalah dua heatmap. Yang pertama menunjukkan kuantil dari distribusi untuk kisi berjalan dari hingga . The -coordinate memberikan probabilitas terkait dengan setiap kuantil. Warna-warna menunjukkan nilai kuantil dengan merah tua menunjukkan nilai sangat besar (positif) dan biru tua menunjukkan nilai negatif besar. Dengan demikian setiap strip vertikal adalah plot kuantil (marginal) yang terkait dengan .P(x) x 0 2π y p P(x)
Peta panas kedua di bawah ini menunjukkan kuantil itu sendiri, diwarnai oleh probabilitas yang sesuai. Misalnya, merah tua berhubungan dengan dan biru tua berhubungan dengan dan . Cyan kira-kira dan . Ini lebih jelas menunjukkan dukungan dari setiap distribusi dan bentuk.p=1/2 p=0 p=1 p=1/4 p=3/4
Beberapa
R
kode contohFungsi diP(x) x
qproc
bawah ini menghitung fungsi kuantil untuk diberikan . Ini menggunakan yang lebih umum untuk menghasilkan kuantil.qtrap
Di bawah ini adalah tes dengan output yang sesuai.
sumber