Recomiendo que veas el paquete lsmeans. Podes especificar comparaciones dentro factores.
Saludos Luciano El 29 de marzo de 2015, 11:36, Miguel Lázaro <lazar...@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-requ...@r-project.org < > r-help-es-requ...@r-project.org> escribió: > > Asunto: Resumen de R-help-es, Vol 73, Envío 42 > Para: r-help-es@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@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-requ...@r-project.org > > Puede contactar con el responsable de la lista escribiendo > a: > r-help-es-ow...@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 <pedroherreropeti...@gmail.com> > To: Ruben Tobalina Ramirez <lagrimaescr...@gmail.com> > Cc: Lista R <r-help-es@r-project.org> > Subject: Re: [R-es] Uso de R a través de Telegram > Message-ID: > > <CAM7H6U_EyWwvWKb4STnXWbNys_Z170zeN2+f6dqvz=me6m+...@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" <lagrimaescr...@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@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.hernandezcastell...@gmail.com> > To: r-help-es@r-project.org > Subject: [R-es] Comparaciones múltiples > Message-ID: > <cafqssgjvun3tssut8gyo0ggaktd8mymszld2stnbwfw6o9f...@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.hernandezcastell...@gmail.com > > [[alternative HTML version deleted]] > > > > ------------------------------ > > Message: 3 > Date: Sat, 28 Mar 2015 23:31:44 +0000 > From: Víctor Granda García <victorgrandagar...@gmail.com> > To: Carlos Hernández-Castellano > <carlos.hernandezcastell...@gmail.com>, > r-help-es@r-project.org > Subject: Re: [R-es] Comparaciones múltiples > Message-ID: > <caopgseatessjj0kju0y7dumtuomd8zr+a3bov3hy3tmppn4...@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.hernandezcastell...@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.hernandezcastell...@gmail.com > > > > [[alternative > HTML version deleted]] > > > > _______________________________________________ > > R-help-es mailing list > > R-help-es@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@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@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es > [[alternative HTML version deleted]] _______________________________________________ R-help-es mailing list R-help-es@r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es