Re: [R] Schmera: a question on R software

From: Douglas Bates <bates_at_stat.wisc.edu>
Date: Wed 18 Oct 2006 - 14:34:52 GMT

On 10/18/06, Torsten Hothorn <Torsten.Hothorn@rzmail.uni-erlangen.de> wrote:
>
> On Wed, 18 Oct 2006, SCHMERA Dénes wrote:
>
> > Dear All,
> >
> > I would like to run a generalized linear mixed model with the software R (one categorical predictor, one random factor, the distribution of the dependent variable is binomial, and the link is logit). Thereafter, I would like to perform multiple comparisons (post hoc test) among the groups of the categorical predictor.
> >
> > Is it possible with the software R?
> >
>
> yes and no. There is the `multcomp' package available but the CRAN version
> is limited to linear models (more or less). However, I'm currently working
> on a re-implementation which has solutions for your problem and much
> more :-)
>
> For the time being, the development version is available from
>
> http://www.imbe.med.uni-erlangen.de/~hothorn/multcomp2
>
> The basic steps would be as follows:
>
> ### fit logistic regression model
> > my.glm <- glm(response ~ myfactor, data = <mydata>, family = binomial())
> ### set up linear hypothesis for coef(my.glm) corresponding to all-pair
> ### comparisons
> > library("multcomp2")
> > gh <- glht(my.glm, linfct = mcp(myfactor = "Tukey"))
> ### simultaneous confidence intervals with plots
> > confint(gh)
> > plot(confint(gh))

The original question was about generalized linear mixed models which is going to be even more difficult. :-(

> > Are traditional methods of multple comparisons (for instance Tukey test) adequate for this purpose or not? Does Tukey test consider the binomial dsitribution of the dependent variable and the logit link?
> >
> > I am waiting for your answer,
> > Denes Schmera
> >
> > ______________________________________________
> > 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
> > and provide commented, minimal, self-contained, reproducible code.
> >
>
> ______________________________________________
> 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
> and provide commented, minimal, self-contained, reproducible code.
>
>
>



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 and provide commented, minimal, self-contained, reproducible code. Received on Thu Oct 19 01:13:33 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Wed 18 Oct 2006 - 15:30:10 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.