एन सामान्य आईआईडी के उत्पाद का अनुमानित वितरण? विशेष मामला μ≈0


12

यह देखते हुए आईआईडी , और , की तलाश में:N≥30Xn≈N(μX,σX2)μX≈0

  1. सटीक बंद फ़ॉर्म वितरण सन्निकटन YN=∏1NXn
  2. असममित ( घातांक ?) एक ही उत्पाद का सन्निकटन

यह अधिक सामान्य प्रश्न का एक विशेष केस है ।μX≈0


1. क्या आपके पास और बारे में कोई जानकारी है ? (यह अच्छा होगा यदि सभी , उदाहरण के लिए।) (2) एक एसिम्प्टोटिक सामान्य सन्निकटन भयानक होगा , क्योंकि एसिम्पोटली दूरस्थ रूप से सामान्य नहीं लगेगा। σ एक्स μ एक्स / σ एक्स » 0 YμXσXμX/σX≫0Y
— whuber

मैंने बस इसके साथ एक त्वरित नाटक किया था। यदि आप रुचि रखते हैं, तो यादृच्छिक चर के उत्पाद के लिए एक सटीक बंद फॉर्म समाधान प्राप्त करना संभव है जो कि iid । गैर-शून्य मामला चीजों को और अधिक जटिल बनाता है। एन ( 0 , σ 2 ) μnN(0,σ2)μ
— वुल्फिज

@whuber (1) कुछ अलग-अलग और साथ कुछ मोंटे कार्लो करने के बाद , मैंने पाया कि वितरण लिए अच्छा व्यवहार करता है और | अब मैं और लिए एक अच्छी अभिव्यक्ति खोजना कि कैसे में कुछ अच्छे सन्निकटन हैं। मैंने टेलर विस्तार के माध्यम से कुछ सन्निकटन बनाए, लेकिन उन्होंने बुरी तरह से दुर्व्यवहार किया। (2) अच्छी तरह से, निश्चित रूप से "दिखता है" ची स्क्वायर के साथ सामान्य की राशि की तरह, इसलिए को सामान्य तक कम किया जा सकता है, अगर सन्निकटन "साबित" होता है। σ एफ एन > 30 | μ एक्स | ≥ 10 σ एक्स μ एफ σ एफ χ 2 एफ एफμσFN>30|μX|≥10σXμFσFχ2FF
— आंद्रेई पॉज़ोलोटिन 18

3
जब , को डिस्ट्रीब्यूशन ( बैरी-एसेन प्रमेय के एक एप्लीकेशन के रूप में to शो) द्वारा अच्छी तरह से अनुमानित किया जाएगा । वाई लॉग ( एक्स )μX≥10σXYlog⁡(X)
— whuber

बैरी-एस्सेन के सीधे प्रत्यक्ष आवेदन से , जो वास्तव में अच्छा है, लेकिन यह कुछ संरचना खो देता है: को नकारात्मक होना चाहिए, पर निर्भर होना चाहिए , आदि शायद, इसे लागू करने के बेहतर तरीके हैं? μएफσएफαFN≈0+1NZμFσFα
— आंद्रेई पॉज़ोलोटिन

जवाबों:


10

शून्य-मीन मामले (भाग बी) में एक सटीक समाधान प्राप्त करना संभव है।

समस्या

Let को iid चर, प्रत्येक को सामान्य पीडीएफ साथ :(X1,…,Xn)nN(0,σ2)f(x)

यहाँ छवि विवरण दर्ज करें

हम की पीडीएफ की तलाश , के लिए∏i=1nXin=2,3,…

उपाय

दो ऐसे नॉर्म्स के उत्पाद की पीडीएफ बस है:

यहाँ छवि विवरण दर्ज करें

... जहाँ मैं उपयोग कर रहा हूँ TransformProductसे समारोह mathStatica के लिए पैकेज मेथेमेटिका । समर्थन का डोमेन है:

यहाँ छवि विवरण दर्ज करें

3, 4, 5 और 6 नॉर्म्स का उत्पाद पुनरावृत्तीय रूप से समान फ़ंक्शन (यहां चार बार) लागू करके प्राप्त किया जाता है:

यहाँ छवि विवरण दर्ज करें

... जहां मीजेर जी फ़ंक्शनMeijerG को दर्शाता है

प्रेरण द्वारा, iid के उत्पाद का pdf यादृच्छिक चर है:nN(0,σ2)


1(2π)n2σnMeijerG[{{},{}},{{01,…,0n},{}},x22nσ2n] for x∈R

त्वरित मोंटे कार्लो की जाँच करें

यहाँ एक त्वरित जाँच है:

  • सैद्धांतिक पीडीएफ सिर्फ प्राप्त (जब और ): लाल रंग का वक्रσ = 3n=6σ=3
  • अनुभवजन्य मोंटे कार्लो पीडीएफ के लिए: squiggly BLUE वक्र

यहाँ छवि विवरण दर्ज करें

ठीक लग रहा है! [ब्लू स्क्विग्ली मोंटे वक्र सटीक लाल-धराशायी वक्र अस्पष्ट है]


उत्कृष्ट, धन्यवाद, कॉलिन। अब मैं देखता हूं कि मुझे आपकी पुस्तक क्यों खरीदनी चाहिए :-) मुझे यह भी आश्चर्य होता है कि क्या किसी भी सरल दिखता है। मेरे वुल्फराम कौशल को धूल चटाने का समय। log(...MeijerG(...))
— आंद्रेई पॉज़ोलोटिन
हमारी साइट का प्रयोग करके, आप स्वीकार करते हैं कि आपने हमारी Cookie Policy और निजता नीति को पढ़ और समझा लिया है।
Licensed under cc by-sa 3.0 with attribution required.