Displaying 1 result from an estimated 1 matches for "bernulli".
Did you mean:
bernoulli
2007 Nov 23
2
rbinom with computed probability
Hello,
I have a loop with probability computed from a logistic model like this:
for (i in 1:300){
p[i]<-exp(-0.834+0.002*x[i]+0.023*z[i])/(1+exp(-0.834+0.002*x[i]+0.023
+z[i]))
x and z generated from normal distribution.
I get 300 different probabilities And I want to generate variables from
bernulli distribution
with P for every observation:
T[i]<-rbinom(1,1,p[i])
But i get missing values for T.
What I'm doing wrong?
Thank you,
Sigalit.
[[alternative HTML version deleted]]