Hola,
Quizas puedas considerar tus contrastes como comparaciones o contrastes
planificados (planned contrasts o creo que también llamados por algunos a-priori
contrasts).
Una búsqueda rápida da algunos buenos ejemplos:
https://www.google.es/search?q=r+anova+planned+contrasts
Por ejemplo
http://faculty.smu.edu/kyler/courses/7311/planned_4up.pdf
Saludos,
Pedro
El 30/03/2015 a las 12:00, r-help-es-request en
r-project.org<mailto:r-help-es-request en r-project.org> escribió:
Message: 2
Date: Sat, 28 Mar 2015 23:10:04 +0100
From: Carlos Hernández-Castellano
<carlos.hernandezcastellano en
gmail.com><mailto:carlos.hernandezcastellano en gmail.com>
To: r-help-es en r-project.org<mailto:r-help-es en r-project.org>
Subject: [R-es] Comparaciones múltiples
Message-ID:
<CAFqSSGjvUN3Tssut8gYO0GGaKtD8MYMSzLd2sTNBwfW6O9foHA en
mail.gmail.com><mailto:CAFqSSGjvUN3Tssut8gYO0GGaKtD8MYMSzLd2sTNBwfW6O9foHA
en mail.gmail.com>
Content-Type: text/plain; charset="UTF-8"
Saludos,
tengo una base de datos con tres variables: especie,
tratamiento y
abundancia.
Quiero hacer comparaciones múltiples con Tukey.
Ya había usado el comando TukeyHSD(aov(x~y)).
La cuestión esque ahora tengo varios niveles del factor
especie y varios
niveles del factor tratamiento.
Si hago TukeyHSD(aov(abundancia~especie:tratamiento)) me
salen todas las
comparaciones posibles, pero yo sólo estoy interesado en
que para cada
especie (i.e. para cada nivel del factor especie) me salgan
las
comparaciones de los datos de abundancia para cada par de
tratamientos
(niveles del factor tratamiento).
Una solución laboriosa sería partir el documento original
en tantos como
especies, pero seguro que alguien me puede sugerir una
solución más
inteligente.
Saludos y gracias,
--
*?*
*Carlos Hernández-Castellano*
Environmental Scientist
Student, MSc Terrestrial Ecology and Biodiversity
Management
Research Collaborator at Centre of Ecological Research and
Forestry
Applications (CREAF)
Ecology Unit - Department of Animal Biology, Plant Biology
and Ecology;
Faculty of Sciences, Autonomous University of Barcelona
(UAB)
08193 Bellaterra, Spain
Email: carlos.hernandezcastellano en
gmail.com<mailto:carlos.hernandezcastellano en gmail.com>
[[alternative HTML version deleted]]
------------------------------
Message: 3
Date: Sat, 28 Mar 2015 23:31:44 +0000
From: Víctor Granda García <victorgrandagarcia en
gmail.com><mailto:victorgrandagarcia en gmail.com>
To: Carlos Hernández-Castellano
<carlos.hernandezcastellano en
gmail.com><mailto:carlos.hernandezcastellano en gmail.com>,
r-help-es en r-project.org<mailto:r-help-es en r-project.org>
Subject: Re: [R-es] Comparaciones múltiples
Message-ID:
<CAOpgSEaTESSJj0kju0y7DuMtuomd8ZR+A3BOv3hY3TmPPn4yWw en
mail.gmail.com><mailto:CAOpgSEaTESSJj0kju0y7DuMtuomd8ZR+A3BOv3hY3TmPPn4yWw
en mail.gmail.com>
Content-Type: text/plain; charset="UTF-8"
Hola Carlos.
La orden TukeyHSD tiene un parámetro "which" que te permite
indicar para
que variable explicativa quieres las comparaciones.
Mira la ayuda de TukeyHSD con '?TukeyHSD' donde tienes un
ejemplo, aunque
en tu caso sería algo parecido a esto:
modelo <- aov(abundancia~especie:tratamiento)
TukeyHSD(modelo, "tratamiento")
Espero que te sirva, un saludo.
El sáb., 28 de marzo de 2015 a las 23:10, Carlos
Hernández-Castellano (<
carlos.hernandezcastellano en gmail.com<mailto:carlos.hernandezcastellano en
gmail.com>>)
escribió:
> Saludos,
>
> tengo una base de datos con tres variables: especie,
tratamiento y
> abundancia.
>
> Quiero hacer comparaciones múltiples con Tukey.
> Ya había usado el comando TukeyHSD(aov(x~y)).
>
> La cuestión esque ahora tengo varios niveles del
factor especie y varios
> niveles del factor tratamiento.
> Si hago TukeyHSD(aov(abundancia~especie:tratamiento))
me salen todas las
> comparaciones posibles, pero yo sólo estoy interesado
en que para cada
> especie (i.e. para cada nivel del factor especie) me
salgan las
> comparaciones de los datos de abundancia para cada par
de tratamientos
> (niveles del factor tratamiento).
>
> Una solución laboriosa sería partir el documento
original en tantos como
> especies, pero seguro que alguien me puede sugerir una
solución más
> inteligente.
>
> Saludos y gracias,
>
>
> --
> **
>
> *Carlos Hernández-Castellano*
>
> Environmental Scientist
>
> Student, MSc Terrestrial Ecology and Biodiversity
Management
>
> Research Collaborator at Centre of Ecological Research
and Forestry
> Applications (CREAF)
>
> Ecology Unit - Department of Animal Biology, Plant
Biology and Ecology;
> Faculty of Sciences, Autonomous University of Barcelona
(UAB)
>
> 08193 Bellaterra, Spain
>
> Email: carlos.hernandezcastellano en
gmail.com<mailto:carlos.hernandezcastellano en gmail.com>
>
> [[alternative
HTML version deleted]]
>
> _______________________________________________
> R-help-es mailing list
> R-help-es en r-project.org<mailto:R-help-es en r-project.org>
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
[[alternative HTML version deleted]]
--
Pedro Concejero
BI & Big Data - Internal Exploitation - Telefónica
I+D<http://www.tid.es>
E-mail: pedro.concejerocerezo en telefonica.com<mailto:pedro.concejerocerezo
en telefonica.com>
skype: pedro.concejero
twitter @ConcejeroPedro<https://twitter.com/ConcejeroPedro>
linkedin pedroconcejero<http://www.linkedin.com/in/pedroconcejero/es>
Entusiasta R, me encontraréis aquí gRupo R madRid
<http://madrid.r-es.org/>
________________________________
Este mensaje y sus adjuntos se dirigen exclusivamente a su destinatario, puede
contener información privilegiada o confidencial y es para uso exclusivo de la
persona o entidad de destino. Si no es usted. el destinatario indicado, queda
notificado de que la lectura, utilización, divulgación y/o copia sin
autorización puede estar prohibida en virtud de la legislación vigente. Si ha
recibido este mensaje por error, le rogamos que nos lo comunique inmediatamente
por esta misma vía y proceda a su destrucción.
The information contained in this transmission is privileged and confidential
information intended only for the use of the individual or entity named above.
If the reader of this message is not the intended recipient, you are hereby
notified that any dissemination, distribution or copying of this communication
is strictly prohibited. If you have received this transmission in error, do not
read it. Please immediately reply to the sender that you have received this
communication in error and then delete it.
Esta mensagem e seus anexos se dirigem exclusivamente ao seu destinatário, pode
conter informação privilegiada ou confidencial e é para uso exclusivo da pessoa
ou entidade de destino. Se não é vossa senhoria o destinatário indicado, fica
notificado de que a leitura, utilização, divulgação e/ou cópia sem autorização
pode estar proibida em virtude da legislação vigente. Se recebeu esta mensagem
por erro, rogamos-lhe que nos o comunique imediatamente por esta mesma via e
proceda a sua destruição
[[alternative HTML version deleted]]