Displaying 1 result from an estimated 1 matches for "sgnbeta".
Did you mean:
pnbeta
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 user specified input
that relates the?data...