First you should create a new binary variable, sex_partners. The 1-1 and 2-2 cases should receive one code (e.g. 1, "Y"), and the 3, 4, 5 should receive the alternate code (e.g. 0, "N"). This is your independent variable. The other variables are nesting variables. Then I would try: proc glimmix data = mydatset ASYCORR ASYCOV; class household pairs sex_partners; model disease = sex_partners /dist = bin link = logit oddsratio solution; random intercept / subject=pairs(household); run;
... View more