Displaying 2 results from an estimated 2 matches for "betasize".
Did you mean:
metasize
2011 Dec 27
1
rbinom
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...
2012 Feb 01
1
Simulation confidence interval
The follwing is a code snippet from a power simulation
program that I'm using:
?
estbeta<-fixef(fitmodel)
?sdebeta<-sqrt(diag(vcov(fitmodel)))
? for(l in 1:betasize)
? {?
?? cibeta<-estbeta[l]-sgnbeta[l]*z1score*sdebeta[l]
??? if(beta[l]*cibeta>0)????????????? powaprox[[l]]<-powaprox[[l]]+1
????? sdepower[l,iter]<-as.numeric(sdebeta[l])
? }
?
Estbeta recovers the fixed effects from a model fitted using lmer.
Beta is defined elsewhere and is a use...