Displaying 1 result from an estimated 1 matches for "tab_psi".
Did you mean:
tab_pro
2012 Nov 12
1
R lmer & SAS glimmix
...have
simplified my code but I am surprised to see I get different results from
the two softwares.
My R code is :
lmer(y~age_cat + (1|cat),data=fic,family=binomial(link = "logit"), NaGQ=1)
My SAS code is :
ods output Glimmix.Glimmix.ParameterEstimates=t_estimates;
proc glimmix data=tab_psi method=laplace;
class age_cat cat;
model psi (event='1') = age_cat / solution dist=B link=logit ;
random intercept / subject=cat;
run;
>From R, I get the following fixed effects
(Intercept) age_cat2. 76-85 ans age_cat3. 66-75 ans age_cat4. 41-65 ans
-3.5766898 -0.01...