Trigonometrische Operationen an Standardabweichungen

14

Addition, Subtraktion, Multiplikation und Division normaler Zufallsvariablen sind gut definiert, aber was ist mit trigonometrischen Operationen?

Nehmen wir zum Beispiel an, ich versuche, den Winkel eines dreieckigen Keils (modelliert als rechtwinkliges Dreieck) mit den beiden Katheten mit den Dimensionen d1 und d2 , die beide als Normalverteilungen beschrieben werden.

Sowohl Intuition als auch Simulation sagen mir, dass die resultierende Verteilung normal ist, mit mittlerem arctan(mean(d1)mean(d2)). Aber gibt es eine Möglichkeit, die Verteilung des resultierenden Winkels zu berechnen? Referenzen, wo ich die Antwort finden würde?

(Für ein bisschen Kontext arbeite ich an der statistischen Toleranz von mechanischen Teilen. Mein erster Impuls wäre, einfach den gesamten Prozess zu simulieren, zu überprüfen, ob das Endergebnis einigermaßen normal ist, und die Standardabweichung zu berechnen. Aber ich frage mich wenn es einen ordentlicheren analytischen Ansatz geben könnte.)

Bossykena
quelle
2
Können Sie bestätigen, dass (a) d1 und d2 die Seitenlängen (und nicht die Winkel) sind? (b) Sie nehmen an, dass der Winkel zwischen ihnen ein rechter Winkel ist (andernfalls ist die atan-Formel verdächtig); und (c) dass Sie an der Verteilung eines der anderen Winkel dieses rechtwinkligen Dreiecks interessiert sind? Vermutlich ist auch die SD jeder Längenverteilung viel kleiner als erwartet, da das Dreieck keine nennenswerte Wahrscheinlichkeit für eine negative Seitenlänge haben sollte :-).
whuber
Genau. Ich habe das Problem umformuliert, um es ein bisschen klarer zu machen. Und ja, die SD wird im Verhältnis zu den Abmessungen klein sein.
Bossykena
Mithilfe von Formeln zur Multiplikation und Addition können Sie die Taylor-Erweiterung ausprobieren.
Vielen Dank für Ihre beiden hervorragenden Antworten, die (soweit ich mit meiner eingeschränkten Statistikkenntnis einschätzen kann) sowohl intuitiv als auch fundiert sind.
Bossykena

Antworten:

15

In dieser Interpretation ist das Dreieck ein rechtwinkliges Dreieck mit den Seitenlängen und Y, die binormal verteilt sind, mit den Erwartungen μ x und μ y , den Standardabweichungen σ x und σ y und der Korrelation ρ . Wir suchen die Verteilung von Arctan ( Y / X ) . Zu diesem Zweck standardisieren Sie X und Y so, dassXYμxμyσxσyρarctan(Y/X)XY

und Y = σ y η + μ y

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

mit und η ändert sich die Normale mit der Korrelation ρ . Sei θ ein Winkel und schreibe der Einfachheit halber q = tan ( θ ) . Dannξηρθq=tan(θ)

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

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

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

Die linke Seite ist eine lineare Kombination von Normals ist, normal ist , mit einem Mittelwert und die Varianz σ 2 y + q 2 σ 2 x - 2 q & rgr; σ x σ y . μyσyqμxσxσy2+q2σx22qρσxσy

Die Differenzierung der Normalen cdf dieser Parameter in Bezug auf ergibt die pdf des Winkels. Der Ausdruck ist ziemlich grausam, aber ein wesentlicher Teil davon ist das Exponentialθ

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

sofort zeigen, dass der Winkel nicht normal verteilt ist. Wie Ihre Simulationen und Ihre Intuition zeigen, sollte dies jedoch annähernd normal sein, vorausgesetzt, die Abweichungen der Seitenlängen sind im Vergleich zu den Längen selbst gering. In diesem Fall sollte eine Sattelpunktnäherung gute Ergebnisse für bestimmte Werte von , μ y , σ x , σ y und ρ liefern , obwohl eine geschlossene allgemeine Lösung nicht verfügbar ist. Die ungefähre Standardabweichung fällt sofort ab, wenn die zweite Ableitung (in Bezug auf θ) gefunden wirdμxμyσxσyρθ) des Logarithmus des pdf (wie in den Gleichungen (2.6) und (3.1) der Referenz gezeigt). Ich empfehle dazu ein Computeralgebrasystem (wie MatLab oder Mathematica)!

whuber
quelle
1
Es gab nie eine Chance, dass es normal verteilt wurde. Es ist ein Winkel! Es werden nur Werte für . [π,π)
Robby McKilliam
1
P (Y / X q) = P (Y qX) ist nicht korrekt, wenn X ein normales rv ist - X kann auch negativ sein.
Ronaf
@ronaf: Da und Y die Seitenlängen eines physischen Dreiecks sind, sollten wir kein negatives X haben ! XYX
Shabbychef
2
@ronaf: Das ist die richtige Idee. Verwendet man vorzeichenbehaftete Seitenlängen und betrachtet den Winkel auch als reellen Wert (statt dessen Wert modulo ), so gibt es in beiden Fällen keine Inkonsistenz mit der Normalität. Ihr Standpunkt, dass die Ungleichung möglicherweise falsch ist, ist ausgezeichnet. Alles, was ich als Antwort tun kann, ist zu behaupten, dass die Gleichung eine ausgezeichnete Annäherung unter den getroffenen Annahmen ist, da die Wahrscheinlichkeit, dass X oder Y negativ sind, vernachlässigbar ist. 2π
Whuber
1
@YBE Ich stimme zu, dass das letzte "+" in meinem Ausdruck nicht dazu gehört - es ist möglicherweise eingedrungen, als ich das TeX-Markup bereinigt habe. Ich habe keine Referenz, weil ich das Derivat selbst berechnet habe.
whuber
12

Es handelt sich um Zirkelstatistiken und insbesondere um eine Zirkelverteilung, die als projizierte Normalverteilung bezeichnet wird .

Aus irgendeinem Grund ist dieses Thema etwas schwer zu googeln, aber die beiden Haupttexte zur Kreisstatistik sind Die statistische Analyse von Kreisdaten nach Fisher und Richtungsstatistik nach Mardia und Jupp.

Eine gründliche Analyse der projizierten Normalverteilung finden Sie auf Seite 46 von Mardia und Jupp. Es gibt geschlossene Ausdrücke (bis auf das Fehlerfunktionsintegral) für die Verteilung, und wie wir bereits vermutet haben, sieht sie ähnlich aus wie normal, wenn ihre "Varianz" (Vorsicht hier, was bedeutet Varianz für eine Zufallsvariable auf einem Kreis? !) ist klein, dh wenn die Verteilung auf einen Punkt (oder Richtung oder Winkel) konzentriert ist.

Robby McKilliam
quelle