Displaying 1 result from an estimated 1 matches for "fixpart".
Did you mean:
fixpar
2011 Dec 27
1
rbinom
...t 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