Operasi trigonometri pada standar deviasi

14

Penambahan, pengurangan, perkalian dan pembagian variabel acak normal didefinisikan dengan baik, tetapi bagaimana dengan operasi trigonometri?

Sebagai contoh, mari kita anggap bahwa saya sedang berusaha mencari sudut baji segitiga (dimodelkan sebagai segitiga sudut kanan) dengan dua dimensi catheti memiliki dan d 2 , baik digambarkan sebagai distribusi normal.d1d2

Baik intuisi maupun simulasi memberi tahu saya bahwa distribusi yang dihasilkan normal, dengan rata-rata arctan(mean(d1)mean(d2)) . Tetapi apakah ada cara untuk menghitung distribusi sudut yang dihasilkan? Referensi di mana saya akan menemukan jawabannya?

(Untuk sedikit konteks, saya sedang mengerjakan toleransi statistik bagian-bagian mekanik. Dorongan pertama saya adalah mensimulasikan seluruh proses, memeriksa apakah hasil akhirnya cukup normal, dan menghitung standar deviasi. Tapi saya bertanya-tanya jika mungkin ada pendekatan analitis yang lebih rapi.)

Bossykena
sumber
2
Bisakah Anda mengonfirmasi bahwa (a) d1 dan d2 adalah panjang sisi (dan bukan sudut); (B) bahwa Anda mengasumsikan sudut di antara mereka adalah sudut yang benar (karena kalau tidak rumus atan diduga); dan (c) Anda tertarik pada distribusi salah satu sudut segitiga siku-siku ini? Juga, mungkin, SD dari setiap distribusi panjang jauh lebih kecil dari yang diharapkan karena segitiga seharusnya tidak memiliki probabilitas yang cukup untuk panjang sisi negatif :-).
whuber
Tepat. Saya telah mengulangi masalahnya agar lebih jelas. Dan ya, SD akan relatif kecil terhadap dimensi.
Bossykena
Menggunakan rumus untuk perkalian dan penambahan, Anda dapat mencoba ekspansi Taylor.
Terima kasih atas kedua jawaban Anda yang luar biasa, yang (sejauh yang dapat saya katakan dengan keahlian statistik saya yang terbatas) bersifat intuitif dan baik.
Bossykena

Jawaban:

15

Dalam interpretasi ini, segitiga adalah segitiga siku-siku panjang sisi dan Y didistribusikan secara biner dengan harapan μ x dan μ y , standar deviasi σ x dan σ y , dan korelasi ρ . Kami mencari distribusi arctan ( Y / X ) . Untuk tujuan ini, standarisasi X dan Y sehinggaXYμxμyσxσyρarctan(Y/X)XY

dan Y = σ y η + μ y

X=σxξ+μx
Y=σyη+μy

dengan dan η varian normal standar dengan korelasi ρ . Biarkan θ menjadi sudut dan untuk kenyamanan tulis q = tan ( θ ) . Kemudianξηρθq=tan(θ)

P[arctan(Y/X)θ]=P[YqX]

=P[σyη+μyq(σxξ+μx)

=P[σyηqσxξqμxμy]

Sisi kiri, yang merupakan kombinasi linear dari Normals, adalah normal, dengan rata-rata dan varians σ 2 y + q 2 σ 2 x - 2 q ρ σ x σ y . μyσyqμxσxσy2+q2σx22qρσxσy

Membedakan cdf Normal dari parameter ini sehubungan dengan θ yields the pdf of the angle. The expression is fairly grisly, but a key part of it is the exponential

exp((μy(σy+1)μx(σx+1)tan(θ))22(2ρσxσytan(θ)+σx2+σy2+tan2(θ))),

showing right away that the angle is not normally distributed. However, as your simulations show and intuition suggests, it should be approximately normal provided the variations of the side lengths are small compared to the lengths themselves. In this case a Saddlepoint approximation ought to yield good results for specific values of μx, μy, σx, σy, and ρ, even though a closed-form general solution is not available. The approximate standard deviation will drop right out upon finding the second derivative (with respect to θ) of the logarithm of the pdf (as shown in equations (2.6) and (3.1) of the reference). I recommend a computer algebra system (like MatLab or Mathematica) for carrying this out!

whuber
sumber
1
There was never any chance of it being normally distributed. It is an angle! It only takes values on [π,π).
Robby McKilliam
1
P(Y/X q) = P(Y qX) is not correct if X is a normal r.v. - X can be negative too.
ronaf
@ronaf: actually, since X and Y are the side lengths of a physical triangle, we should not have negative X!
shabbychef
2
@ronaf: That's the right idea. If one uses signed side lengths and also considers the angle as a real value (rather than its value modulo 2π), there is no inconsistency with normality in either case. Your point about the inequality possibly being wrong is excellent. All I can do in response is to claim that the equation is an excellent approximation under the assumptions made because the chance of X or Y being negative is negligible.
whuber
1
@YBE I agree that the last "+" in my expression looks like it doesn't belong--it might have slipped in when I was cleaning up the TeX markup. I don't have a reference because I computed the derivative myself.
whuber
12

You are looking at circular statistics and in particular a circular distribution called the projected normal distribution.

For some reason this topic can be a little hard to google, but the two major texts on circular statistics are The Statistical Analysis of Circular Data by Fisher and Directional Statistics by Mardia and Jupp.

For a thorough analysis of the projected normal distribution see page 46 of Mardia and Jupp. There are closed form expressions (up to the error function integral) for the distribution, and as whuber has suggested, it looks similar to the normal when its `variance' (careful here, what does variance mean for a random variable on a circle?!) is small, i.e. when the distribution is quite concentrated at one point (or direction or angle).

Robby McKilliam
sumber