Jim Silverton <jim.silverton <at> gmail.com> writes:
>
> Hello all,
> I am doing glm with a negative binomial link.
> I have two treatments and 3 replicates in each treatment. My question is
> this, how can I simulate data for the the columns from the null and
[alternative?]
> distribution.
>
Simulate from the null:
fit the null model and use the simulate method:
nullmodel <- glm.nb(response~1,data=yourdata)
simulate(nullmodel)
Simulate from the fitted distribution: as above, but
fit the full model.