小弟现在在用工业企业数据库，需要计算同一年份、同一地区（省份代码前两位11-15、21-23、31-37、41-46、50-54、61-65共31个地区）、同一二级行业（2位代码11-43其中39缺失）下的某一变量hhi（各公司占全行业比重平方和）及CR4（前四大公司之和/全行业）指数
无奈刚学stata不久，循环语句更是看的云里雾里。
————————————————————————————————————————
下面是我的笨方法，但如果真的这样操作岂不是得重复劳动上万次。。。
*处理07年、代码为11地区、次级行业代码为13的公司数据
egen rd071113=sum(syrd) if dq_sub==11& hydm_sub==13& year=2007              //求07年11地区，行业代码为13的syrd和
gen h071113=(syrd/(rd071113))^2 if dq_sub==11& hydm_sub==13& year=2007//求07年11地区，行业代码为13的各公司syrd 占行业syrd比重的平方
egen HHI071113=sum(h071113) if dq_sub==11& hydm_sub==13& year=2007   //求出syrd的HHI指数
egen cr071113=sum(syrd) in 1/4 if dq_sub==11& hydm_sub==13& year=2007   //求出07年11地区，13行业syrd最大的前4家公司syrd和（之前已对syrd排序）
gen CR=cr071113/rd071113 if dq_sub==11& hydm_sub==13& year=2007          //07年11地区，13行业syrd最大的前4家公司syrd和占行业syrd的比重 ————————————————————————————————————————————
所以简单地说就是如何能让stata自行计算不同年份不同地区不同行业下某一变量的HHI指数及CR4甚至GINI指数？