Hi, I am trying to compare negative binomial models for the prediction of sports games (I know that Poisson models would be better but I'm just trying Negative Binomial at the moment). But, to compare the models I need them to have the same theta value. How can I change the explanatory variables while maintaining the theta value? Thanks for the help. David -- View this message in context: http://r.789695.n4.nabble.com/Comparing-negative-binomial-models-tp3985353p3985353.html Sent from the R help mailing list archive at Nabble.com.
Comparing negative binomial models
2 messages · Davg, Ben Bolker
Davg <davidgrimsey <at> hotmail.com> writes:
I am trying to compare negative binomial models for the prediction of sports games (I know that Poisson models would be better but I'm just trying Negative Binomial at the moment). But, to compare the models I need them to have the same theta value. How can I change the explanatory variables while maintaining the theta value?
library(MASS) modelfit <- glm(..., family=negative.binomial(theta=theta_value)) I believe that drop1() applied to a glm.nb() fit does this correctly (i.e. holds the theta parameter fixed at the estimate from the most complex model).