楼主是要做分组回归把，那就是三个地区要做三个回归，都要用到 if 进行筛选，还要做一个总的全国回归，也就是四个回归。我的命令是这样的：
假设地区region有东中西三个赋值，

则：
xtreg lny lnx1 lnx2 lnx3 lnx4 b1.region
eststo m1
xtreg lny lnx1 lnx2 lnx3 lnx4 if  region==1
eststo m2
xtreg lny lnx1 lnx2 lnx3 lnx4 if region==2
eststo m3
xtreg lny lnx1 lnx2 lnx3 lnx4 if region==3
eststo m4

esttab m1 m2 m3 4,se(3) star(* 0.1 ** 0.05 *** 0.01) b(%6.3f) t(%6.3f) compress nogap ///
varwidth(10) nor pr2 replace

非常感谢，按你说的做出来了，不过我是分了四个区。但是不知道b1是什么意思呢。具体结果不太会解读，1.region 没有数据是什么原因呢，下面括号里的含义是什么啊，我这个结果还行吗

---------------------------------------------------------------------------
                 (1)          (2)          (3)          (4)          (5)   
               lngdp        lngdp        lngdp        lngdp        lngdp   
---------------------------------------------------------------------------
lntg          -0.006***    -0.011***    -0.004       -0.002       -0.013*  
             (0.002)      (0.003)      (0.003)      (0.003)      (0.007)   
lncy           0.254***     0.324***     0.349***     0.094***     0.685***
             (0.024)      (0.048)      (0.049)      (0.029)      (0.066)   
lnzx           0.136***     0.107***     0.140***     0.101***     0.389***
             (0.015)      (0.024)      (0.030)      (0.027)      (0.083)   
lnzb           0.599***     0.636***     0.535***     0.690***     0.337***
             (0.015)      (0.031)      (0.024)      (0.028)      (0.048)   
0.region       0.164**                                                     
             (0.064)                                                      
1.region       0.000                                                      
                 (.)                                                      
2.region      -0.134**                                                     
             (0.067)                                                      
3.region       0.213**                                                     
             (0.091)                                                      
_cons         -0.262*      -0.520**     -0.442        0.252       -4.005***
             (0.146)      (0.213)      (0.280)      (0.186)      (0.707)   
---------------------------------------------------------------------------
N                426          138          112          134           42   
pseudo R~q                                                                 
---------------------------------------------------------------------------
Standard errors in parentheses
* p<0.1, ** p<0.05, *** p<0.01

请问楼主如何31省生成3个地区的虚拟变量呢

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 "东北"
  之前参考的论坛答案，

十分感谢