Rasgele bir değişkenin uzun dikdörtgen bir matrisle doğrusal dönüşümü


12

Diyelim ki , olasılık yoğunluk fonksiyonuna sahip bir dağılımdan alınan rastgele bir . Biz doğrusal bir tam-sıraya göre dönüştürmek ise matris almak , daha sonra yoğunluğu verilirXRnfX(x)n×nAY=AXY

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

Şimdi dönüştürmek ki X yerine göre m×n matris B ile, m>n veren Z=BX . Açıkça ZRm , ancak n boyutlu bir altuzay G \ alt kümesi \ mathbb {R} ^ m "yaşıyor" GRm. G'de olduğunu bildiğimiz için \ vec {Z} ' nin koşullu yoğunluğu nedir ?ZG

İlk içgüdüm B'nin sahte tersini kullanmaktı B. Eğer B=USVT tekil değer ayrışımı olduğunu B , daha sonra B+=VS+UT yalancı ters olduğu S+ diyagonal matris sıfır olmayan girişleri çevrilmesi ile oluşturulur S . Bunun

fZ(z)=1|det+S|fX(B+z),
burada det+S sıfır olmayan tekil değerlerin ürünü anlamına gelir.

Bu akıl yürütme, burada verilen ve burada ve bu CrossValidated mesajında belirtilen tek bir normalin (değişkenin uygun alt alanda yaşadığı bilgisine bağlı olarak) yoğunluğu ile uyumludur .

Ama bu doğru değil! Normalleştirme sabiti kapalı. Aşağıdaki durum dikkate alınarak bir (önemsiz) karşı örnek verilir: XN(0,1) ,

Y=(11)X=(XX).
Burada B matrisi Byukarıdakiler sadece vektörlerdir. Yalancı tersi
B+=(1/21/2)
ve det+B=2 . Yukarıdaki mantık
fY(y)=12π2exp(12yT(B+)TB+y),
ancak bu aslında ( y = x satırında y=x) \ frac {1} {\ sqrt ile bütünleşir 12. Bu durumda, \ vec {Y} girişlerinden birini bırakabileceğinizi Yanlıyorum, ancak B çok daha büyük olduğunda , bırakılacak giriş kümesini tanımlamak sinir bozucu. Sahte ters akıl yürütme neden çalışmıyor? Bir dizi rasgele değişkenin "uzun" bir matrisle doğrusal dönüşümünün yoğunluk fonksiyonu için genel bir formül var mı? Herhangi bir referans da büyük takdir edilecektir.

Yanıtlar:


2

Gelecekte bununla karşılaşabilecek olanlar için ... hatanın kaynağı aslında entegrasyondan kaynaklanmaktadır. Yukarıdaki örnekte, entegrasyon çizgisi üzerinden gerçekleşir . Bu nedenle, integrali alırken çizgiyi "parametrelendirmek" ve parametreleştirmenin Jacobian'ını düşünmek gerekir, çünkü eksenindeki her birim adımı satırdaki uzunluk adımlarına karşılık gelir . Örtük olarak kullandığım parametreleştirme , diğer bir deyişle nin her iki özdeş girişini de değere göre belirterek verildi . Bu Jacobi sahip düzgün ile iptal,y=xx2x(x,x)y22 (tam olarak aynı Jacobian'dan).

Örnek yapay olarak basitti - genel bir dönüşüm için, problem bağlamında doğal olan çıktı için başka bir parametrelendirme olabilir. Parametrizasyonu aynı alt uzay kapsamalıdır yana olarak ve bu alt uzay bir hiper olup, parametreleştirme kendisi muhtemelen doğrusal olmaktır. Parametrelendirme nin matris temsili çağrıldığında , gereksinim basitçe aynı sütun alanına sahip olmasıdır (aynı hiper düzlemi kapsar). Sonra nihai yoğunlukBGBm×nLB

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

Genel olarak, bu kurulum biraz garip ve bence yapılacak en doğru şey, doğrusal olarak bağımsız bir satırı kümesi bulmak ve satırların geri kalanını kaldırmak (dönüştürülmüş değişkenin karşılık gelen bileşenleri ile birlikte ) bir kare matris elde etmek için . Sonra sorun tam sıralı durumuna düşer ( tam sütun sıralamasına sahip olduğu varsayılarak ).BzB^n×nB

Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.