Hi there, I had the similar problem, when I removed the ddfm=kr from the model line it gets all solved. but I don't understand why it would like that, does it affect the bias? here is my origional code. thanks, data proximate; input location$ rep GE ; response=GE;var='GE';output; cards; Kansas 1 4255.79 Kansas 2 4338.34 Kansas 3 4268.54 Kansas 4 4239.63 Kansas 5 4201.55 Kansas 6 4230.07 LA 1 3986.36 LA 2 4274.80 LA 3 4024.32 LA 4 4239.64 LA 5 4165.36 LA 6 4243.31 Dallas 1 4239.02 Dallas 2 4195.14 Dallas 3 4151.11 Dallas 4 4097.20 Dallas 5 4127.85 Dallas 6 4177.10 GA 1 4344.53 GA 2 4269.53 GA 3 4336.44 GA 4 4514.79 GA 5 4388.69 GA 6 4237.09 Stockton 1 4099.26 Stockton 2 4010.61 Stockton 3 3960.10 Stockton 4 4024.82 Stockton 5 4238.10 Stockton 6 4104.96 Baltimore 1 4086.10 Baltimore 2 4168.86 Baltimore 3 4064.65 Baltimore 4 4460.57 Baltimore 5 4151.72 Baltimore 6 3982.83 ; RUN; proc print data=proximate; run; proc sort data=proximate; by var; run; proc glimmix data=proximate plots=RESIDUALPANEL; by var; class location rep; model response = location/ddfm=kr; random rep(location); lsmeans location / pdiff lines adjust=tukey; run;
... View more