我大概要run以下这两个方程：
1. 第一步，Loss=a+bx+u 是一个panel data，但是我已经生成好了dummy variable， 所以直接用reg
所以我的code大致是： reg loss x, robust (robust 加不加另讨论)

2. 第二步，我要run一个多次方程：
y=a+b*Loss+c*Loss^2+d*Loss^3+u
由于上一步的Loss是genenrate生成的，并不能保证满足OLS的要求，所以我必须要做一个bootstrap

那这个程序在stata应该大致如何实现呢？给出类似的例子也好，我Cameron的那本书bootstrap是用heckman two stage的讲了老半天，我没看懂…………T_T