Bayesian ऑनलाइन बदलाव का पता लगाने (सीमांत भविष्य कहनेवाला वितरण)


9

मैं एडम्स और मैकके ( लिंक ) द्वारा बायेसियन ऑनलाइन चेंजप्वाइंट डिटेक्शन पेपर पढ़ रहा हूं ।

लेखक सीमांत भविष्य कहनेवाला वितरण लिखकर शुरू करते हैं: जहां

P(xt+1|x1:t)=rtP(xt+1|rt,xt(r))P(rt|x1:t)(1)
  • xt समय पर अवलोकन ;t
  • x1:t अवलोकन के सेट को दर्शाता है जब तक कि समय ;t
  • rtN वर्तमान रननलोफ्ट (अंतिम बदलाव के बाद का समय, 0 हो सकता है); तथा
  • xt(r) रन जुड़ी टिप्पणियों का ।rt

Eq। 1 औपचारिक रूप से सही है (@JuhoKokkala द्वारा नीचे दिए गए उत्तर देखें), लेकिन मेरी समझ यह है कि यदि आप वास्तव में बारे में एक भविष्यवाणी करना चाहते हैं, तो आपको इसे निम्नानुसार विस्तारित करना होगा:xt+1

P(xt+1|x1:t)=rt,rt+1P(xt+1|rt+1,xt(r))P(rt|x1:t)P(rt+1|rt)(1b)

मेरा तर्क यह है कि भविष्य में ( + ) समय t + 1 में एक बदलाव हो सकता है t+1, लेकिन पीछे P(rt|x1:t) केवल t तक शामिल हैt

मुद्दा यह है, कागज में लेखक हमें ईक बनाते हैं। 1 जैसा है (Eqs देखें। 3 और 11 पेपर में), और 1 बी नहीं । तो, वे उचित रूप में समय में एक changepoint की संभावना की उपेक्षा t+1 जब भविष्यवाणी xt+1 समय में उपलब्ध आंकड़ों से t । धारा 2 की शुरुआत में वे कहते हैं कि एन पास

हम मानते हैं कि हम एक दिए गए रन लंबाई r_t पर [वितरण के लिए [ xt+1 ] के लिए भविष्य कहनेवाला गणना कर सकते हैं ।rt

जो शायद जहां चाल है। लेकिन सामान्य तौर पर, यह भविष्य कहनेवाला वितरण ईक की तरह कुछ दिखना चाहिए। 1b; वह नहीं है जो वे करते हैं (Eq 11)।

इसलिए, मुझे यकीन नहीं है कि मैं समझता हूं कि क्या हो रहा है। शायद अंकन के साथ कुछ मजेदार चल रहा है।


संदर्भ

  • एडम्स, आरपी, और मैकके, डीजे (2007)। बायेसियन ऑनलाइन चेंजप्वाइंट डिटेक्शन। arXiv प्रीप्रिंट arXiv: 0710.3742।

एक संभावित व्याख्या यह है कि टाइम के अंत में रन लंबाई का प्रतिनिधित्व करता , जो कि समय में बदलाव के बाद । इसके साथ ही, Eq। 1 समझ में आता है। वास्तव में, एल्गोरिथ्म का एक आरंभीकरण सेट करके है जो मानता है कि पर शुरू होने से ठीक पहले एक बदलाव है । हालाँकि, चित्र 1 गलत है (या कम से कम भ्रामक) है कि यदि और बीच कोई परिवर्तन बिंदु है , और और बीच चित्र 1a में दर्शाया गया है, तो औरrtttP(r0=0)=1t=1t=4t=5t=10t=11r4r10इस संकेतन के अनुसार 0 होना चाहिए, और छवि 1 बी के अनुसार और नहीं होना चाहिए । r5r11
लैकरबी

1
इक में कुछ अजीब सा चल रहा है। 3 अंतिम पंक्ति में सारांश में मध्य कारक के रूप में जबकि मैंने सोचा कि में शामिल है । मुझे संदेह है कि और ने रूप में स्थान बदल दिया है। Eq में। 11, दाहिने हाथ की ओर निर्भर करता है जो बाएं हाथ की ओर बिल्कुल भी दिखाई नहीं देता है, इसलिए या तो कुछ गड़बड़ है या मैं बिल्कुल भी संकेतन नहीं समझता। P(xtrt1,xt(r))xt(r)xttt1P(xtrt,xt1(r))xt(r)
जुहो कोक्कला

@JuhoKokkala: मुझे खुशी है कि मैं उस भावना के साथ अकेला नहीं हूं ...
लैकरबी

1
@lacerbi, मेरे पास इस पेपर के बारे में एक और सवाल है, और आपको लगता है कि आप इसका जवाब देने में सक्षम हो सकते हैं क्योंकि आप काम से परिचित हैं: आंकड़े . stackexchange.com/questions/419988
gwg

जवाबों:


5

दोनों (1) और (1 बी) सही हैं। ओपी के पास यह अधिकार है कि (इस मॉडल में) पर एक बदलाव हो सकता है , और बात पर निर्भर करता है कि क्या कोई बदलाव है। यह (1) के संभावित मानों के साथ कोई समस्या नहीं है, द्वारा पूरी तरह से "कवर" किया गया है । अर्थ है सशर्त वितरण । यह सशर्त वितरण औसत "सब कुछ" पर है, जिसमें , सशर्त पर । जैसे कोई लिख सकता है, कह सकता है,t+1xt+1rt+1P(xt+1rt,x1:t)P(xt+1|rt,x1:t)xt+1(rt,x1:t)rt+1(rt,x1:t)P(xt+1000|xt), जो बदलावों के सभी संभावित विन्यासों के साथ-साथ और बीच होने वाले मूल्यों को ध्यान में रखेगा ।xitt+1000

शेष में, मैं पहले (1) और फिर (1 बी) पर आधारित (1) प्राप्त करता हूं।

(1) की व्युत्पत्ति

किसी भी यादृच्छिक चर , हमारे पास जब तक असतत है तब तक (अन्यथा योग को एक अभिन्न द्वारा प्रतिस्थापित किया जाना चाहिए)। इसे : पर लागू किया जा रहा हैA,B,C

P(AB)=cP(AB,C=c)P(C=cB),
Cxt+1,x1:t,rt

P(xt+1x1:t)=rtP(xt+1rt,x1:t)P(rtx1:t),
जो कोई बात नहीं रखता है जो , , बीच निर्भरताएं हैं, यानी कोई मॉडल मान्यताओं को अभी तक नहीं किया गया है इस्तेमाल किया गया। वर्तमान मॉडल में, दिया माना जाता है * के मूल्यों की सशर्त स्वतंत्र होने के लिए रन से पहले से । इसका तात्पर्य । इसे पिछले समीकरण में शामिल करते हुए, हम प्राप्त करते हैंrtx1:txt+1xt+1rt,xt(r)xxt(r)P(xt+1rt,x1:t)=P(xt+1rt,xt(r))

P(xt+1x1:t)=rtP(xt+1rt,xt(r))P(rtx1:t),(1)
जो ओपी में (1) है।

(1 बी) की व्युत्पत्ति

हम में से अपघटन पर विचार करें के संभावित मूल्यों पर : P(xt+1rt,xt(r))rt+1

P(xt+1rt,xt(r))=rt+1P(xt+1rt+1,rt,xt(r))P(rt+1rt,xt(r)).

चूँकि यह माना जाता है कि क्या कोई परिवर्तन बिंदु पर होता है ( और ) के इतिहास पर निर्भर नहीं होता है , हमारे पास । इसके अलावा, चूंकि यह निर्धारित करता है कि के समान रन में है , हमारे पास । इन दो सरलीकरणों को उपर्युक्त कारक में प्रतिस्थापित करते हुए, हम t+1xtxt+1xP(rt+1rt,xt(r))=P(rt+1rt)rt+1xt+1xtP(xt+1rt+1,rt,xt(r))=P(xt+1rt+1,xt(r))

P(xt+1rt,xt(r))=rt+1P(xt+1rt+1,xt(r))P(rt+1rt).
इसे (1) में प्रतिस्थापित करते हुए, हम जो OP (1b) है।
P(xt+1x1:t)=rt(rt+1P(xt+1rt+1,xt(r))P(rt+1rt))P(rtx1:t),(1b)

* मॉडल की सशर्त स्वतंत्रता मान्यताओं पर टिप्पणी

कागज को जल्दी से ब्राउज़ करने के आधार पर, मैं व्यक्तिगत रूप से सशर्त स्वतंत्रता संपत्तियों को अधिक स्पष्ट रूप से कहीं और कहा जाना पसंद करूंगा, लेकिन मुझे लगता है कि इरादा यह है कि मार्कोवियन है और : s अलग-अलग रन से जुड़ा हुआ है स्वतंत्र (रन दिए गए) हैं।rx


1
(+1) धन्यवाद। हां, निश्चित रूप से, मैं समझता हूं कि Eq। 1 औपचारिक रूप से सही है यदि कोई मानता है कि पर निहित हाशिए पर है । समस्या यह है कि बाद में लेखकों ने भविष्यवाणियाँ कीं (Eq। 11 कागज में, और स्पष्ट रूप से Eq 3 में) और वे जब वे ले रहे हैं तो से अधिक हाशिए पर नहीं रहे हैं। rt+1rt+1
लैकरबी

1
ओह। तब ऐसा लगता है कि मैंने इस प्रश्न को गलत समझा - क्या मुझे इसे हटाना चाहिए? आप प्रश्न को स्पष्ट करना चाह सकते हैं, वर्तमान में ऐसा लगता है (1) किसी तरह गलत है (शायद उपयोगी नहीं है)
जुहो कोक्कला

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