एक लंबा आयताकार मैट्रिक्स द्वारा एक यादृच्छिक चर का रैखिक परिवर्तन


12

मान लें कि हमारे पास एक यादृच्छिक वेक्टर , जो वितरण घनत्व से संभाव्यता घनत्व फ़ंक्शन । यदि हम इसे प्राप्त करने के लिए इसे पूर्ण-श्रेणी मैट्रिक्स से बदलते हैं , तो का घनत्व द्वारा दिया जाता है।XRnfX(x)n×nAY=AXY

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

अब कहते हैं कि हम बदलना X बजाय एक से m×n मैट्रिक्स B , के साथ m>n , दे रही है Z=BX । स्पष्ट रूप से ZRm , लेकिन यह "एक n -dimensional उप-स्थान GRm\ Vec {Z} का सशर्त घनत्व क्या है Z, यह देखते हुए कि हम जानते हैं कि यह G ?

मेरी पहली वृत्ति बी के छद्म-विलोम का उपयोग करना था B। यदि B = USV ^ T , BB=USVT का एकवचन मान अपघटन है , तो B ^ + = VS ^ + U ^ T छद्म व्युत्क्रम है, जहां S ^ + का निर्माण द्विवार्षिक मैट्रिक्स S की गैर-शून्य प्रविष्टियों को सम्मिलित करके किया जाता है । मैंने अनुमान लगाया कि यह f_ \ vec {Z} (\ vec {z}) = \ frac {1} {\ left। \ Det ^ + S \ right |} f_ \ vec {X} (B ^ + vec) को देगा। {z}), जहां by \ det ^ + S का अर्थ है गैर-शून्य एकवचन मानों का गुणनफल।BB+=VS+UTS+S

fZ(z)=1|det+S|fX(B+z),
det+S

इस तर्क के लिए घनत्व के साथ सहमत हैं एक विलक्षण सामान्य (ज्ञान पर वातानुकूलित है कि उपयुक्त उपस्पेस पर चर जीवन) दिया यहाँ और उल्लेख यहाँ भी और में इस CrossValidated पोस्ट

लेकिन यह सही नहीं है! सामान्यीकरण स्थिरांक बंद है। A (ट्रिवियल) प्रतिधारण निम्नलिखित मामले पर विचार करके दिया जाता है: XN(0,1) , let

Y=(11)X=(XX).
यहां ऊपर से मैट्रिक्स B सिर्फ वेक्टर है। इसका छद्म उलटा
B+=(1/21/2)
और det+B=2 । ऊपर से तर्क
fY(y)=12π2exp(12yT(B+)TB+y),
लेकिन यह वास्तव में एकीकृत करता है (लाइन y=x ) 12। मुझे लगता है कि इस मामले में आप बस अपने किए गए \ vec {Y} की प्रविष्टियों में से एक को छोड़ Yसकते हैं, लेकिन जब B बहुत बड़ी हो तो ड्रॉप करने के लिए प्रविष्टियों के सेट की पहचान करना कष्टप्रद है। छद्म उलटा तर्क कार्य क्यों नहीं करता है? क्या "लंबा" मैट्रिक्स द्वारा यादृच्छिक चर के सेट के रैखिक परिवर्तन के घनत्व समारोह के लिए एक सामान्य सूत्र है? किसी भी संदर्भ की बहुत सराहना की जाएगी।

जवाबों:


2

उन लोगों के लिए जो भविष्य में इस पर चल सकते हैं ... त्रुटि का स्रोत वास्तव में एकीकरण से उपजा है। ऊपर के उदाहरण में, लाइन पर एकीकरण होता है । इसलिए यह आवश्यक है कि लाइन को "पैरामीरिज़" करें और इंटीग्रल लेते समय पैराब्रिज़ेशन के जैकबियन पर विचार करें, क्योंकि -axis में प्रत्येक इकाई चरण लाइन पर लंबाई चरणों से मेल खाती है । मैं जिस का उपयोग कर रहा था, वह द्वारा दिया गया था , दूसरे शब्दों में मान द्वारा की दोनों समान प्रविष्टियों को निर्दिष्ट करता है । इसमें जैकबियन , जो बड़े करीने से रद्द करता हैy=xx2x(x,x)y22 (वास्तव में एक ही याकूब से आ रहा है) हर में।

उदाहरण कृत्रिम रूप से सरल था - सामान्य परिवर्तन , किसी के पास आउटपुट के लिए एक और पैरामीरिजेशन हो सकता है जो समस्या के संदर्भ में स्वाभाविक है। चूंकि पैरामीरीज़ेशन को रूप में एक ही उप-क्षेत्र को कवर करने की आवश्यकता होती है , और यह उप-क्षेत्र एक हाइपरप्लेन है, पैरामीटर के रैखिक होने की संभावना है। पैरामीरिजेशन के मैट्रिक्स प्रतिनिधित्व को बुलाते हुए , आवश्यकता बस इतनी है कि इसमें के समान स्तंभ स्थान है (समान हाइपरप्लेन को कवर करें)। फिर अंतिम घनत्व हो जाता हैBGBm×nLB

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

सामान्य तौर पर, यह सेटअप एक प्रकार का अजीब है, और मुझे लगता है कि सही काम करना की पंक्तियों का एक अधिकतम रैखिक रूप से स्वतंत्र सेट ढूंढना है , और बाकी पंक्तियों को हटा दें (साथ ही परिवर्तित चर के संगत घटकों के साथ) ) एक वर्ग मैट्रिक्स प्राप्त करने के लिए । फिर समस्या पूर्ण-रैंक केस तक कम हो जाती है (मान लें कि में पूर्ण स्तंभ रैंक है)।BzB^n×nB

हमारी साइट का प्रयोग करके, आप स्वीकार करते हैं कि आपने हमारी Cookie Policy और निजता नीति को पढ़ और समझा लिया है।
Licensed under cc by-sa 3.0 with attribution required.