# Adding a square root link function to an overdispersed negative binomial GLM

I'm analyzing nematode count data (80 data points) from a randomized block design in which I have two factors with both four levels (Plant and Inoc). The data show heavy overdispersion when analyzed with a glm with poisson-distribution. For one of the plant levels counts are very low (below 100), while for others they almost reach 100000. Therefore I analyzed the data with a Negative-binomial distribution (glm.nb in R MASS-package) which results in limited overdispersion of ~1.5 (Res. def. 89.12 over 60DF). The predicted versus residual plots look fine, but I'm still worried about the amount of overdispersion. I thought it might be worth it to try to change the link function from log to sqrt to see whether that helps. To do this you have to give starting values for your coefficients, like model<-glm.nb(Nem_gramroot ~ Block + Plant + Plant*Inoc, data=data, link = sqrt, start = c(10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10)). I've tried putting in 0's, 1's, 10's an the coefficients of the model based on the log link function. However, I always get the message step size truncated: out of bounds. Any ideas on this? And moreover, any thoughts on using sqrt link functions to reduce overdispersion in negative binomial glms?