Dear R Community, I have encountered a problem while using the R function glm.nb. The code that produce the error was following two lines:
group=c(1,1,1,1,0,0,0,0) fit=glm.nb(y~group) While the y contains 8 sets of number like: gene275 0 1 0 0 1 5 1 0 Error message: Error in while ((it <- it + 1) < limit && abs(del) > eps) { : missing value where TRUE/FALSE needed Calls: glm.nb -> as.vector -> theta.ml In addition: There were 50 or more warnings (use warnings() to see the first 50) Execution halted Information of my system: > sessionInfo() R version 3.0.1 (2013-05-16) Platform: x86_64-unknown-linux-gnu (64-bit) locale: [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C [3] LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8 [5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8 [7] LC_PAPER=C LC_NAME=C [9] LC_ADDRESS=C LC_TELEPHONE=C [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C attached base packages: [1] stats graphics grDevices utils datasets methods base Does anyone happen to have some hit on how to solve this? Appreciate for any response. Thanks in advance, Daofeng [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.