Hi everybody!

I have a lingitudinal dataset (individuals x years, individuals live in different areas) and I am using the following command:
reghdfe y x, absorb(area=FE year) vce(cluster pro)
where FE are the estimated fixed effects for the N areas. How can I save the cov matrix for the FE parameters? I'd like to have a NxN matrix in the .dta file.

Thanks!

Giorgia