आप दो गाऊसी प्रक्रियाओं की तुलना कैसे करते हैं?


14

कुल्बैक-लिबलर विचलन दो संभावना घनत्व कार्यों की तुलना करने के लिए एक मीट्रिक है, लेकिन दो जीपी के और की तुलना करने के लिए किस मीट्रिक का उपयोग किया जाता है ?वाईXY


d(X,Y)=E[supt|X(t)Y(t)|]
ज़ेन

@ ज़ेन: यदि आपके पास समय है, तो मुझे इस दूरी मीट्रिक के बारे में अधिक जानने की दिलचस्पी है।
नील जी

हाय, नील। मैं इसके बारे में ज्यादा नहीं जानता। कृपया, मेरे उत्तर को देखें।
ज़ेन

जवाबों:


8

यह याद रखें कि गौसियन प्रक्रियाओं का वितरण संभवतया अनंत लिए मल्टीवेरेट गौसियन का विस्तार है । इस प्रकार, आप जीपी प्रायिकता वितरणों के बीच केएल विचलन का उपयोग पर एकीकृत करके कर सकते हैं :एक्स एक्स आर एक्सXRXRX

DKL(P|Q)=RXlogdPdQdP.

आप अपने जीपी वितरण के अनुसार बार-बार नमूने प्रक्रियाओं द्वारा एक विखंडित पर संख्यात्मक रूप से इस मात्रा को अनुमानित करने के लिए MC विधियों का उपयोग कर सकते हैं । मुझे नहीं पता कि क्या अभिसरण गति पर्याप्त रूप से अच्छी है ...X

इस बात पर ध्यान दें कि यदि साथ परिमित है , तो यदि आप बहुभिन्नरूपी सामान्य वितरण के लिए केएल विचलन पर वापस आते हैं: | एक्स | = एन डी के एल ( जी पी ( μ 1 , के 1 ) , जी पी ( μ 2 , के 2 ) ) = 1X|X|=n

DKL(GP(μ1,K1),GP(μ2,K2))=12(tr(K21K1)+(μ2μ1)K21(μ2μ1)n+log|K2||K1|)

मैं आपके द्वारा उल्लिखित दो साधनों (mu1 और mu2) की गणना कैसे कर सकता हूं। या मुझे उन्हें गॉसियन प्रक्रिया के लिए हमेशा की तरह शून्य के बराबर लेना चाहिए?
मराट जकीरोव

4

याद रखें कि अगर एक गॉसियन प्रक्रिया है जिसका मतलब है फ़ंक्शन और सहसंयोजक फ़ंक्शन , तो, हर , यादृच्छिक वेक्टर का अर्थ वेक्टर और सहसंयोजक मैट्रिक्स साथ एक बहुभिन्नरूपी सामान्य वितरण है। , जहां हमने सामान्य संक्षिप्त नाम । एम कश्मीर टी 1 , ... , टी कश्मीरटी ( एक्स ( टी 1 ) , ... , एक्स ( टी कश्मीर ) ) ( मीटर ( टी 1 ) , ... , मीटर ( टी कश्मीर ) ) Σ = ( σ मैं जे ) = ( कश्मीर ( टी मैंX:T×ΩRmKt1,,tkT(X(t1),,X(tk))(m(t1),,m(tk))X ( t ) = X ( t ,Σ=(σij)=(K(ti,tj))X(t)=X(t,)

प्रत्येक बोधन एक वास्तविक कार्य है जिसका डोमेन सूचकांक । मान लीजिए कि । दो को देखते हुए गाऊसी प्रक्रियाओं और , दो प्रतीति के बीच एक आम दूरी और है। इसलिए, यह दो प्रक्रियाओं के बीच की दूरी को परिभाषित करने के प्राकृतिक लगता है और के रूप में टी टी = [ 0 , 1 ] एक्स वाई एक्स (X(,ω)TT=[0,1]XYवाई (X(,ω)Y(,ω)supt[0,1]|X(t,ω)Y(t,ω)|XY

d(X,Y)=E[supt[0,1]|X(t)Y(t)|].()
मुझे नहीं पता कि इस दूरी के लिए एक विश्लेषणात्मक अभिव्यक्ति है, लेकिन मेरा मानना ​​है कि आप एक मोंटे कार्लो सन्निकटन की गणना निम्नानुसार कर सकते हैं। कुछ ठीक ग्रिड को ठीक करें , और नमूने आकर्षित और से सामान्य यादृच्छिक वैक्टर और , क्रमशः, । लगभग by 0t1<<tk1(xi1,,xik)(yi1,,yik)(X(t1),,X(tk))(Y(t1),,Y(tk))i=1,,Nd(X,Y)
1Ni=1Nmax1jk|xijyij|.

आप प्रत्येक वेक्टर से नमूना कैसे लेते हैं? यदि आप प्रत्येक जीपी में केवल उन साधनों का नमूना लेते हैं जिन्हें आप भिन्नताओं को ध्यान में नहीं रखते हैं। अन्यथा आपको एक नमूना तकनीक विकसित करनी होगी जो सुसंगत हो।
पुष्कर

यह एक उत्कृष्ट संसाधन है: gaussianprocess.org/gpml/chapters
Zen

आप इस प्रश्न के सभी उत्तर भी पढ़ सकते हैं: आंकड़े.stackexchange.com/questions/30652/…
Zen

वेतन ध्यान है कि इस के बाद से एक दूरी नहीं है । केएल दो वितरण और नहीं दो प्रतीति तुलना के रूप में, दो ग्राम पंचायतों के बीच जेन की दूरी के रूप में परिभाषित किया जाना चाहिए ( जी 1 , जी 2 ) = एक्स ~ जी 1 , वाई ~ जी 2 [ sup टी | X ( t ) - Y ( t ) | ] , और हम है कि एक्स ~d(X,X)0d(G1,G2)=EXG1,YG2[supt|X(t)Y(t)|]गैर के लिए गाऊसी प्रक्रिया विकृतजीEXG,YGsupt|X(t)Y(t)|>0G
एमिल

@Emile: यह कैसे है कि का उपयोग कर परिभाषा ( * ) ? d(X,X)0()
ज़ेन
हमारी साइट का प्रयोग करके, आप स्वीकार करते हैं कि आपने हमारी Cookie Policy और निजता नीति को पढ़ और समझा लिया है।
Licensed under cc by-sa 3.0 with attribution required.