I don't run into problems doing this:
x=rnorm(300)
z=rnorm(300)
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])); T[i]<-rbinom(1,1,p[i])
}
hth, Ingmar
On 23 Nov 2007, at 10:25, sigalit mangut-leiba wrote:
> 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]]
>
> ______________________________________________
> R-help@r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-
> guide.html
> and provide commented, minimal, self-contained, reproducible code.
Ingmar Visser
Department of Psychology, University of Amsterdam
Roetersstraat 15
1018 WB Amsterdam
The Netherlands
t: +31-20-5256723
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
