[R-es] Comparaciones múltiples

Lun Mar 30 12:16:34 CEST 2015

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]]