I am not real sure, but would some sort of hierarchical model, such as GLIMMIX work for what you want? proc glimmix data=have; class id1 id2; model depvar=indvar; random intercept/subject=id1; random _residual_/subject=id2(id1) type=vc; run; Steve Denham
... View more