জেএজেএসে সেন্সরিংয়ের সমস্যাটি কীভাবে ফিট করতে হয় সে সম্পর্কে আমার একটি প্রশ্ন রয়েছে।
আমি একটি দ্বিবিড়ীয় মিশ্রণটি সাধারণভাবে পর্যবেক্ষণ করি যেখানে এক্স মানগুলিতে পরিমাপের ত্রুটি রয়েছে। আমি পর্যবেক্ষণ করা সেন্সরকৃত মানগুলির সত্যিকারের অন্তর্নিহিত 'উপায়গুলি' মডেল করতে চাই।
আমার এখন যা আছে তা এখানে:
for (i in 1:n){
x[i,1:2]~dmnorm(mu[z[i],1:2], tau[z[i],1:2,1:2])
z[i]~dcat(prob[ ])
}
ওয়াইতে পরিমাপের ত্রুটিও রয়েছে। আমি যা করতে চাই তা হ'ল:
for (i in 1:n){
x_obs[i] ~ dnorm(x_true[i],prec_x)I(x_true[i],)
y_obs[i] ~ dnorm(y_true[i],prec_y)
c(x_true[i]:y_true[i])~dmnorm(mu[ z [ i ],1:2], tau[z[i],1:2,1:2])
z[i]~dcat(prob[ ])
}
#priors for measurement error
e_x~dunif(.1,.9)
prec_x<-1/pow(e_x,2)
e_y~dunif(2,4)
prec_y<-1/pow(e_y,2)
স্পষ্টতই সি কমান্ডটি জেজিএসে বৈধ নয়।
আগাম ধন্যবাদ.