কিভাবে rjags দিয়ে পূর্বাভাস উত্পন্ন?

Question

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

আমি বিশ্বাস করি http://sourceforge.net/p/mcmc-jags/discussion/610037/thread/0ecab41c একই প্রশ্ন জিজ্ঞাসা করছে।

Comments

লিঙ্কটি নষ্ট হয়ে গেছে। আপনি এটি আপডেট করতে পারেন, দয়া করে?

— chl

সম্পন্ন. দেখে মনে হচ্ছে সোর্সফোজের পুনরায় সাজানো

— কোয়ানটিটিভেটিভ

ধন্যবাদ! (এসএফ তাদের আপেক্ষিক পথ পরিবর্তন করার পরে আমি নিজেই থ্রেডটি পাইনি))

— chl

Answer 1

সাধারণত আপনি জ্যাজিএস-এ ভবিষ্যদ্বাণী করতে পারেন। নীচে নির্ভরশীল পরিবর্তনশীল এবং বয়স এবং ভবিষ্যদ্বাণীকারী হিসাবে ধূমপান সূচক হিসাবে FEV (ফুসফুসের ক্ষমতার সাথে কিছু করার জন্য) সহ একটি রিগ্রেশন উদাহরণ রয়েছে।

FEV20s এবং FEV20ns হ'ল 20 বছর বয়সী ধূমপায়ী এবং 20 বছর বয়সী ধূমপায়ী নানার জন্য পূর্বাভাস হওয়া FEV মান।

model
{
for(i in 1:n){
    FEV[i] ~ dnorm(mu[i],tau)
    mu[i] <- beta[1] + beta[2]*Age[i] + beta[3]*Smoke[i]  + beta[4]*Age[i]*Smoke[i]
}

#priors
beta[1] ~ dnorm(0,0.001)
beta[2] ~ dnorm(0,0.001)
beta[3] ~ dnorm(0,0.001)
beta[4] ~ dnorm(0,0.001)
tau ~ dgamma(0.001,0.001)
sigma<-1/sqrt(tau) 

## Predict the FEV for a 20 year old smoker and for a 20 year old nonsmoker
mu20s <-  beta[1] + (beta[2]+beta[4])*20 + beta[3]
mu20ns <-  beta[1] + beta[2]*20 
FEV20s ~ dnorm(mu20s,tau)
FEV20ns ~ dnorm(mu20ns,tau)
}

এর উদাহরণ: বায়েশিয়ান আইডিয়া এবং ডেটা বিশ্লেষণ

Comments

পয়েন্টারের জন্য ধন্যবাদ - আমি কেবল আমার পরীক্ষার ডেটা জেএজেজে প্রেরণের বিষয়ে চিন্তা করিনি, তবে এটি করা উচিত।

— পরিমাণগত

পুরো মডেলটি রিফিট না করে এই পূর্বাভাসগুলি উত্পন্ন করার কোনও উপায় আছে কি? যদি ভবিষ্যতটি তৈরির পূর্বাভাসকে ব্যাপকভাবে সমান্তরাল করা যথেষ্ট সহজ হত তবে, যদি পুরো মডেলটিকে পুনর্বিবেচনার প্রয়োজন হয় তবে এটি সম্ভব নয়।

— কলিন