Recomiendo que veas el paquete lsmeans. Podes especificar comparaciones
dentro factores.
Saludos
Luciano
El 29 de marzo de 2015, 11:36, Miguel Lázaro <lazarolv en yahoo.es>
escribió:
> 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
>
> _______________________________________________
> 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]]