Hi all,
I am using the negative binomila glm in MASS. This is my data alled data1:


     Reps Treats counts
HSM1    1      0     21
HSM2    2      0     34
HSM3    3      0     27
PTM1    1      1     32
PTM2    2      1     20
PTM3    3      1     23


I goal is to do a GLM to extract the p-values from Treatment. row.
However, I can access the coefficients but not the p-values.

On another note, I want to simulate from the null. Is it more appropriate to
use:

nullmodel <- glm.nb(counts~1,data=data.use)

OR

should I use:  nullmodel <- glm.nb(counts~ Reps,data=data.use)  instead?

FInally, if I simulate 10 data simulations of the counts, Is this the way to
store the p-values?

response.sim = simulate(nullmodel,10)
for (j in 1:10)
{
counts.null = response.sim[,j]
new.data.mat = data.frame(Reps,Treats,counts.null)
new.glm <- glm.nb(counts.null ~ Reps, data = new.data.mat)
storage.ps[j] <- 1-pchisq(new.glm$null.deviance,df=new.glm$df.null)
 }


The null is there is no differences between the 2 treatments
The alternative is there is a difference between the 2 treatments.

Any help is greatly appreciated.

-- 
Thanks,
Jim.

	[[alternative HTML version deleted]]

