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.
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;
Don't miss out on SAS Innovate - Register now for the FREE Livestream!
Can't make it to Vegas? No problem! Watch our general sessions LIVE or on-demand starting April 17th. Hear from SAS execs, best-selling author Adam Grant, Hot Ones host Sean Evans, top tech journalist Kara Swisher, AI expert Cassie Kozyrkov, and the mind-blowing dance crew iLuminate! Plus, get access to over 20 breakout sessions.
ANOVA, or Analysis Of Variance, is used to compare the averages or means of two or more populations to better understand how they differ. Watch this tutorial for more.
Find more tutorials on the SAS Users YouTube channel.