G�ran Brostr�m <[EMAIL PROTECTED]> writes: > > Try it with > > > > b=c(rep(0,8),rmultinom(1,24,rep(1/8,8))) > > This is the conditional distribution, given the sufficient statistic.
Yes (and that was of course the point). I was surprised that there is a difference at all (between this and rep(3,8) and also between different calls to rmultinom). Theoretically, the MLE should be the same, but of course the "saturated models" have a different likelihoods, so the deviances vary and there may be variations in the number of iterations taken in the (divergent) fitting process. -- O__ ---- Peter Dalgaard Blegdamsvej 3 c/ /'_ --- Dept. of Biostatistics 2200 Cph. N (*) \(*) -- University of Copenhagen Denmark Ph: (+45) 35327918 ~~~~~~~~~~ - ([EMAIL PROTECTED]) FAX: (+45) 35327907 ______________________________________________ [email protected] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
