आर में मिक्सचर बायेसियन रिग्रेशन का गिब्स नमूना कैसे बनाया जाए?

Question

मैं तीन मापदंडों के गिब्स नमूने पर काम कर रहा हूं और हम तीन मापदंडों के पूर्ण सशर्त वितरण को जानते हैं।

Answer 1

नोट: यह सिर्फ एक टिप्पणी है लेकिन एक लंबे कोड को उद्धृत करने के लिए, मैंने इसे यहां रखा है।

for (ite in 2:NSim){
#Full conditional for pi

 pi[ite]=rbeta(1, sum(delta[ite-1,])+0.5, sum(1-delta[ite-1,])+0.5)

 #Full conditional for delta
 for(j in 1:4){
  p1=pi[ite]*exp(-beta[ite-1,j]^2/(20))
  p0=((1-pi[ite])*10^3)*exp(-500*beta[ite-1,j]^2)
  cat('\n',ite,j,(p1/(p0+p1)))
  delta[ite,j]=rbinom(1, 1,prob=(p1/(p0+p1)))
  }

त्रुटि कहती है कि हो सकता है $p1$ तथा $p0$NA हैं। इस प्रकार की त्रुटि की जाँच करते समय मेरा अनुभव यह है कि लूप के बजाय, बस दें$ite=2$ तथा $j=1$। अपने सूत्र के रूप में p1 और p0 की गणना करें। ध्यान से देखें कि वे NA हैं या नहीं। यदि कोड सममित हैं$iter$ तथा $j$, और यदि आप इस मामले के लिए त्रुटि को ठीक कर सकते हैं, तो यह इसे पारित करेगा।