BookmarkSubscribeRSS Feed
econjoe
Calcite | Level 5

Is there a SAS equivalent to Stata's areg command? Specifically, I'm looking for a procedure that will replicate the following Stata command:

 

areg depvar indvar, absorb(id1) cluster(id2)

 

In this case id1 is nested within id2.

 

The proc genmod below clusters the standard errors at the id2 level, but is not able to absorb id1.
 
proc genmod;
  class id2;
  model depvar = indvar;
  repeated subject=id2 / type=ind;
run;
 
The proc glm below absorbs id1 but cannot cluster the standard errors at the id2 level.
 
proc glm;
  absorb id1;
  model depvar = indvar / solution noint;
run;
 
Is there a simple way to both absorb id1 and cluster on id2?
1 REPLY 1
SteveDenham
Jade | Level 19

I am not real sure, but would some sort of hierarchical model, such as GLIMMIX work for what you want?

 

proc glimmix data=have;

class id1 id2;

model depvar=indvar;

random intercept/subject=id1;

random _residual_/subject=id2(id1) type=vc;

run;

 

Steve Denham

sas-innovate-wordmark-2025-midnight.png

Register Today!

Join us for SAS Innovate 2025, our biggest and most exciting global event of the year, in Orlando, FL, from May 6-9. Sign up by March 14 for just $795.


Register now!

What is Bayesian Analysis?

Learn the difference between classical and Bayesian statistical approaches and see a few PROC examples to perform Bayesian analysis in this video.

Find more tutorials on the SAS Users YouTube channel.

SAS Training: Just a Click Away

 Ready to level-up your skills? Choose your own adventure.

Browse our catalog!

Discussion stats
  • 1 reply
  • 3200 views
  • 0 likes
  • 2 in conversation