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
- /
- PROC GLM/PROC REG "Strata statement" equivalent

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

01-27-2016 03:06 PM

I have a regression analysis in which I have a continous outcome but I also have data which is matched on certain variables. Originally this analysis was completed by using PROC GLM (with a different sample not case-control matched) with a simple code like the following:

proc glm data=ds;

model outcome = clinical_trait_1;

weight iwm;

*also had an inverse weighting method;

run;

If I had a binary outcome I would use PROC LOGISTIC with a strata statement for a conditional logistic regression method- but in this instance I'm attempting to use the PROC GLM framework: there is no strata statement in this or PROC REG. Is there an equivalant?

A few ideas I've had so far but haven't had much success:

*PROC SURVEYREG -- Doesn't seem to have same output to match PROC REG/GLM prior to accounting for strata statement.

*PROC GENMOD with dist=normal link=identity -- Doesn't seem to match output found in PROC REG/GLM again.. Seems closer, but still differences in standard error.

*PROC GLM/PROC REG with a random statement-- would this work for a matched ID?

Any insight on this would be greatly appreciated. (on SAS 9.4)

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

01-27-2016 04:45 PM

Yes, I agree it would be nice to have strata statements in proc glm/proc reg. What it should do is to allow to include a group variable in the model, and it should be interpreted as conditioning with the sum within each level of the strata variable.

I tried apply some semiparametric theory, which in principle work like projecting the score functions of the parameters of interest into the orthogonal space to the score functions of the nuisance parameters belonging the strata variable. This give estimating equations only for the parameters of interest. (same princip as how strata works in cox regression). What a strata statement should do is then to tell the procedure to instead only solve the estimating equations for parameters in the model line, instead of the "raw" score equations.

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

01-28-2016 10:14 AM

Interesting. Recently we were also thinking about using the cluster statement in place also.