我用的工业企业数据库，想计算每年每个城市的行业HHI，用就业人数算的，目前我可以算到某一城市中某一行业总人数占这个城市所有行业总人数的平方，但是怎么用命令将这一平方按照行业加总呢？我想得到的结果是1101的hyHHI就是sic为13、14、15、16、17的rate^2的加总（五个数字的加总），我在编程中的结果总是所有rate^2的加总（9个数字的加总），不知道怎么实现只加总相同城市不同行业代码的rate^2。。。

year	sic	fipss	emp	hyemp	tot_hyemp	rate	rate^2	hyHHI
2005	13	1101	12	16	70	0.228571	0.0522449
2005	14	1101	8	22	70	0.314286	0.0987755
2005	13	1101	4	16	70	0.228571	0.0522449
2005	16	1101	6	14	70	0.2	0.04
2005	14	1101	4	22	70	0.314286	0.0987755
2005	15	1101	7	7	70	0.1	0.01
2005	16	1101	8	14	70	0.2	0.04
2005	17	1101	11	11	70	0.157143	0.0246939
2005	14	1101	10	22	70	0.314286	0.0987755
2005	13	1201	13	41	88	0.465909	0.2170713
2005	16	1201	16	23	88	0.261364	0.068311
2005	13	1201	5	41	88	0.465909	0.2170713
2005	14	1201	7	7	88	0.079545	0.0063275
2005	15	1201	9	17	88	0.193182	0.0373192
2005	13	1201	10	41	88	0.465909	0.2170713
2005	15	1201	8	17	88	0.193182	0.0373192
2005	16	1201	7	23	88	0.261364	0.068311
2005	13	1201	13	41	88	0.465909	0.2170713