জেজিএসে সেন্সরিং / কাটাকাটি


10

জেএজেএসে সেন্সরিংয়ের সমস্যাটি কীভাবে ফিট করতে হয় সে সম্পর্কে আমার একটি প্রশ্ন রয়েছে।

আমি একটি দ্বিবিড়ীয় মিশ্রণটি সাধারণভাবে পর্যবেক্ষণ করি যেখানে এক্স মানগুলিতে পরিমাপের ত্রুটি রয়েছে। আমি পর্যবেক্ষণ করা সেন্সরকৃত মানগুলির সত্যিকারের অন্তর্নিহিত 'উপায়গুলি' মডেল করতে চাই।

xtrue+ϵ=xobserved ϵN(0,sd=.5)

আমার এখন যা আছে তা এখানে:

 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)

স্পষ্টতই সি কমান্ডটি জেজিএসে বৈধ নয়।

আগাম ধন্যবাদ.


3
ছাঁটাই করতে টি (-, -) ব্যবহার করুন, তবে
সেন্সারিং

উত্তর:


9

আপনি এটি যা খুঁজছেন সম্ভবত:

x_obs[i] ~ dnorm(x_true[i],prec_x)T(x_true[i], )

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

জাগগুলি কীভাবে ছাঁটাই এবং সেন্সর ব্যবহার করে সে সম্পর্কে আরও তথ্যের জন্য ব্যবহারকারীর ম্যানুয়ালটি পড়ুন ।


3

টিপসের জন্য ধন্যবাদ ডেভিড। আমি জাগস সমর্থন ফোরামে এই প্রশ্নটি পোস্ট করেছি এবং একটি দরকারী উত্তর পেয়েছি। মূলটি ছিল 'সত্য' মানগুলির জন্য একটি দ্বিমাত্রিক অ্যারে ব্যবহার করা to

for (j in 1:n){ 
  x_obs[j] ~ dnorm(xy_true[j,1], prec_x)T(xy_true[j,1],) 
  y_obs[j] ~ dnorm(xy_true[j,2], prec_y)
  xy_true[j, ] ~ dmnorm(mu[ z [j],1:2], tau[z[j],1:2,1:2]) 
  z[j]~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) 
আমাদের সাইট ব্যবহার করে, আপনি স্বীকার করেছেন যে আপনি আমাদের কুকি নীতি এবং গোপনীয়তা নীতিটি পড়েছেন এবং বুঝতে পেরেছেন ।
Licensed under cc by-sa 3.0 with attribution required.