06-22-2011 05:06 PM
06-23-2011 07:15 AM
06-23-2011 08:42 PM
Yes, a subject specification is required. If there are no subject blocks and all data should be regarded as a single subject, then specify subject=intercept.
06-24-2011 05:18 PM
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;
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;