হিট এন্ড এমসিএমসি চালান

Question

আমি হিট বাস্তবায়ন এবং এমসিসিএম অ্যালগরিদম চালানোর চেষ্টা করছি, তবে কীভাবে এটি করা যায় তা বুঝতে আমার কিছুটা সমস্যা হচ্ছে। সাধারণ ধারণাটি নিম্নরূপ:

এমএইচে একটি প্রস্তাব লাফ উত্পন্ন করতে, আমরা:

• ইউনিট গোলকের surface of পৃষ্ঠের বিতরণ থেকে একটি দিক তৈরি করুন$d$$\mathcal{O}$
• সীমাবদ্ধ স্থান সহ একটি স্বাক্ষরিত দূরত্ব উত্পন্ন করুন ।$\lambda$

তবে এটিকে আর (বা অন্য কোনও ভাষায়) বাস্তবায়নের বিষয়ে আমার কী ধারণা করা উচিত তা সম্পর্কে আমার কোনও ধারণা নেই।

কারও কি কোডের স্নিপেট রয়েছে যা আমাকে সঠিক দিকে নির্দেশ করবে?

বিটিডাব্লু, আমি যে লাইব্রেরিতে এই পদ্ধতিটি করি সে বিষয়ে আগ্রহী নই, আমি নিজে চেষ্টা করে কোড করতে চাই।

অনেক ধন্যবাদ.

Comments

এই অ্যালগরিদমের আগে কখনও শুনিনি, তবে এটি বেশ আকর্ষণীয় বলে মনে হচ্ছে। আপনি একটি ব্যাখ্যামূলক উত্স একটি লিঙ্ক প্রদান করতে পারেন? ধন্যবাদ !

— স্টেফেন

@ স্টেফেন এখানে মূল কাগজের একটি লিঙ্ক - ভাল আসলে প্রযুক্তিগত প্রতিবেদন। আমি কখনও এই পদ্ধতি ব্যবহার করি নি।

— csgillespie

অন্যান্য পদ্ধতিতে হিট অ্যান্ড এমসিএমসি চালানোর সুবিধা কী? অভিযানের গতি?

— রকসায়েন্স

@ ফ্রেড: এজন্য আমি অনুদানের প্রস্তাব দিয়েছি। আমি হিট অ্যান্ড রান সম্পর্কে পড়েছি কিন্তু কোন পরিস্থিতিতে এটি কার্যকর হবে তা জানি না।

— csgillespie

লিঙ্কটি এখন ভেঙে গেছে।

— ডাকোনেসস

Answer 1

আপনার সরবরাহিত কাগজের দিকে আমি নজর দিইনি, তবে যাই হোক আমাকে যেতে দিন:

আপনি যদি একটি থাকে তাহলে -dimensional প্যারামিটার স্থান একটি র্যান্ডম অভিমুখ তৈরি করতে পারেন অবিশেষে সঙ্গে ইউনিট গোলক পৃষ্ঠতলে বিতরণ$p$$d$

x <- rnorm(p)
d <- x/sqrt(sum(x^2))

(সিএফ উইকি )

তারপরে, প্রত্যাখ্যান স্যাম্পলিংয়ের জন্য প্রস্তাব উত্পন্ন করার জন্য এটি ব্যবহার করুন (ধরে নিবেন যে আপনি আসলে এর বিতরণটি মূল্যায়ন করতে পারেন )।$d$$d$

ধরে নেওয়া যাক আপনি অবস্থান শুরু এবং গ্রহণ করেছি , একটি প্রস্তাব উৎপন্ন সঙ্গে$x$$d$$y$

 lambda <- r<SOMEDISTRIBUTION>(foo, bar)
 y <- x + lambda * d

এবং তে স্থানান্তরিত হবে কি না তা সিদ্ধান্ত নিতে একটি মহানগর-হেস্টিংস-পদক্ষেপ করুন ।$y$

অবশ্যই, এটি কতটা ভাল কাজ করতে পারে তা এর বিতরণের উপর নির্ভর করবে এবং প্রত্যাখ্যানের নমুনা পদক্ষেপে এর ঘনত্ব মূল্যায়ন (বারবার) করা কত ব্যয়বহুল, তবে যেহেতু জন্য প্রস্তাবনা উত্পন্ন করা সস্তা, আপনি এটি থেকে দূরে সরে যেতে পারেন।$d$$d$

---

@ সিগিল্লেস্পির সুবিধার জন্য যুক্ত করা হয়েছে:

কিছু গুগলিংয়ের মাধ্যমে আমি যেটি সংগ্রহ করতে পেরেছিলাম, সেগুলি থেকে হিট-এন্ড-রান এমসিএমসি প্রাথমিকভাবে দ্রুত মিক্সিংয়ের জন্য দরকারী যদি আপনার একটি (মাল্টিভারিয়েট) লক্ষ্য থাকে যা নির্বিচারে আবদ্ধ থাকে তবে প্রয়োজনীয় সংযুক্ত সমর্থন নেই, কারণ এটি আপনাকে যে কোনও বিন্দু থেকে সরিয়ে নিতে সক্ষম করে তোলে এক পদক্ষেপে অন্য যে কোনও ব্যক্তির সমর্থনে। আরও এখানে এবং এখানে ।

Comments

যদি এই উত্তর সন্তোষজনক না হয়, আপনি কেন ব্যাখ্যা করতে পারেন?

— জন সালভাটিয়ার

আমি @ ফ্রেড সম্পর্কে নিশ্চিত নই, তবে আমি যখন অনুগ্রহের প্রস্তাব দিই তখন আমি মনে করি আমি হিট-এন্ড-রান এমসিসিএম সম্পর্কে আরও কিছুটা অন্তর্দৃষ্টি চেয়েছি। উদাহরণস্বরূপ, কোন ধরণের সমস্যার জন্য এটি সবচেয়ে উপযুক্ত হবে। অবশ্যই, যদি অন্য কোনও উত্তর না থাকে, তবে এই প্রশ্নটি অনুগ্রহ অর্জন করবে।

— csgillespie

@ সিএসগিলিসপি: আমি আপনার আগ্রহ আরও ভালভাবে যুক্ত করতে আমার উত্তর সম্পাদনা করেছি। এটি যেন বলা না হয় যে আমি অনুগ্রহের অযোগ্য ছিলাম। ;)

— ফেবিয়ানরা

লিঙ্কটির জন্য অনেক ধন্যবাদ। আমি অনুগ্রহটি রাখার একটি কারণ হ'ল আমার গুগল অনুসন্ধানগুলি পদ্ধতিটির কয়েকটি গাণিতিক আলোচনা সঞ্চার করেছিল, তবে ব্যবহারিক প্রয়োগগুলির ক্ষেত্রে সামান্য। অনুগ্রহ প্রদানের আগে যদি আমি আরও 48 ঘন্টা অপেক্ষা করি তবে দয়া করে এটিকে সামান্য হিসাবে নেবেন না (এটি একটি বিশেষ উদার অনুগ্রহ!)

— csgillespie

Answer 2

আমি যখন হিট-এন্ড-রানের মূল সূত্রটি খুঁজছিলাম তখন আমি আপনার প্রশ্নটি জুড়ে এসেছি। তার জন্য ধন্যবাদ! আমি সাম্প্রতিক এই ব্লগের শেষে পাইএমসির জন্য হিট-ও-রানের একটি প্রমাণ-ধারণা ধারণাকে একসাথে রেখেছি ।