Heterogeneous difference in differences for panel data Use xthdidregress With STATA 18