[R-es] comparaciones múltiples

[Miguel Lázaro]

Hola a todos,
yo he tenido el mismo problema y después de hablar con mucha gente que a su vez habló con mucha gente, he optado por la solución laboriosa y parto la matriz general en tantas matrices como niveles tenga el factor -una vez la interacción es significativa. Es un procedimiento laborioso pero incontestable para cualquier revisor, creo. 
Saludos
Miguel
--------------------------------------------
El dom, 29/3/15, r-help-es-request en r-project.org <r-help-es-request en r-project.org> escribió:

 Asunto: Resumen de R-help-es, Vol 73, Envío 42
 Para: r-help-es en r-project.org
 Fecha: domingo, 29 de marzo, 2015 11:00

 Envíe los mensajes para la lista
 R-help-es a
     r-help-es en r-project.org

 Para subscribirse o anular su subscripción a través de la
 WEB
     https://stat.ethz.ch/mailman/listinfo/r-help-es

 O por correo electrónico, enviando un mensaje con el texto
 "help" en
 el asunto (subject) o en el cuerpo a:
     r-help-es-request en r-project.org

 Puede contactar con el responsable de la lista escribiendo
 a:
     r-help-es-owner en r-project.org

 Si responde a algún contenido de este mensaje, por favor,
 edite la
 linea del asunto (subject) para que el texto sea mas
 especifico que:
 "Re: Contents of R-help-es digest...". Además, por favor,
 incluya en
 la respuesta sólo aquellas partes del mensaje a las que
 está
 respondiendo.


 Asuntos del día:

    1. Re: Uso de R a través de Telegram
 (Pedro Herrero Petisco)
    2. Comparaciones múltiples (Carlos
 Hernández-Castellano)
    3. Re: Comparaciones múltiples (Víctor
 Granda García)


 ----------------------------------------------------------------------

 Message: 1
 Date: Sat, 28 Mar 2015 12:06:35 +0100
 From: Pedro Herrero Petisco <pedroherreropetisco en gmail.com>
 To: Ruben Tobalina Ramirez <lagrimaescrita en gmail.com>
 Cc: Lista R <r-help-es en r-project.org>
 Subject: Re: [R-es] Uso de R a través de Telegram
 Message-ID:
    
 <CAM7H6U_EyWwvWKb4STnXWbNys_Z170zeN2+f6dqvz=ME6m+qmg en mail.gmail.com>
 Content-Type: text/plain; charset="UTF-8"

 Muchas gracias por compartir, la verdad es que está muy
 chulo :-)
 El 28/03/2015 09:29, "Ruben Tobalina Ramirez" <lagrimaescrita en gmail.com>
 escribió:

 > Buenas,
 >
 > He descubierto y queria compartirla con vosotros. La
 verdad es que
 > tiene buena pinta, es curiosa, pero no sé si tiene
 mucha utilidad.
 > Tele R es una "app" para usar R a través de Telegram.
 Básicamente
 > envias el código de R a una cuenta de Telegram que
 esta conectada a un
 > servidor en la nube con R y te envia a tu telefono los
 resultados,
 > vamos como si tuvieses R en tu móvil. Os envio su werb
 por si os
 > interesa: http://telemath.altervista.org/TeleR.html
 >
 > y existe lo mismo para Octave:
 >
 > http://nbviewer.ipython.org/github/CAChemE/lightning-talks/blob/master/2015-02/TeleOctave.ipynb
 >
 > Un saludo
 >
 > Rubén
 >
 > _______________________________________________
 > R-help-es mailing list
 > R-help-es en r-project.org
 > https://stat.ethz.ch/mailman/listinfo/r-help-es
 >

     [[alternative HTML version deleted]]



 ------------------------------

 Message: 2
 Date: Sat, 28 Mar 2015 23:10:04 +0100
 From: Carlos Hernández-Castellano
     <carlos.hernandezcastellano en gmail.com>
 To: r-help-es en r-project.org
 Subject: [R-es] Comparaciones múltiples
 Message-ID:
     <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

     [[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>
 To: Carlos Hernández-Castellano
     <carlos.hernandezcastellano en gmail.com>,
 r-help-es en r-project.org
 Subject: Re: [R-es] Comparaciones múltiples
 Message-ID:
     <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>)
 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
 >
 >         [[alternative
 HTML version deleted]]
 >
 > _______________________________________________
 > R-help-es mailing list
 > R-help-es en r-project.org
 > https://stat.ethz.ch/mailman/listinfo/r-help-es
 >

     [[alternative HTML version deleted]]



 ------------------------------

 Subject: Pié de página del digest

 _______________________________________________
 R-help-es mailing list
 R-help-es en r-project.org
 https://stat.ethz.ch/mailman/listinfo/r-help-es


 ------------------------------

 Fin de Resumen de R-help-es, Vol 73, Envío 42