Addition, subtraction, multiplication and division of normal random variables are well defined, but what about trigonometric operations?

For instance, let us suppose that I’m trying to find the angle of a triangular wedge (modelled as a right-angle triangle) with the two catheti having dimensions d1 and d2, both described as normal distributions.

Both intuition and simulation tell me that the resulting distribution is normal, with mean arctan(mean(d1)mean(d2)). But is there a way to compute the distribution of the resulting angle? References on where I’d find the answer?

(For a bit of context, I’m working on statistical tolerance of mechanical parts. My first impulse would be to simply simulate the whole process, check if the end result is reasonably normal, and compute the standard deviation. But I’m wondering if there might be a neater analytical approach.)

**Answer**

In this interpretation, the triangle is a right triangle of side lengths X and Y distributed binormally with expectations μx and μy, standard deviations σx and σy, and correlation ρ. We seek the distribution of arctan(Y/X). To this end, standardize X and Y so that

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

with ξ and η standard normal variates with correlation ρ. Let θ be an angle and for convenience write q=tan(θ). Then

P[arctan(Y/X)≤θ]=P[Y≤qX]

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

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

The left hand side, being a linear combination of Normals, is normal, with mean μyσy−qμxσx and variance σ2y+q2σ2x−2qρσxσy.

Differentiating the Normal cdf of these parameters with respect to θ 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(θ)+σ2x+σ2y+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!

**Attribution***Source : Link , Question Author : Bossykena , Answer Author : whuber*