Displaying 1 result from an estimated 1 matches for "4v2".
Did you mean:
42
2010 Oct 22
1
getting all contrasts from glm
...glm(cbind(events,patnums-events) ~ as.factor(treat) + as.factor(numtrial), family = binomial )
Is this example there are 4 treatments . In the glm object I can find the contrasts of the main treats vs the first i.e. 2v1, 3v1 and 4v1 ... however I would like to get the complete set including 3v2, 4v2, and 4v3 ... along with the Std. Errors of all contrasts. I've tried the contrasts argument of glm and attempted to use model.matrix.default but with no luck.
Any suggestions how to get all treat contrasts with SE's?
Thanks
J
===============================
Dr. Jim Maas
University of...