Proc GLIMMIX - trying to program 3 random intercepts

[KPCklebspn]

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!

[KPCklebspn]

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

[PGStats]

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