以下代码不是mata，但运行一遍只需55秒左右，我觉得时间还是可以接受的（最长的时间花在1000次回归上，bsample只用了4秒）。

set rmsg on
clear
set obs 900   // T
gen x1 = runiform()
gen x2 = runiform()
gen y = 2*x1 + 3*x2 + rnormal()
reg y x1 x2
predict yhat
predict res, residual

expand 1000
bysort res: gen n = _n
sort n

preserve
keep res
bsample 900000  // T*1000
save res, replace
restore

drop res
merge 1:1 _n using res, nogen

gen ynew = yhat + res
statsby _b, by(n) clear: reg ynew x1 x2

你好，看到你的回答，你应该是比较熟悉用stata来抽样的了，我现在有个数据是640万个样本，要抽取30% 的样本，stata打开总体后，写什么语句可以抽取出来这30%的样本呢？希望能够指导下，万分感谢！

抽取30%的样本有很简单的命令：
sample 30
就可以解决问题了

其实，楼主所说的问题就是bootstrap吧，可以在回归命令前面加上bootstrap就可以实现了：
bootstrap,reps(1000):reg y x