Displaying 1 result from an estimated 1 matches for "currphi2".
Did you mean:
currphi1
2009 Mar 28
1
Error in R??
...code is part of a bigger program.
##sample from full conditional distribution of Si
#Prob(Si = 1)
for (j in 1:N)
{
numerat =
currphi1*exp((-1/(2*currvar))*t(Y1-(t(XX1)%*%currbeta1))%*%(Y1-(t(XX1)%*%currbeta1)))
denomin =
currphi2*exp((-1/(2*currvar))*t(Y2-(t(XX2)%*%currbeta2))%*%(Y2-(t(XX2)%*%currbeta2)))
sum=denomin + numerat
ProbSi = numerat/sum
arunofSi[j]=rbinom(1,1,ProbSi) #Generate 50 Bernoulli rvs and
assign them to arunofSi array
}
N0=sum(arunofSi==0)...