Hello!
I am using glmmML for a logitic regression with random effect.
I use the posterior.mode as an estimate for the random effects.
These can be very different from the estimates obtained using SAS , NLMIXED
in the random with out= option. (all the fixed and standard error of random
effect estimators are almost identical)
Can someone explain to me why is that.

The codes I use:
R:
glmm1<-glmmML(mort30 ~ x , data=dat2,cluster=hospital,family=binomial)
print(sort(glmm1$posterior.mode))

SAS:
*

proc* *nlmixed* data*=*dat*;*

eta = b0 + b1*x+ u;

expeta = exp(eta);

p = expeta/(*1*+expeta);

model mort30 ~ binomial(*1*,p);

random u ~ normal(*0*,s2) subject=hospital out=blue;

*run*;
*

proc* *sort* data=blue;by estimate;
*

proc* *print* data=blue;*run*;

**

*THANKS FOR THE HELP*

*RON*

	[[alternative HTML version deleted]]

