आप एक परिवर्तित रैंडम वेरिएबल के निम्न क्रम के क्षणों का अनुमान प्राप्त करने के लिए टेलर सीरीज़ का उपयोग कर सकते हैं। यदि वितरण माध्य (एक विशेष अर्थ में) के आसपास काफी 'तंग' है, तो अनुमान बहुत अच्छा हो सकता है।
इसलिए उदाहरण के लिए
g(X)=g(μ)+(X−μ)g′(μ)+(X−μ)22g′′(μ)+…
इसलिए
Var[g(X)]===Var[g(μ)+(X−μ)g′(μ)+(X−μ)22g′′(μ)+…]Var[(X−μ)g′(μ)+(X−μ)22g′′(μ)+…]g′(μ)2Var[(X−μ)]+2g′(μ)Cov[(X−μ),(X−μ)22g′′(μ)+…]+Var[(X−μ)22g′′(μ)+…]
अक्सर केवल पहला शब्द लिया जाता है
Var[g(X)]≈g′(μ)2Var(X)
इस मामले में (मान लें कि मैंने कोई गलती नहीं की), , ।g(X)=1XVar[1X]≈1μ4Var(X)
विकिपीडिया: बेतरतीब चर के कार्यों के लिए टेलर विस्तार
---
इसे उदाहरण देने के लिए कुछ उदाहरण। मैं आर में दो (गामा-वितरित) नमूने उत्पन्न करूंगा, एक माध्य के बारे में 'नॉट-सो-टाइट' वितरण के साथ और थोड़ा हल्का।
a <- rgamma(1000,10,1) # mean and variance 10; the mean is not many sds from 0
var(a)
[1] 10.20819 # reasonably close to the population variance
सन्निकटन से पता चलता है कि करीब होना चाहिए1/a(1/10)4×10=0.001
var(1/a)
[1] 0.00147171
बीजगणितीय गणना में कहा गया है कि वास्तविक जनसंख्या विचरण1/648≈0.00154
अब तंग एक के लिए:
a <- rgamma(1000,100,10) # should have mean 10 and variance 1
var(a)
[1] 1.069147
सन्निकटन से पता चलता है कि करीब होना चाहिए1/a(1/10)4×1=0.0001
var(1/a)
[1] 0.0001122586
बीजीय गणना से पता चलता है कि पारस्परिक की जनसंख्या भिन्नता ।102992×98≈0.000104