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, [email protected]<mailto:[email protected]> escribi�: Message: 2 Date: Sat, 28 Mar 2015 23:10:04 +0100 From: Carlos Hern�ndez-Castellano <[email protected]><mailto:[email protected]> To: [email protected]<mailto:[email protected]> Subject: [R-es] Comparaciones m�ltiples Message-ID: <cafqssgjvun3tssut8gyo0ggaktd8mymszld2stnbwfw6o9f...@mail.gmail.com><mailto: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: [email protected]<mailto:[email protected]> [[alternative HTML version deleted]] ------------------------------ Message: 3 Date: Sat, 28 Mar 2015 23:31:44 +0000 From: V�ctor Granda Garc�a <[email protected]><mailto:[email protected]> To: Carlos Hern�ndez-Castellano <[email protected]><mailto:[email protected]>, [email protected]<mailto:[email protected]> Subject: Re: [R-es] Comparaciones m�ltiples Message-ID: <caopgseatessjj0kju0y7dumtuomd8zr+a3bov3hy3tmppn4...@mail.gmail.com><mailto: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 (< [email protected]<mailto:[email protected]>>) 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: > [email protected]<mailto:[email protected]> > > [[alternative HTML version deleted]] > > _______________________________________________ > R-help-es mailing list > [email protected]<mailto:[email protected]> > 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: [email protected]<mailto:[email protected]> 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]]
_______________________________________________ R-help-es mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/r-help-es
