আমি এমএমসিসি.লিস্ট থেকে প্লট.বাগস এবং প্লট.জ্যাগগুলির অনুরূপ প্লট কীভাবে তৈরি করতে পারি? [বন্ধ]

Question

আর থেকে আউটপুট চমৎকার সারসংক্ষেপ প্লট পাবে বলে মনে হয় bugsএবং jagsফাংশন দ্বারা উত্পন্ন বস্তু R2WinBUGS :: বাগ এবং R2jags: Jags ।

তবে আমি rjagsপ্যাকেজটি ব্যবহার করছি । আমি যখন ফলাফলগুলি rjags::coda.samplesব্যবহার করে ফাংশনের ফলাফলগুলি প্লট করার চেষ্টা করি তখন R2WinBUGS::plot.mcmc.listপ্রতিটি প্যারামিটারের জন্য ডায়গনিস্টিক প্লট (প্যারামিটারের ঘনত্ব, চেইন টাইম সিরিজ, স্বতঃসংশোধন) হয়।

অ্যান্ড্রু গেলম্যানের টিউটোরিয়াল "রানিং উইনবুগস এবং ওপেনব্যাগস অফ আর থেকে" থেকে নীচে আমি যে ধরণের প্লট তৈরি করতে চাই তার নীচে দেওয়া হল । এগুলি ব্যবহার করে উত্পাদিত হয়েছিল plot.pugs।

সমস্যাটি হ'ল plot.bugsকোনও bugsবস্তুকে আর্গুমেন্ট হিসাবে plot.mcmc.listগ্রহণ করে , যখন আউটপুট নেয় coda.samples।

এখানে একটি উদাহরণ (থেকে coda.samples):

 library(rjags)
 data(LINE)
 LINE$recompile()
 LINE.out <- coda.samples(LINE, c("alpha","beta","sigma"), n.iter=1000)
 plot(LINE.out)

আমার যা দরকার তা হয়

• উত্পাদিতর মতো একই, তথ্য সমৃদ্ধ, এক পৃষ্ঠার সংক্ষিপ্ত প্লট তৈরি করার একটি উপায় plot.bugs
• একটি ফাংশন যা LINE.outবাগের বস্তুতে রূপান্তরিত করে বা

Answer 1

যেহেতু কোনও উত্তর নেই, তাই আমি এখন পর্যন্ত যা অর্জন করেছি তা অন্তত পোস্ট করব:

as.bugs.arrayফাংশন R2WinBUGSপ্যাকেজ এই কাজের জন্য তৈরি করা হয়েছে। ডকুমেন্টেশন অনুযায়ী ( ?as.bugs.array):

মার্কভ চেইন সিমুলেশনগুলি থেকে ফলাফলগুলিকে রূপান্তরিত করা ফাংশন, এটি বাগের বিষয়বস্তুতে BUGS থেকে নাও আসতে পারে। প্লট.বাগগুলি সহ ফলাফলগুলি প্রদর্শন করতে সাধারণত ব্যবহৃত হয়।

সুতরাং, LINE.outআপনার উদাহরণ থেকে একটি প্লট পাওয়া সম্ভব , যদিও এটি সঠিক ভেরিয়েবলগুলি প্লট করে না:

plot(as.bugs.array(sims.array = as.array(LINE.out)))        

রূপান্তর করার সঠিক উপায়টি নির্ধারণ করতে আরও কিছুটা সময় লাগবে LINE.outএবং LINE.samplesঅবজেক্টটি থেকে example(jags.samples)শুরু করার সহজ জায়গা হতে পারে।

Answer 2

নিম্নলিখিতটি আমার পক্ষে কাজ করে বলে মনে হচ্ছে:

require(R2jags)
m <-jags(data=d,inits=i,pars,n.iter=1000,n.chains=3,model.file="foo.txt",DIC=F)
m <- autojags(m)
plot(m)

এখানে একটি পুনরুত্পাদনযোগ্য উদাহরণ:

example(jags)
plot(jagsfit)

Comments

এটি একটি সহায়ক ক্লু, তবে একটি mcmc.list(যতদূর আমি বলতে পারি) দিয়ে শুরু করার সমস্যাটি সমাধান করে না ।

— ডেভিড লেবাউর