[R] glm with negative binomial family

2006-06-09 Thread Elizabeth Rainwater
I am analysing parasite egg count data and am having trouble with glm with a negative binomial family. In my first data set, 55% of the 3000 cases have a zero count, and the non-zero counts range from 94 to 145,781. Eventually, I want to run bic.glm, so I need to be able to use glm(family=

Re: [R] glm with negative binomial family

2006-06-09 Thread Prof Brian Ripley
Did you consult the help page for neg.bin? Details: These are not intended to be called by the user. Some are for compatibility with earlier versions of MASS (the book). The book discusses the family negative.binomial() on p.206, and it is that which glm.nb uses. On Fri, 9 Jun

Re: [R] glm with negative binomial family

2006-06-09 Thread Elizabeth Rainwater
Many thanks to Dr Ripley, who told me to use negative.binomial instead of neg.bin, and to Dr Lancelot, who pointed toward the excellent aod. My models are now converging and the output I am getting seems reasonable. Cheers, Elizabeth On 6/9/06, Prof Brian Ripley [EMAIL PROTECTED] wrote: Did you