用stata做DF检验，写了以下的命令，想重复1000次蒙特卡洛模拟，并计算拒绝原假设的概率，要怎么编写命令啊！

set obs 100

gen t=_n

tsset t

gen u=invnormal(uniform())

gen y=0 in 1

replace y=0.1+ 0.01*t+ l.y + u if t >1  

dfuller y, trend   

各位大神求助！