对于零膨胀泊松模型，模拟产生协变量x1,x2服从[0,1]均匀分布，然后产生服从ZIP模型的数据yi,其中yi的概率密度函数为

  yi=0时P（Yi=yi）=p(i)+(1-p(i))*exp(-lambda(i)),yi>0时P(Yi=yi)=(1-p(i))*lambda(i)^yi*1/(yi!)*exp(-lambda(i)),p(i),lambda(i)满足关系：

lambda(i)=exp(1+x1+x2),

p(i)=exp(0.5+0.5*x1+0.5*x2)/[1+exp(0.5+0.5*x1+0.5*x2))]

请问哪位大侠会写SAS程序啊？？？