[R] posthoc tests with ANCOVA
Torsten Hothorn
Torsten.Hothorn at rzmail.uni-erlangen.de
Wed Jan 10 11:02:21 CET 2007
On Tue, 9 Jan 2007, Walter Durka wrote:
> dear all,
>
Walter,
_please_ cc questions on contributed packages to the maintainer, since not
everybody follows r-help closely...
> I want to perform a posthoc test for my ANCOVA:
> a1<-aov(seeds~treatment*length)
>
> With
> summary(glht(a1, linfct = mcp(treatment = "Tukey")))
> R tells me: "covariate interactions found -- please choose appropriate
> contrast"
>
one needs to specify a certain value for `length' which, I assume, is a
numeric covariate, right? The current interface doesn't support this (this
on my to-do-list), however, you can set up the matrix of linear functions
by yourself (contact me privately if you have problems to do that).
> How do I build these contrasts?
>
> Ideally, I would like to have the posthoc test for the ANCOVA including
> a block-effect
> a2<-aov(seeds~treatment*length+Error(site))
>
> How do I make a posthoc test here?
>
its on the to-do-list as well :-(
Torsten
> Thanks for any comments
> Walter
>
>
> --
>
> *****
> Dr. Walter Durka
> Department Biozönoseforschung
> Department of community ecology
>
> Helmholtz-Zentrum für Umweltforschung GmbH - UFZ
> Helmholtz Centre for Environmental Research - UFZ
> Theodor-Lieser-Str. 4 / 06120 Halle / Germany
>
> walter.durka at ufz.de / http://www.ufz.de/index.php?en=798
> phone +49 345 558 5314 / Fax +49 345 558 5329
>
> +++++++++++++++++++++++++++++++++++++++++++++++++
>
> Das UFZ hat einen neuen Namen:
> Helmholtz-Zentrum für Umweltforschung GmbH - UFZ
>
> The UFZ has a new name:
> Helmholtz Centre for Environmental Research - UFZ
>
> +++++++++++++++++++++++++++++++++++++++++++++++++
>
> ______________________________________________
> R-help at 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.
>
>
More information about the R-help
mailing list