Thanks for your suggestion! When I try the code below, it works. However, the standard error of SP(SPH) was showed "." in the covariance parameter estimates. The odds ratios do not change no matter I change the type of covariance to vc or cs. Is it correct? The ID is the unique ID for case and control and the "x_cor" and "y_cor" are the coordinate information. The geographic distribution of case and control points are very similar. PROC GLIMMIX data=BROWN; CLASS ID; MODEL DISEASE= ode200m grass200m EM1/ dist=binary solution oddsratio; random _residual_ /subject=ID type=sp(sph) (x_cor y_cor); output out=abd pred=p resid=r; RUN;
... View more