दो मापदंडों के उत्पाद के लिए आत्मविश्वास अंतराल


11

मान लें कि हमारे पास दो पैरामीटर हैं, और । हमारे पास इन मापदंडों के लिए दो अधिकतम संभावना अनुमानक और और दो आत्मविश्वास अंतराल हैं। क्या लिए एक विश्वास अंतराल बनाने का एक तरीका है ?पी 2 ^ पी 1 ^ पी 2 पी 1 पी 2p1p2p1^p2^p1p2

जवाबों:


13

आप की मानक त्रुटि की गणना करने के लिए डेल्टा विधि का उपयोग कर सकते हैं । डेल्टा विधि बताती है कि एक फ़ंक्शन के विचरण का एक अनुमान इस प्रकार दिया गया है: दूसरी ओर की अपेक्षा का सन्निकटन इसके द्वारा दिया गया है: तो अपेक्षा केवल कार्य है। आपका फ़ंक्शन है: । की उम्मीद बस होगी:p1^p2^g(t)

Var(g(t))i=1kgi(θ)2Var(ti)+2i>jgi(θ)gj(θ)Cov(ti,tj)
g(t)
E(g(t))g(θ)
g(t)g(p1,p2)=p1p2g(p1,p2)=p1p2p1p2 । विचरण के लिए, हमें के आंशिक व्युत्पत्ति की आवश्यकता है : g(p1,p2)
p1g(p1p2)=p2p2g(p1p2)=p1

ऊपर दिए गए विचरण के लिए फ़ंक्शन का उपयोग करना, हमें मिलता है:

Var(p1^p2^)=p2^2Var(p1^)+p1^2Var(p2^)+2p1^p2^Cov(p1^,p2^)
मानक त्रुटि तो बस उपरोक्त अभिव्यक्ति का वर्गमूल होगा। एक बार जब आप मानक त्रुटि प्राप्त कर लेते हैं, तो यह सीधे सीधे लिए 95% विश्वास अंतराल की गणना करने के लिए आगे होता है :p1^p2^p1^p2^±1.96SE^(p1^p2^)

की मानक त्रुटि को शांत करने के लिए , आपको और के विचरण की आवश्यकता होती है, जिसे आप आमतौर पर प्राप्त कर सकते हैं द्वारा विचरण-सहप्रसरण मैट्रिक्स जो आपके मामले में एक 2x2-मैट्रिक्स होगा क्योंकि आप दो अनुमान हैं। विचरण-सहसंयोजक मैट्रिक्स में विकर्ण तत्व और के भिन्न रूप हैं, जबकि ऑफ-विकर्ण तत्व covariance of और (मैट्रिक्स सममित है)। जैसा कि @gung टिप्पणियों में उल्लेख करते हैं, अधिकांश सांख्यिकीय सॉफ्टवेयर्स द्वारा विचरण-कोविरियन मैट्रिक्स को निकाला जा सकता है। कभी-कभी, अनुमान एल्गोरिदम प्रदान करते हैंp1^p2^p1^p2^ Σ ^ पी 1 ^ पी 2 ^ पी 1 ^ पी 2 Σp1^p2^p1^p2^हेसियन मैट्रिक्स (मैं यहाँ उस बारे में विवरण में नहीं जाऊंगा), और विचरण-सहसंयोजक मैट्रिक्स का अनुमान नकारात्मक हेसियन के व्युत्क्रम से लगाया जा सकता है (लेकिन केवल तभी जब आप लॉग-लाइबिलिटी को अधिकतम कर रहे हैं! ( इस पोस्ट को देखें )। फिर से, अपने सांख्यिकीय सॉफ़्टवेयर और / या वेब को हेसियन निकालने के तरीके पर और मैट्रिक्स के व्युत्क्रम की गणना करने के तरीके के बारे में सलाह लें।

वैकल्पिक रूप से, आप निम्न अंतराल में आत्मविश्वास अंतरालों से और संस्करण प्राप्त कर सकते हैं (यह 95% के लिए मान्य है): । के लिए एक -CI, अनुमानित मानक त्रुटि है: , जहां है (के लिए मानक सामान्य वितरण की quantile , )। फिर,p1^p2^SE(p1^)=(upper limitlower limit)/3.92100(1α)%SE(p1^)=(upper limitlower limit)/(2z1α/2)z1α/2(1α/2)α=0.05z0.9751.96Var(p1^)=SE(p1^)2। वही के विचरण के लिए सही है । हमें और भी देखने की आवश्यकता है (पैराग्राफ देखें)। यदि और स्वतंत्र हैं, तो सहसंयोजक शून्य है और हम पद छोड़ सकते हैं।p2^p1^p2^p1^p2^

यह पत्र अतिरिक्त जानकारी प्रदान कर सकता है।


4
+1। मापदंडों और उनके - मैट्रिक्स के जांच करके पाया जा सकता है , जो कि अधिकांश सांख्यिकीय सॉफ़्टवेयर प्रदान कर सकते हैं। जैसे, R में, यह है ; vcov ; & SAS में, PROC REG में मॉडल स्टेटमेंट के विकल्प के रूप में जोड़ा जाता हैβcovb
गंग - मोनिका

1
@ गूँग पेन्ट्री के एक बिंदु पर यह इंगित करने के लायक हो सकता है (क्योंकि मुझे पता है कि यह कुछ लोगों को भ्रमित करता है) कि यह वास्तव में बजाय ( ) का प्रसरण-सहसंयोजक मैट्रिक्स है (और वास्तव में यह वास्तव में ऐसा नहीं है , क्योंकि मानक विचलन का नमूना से अनुमान लगाया जाना है, इसलिए यह वास्तव में अनुमानित विचरण-सहसंयोजक मैट्रिक्स है ..) बीटाβ^β
सिल्वरफ़िश

3
@ सिल्वरफ़िश, विधिवत पीछा किया। अगली बार मैं कहूँगा " - मैट्रिक्स ।" β^
गूँज - मोनिका

1
आप एक प्रोफ़ाइल संभावना समारोह बनाने की कोशिश कर सकते हैं! और उस से विश्वास अंतराल का निर्माण।
kjetil b halvorsen

ना के बाद से यह एक पैरामीटर है? var(p1)=0
user0

1

मुझे उत्पाद के विचरण की गणना के लिए एक अलग समीकरण मिला।

यदि x और y स्वतंत्र रूप से वितरित किए जाते हैं, तो उत्पाद का विचरण अपेक्षाकृत सीधा होता है: V (x * y) = V (y) * E (x) ^ 2 + V (x) * E (y) ^ 2 + V ( x) * V (y) ये परिणाम तीन या अधिक चर (गुडमैन 1960) से जुड़े मामलों को भी सामान्य करते हैं। स्रोत: रेगुलेटिंग पेस्टीसाइड्स (1980), परिशिष्ट एफ

Coolserdash: अंतिम घटक V (x) * V (y) आपके समीकरण में गायब है। क्या संदर्भित पुस्तक (कीटनाशकों का विनियमन) गलत है?

इसके अलावा, दोनों समीकरण सही नहीं हो सकते हैं। " ... हम दिखाते हैं कि तीन स्वतंत्र सामान्य चर के उत्पाद का वितरण सामान्य नहीं है ।" ( स्रोत )। मैं दो सामान्य रूप से वितरित चर के उत्पाद में भी कुछ सकारात्मक तिरछा होने की उम्मीद करूंगा ।


0
  1. CI / 2 / 1.96 = se की लंबाई, यानी A या B की मानक त्रुटि
  2. se ^ 2 = var, यानी अनुमान A या B का प्रसरण
  3. अनुमानित A या B का उपयोग A या B, अर्थात E (A) या E (B) के साधन के रूप में करें
  4. इस पृष्ठ का अनुसरण करें http://falkenblog.blogspot.se/2008/07/formula-for-varxy.html var (A * B), अर्थात var (C) प्राप्त करने के लिए
  5. वर्जन (C) का वर्गमूल C का से है
  6. (C - 1.96 * se (C), C + 1.96 * se (C)) C का 95% CI है

ध्यान दें कि यदि आपके ए और बी सहसंबद्ध हैं, तो आपको उनके सहवास पर भी विचार करने की आवश्यकता है।

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