Re: [R] Something changed and glm(..., family=binomial) doesn't work now

From: White, Charles E WRAIR-Wash DC <charles.edwin.white_at_us.army.mil>
Date: Sun 19 Feb 2006 - 00:52:52 EST


"Are you using data=myframe in the glm call? Otherwise, you're open to masking by variables in your workspace."

FYI, the data statement doesn't seem to be the shortcut it is elsewhere when the form of the argument is glm(cbind(success,failure)~predictor, family=binomial, data=myframe). The data statement doesn't seem to put 'myframe' in the environment where cbind operates. Choices for the user are to fully quality names inside of cbind [i.e., cbind(myframe$success,myframe$failure)] or to attach 'myframe' to the working environment. The same thing happens in lmer.

Chuck

Charles E. White, Senior Biostatistician, MS Walter Reed Army Institute of Research
503 Robert Grant Ave., Room 1w102
Silver Spring, MD 20910-1557
301 319-9781
Personal/Professional Site: http://users.starpower.net/cwhite571/professional/



R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Sun Feb 19 00:59:29 2006

This archive was generated by hypermail 2.1.8 : Mon 20 Feb 2006 - 14:08:41 EST