PROC MODEL can correct the standard errors for heteroscedasticity using the HCCME= option on the FIT statement. This option does not allow you to specify the structure of the clustering of the data, but it may give you results closer to what you're looking for.
... View more