Bagaimana , koordinat kutub, didistribusikan ketika dan kapan ?

19

Biarkan koordinat Cartesian dari titik acak dipilih st .x,y(x,y)U(10,10)×U(10,10)

Dengan demikian, jari-jari, , tidak merata seperti yang tersirat oleh 's pdf .ρ=x2+y2ρ

Meskipun demikian saya berharap hampir seragam, tidak termasuk artefak karena sisa 4 di tepi:θ=arctanyx

masukkan deskripsi gambar di sini

Berikut ini adalah fungsi kepadatan probabilitas yang dihitung secara grafis dari dan : θρmasukkan deskripsi gambar di sini

Sekarang jika saya membiarkan didistribusikan st maka tampaknya terdistribusi secara merata:x , y N ( 0 , 20 2 ) × N ( 0 , 20 2 ) θx,yx,yN(0,202)×N(0,202)θ

masukkan deskripsi gambar di sini

Mengapa tidak seragam ketika dan seragam ketika ?( x , y ) U ( - 10 , 10 ) × U ( - 10 , 10 ) x , y N ( 0 , 20 2 ) × N ( 0 , 20 2 )θ(x,y)U(10,10)×U(10,10)x,yN(0,202)×N(0,202)

Kode Matlab yang saya gunakan:

number_of_points = 100000;
rng('shuffle')

a = -10;
b = 10;
r = (b-a).*randn(2,number_of_points);
r = reshape(r, [2,number_of_points]);
I = eye(2);
e1 = I(:,1); e2 = I(:,2);
theta = inf*ones(1,number_of_points);
rho = inf*ones(1,number_of_points);

for i=1:length(r(1,:))
    x = r(:,i);
    [theta(i),rho(i)] = cart2pol(x(1),x(2));        
end

figure
M=3;N=1; bins = 360;
subplot(M,N,1); 
histogram(rad2deg(theta), bins)
title('Polar angle coordinate p.d.f');

subplot(M,N,2); 
histogram(rho, bins);
title('Polar radius coordinate p.d.f');

subplot(M,N,3); 
histogram(r(:));
title('The x-y cooridnates distrbution (p.d.f)');

Mengganti baris ke-3: r = (b-a).*randn(2,number_of_points);dengan r = (b-a).*randn(2,number_of_points) +a ;akan mengubah distribusi dari normal ke seragam.(x,y)

0x90
sumber
5
Pertanyaannya terlihat lebih cantik dan lebih cantik dengan setiap pengeditan dan judul pertanyaannya lebih jelas dan lebih ringkas. Bagus @ 0x90.
Michael R. Chernick
3
+1. Sangat menarik bahwa distribusi normal adalah satu-satunya yang mengarah ke sudut terdistribusi secara seragam (yaitu distribusi 2D simetris rotasi), lihat stats.stackexchange.com/a/255417/28666 .
Amoeba berkata Reinstate Monica

Jawaban:

13

Anda mengacu pada transformasi dari sepasang varian independen ke representasi polar ( R , θ ) (jari-jari dan sudut), dan kemudian melihat distribusi marginal θ .(X,Y)(R,θ)θ

Saya akan menawarkan penjelasan yang agak intuitif (meskipun derivasi matematis dari kepadatan pada dasarnya melakukan apa yang saya uraikan secara informal).

Perhatikan bahwa jika Anda skala dua variabel, X dan Y oleh beberapa skala umum (misalnya pergi dari U (-1,1) ke U (-10,10) atau dari N (0,1) ke N (0,20) pada kedua variabel secara bersamaan) yang tidak membuat perbedaan pada distribusi sudut (hanya mempengaruhi skala distribusi jari-jari). Jadi mari kita pertimbangkan kasus unit saja.

Pertama-tama pertimbangkan apa yang terjadi dengan kasing seragam. Perhatikan bahwa distribusinya seragam di atas satuan kuadrat, sehingga kepadatan probabilitas di wilayah yang terkandung dalam sebanding dengan luas wilayah. Secara khusus, lihat kepadatan yang terkait dengan elemen sudut, d θ dekat horizontal (dekat sudut ) dan pada diagonal (dekat sudut ):[1,1]2dθθ = π / 4θ=0θ=π/4

masukkan deskripsi gambar di sini

Jelas elemen probabilitas (yaitu area) yang sesuai dengan elemen sudut ( ) lebih besar ketika sudutnya dekat dengan salah satu diagonal. Memang pertimbangkan untuk menuliskan lingkaran di dalam kotak; area yang direntang oleh sudut kecil yang diberikan di dalam lingkaran adalah konstan, dan kemudian bagian di luar lingkaran tumbuh ketika kita mendekati diagonal, di mana ia pada titik maksimum. d θdfθdθ

Ini sepenuhnya menjelaskan pola yang Anda lihat dalam simulasi.

Memang, kita dapat melihat bahwa kerapatan harus proporsional dengan panjang segmen dari pusat persegi ke ujungnya; trigonometri sederhana sudah cukup untuk memperoleh kerapatan dari sana dan kemudian mudah untuk menemukan konstanta yang diperlukan untuk membuat kerapatan berintegrasi ke 1.

[Sunting: menambahkan bit ini selanjutnya untuk membahas radius, karena pertanyaan telah berubah sejak jawaban awal saya.]

Perhatikan bahwa jika kita memiliki distribusi yang seragam di atas lingkaran satuan (yaitu yang kita tulis di alun-alun sebelumnya) maka kepadatan jari-jari untuk itu akan sebanding dengan jari-jari (pertimbangkan luas elemen annular kecil lebar di radius - yaitu antara dan - memiliki luas proporsional dengan ). Kemudian ketika kita melewati di luar lingkaran, daerah annular baru dengan jari-jari lebih besar hanya mendapatkan kontribusi kepadatan dari bagian di alun-alun, sehingga kepadatan berkurang (awalnya cukup cepat, kemudian lebih lambat) antara dan . (Sekali lagi, pengertian geometris yang cukup sederhana cukup untuk mendapatkan bentuk fungsional kepadatan jika diperlukan.)drrrr+drr12


Sebaliknya, jika distribusi sambungan simetris secara rotasi tentang asal maka elemen probabilitas pada beberapa sudut tidak tergantung pada sudut (ini pada dasarnya adalah sebuah tautologi!). Distribusi bivariat dari dua standar Gaussians secara simetris berputar tentang asal:

masukkan deskripsi gambar di sini

(kode untuk gambar ini berdasarkan kode Elan Cohen di sini tetapi ada alternatif yang bagus di sini , dan sesuatu di antara keduanya di sini )

Akibatnya volume yang terkandung dalam beberapa sudut adalah sama untuk setiap , sehingga kepadatan yang terkait dengan sudut adalah seragam pada .dθθ[0,2π)

[Trik polar biasanya digunakan untuk mengintegrasikan kerapatan normal di atas garis nyata dapat digunakan untuk mengetahui bahwa kerapatan jari-jari kuadrat adalah eksponensial negatif, dan dari sana kerapatan jari-jari mudah diidentifikasi dengan argumen transformasi sederhana dari fungsi distribusi]

Glen_b -Reinstate Monica
sumber
4
Keempat paku dalam distribusi memang karena empat sudut alun-alun . Perhatikan bahwa setiap distribusi simetris bola akan mengarah ke distribusi Seragam pada , dimulai dengan Seragam pada bola dan lingkaran yang berpusat pada . θ(10,10)2θ(0,0)
Xi'an
2
+1. Sangat menarik bahwa distribusi normal adalah satu-satunya yang mengarah ke distribusi 2D simetris rotasi, lihat stats.stackexchange.com/a/255417/28666 . Ini mengejutkan saya.
Amoeba berkata Reinstate Monica
3
@amoeba Ya, itu satu-satunya distribusi simetris bundar yang merupakan produk dari margin independen.
Glen_b -Reinstate Monica
2
Saya pikir itu luar biasa. Pertimbangkan menyebutkannya dalam jawaban Anda!
Amuba mengatakan Reinstate Monica
6

Saya akan menjawab pertanyaan tentang kasus normal yang mengarah ke distribusi seragam. Diketahui bahwa jika dan Y adalah independen dan terdistribusi normal, kontur kerapatan probabilitas konstan adalah sebuah lingkaran pada bidang x - y . Jari-jari R = XYxy memilikidistribusi Rayleigh. Untuk diskusi yang baik tentang ini, artikel wikipedia berjudul Distribusi Rayleigh.R=X2+Y2

Sekarang mari kita lihat variabel acak dan Y menggunakan koordinat kutub.XY

X=rcos(θ)Y=rsin(θ)X2+Y2=r2θ(0,2π)rXY0

XN(0,1)YN(0,1)

f(x,y)=(1/2π)exp[([x2+y2])/2]g(r,θ)x=rsin(θ)y=rcos(θ)r=x2+y2θ=arctan(x/y)f(x,y)g(r,θ)rexp[(r2)/(2π)]r00θ2πrr1/(2π)

Michael R. Chernick
sumber
Apa artinya adalah bahwa jika Anda melihat ketinggian kepadatan bivariat pada jarak radial tetap dari pusat (dalam hal ini asalnya) itu akan menjadi nilai yang sama pada semua titik pada lingkaran itu.
Michael R. Chernick
@ 0x90 Ya, tautan Anda menunjukkan salah satu cara untuk melihatnya adalah dengan melihat bentuk kuadratik dalam eksponen kerapatan. Jadi secara umum untuk pengaturan normal bivariat yang eksponen ke konstanta menentukan kontur kerapatan konstan dan persamaan itu adalah salah satu elips. dalam kasus khusus ketika matriks kovarians adalah matirix identitas skala, elips menyederhanakan menjadi lingkaran.
Michael R. Chernick
2
X,Y0Cauchy(0,1)arctan arctan(X/Y)
1
@ Francis Sebagian besar saya menghargai editing menyeluruh Anda dari semua persamaan saya. Saya juga ingin mengatakan bahwa komentar Anda di atas jelas menunjukkan pendekatan imajinatif untuk menyelesaikan masalah keseragaman dengan theta. Saya yakin beberapa orang akan setuju bahwa ini lebih mudah.
Michael R. Chernick
6

θ(X,Y)[1,1]×[1,1]14014

Wilayah yang menarik untuk pertanyaan kami adalah sektor merah pada gambar ini: persegi dengan sektor yang diarsir

θθ+dθθθ+dθθ

θ[π4,π4]π2

1cosθ

1cos(θ+dθ)=1cosθ+sinθcos2θdθ.

abα12absinα

12(1cosθ)(1cosθ+sinθcos2θdθ)sindθ=dθ2cos2θ
dθsindθ=dθ

θ

18cos2θ
θ[π4,π4]π2

Verifikasi:

x <- runif(1e6, -1, 1)
y <- runif(1e6, -1, 1)
hist( atan2(y,x), freq=FALSE, breaks=100)
theta <- seq(-pi, pi, length=500)
lines(theta, 0.125/cos((theta + pi/4)%%(pi/2) - pi/4)**2, col="red" )

histogram + kepadatan

Elvis
sumber