destring stkcd, generate(yy)
xtset yy year,yearly
egen ind_year=group(ind year),label lname(ind_year)
sum ind_year
global N=r(max)
reg y x1 x2
xtreg y x1 x2, fe