Lineare Transformation einer Zufallsvariablen durch eine hohe rechteckige Matrix

12

Nehmen wir an, wir haben einen Zufallsvektor , der aus einer Verteilung mit der Wahrscheinlichkeitsdichtefunktion f X ( x ) gezogen wird . Wenn wir es linear durch eine n × n- Matrix A mit vollem Rang transformieren , um Y = A X zu erhalten , ist die Dichte von Y gegeben durch f Y ( y ) = 1XRnfX(x)n×nAY=AXY

fY(y)=1|detA|fX(A1y).

Nehmen wir nun an, wir transformieren stattdessen durch eine m × n- Matrix B mit m > n , was Z = B X ergibt . Klar Z R m , aber es "lebt" auf einem n- dimensionalen Unterraum G R m . Was ist die bedingte Dichte von Z , wenn wir wissen, dass es in G liegt ?Xm×nBm>nZ=BXZRmnGRmZG

Mein erster Instinkt war, die Pseudo-Inverse von . Wenn B = U S V T die Singularwertzerlegung von B ist , dann ist B + = V S + U T die Pseudo-Inverse, wobei S + durch Invertieren der Nicht-Null-Einträge der Diagonalmatrix S gebildet wird . Ich vermutete, dass dies f Z ( z ) = 1 ergeben würdeBB=USVTBB+=VS+UTS+S wobei durchdet+

fZ(z)=1|det+S|fX(B+z),
das Produkt der Singularwerte ungleich Null gemeint ist.det+S

Diese Argumentation stimmt mit der hier angegebenen und erwähnten Dichte für eine singuläre Normale überein (abhängig von dem Wissen, dass die Variable auf dem entsprechenden Unterraum lebt) auch hier und in diesem CrossValidated-Beitrag erwähnt wird .

Aber es ist nicht richtig! Die Normalisierungskonstante ist ausgeschaltet. Ein (triviales) Gegenbeispiel ergibt sich aus folgendem Fall: Mit sei Y = ( 1 1 ) X = ( X.XN(0,1) Hier ist die MatrixBvon oben nur der Ein-Vektor. Seine pseudoinversen istB+=( 1 / 2 1 / 2 ) unddet+B=

Y=(11)X=(XX).
B
B+=(1/21/2)
. Die Argumentation von oben würdefY (y )= 1 vorschlagendet+B=2 aber dies integriert tatsächlich (auf der Liniey=x) zu1
fY(y)=12π2exp(12yT(B+)TB+y),
y=x . Mir ist klar, dass Sie in diesem Fall einfach einen der Einträge vonYlöschen können, aber wennBviel größer ist, ist es ärgerlich, die zu löschenden Einträge zu identifizieren. Warum funktioniert das pseudo-inverse Denken nicht? Gibt es eine allgemeine Formel für die Dichtefunktion einer linearen Transformation einer Menge von Zufallsvariablen durch eine "große" Matrix? Alle Referenzen wäre auch sehr dankbar.12YB
Dan
quelle

Antworten:

2

y=xx2x(x,x)y22 (kommt aus genau dem gleichen Jacobian) im Nenner.

BGBm×nLB

fZ(z)=|det+L||det+B|fX(B+z).

BzB^n×nB

Dan
quelle