[R] rbinom

Scott Raynaud scott.raynaud at yahoo.com
Tue Dec 27 15:47:45 CET 2011


I have the following code (which I did not write) that generates 
data based on a logistic model.  I'm only getting a single record 
with y=1.  It seems implausible that in 50k cases that have a 
single y=1.  Does that ring alarm bells for anyone else?
 
beta<-c(-1.585600,-0.246900)
betasize<-length(beta)
meanpred<-c(0,35.900000)
varpred<-c(0,1.000000)
#loop code
x<-matrix(1,length,betasize) #length set to 50k
#loop code
  x[,2]<-rnorm(length,meanpred[2],sqrt(varpred[2])) #length set to 50k
   fixpart<-x%*%beta
    binomprob<-exp(fixpart)/(1+exp(fixpart))
     data$y<-rbinom(n1,1,binomprob)
#more loop code



More information about the R-help mailing list