Hi Peter, On Tue, Nov 27, 2012 at 8:05 PM, Peter Ehlers <ehl...@ucalgary.ca> wrote: > > On 2012-11-27 14:34, Steve Lianoglou wrote: [snip] > Steve: > re a matrix response: see MASS (the book, 4ed) page 191; also found > in the ch07.R file in the /library/MASS/scripts folder. I seem to > recall that this is mentioned somewhere in the docs, but put my finger > on it now.
Indeed -- thanks for the pointer. Well then ... let me try to recover and at least offer some decent advice to Craig's original post. I'd still recommend avoiding the use of `attach` and favor passing in your data data.fame using the `data` param of the call to `glm`. I'd also call the `avoid` data.frame something else, to avoid potential confusion (if only in your mind) with the column `avoid` inside the `avoid` data.frame. Anyhow, when you pass in your data as param, things work ok: R> model1 <- glm(cbind(avoid, noavoid) ~ treatment, binomial, avoid) R> summary(model1) Call: glm(formula = cbind(avoid, noavoid) ~ treatment, family = binomial, data = avoid) Deviance Residuals: Min 1Q Median 3Q Max -3.6480 -1.3350 0.4297 1.5706 2.6200 ... OK ... well, hope that helped somewhat. -steve -- Steve Lianoglou Graduate Student: Computational Systems Biology | Memorial Sloan-Kettering Cancer Center | Weill Medical College of Cornell University Contact Info: http://cbio.mskcc.org/~lianos/contact ______________________________________________ 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.