连老师：
您好！
我在学习Monte Carlo Simulation的编程，碰到一个小问题，不知道怎么回事。
运行下面的程序，无法得出结果。return list命令也没有东西。 如果最后用simulate  _b _se, reps(10) nodots: ppmlsim_k1 x1 x2
则错误信息为“no variables defined r(111);”。请指点一下哪里出了问题。（另外，我用OLS回归模拟没有出现这个问题。模拟结束会返回关于_b _se的信息，不知道这个该怎么设置？）

谢谢！

*--------------------------------------------------------------------------------
program define ppmlsim_k1,rclass
version 10.0
clear
set obs 1000
gen x1=rnormal()
gen x2 = uniform()
replace x2 = (x2<=0.40)
gen v=1/(exp(1*x1+1*x2))
gen f=rnormal(-(ln(v+1))/2,ln(v+1))
gen y=exp(1*x1+1*x2+f)
syntax varlist
tempvar y1 k
       gen `k'=1
        gen `y1' = y-`k'
        replace `y1'=0 if `y1'<0
        poisson `y1' `varlist',robust
        
end
*------------------------------------------
      simulate , reps(10) nodots: ppmlsim_k1 x1 x2
      sum