想解决这样一个问题 从回归y=a+bx+e中得到残差e,假设样本容量是T，然后从e中随机抽取组成一个新的残差项e1，这个新的残差项的大小也是T，这样重复1000次，就组成1000组容量是T的新的残差项e1, 然后将e1加到原回归方程中得到1000组新的y, 令其为y1, 即y1=bx+e1,b 的顺序不变就是扩大为1000组b同时假设a=0， 最后再用新得到的y1和原b 做回归。

之前用非mata的语言写过，不过因为有上百组y，光进行一次这样的运算就耗费很长时间，不知道用mata是不是可以更便捷，特别是如果用mata怎么写这个随机1000次抽样的问题。。万分感谢。。。