similar to: R-beta: More glm bugs

A non-text attachment was scrubbed... Name: not available Type: text Size: 686 bytes Desc: not available Url : https://stat.ethz.ch/pipermail/r-help/attachments/19971216/6e13087c/attachment.pl

G'day all, I had a look at the GLM code of R (1.4.1) and I believe that there are problems with the function "glm.fit" that may bite in rare circumstances. Note, I have no data set with which I ran into trouble. This report is solely based on having a look at the code. Below I append a listing of the glm.fit function as produced by my system. I have added line numbers so that I

A non-text attachment was scrubbed... Name: not available Type: text Size: 997 bytes Desc: not available Url : https://stat.ethz.ch/pipermail/r-help/attachments/19980618/ee08ba8d/attachment.pl

A non-text attachment was scrubbed... Name: not available Type: text Size: 750 bytes Desc: not available Url : https://stat.ethz.ch/pipermail/r-help/attachments/19980403/5ef52eac/attachment.pl

A non-text attachment was scrubbed... Name: not available Type: text Size: 668 bytes Desc: not available Url : https://stat.ethz.ch/pipermail/r-help/attachments/19980514/af07cf17/attachment.pl

A non-text attachment was scrubbed... Name: not available Type: text Size: 1262 bytes Desc: not available Url : https://stat.ethz.ch/pipermail/r-help/attachments/19980602/3cf08662/attachment.pl

A non-text attachment was scrubbed... Name: not available Type: text Size: 1262 bytes Desc: not available Url : https://stat.ethz.ch/pipermail/r-help/attachments/19980602/3cf08662/attachment.pl

Hi R-help, Is there such a thing as a function in R for fitting a GLM where the response is distributed as a Beta distribution? In my case, the response variable is a percentage ([0,1] and continuous). The current glm() function in R doesn't include the Beta distribution. Thank you for any help on this topic. Sincerely, Sharon K?hlmann

A non-text attachment was scrubbed... Name: not available Type: text Size: 1659 bytes Desc: not available Url : https://stat.ethz.ch/pipermail/r-help/attachments/19980107/3a4902ef/attachment.pl

A non-text attachment was scrubbed... Name: not available Type: text Size: 1659 bytes Desc: not available Url : https://stat.ethz.ch/pipermail/r-help/attachments/19980107/3a4902ef/attachment.pl

I have been implementing a proposal of Jim Lindsey for glm(.) to return AIC values, and print.glm(.) and print.summary.glm(.) printing them.... however: >>>>> "Jim" == Jim Lindsey <jlindsey@luc.ac.be> writes: Jim> The problem still remains of getting the correct AIC when the user Jim> wants the scale parameter to be fixed. (The calculation should

Is there any logical reason why glm prints out the labels of factor levels after variable names when baseline contrasts (contr.treatment) are used but the codes for the levels when mean contrasts (contr.sum) are used? Jim -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info",

A non-text attachment was scrubbed... Name: not available Type: text Size: 478 bytes Desc: not available Url : https://stat.ethz.ch/pipermail/r-help/attachments/19971103/ae2f1119/attachment.pl

A non-text attachment was scrubbed... Name: not available Type: text Size: 1287 bytes Desc: not available Url : https://stat.ethz.ch/pipermail/r-help/attachments/19980406/39dc7891/attachment.pl

A non-text attachment was scrubbed... Name: not available Type: text Size: 645 bytes Desc: not available Url : https://stat.ethz.ch/pipermail/r-help/attachments/19980424/94b87a18/attachment.pl

A non-text attachment was scrubbed... Name: not available Type: text Size: 520 bytes Desc: not available Url : https://stat.ethz.ch/pipermail/r-help/attachments/19980424/b1195857/attachment.pl

A non-text attachment was scrubbed... Name: not available Type: text Size: 1186 bytes Desc: not available Url : https://stat.ethz.ch/pipermail/r-help/attachments/19980814/d985de5f/attachment.pl

A non-text attachment was scrubbed... Name: not available Type: text Size: 1215 bytes Desc: not available Url : https://stat.ethz.ch/pipermail/r-help/attachments/19980225/6146ab11/attachment.pl

>>>>> "Jim" == Jim Lindsey <jlindsey@alpha.luc.ac.be> writes in a private mail to me Jim> That bug in seq() really has shaken my faith in R. How Jim> many other known bugs are there that we on r-help are not informed Jim> of and which can falsify our statistical analyses??? It is Jim> imperative that patches for such things be made

Dear all, I am taking a course in GLM given by a devoted SAS user. He has given us a homework where a mixed GLM with a logistic link and binomially distributed observations should be fitted. I know of the library nlme for mixed effect models but as I understand it, one cannot choose between different links and distributions in the functions provided there. I have so far managed very well with