turn on suggestions

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

Find a Community

- Home
- /
- Analytics
- /
- Stat Procs
- /
- How to adjust for spatial autocorrelation in the l...

Topic Options

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

06-22-2011 05:06 PM

Hi!

I conducted a case control study about the disease incidence and land cover types. The contro points were random sampled based on the population density.

How to adjust for the spatial autocorrelation in thr logistic model?

I have tried to do it by using GLIMMIX

PROC GLIMMIX data=test;

MODEL disease= grass200m

/ dist=binary solution oddsratio;

random _residual_ / type=sp(exp) (x_cor y_cor);

run;

It does not work. There is a warning message: "Obtaining minimum variance quadratic unbiased estimates as starting values for the covariance parameters failed."

How do I modify my code? Thanks a lot!

I conducted a case control study about the disease incidence and land cover types. The contro points were random sampled based on the population density.

How to adjust for the spatial autocorrelation in thr logistic model?

I have tried to do it by using GLIMMIX

PROC GLIMMIX data=test;

MODEL disease= grass200m

/ dist=binary solution oddsratio;

random _residual_ / type=sp(exp) (x_cor y_cor);

run;

It does not work. There is a warning message: "Obtaining minimum variance quadratic unbiased estimates as starting values for the covariance parameters failed."

How do I modify my code? Thanks a lot!

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

06-23-2011 07:15 AM

I think you will need a subject= specification in the random statement. Try:

PROC GLIMMIX data=test;

CLASS subjid;

MODEL disease= grass200m

/ dist=binary solution oddsratio;

random _residual_ / type=sp(exp) (x_cor y_cor) subject=subjid;

run;

Each subject should have its own set of x and y coordinates.

SteveDenham

PROC GLIMMIX data=test;

CLASS subjid;

MODEL disease= grass200m

/ dist=binary solution oddsratio;

random _residual_ / type=sp(exp) (x_cor y_cor) subject=subjid;

run;

Each subject should have its own set of x and y coordinates.

SteveDenham

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

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.

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

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;

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;