dear listers, i tried to use glm.nb to estimate a nega. binomial but have no luck to get the result. here is the code: model <- glm.nb(Y ~ ., data = mydata) I am not sure if I have missed anything. thanks for your insight! wensui