topic Proc GLIMMIX - trying to program 3 random intercepts in Statistical Procedures

Proc GLIMMIX - trying to program 3 random intercepts

KPCklebspn — Tue, 04 Dec 2018 15:50:46 GMT

Hi all,

I want 3 random intercepts (hospital, unit, room). When I program the intercept only model, I get 0 as the room level variance estimate which seems wrong. Did I program it correctly?

proc glimmix data = master;
class Hospital UniqueUnit UniqueRm;
model Result1 = / dist=bin solution cl link=logit oddsratio;
random int / subject=hospital type=un s;
random int / subject=UniqueUnit(hospital) type=un s;
random int / subject=UniqueRm (UniqueUnit hospital) type=un s;
run;

Thanks!

Re: Proc GLIMMIX - trying to program 3 random intercepts

KPCklebspn — Tue, 04 Dec 2018 16:05:28 GMT

To clarify the data has 4 levels: patient, room, unit, hospital

Re: Proc GLIMMIX - trying to program 3 random intercepts

PGStats — Tue, 04 Dec 2018 20:32:44 GMT

You would need more than one patient per room to estimate the room level variance.