Operações trigonométricas em desvios padrão

14

Adição, subtração, multiplicação e divisão de variáveis ​​aleatórias normais estão bem definidas, mas e as operações trigonométricas?

Por exemplo, suponhamos que estou tentando encontrar o ângulo de uma cunha triangular (modelada como um triângulo de ângulo reto) com os dois catetos com as dimensões d1 e , ambos descritos como distribuições normais.d2

Tanto a intuição quanto a simulação me dizem que a distribuição resultante é normal, com média . Mas existe uma maneira de calcular a distribuição do ângulo resultante? Referências sobre onde eu encontraria a resposta?arctan(mean(d1)mean(d2))

(Para um pouco de contexto, estou trabalhando na tolerância estatística de peças mecânicas. Meu primeiro impulso seria simplesmente simular todo o processo, verificar se o resultado final é razoavelmente normal e calcular o desvio padrão. Mas estou pensando se houver uma abordagem analítica mais clara.)

Bossykena
fonte
2
Você poderia confirmar que (a) d1 e d2 são os comprimentos laterais (e não os ângulos); (b) que você está assumindo que o ângulo entre eles é um ângulo reto (caso contrário, a fórmula atan é suspeita); e (c) que você está interessado na distribuição de um dos outros ângulos deste triângulo retângulo? Além disso, presumivelmente, o DP de cada distribuição de comprimento é muito menor que a expectativa, porque o triângulo não deve ter nenhuma probabilidade apreciável de um comprimento lateral negativo :-).
whuber
Exato. Eu reformulei o problema para torná-lo um pouco mais claro. E sim, o SD será pequeno em relação às dimensões.
Bossykena 30/09/10
Usando fórmulas para multiplicação e adição, você pode tentar a expansão de Taylor.
Obrigado por suas excelentes respostas, que (até onde posso dizer com meu conhecimento limitado em estatísticas) são intuitivas e válidas.
Bossykena 1/10/10

Respostas:

15

Nesta interpretação, o triângulo é um triângulo retângulo de comprimentos laterais e Y distribuídos binormalmente com expectativas μ x e μ y , desvios padrão σ x e σ y e correlação ρ . Buscamos a distribuição de arctan ( Y / X ) . Para esse fim, padronize X e Y para queXYμxμyσxσyρarctan(Y/X)XY

e Y = σ y η + μ y

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

com e η variáveis ​​normais padrão com correlação ρ . Seja θ um ângulo e, por conveniência, escreva q = tan ( θ ) . Entãoξηρθq=tan(θ)

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

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

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

O lado esquerdo, sendo uma combinação linear de normais, é normal, com média e variação σ 2 y + q 2 σ 2 x - 2 q ρ σ x σ y . μyσyqμxσxσy2+q2σx22qρσxσy

A diferenciação do cdf normal desses parâmetros em relação a produz o pdf do ângulo. A expressão é bastante pavorosa, mas uma parte importante é a exponencialθ

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

mostrando imediatamente que o ângulo não é normalmente distribuído. No entanto, como as simulações mostram e a intuição sugere, deve ser aproximadamente normal, desde que as variações dos comprimentos laterais sejam pequenas em comparação com os comprimentos propriamente ditos. Nesse caso, uma aproximação do ponto de sela deve produzir bons resultados para valores específicos de , μ y , σ x , σ y e ρ , mesmo que uma solução geral de forma fechada não esteja disponível. O desvio padrão aproximado será eliminado logo após encontrar a segunda derivada (com relação a θμxμyσxσyρθ) do logaritmo do pdf (como mostrado nas equações (2.6) e (3.1) da referência). Eu recomendo um sistema de álgebra computacional (como MatLab ou Mathematica) para realizar isso!

whuber
fonte
1
Nunca houve chance de distribuição normal. É um ângulo! Leva apenas valores em . [π,π)
Robby McKilliam
1
P (Y / X q) = P (Y qX) não está correto se X é um rv normal - X também pode ser negativo.
Ronaf 1/10/10
@ronaf: na verdade, já que e Y são os comprimentos laterais de um triângulo física, devemos não têm negativo X ! XYX
Shabbychef
2
@ronaf: Essa é a ideia certa. Se alguém usa comprimentos laterais assinados e também considera o ângulo como um valor real (em vez de seu valor módulo ), não há inconsistência com a normalidade em nenhum dos casos. Seu ponto de vista sobre a desigualdade possivelmente estar errada é excelente. Tudo o que posso fazer em resposta é afirmar que a equação é uma excelente aproximação sob as suposições feitas, porque a chance de X ou Y ser negativo é desprezível. 2π
whuber
1
@YBE Concordo que o último "+" na minha expressão parece não pertencer - ele pode ter surgido quando eu estava limpando a marcação TeX. Eu não tenho uma referência porque eu mesmo calculei a derivada.
whuber
12

Você está olhando para estatísticas circulares e, em particular, uma distribuição circular chamada distribuição normal projetada .

Por alguma razão, esse tópico pode ser um pouco difícil de pesquisar no Google, mas os dois principais textos sobre estatística circular são The Statistical Analysis of Circular Data by Fisher e Directional Statistics por Mardia e Jupp.

Para uma análise completa da distribuição normal projetada, consulte a página 46 de Mardia e Jupp. Existem expressões de formulário fechado (até a integral da função de erro) para a distribuição e, como sugeriu o whuber, ele se parece com o normal quando sua `variância '(cuidado aqui, o que significa variância para uma variável aleatória em um círculo? !) é pequeno, ou seja, quando a distribuição está bastante concentrada em um ponto (ou direção ou ângulo).

Robby McKilliam
fonte