I have the following code for my two-way fixed effect model. In my dataset "data1", there are a few "units" in each "firm". How I can get the heteroskadasticity-robust standard errors clustered by "firm". proc sort data=data1; by unit time; run; proc panel data=data1; id unit time; model DV = x /FIXTWO; run;
... View more