Displaying 1 result from an estimated 1 matches for "currbeta1".
2009 Mar 28
1
Error in R??
...50 or more warnings (use warnings() to see the first
50)
The 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 aru...