最近看了一篇文献，想复现该篇文献中的模拟，其中有一句是”To produce count outcomes, we draw y from anegative binomial distribution with conditional mean E[ y  | x ] = exp (β1 *X1  + x 2  ) and overdispersion parameter αNB , which captures deviations from the conditional mean-variance equality restriction imposed by Poisson re- gression.“

求助怎么生成这种条件均值的负二项分布随机数呢

我现在写的code如下：
g Ey=exp(2*x1+x2)
g p=Ey/(Ey+2)
g y=rnbinomial(5000,p)


按照我现在的理解是我用X1和X2计算一个条件均值，然后overdispersion parameter αNB我设置的是0.5，所以r就等于2，根据负二项分布的均值计算出P，我一共要生成5000个观测值

请问我这个做法是正确的吗？因为我发现我这样生成的数据好像不太对，但我也不知道错在哪儿，求各位大神帮我看看