मैं वर्तमान में R से JAGS का उपयोग करके पदानुक्रमित बायेसियन मॉडल सीख रहा हूं, और पायथन ( "हैकर्स के लिए बायेसियन मेथड्स" ) का उपयोग करके भी pymc ।
मैं इस पोस्ट से कुछ अंतर्ज्ञान प्राप्त कर सकता हूं : "आप संख्याओं के ढेर के साथ समाप्त हो जाएंगे जो दिखता है" जैसे कि "आप किसी भी तरह से उस जटिल वितरण से स्वतंत्र नमूने लेने में कामयाब रहे जिसके बारे में आप जानना चाहते थे।" यह कुछ ऐसा है जैसे मैं सशर्त संभावना दे सकता हूं, फिर मैं सशर्त संभाव्यता के आधार पर एक स्मृतिहीन प्रक्रिया उत्पन्न कर सकता हूं। जब मैं प्रक्रिया को लंबे समय तक उत्पन्न करता हूं, तो संयुक्त संभावना परिवर्तित हो सकती है। और फिर मैं उत्पन्न अनुक्रम के अंत में संख्याओं का ढेर ले सकता हूं। यह वैसे ही है जैसे मैं जटिल संयुक्त वितरण से स्वतंत्र नमूने लेता हूं। उदाहरण के लिए, मैं हिस्टोग्राम बना सकता हूं और यह वितरण फ़ंक्शन को अनुमानित कर सकता है।
फिर मेरी समस्या यह है कि क्या मुझे यह साबित करने की जरूरत है कि क्या एमसीएमसी एक निश्चित मॉडल के लिए अभिसरण करती है? मैं यह जानने के लिए प्रेरित हूं क्योंकि मैंने पहले GMM और LDA (ग्राफिकल मॉडल) के लिए EM एल्गोरिथ्म सीखा था। अगर मैं सिर्फ यह साबित किए बिना एमसीएमसी एल्गोरिथ्म का उपयोग कर सकता हूं कि क्या यह परिवर्तित होता है, तो यह ईएम की तुलना में बहुत अधिक समय बचा सकता है। चूँकि मुझे अपेक्षित लॉग लाइबिलिटी फंक्शन की गणना करनी होगी (पोस्टीरियर संभावना की गणना करनी होगी), और उसके बाद अपेक्षित लॉग लाइबिलिटी को अधिकतम करें। यह स्पष्ट रूप से एमसीएमसी की तुलना में अधिक बोझिल है (मुझे बस सशर्त संभावना तैयार करने की आवश्यकता है)।
मैं यह भी सोच रहा हूं कि क्या संभावना समारोह और पूर्व वितरण संयुग्म हैं। क्या इसका मतलब यह है कि MCMC को अभिसरण करना चाहिए? मैं MCMC और EM की सीमाओं के बारे में सोच रहा हूँ।