如果我要计算河北省每个城市的平均受教育水平，需要如下命令：
mean edu if region>130099 & region<130101
mean edu if region>130199 & region<130201
mean edu if region>130299 & region<130301
……
需要至少重复11次，如果我使用如下语句：
forv region=130100(100)131100 {
mean edu
}
结果stata把全国的平均受教育水平重复算了11遍，说明上述语句有问题。有没有简单方法能一次性分别算出11个城市的平均受教育水平（最好用循环语句，其中石家庄市所有个体的region代码均为130100，唐山130200，秦皇岛130300，其他城市以此类推）？