I know one method can be:
Code:
regress y x1 x2 i.id, vce(robust)
But the problem is in my data set, I have a lot of id. And it kills my Stata.