Dans cette interprétation, le triangle est un triangle rectangle de longueurs latérales et Y distribuées binormalement avec des attentes μ x et μ y , des écarts types σ x et σ y et une corrélation ρ . Nous recherchons la distribution de l' arctan ( Y / X ) . À cette fin, normalisez X et Y afin queXYμxμyσxσyρarctan(Y/X)XY
et Y = σ y η + μ y
X=σxξ+μx
Y=σyη+μy
avec et η la normale normale varie avec la corrélation ρ . Soit θ un angle et pour plus de commodité, écrivez q = tan ( θ ) . alorsξηρθq=tan(θ)
P[arctan(Y/X)≤θ]=P[Y≤qX]
=P[σyη+μy≤q(σxξ+μx)
=P[σyη−qσxξ≤qμx−μy]
Le côté gauche, étant une combinaison linéaire de normales, est normal, avec une moyenne et une variance σ 2 y + q 2 σ 2 x - 2 q ρ σ x σ y . μyσy−qμxσxσ2y+q2σ2x−2qρσxσy
La différenciation du cdf normal de ces paramètres par rapport à donne le pdf de l'angle. L'expression est assez macabre, mais une partie clé de celle-ci est l'exponentielleθ
exp(−(μy(σy+1)−μx(σx+1)tan(θ))22(−2ρσxσytan(θ)+σ2x+σ2y+tan2(θ))),
montrant tout de suite que l'angle n'est pas normalement distribué. Cependant, comme le montrent vos simulations et votre intuition, cela devrait être à peu près normal à condition que les variations des longueurs des côtés soient faibles par rapport aux longueurs elles-mêmes. Dans ce cas, une approximation Saddlepoint devrait donner de bons résultats pour des valeurs spécifiques de , μ y , σ x , σ y et ρ , même si une solution générale de forme fermée n'est pas disponible. L'écart type approximatif disparaîtra dès la découverte de la dérivée seconde (par rapport à θμxμyσxσyρθ) du logarithme du pdf (comme indiqué dans les équations (2.6) et (3.1) de la référence). Je recommande un système d'algèbre informatique (comme MatLab ou Mathematica) pour réaliser cela!