اقترانات من الماضي

من بين خوارزميات مونت كارلو لسلاسل ماركوف (MCMC) ، يُعدّ الاقتران من الماضي طريقةً لأخذ عينات من التوزيع الثابت لسلسلة ماركوف . وعلى عكس العديد من خوارزميات MCMC، يُعطي الاقتران من الماضي، من حيث المبدأ، عينةً مثاليةً من التوزيع الثابت . وقد ابتكره جيمس بروب وديفيد ويلسون عام 1996.

الفكرة الأساسية

لنفترض سلسلة ماركوف غير دورية غير قابلة للاختزال ذات حالات محدودةم{\displaystyle M}مع فضاء الحالةS{\displaystyle S}والتوزيع الثابت (الفريد)π{\displaystyle \pi }(π{\displaystyle \pi }(هو متجه احتمالي ). لنفترض أننا توصلنا إلى توزيع احتماليμ{\displaystyle \mu }على مجموعة الخرائطو:SS{\displaystyle f:S\to S}مع الخاصية التي لكل ثابتsS{\displaystyle s\in S}صورتهاو(s){\displaystyle f(s)}يتم توزيعها وفقًا لاحتمالية الانتقال لـم{\displaystyle M}من الولايةs{\displaystyle s}ومن أمثلة هذا التوزيع الاحتمالي التوزيع الذيو(s){\displaystyle f(s)}مستقل عنو(s){\displaystyle f(s')}حينماss{\displaystyle s\neq s'}لكن من المفيد غالبًا النظر في توزيعات أخرى. لنفترض الآنوج{\displaystyle f_{j}}لجZ{\displaystyle j\in \mathbb {Z} }كن عينات مستقلة منμ{\displaystyle \mu }.

لنفترض أنx{\displaystyle x}يتم اختيارها عشوائياً وفقاً لـπ{\displaystyle \pi }وهو مستقل عن التسلسلوج{\displaystyle f_{j}}(لا داعي للقلق في الوقت الحالي بشأن مكان هذاx{\displaystyle x}(يأتي من.) ثمو-1(x){\displaystyle f_{-1}(x)}يتم توزيعها أيضًا وفقًا لـπ{\displaystyle \pi }، لأنπ{\displaystyle \pi }يكونم{\displaystyle M}-ثابتة وافتراضنا بشأن قانونو{\displaystyle f}. يُعرِّف

Fج:=و-1و-2و-ج.{\displaystyle F_{j}:=f_{-1}\circ f_{-2}\circ \cdots \circ f_{-j}.}

ثم يترتب على ذلك بالاستقراء أنFج(x){\displaystyle F_{j}(x)}يتم توزيعها أيضًا وفقًا لـπ{\displaystyle \pi }لكلجشمال{\displaystyle j\in \mathbb {N} }ومع ذلك، قد يحدث ذلك بالنسبة للبعضنشمال{\displaystyle n\in \mathbb {N} }صورة الخريطةFن{\displaystyle F_{n}}هو عنصر واحد منS{\displaystyle S}. بعبارة أخرى،Fن(x)=Fن(y){\displaystyle F_{n}(x)=F_{n}(y)}لكلyS{\displaystyle y\in S}لذلك، لسنا بحاجة إلى الوصول إلىx{\displaystyle x}من أجل الحسابFن(x){\displaystyle F_{n}(x)}ثم تتضمن الخوارزمية إيجاد بعضنشمال{\displaystyle n\in \mathbb {N} }بحيثFن(S){\displaystyle F_{n}(S)}هو عنصر منفرد ، وإخراج عنصر هذا العنصر المنفرد. تصميم توزيع جيدμ{\displaystyle \mu }والتي تتمثل مهمة إيجاد مثل هذان{\displaystyle n}والحوسبةFن{\displaystyle F_{n}}ليس من الواضح دائمًا أن التكلفة ليست باهظة، ولكن تم تحقيق ذلك بنجاح في العديد من الحالات المهمة. [ 1 ]

الحالة الرتيبة

توجد فئة خاصة من سلاسل ماركوف تتميز بخيارات جيدة بشكل خاص لـμ{\displaystyle \mu }وأداة لتحديد ما إذا|Fن(S)|=1{\displaystyle |F_{n}(S)|=1}. (هنا||{\displaystyle |\cdot |}(يشير إلى العدد الأصلي .) لنفترض أنS{\displaystyle S}هي مجموعة مرتبة جزئياً بترتيب{\displaystyle \leq }، والذي يحتوي على عنصر بسيط فريدs0{\displaystyle s_{0}}وعنصر أقصى فريدs1{\displaystyle s_{1}}أي كلsS{\displaystyle s\in S}يرضيs0ss1{\displaystyle s_{0}\leq s\leq s_{1}}. افترض أيضاً أنμ{\displaystyle \mu }يمكن اختيارها لتكون مدعومة على مجموعة الخرائط الرتيبةو:SS{\displaystyle f:S\to S}ومن ثم يسهل أن نرى ذلك|Fن(S)|=1{\displaystyle |F_{n}(S)|=1}إذا وفقط إذاFن(s0)=Fن(s1){\displaystyle F_{n}(s_{0})=F_{n}(s_{1})}، منذFن{\displaystyle F_{n}}دالة رتيبة. لذا، يصبح التحقق من ذلك سهلاً للغاية. يمكن للخوارزمية أن تتابع باختيارن:=ن0{\displaystyle n:=n_{0}}لبعض الثوابتن0{\displaystyle n_{0}}أخذ عينات من الخرائطو-1،...،و-ن{\displaystyle f_{-1},\dots ,f_{-n}}، وإخراجFن(s0){\displaystyle F_{n}(s_{0})}لوFن(s0)=Fن(s1){\displaystyle F_{n}(s_{0})=F_{n}(s_{1})}. لوFن(s0)Fن(s1){\displaystyle F_{n}(s_{0})\neq F_{n}(s_{1})}تتم الخوارزمية عن طريق المضاعفةن{\displaystyle n}وتكرار ذلك حسب الحاجة حتى يتم الحصول على مخرجات. (لكن الخوارزمية لا تعيد أخذ عينات من الخرائط).و-ج{\displaystyle f_{-j}}والتي تم أخذ عينات منها بالفعل؛ ويستخدم الخرائط التي تم أخذ عينات منها مسبقًا عند الحاجة.)

مراجع

  • بروب، جيمس غاري؛ ويلسون، ديفيد بروس (1996)، وقائع المؤتمر الدولي السابع حول الهياكل والخوارزميات العشوائية (أتلانتا، جورجيا، 1995) ، الصفحات 223-252 ، MR 1611693  
  • بروب، جيمس؛ ويلسون، ديفيد (1998)، "الاقتران من الماضي: دليل المستخدم"، دراسات استقصائية مصغرة في الاحتمالات المنفصلة (برينستون، نيوجيرسي، 1997) ، سلسلة DIMACS في الرياضيات المنفصلة، ​​النظرية، وعلوم الحاسوب، المجلد  41، بروفيدنس، رود آيلاند: الجمعية الرياضية الأمريكية ، الصفحات 181-192 ، doi : 10.1090/dimacs/041/09 ، ISBN  9780821808276، MR 1630414 ، S2CID 2781385