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]]
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]]
Muchas gracias! 2015-03-29 0:31 GMT+01:00 Víctor Granda García <victorgrandagarcia en gmail.com>:> 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 >> >-- *?* *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]]