请教一下各位前辈！我现在要把地级市分成三个地区（中西东部）的虚拟变量，进行地区异质性分析；我找了很多帖子，都是根据数量较少的省份划分，但是31个省份相对于我的地级市真的少太多了，我不知道该咋分，感觉照上面的方法工作量会很大，他们的划分方法是这样：gen east = inlist(province,2,4,6,9,10,15,23,24,27,31)
gen west = inlist(province,3,5,8,20,21,22,26,28,29,30)
gen mid = inlist(province,1,7,12,13,16,19,25,14)
gen northeast = inlist(province,11,17,18)

gen region=0 // 东部
replace region=1 if province==1 |  province==7 |  province==12 |  province==13 |  province==16 | ///
province==19 |  province==25 |  province==14  //中部
replace region=2 if  province==3 |  province==5 |  province==8 |  province==20 |  province==21 | ///
province==22 |  province==26 |  province==28 |  province==29 |  province==30  //西部
replace region=3 if  province==11 |  province==17 |  province==18  //东北

label define region_lab 0"东部" 1 "中部" 2 "西部" 3 "东北"
怕格式乱，附上这个方法的截图