抱歉，你写的程序有点没看懂。
这是我写的程序，只有前半部分。
local y=1
while `y'<=100   {
set seed 100
preserve
insheet using C:\Users\dell\Desktop\12345.csv
reg  taccta ta revrecta ppea roa1,noc
predict e,r
gen name=_n
keep  v1 e name
sort name
save 1,replace
sample 20
sort name
save 2,replace
use 1
sort name
merge name using 2
gen part=1 if _merge==3
replace part=0 if part==.
reg  e part   ——————————————————————主要的回归
gen ncoeff=_b[part] ——————————————————之前你告诉我的那个语句   
【这里该怎么写】
restore
local y=`y'+1
}
我需要把这100次的part的回归系数保存到一个文件中，然后在对这100个系数做T检验或是别的检验。
同时还想问你一下STATA怎么做二项式检验。

试试这个

不好意思，请问你给我的那两句加在哪里？可以通过给你发消息想你请教吗？

就是把你标注的那两句换为这两句命令
你可以把自己的数据修改，然后传一个样本数据到论坛，可以方便调试。

insheet using C:\Users\dell\Desktop\12345.csv
reg  taccta ta revrecta ppea roa1,noc
predict e,r
gen name=_n
keep  v1 e name
sort name
save 1,replace
sample 20
sort name
save 2,replace
use 1
sort name
merge name using 2
gen part=1 if _merge==3
replace part=0 if part==.
reg  e part
gen ncoeff=_b[part]
以上这步会得到这样的结果
      |     v1           e   name   _merge   part     ncoeff |
      |------------------------------------------------------|
   1. |      2   -.0596653      1        1      0   .0348098 |
   2. |      4     .132309      2        1      0   .0348098 |
   3. |      5    .1781067      3        1      0   .0348098 |
   4. |      6   -.1426554      4        1      0   .0348098 |
   5. |      7    .4624781      5        1      0   .0348098 |
      |------------------------------------------------------|
   6. |      8    .0878844      6        1      0   .0348098 |
   7. |      9    .2980676      7        1      0   .0348098 |
   8. |     11   -.2907471      8        1      0   .0348098 |
   9. |     12    2.144006      9        1      0   .0348098 |
  10. |     14   -.3491684     10        1      0   .0348098 |
      |------------------------------------------------------|
  11. |     16    .2977334     11        1      0   .0348098 |
  12. |     18    .4632861     12        3      1   .0348098 |
  13. |     19    .5795037     13        1      0   .0348098 |
  14. |     20    1.202998     14        3      1   .0348098 |
  15. |     21    .9889596     15        1      0   .0348098 |
最有一列不是那个系数吗？我需要将它带入到另一个数据文件中。这样用循环语句做下来不就可以得到100次随机抽样回归得到的part的系数了吗？可是我不知道怎么编这段程序。

没有你的数据，你按照自己做小改动试试

非常感谢您。实在不好意思再麻烦您，可是我还不知道能找谁帮忙，您给我写的程序我执行了，可是执行不了。我把实验的数据传上来，麻烦您再帮我调试一下。

不需要分年抽样了，对总样本随机抽样20%，重复100次就行。辛苦您了。

多谢高人，在你的帮助下，刚才我试出来了，谢谢您。

模型大体的程序就是这样
其实用我刚开始给的simulate那个程序会好些
因为Stata自带的simulate比较方便。

您好，还需要再麻烦你一下。通过您之前的帮助我已经知道如何将part的回归系数提取出来了，现在还有一个问题我还需要提取出part的p-value,我在网上找了一下有以下这两段命令：
local t = _b[part`i']/_se[part`i']
local p=2*ttail(e(df_r),abs(`t'))
这是算p-value的，我现在不知道如何将它增加到矩阵中，按照您先前叫我的，这样做吗？
mat b[`i',1] = p[part`i']  可是我没做出来，希望您还能帮助我。

你好，请问如果是1000次重复回归那，那两个10应该怎么改？

您好，我您给我的代码试验了一次，有一个错误提示conformability error，请问这个怎么办？可不可以分两次做，就是说一套代码专门获得系数的，一套代码专门用来获得P值得

这个是可以的吧，如果不行，把你的数据上传。

太好了，您真强，成功了。请问您是怎么学的STATA，感觉您逻辑太清晰了。

求问，您试出来最终是代码是怎么编的呢？还想请教下您后来二项检验的代码是怎么编的？以及文献中说对jones模型进行分年度分行业回归，最终得出的系数是一个，怎么理解您解决了吗？