2013-04-17 6 views
7

में यादृच्छिक नमूने आकर्षित करने के लिए एक कस्टम संभावना वितरण बनाना मैं एक मोंटेकारलो प्रकार सिमुलेशन का उपयोग करके चीजों की संभाव्यता वितरण की मनमानी संख्या को जोड़ना चाहता हूं। मैं यादृच्छिक रूप से कुछ के निरंतर वितरण का नमूना देना चाहता हूं और उन्हें अन्य निरंतर वितरण के अन्य यादृच्छिक नमूने में जोड़ना चाहता हूं, अंततः उनके संयोजन के लिए संभावित वितरण प्राप्त करना। वितरण स्वयं अनुभवजन्य हैं - वे एक कार्य नहीं हैं लेकिन पी 99 = 2.4, पी 0 9 = 7.12, पी 50 = 24.53, पी 10 = 82.14 के रूप में और इसी तरह (वास्तव में उन बिंदुओं का एक गुच्छा है)। वितरण कम या ज्यादा असामान्य हैं, इसलिए यदि आवश्यक हो तो उन्हें लॉगनोर्मल के रूप में अनुमानित करना शायद ठीक होगा। लेकिन मैं इसे साइपी के lognorm function में कैसे दर्ज कर सकता हूं? या इसे साइपी, या सामान्य रूप से अजगर में कुछ और तरीका है?साइपी

मुझे उम्मीद है कि यह स्पष्ट है कि मैं क्या करने की कोशिश कर रहा हूं। बहुत बहुत धन्यवाद, एलेक्स

उत्तर

2

ऐसा लगता है कि आपके पास अनिवार्य रूप से संभावना घनत्व का हिस्टोग्राम है। एक चीज जो आप कर सकते हैं उसके बाद अपने अनुभवजन्य वितरण के साथ inverse transform sampling का उपयोग करना है।

एक विकल्प के रूप में, यदि आप एक वितरण (लॉगऑनम या किसी अन्य एक) के एक निश्चित कार्यात्मक रूप की अपेक्षा करते हैं, तो आप संबंधित कार्यात्मक रूप से डेटा को फ़िट करने का प्रयास कर सकते हैं।

+0

हम्म अच्छी तरह से मैं उन्हें लॉगऑनॉर्मल के करीब होने की उम्मीद करता हूं, क्या आप जानते हैं कि मैं डेटा को फ़िट करने के बारे में कैसे जाऊंगा? मुझे [विवरण पृष्ठ] (http://docs.scipy.org/doc/scipy/reference/generated/scipy.stats.lognorm.html#scipy पर दो दिए गए बिंदुओं के माध्यम से एक असामान्य वितरण को फिट करने के तरीके के बारे में उलझन में छोड़ दिया गया था .stats.lognorm)। अगर मेरे पास पी 0 9 था जो एक्स और पी 10 था जो वाई था, तो मैं दोबारा फिट करने वाला एक असामान्य वितरण कैसे करूं? –

+0

मुझे इस दो-बिंदु वाली चीज़ के बारे में निश्चित नहीं है जिसके बारे में आप बात कर रहे हैं। मैं क्या करूँगा, मैं पहले जांच करता हूं कि वितरण वास्तव में * लॉग-सामान्य है या नहीं। इसके लिए, मैं इस तथ्य का उपयोग करता हूं कि सामान्य वितरण के लिए क्षणों को एक बहुत ही सटीक तरीके से (सीएफ विकिपीडिया या मैथवर्ल्ड या अन्यत्र) से संबंधित होना चाहिए, और यदि 'एक्स' लॉग-सामान्य है, तो 'लॉग एक्स' है आम तौर पर वितरित: 'लॉग एक्स' के कई पहले पलों की गणना करें। –

+0

ठीक है, लेकिन मुझे इसे scipy ऑब्जेक्ट में लोड नहीं करना है ताकि मैं ऊपर से लॉगऑर्म फ़ंक्शन पेज में उल्लिखित 'R = lognorm.rvs (s, size = 100) 'का उपयोग करके यादृच्छिक नमूने प्राप्त कर सकूं? –

संबंधित मुद्दे