Hola a tod en s,
Tengo tres variables medidas sobre 300 individuos, con un factor
principal, Genotipo. Lo que acostumbraba hacer eran tres anovas
independientes para cada una de las variables y despues hacer un test
de comparacion multiple (LSD,Tukey,..). Mi interes reacae en hacer lo
mismo pero con el procediento manova, el interes es del Referee
basicamente. Lo que he hecho es hacer el manova:
manova(cbind(A1, A2, A3) ~ Genotipo, data = a_Amil)
Terms:
Genotip Residuals
resp 1 18.270 4512.025
resp 2 3.591 2117.572
resp 3 652.406 7632.530
Deg. of Freedom 2 665
Residual standard error: 2.604803 1.784466 3.387844
testar que hay dif.Significativas y con el residual estandard del
error hacer el LSD. Mis dudas son,
- Creo que no esta bien echo el procedimento, por que si hago el anova
individual para cada variable los MSE son los mismos.
- Hay algun procedimento para hacer comparacion multiple de medianas
utilizando manova en R?
- Hay alguna forma para hacer una representacion grafica de los
resultados, he probado T2Hotteling, si mucho exito.
- Puedo de alguna forma argumentar que el test anova es
suficientemente valido para este proposito?
Muchas gracias por adelantado.
Marcal
Hola, Mira el ejemplo que aparece en "Handbook of Statistical Analysis using R" [1] en la sección "5.3.3 Water Hardness and Mortality". [1]: http://www.amazon.es/Handbook-Statistical-Analyses-Using/dp/1420079336/ref=sr_1_fkmr2_1?ie=UTF8&qid=1345457394&sr=8-1-fkmr2 Saludos, Carlos Ortega www.qualityexcellence.es El 20 de agosto de 2012 01:29, Marçal Plans <marcal.plans@upc.edu> escribió:> Hola a tod@s, > > Tengo tres variables medidas sobre 300 individuos, con un factor > principal, Genotipo. Lo que acostumbraba hacer eran tres anovas > independientes para cada una de las variables y despues hacer un test de > comparacion multiple (LSD,Tukey,..). Mi interes reacae en hacer lo mismo > pero con el procediento manova, el interes es del Referee basicamente. Lo > que he hecho es hacer el manova: > > manova(cbind(A1, A2, A3) ~ Genotipo, data = a_Amil) > > Terms: > Genotip Residuals > resp 1 18.270 4512.025 > resp 2 3.591 2117.572 > resp 3 652.406 7632.530 > Deg. of Freedom 2 665 > > Residual standard error: 2.604803 1.784466 3.387844 > > testar que hay dif.Significativas y con el residual estandard del error > hacer el LSD. Mis dudas son, > > - Creo que no esta bien echo el procedimento, por que si hago el anova > individual para cada variable los MSE son los mismos. > > - Hay algun procedimento para hacer comparacion multiple de medianas > utilizando manova en R? > > - Hay alguna forma para hacer una representacion grafica de los > resultados, he probado T2Hotteling, si mucho exito. > > - Puedo de alguna forma argumentar que el test anova es suficientemente > valido para este proposito? > > Muchas gracias por adelantado. > > Marcal > > ______________________________**_________________ > R-help-es mailing list > R-help-es@r-project.org > https://stat.ethz.ch/mailman/**listinfo/r-help-es<https://stat.ethz.ch/mailman/listinfo/r-help-es> >-- Saludos, Carlos Ortega www.qualityexcellence.es [[alternative HTML version deleted]]
Aunque creo que corresponde a la primera edición, el capítulo está disponible libremente aquí: http://cran.r-project.org/web/packages/HSAUR/vignettes/Ch_analysis_of_variance.pdf e incluye el detalle al que me refería en el correo anterior. Saludos, Carlos Ortega www.qualityexcellence.es El 20 de agosto de 2012 12:11, Carlos Ortega <cof@qualityexcellence.es>escribió:> Hola, > > Mira el ejemplo que aparece en "Handbook of Statistical Analysis using R" > [1] en la sección "5.3.3 Water Hardness and Mortality". > > [1]: > http://www.amazon.es/Handbook-Statistical-Analyses-Using/dp/1420079336/ref=sr_1_fkmr2_1?ie=UTF8&qid=1345457394&sr=8-1-fkmr2 > > Saludos, > Carlos Ortega > www.qualityexcellence.es > > El 20 de agosto de 2012 01:29, Marçal Plans <marcal.plans@upc.edu>escribió: > > Hola a tod@s, >> >> Tengo tres variables medidas sobre 300 individuos, con un factor >> principal, Genotipo. Lo que acostumbraba hacer eran tres anovas >> independientes para cada una de las variables y despues hacer un test de >> comparacion multiple (LSD,Tukey,..). Mi interes reacae en hacer lo mismo >> pero con el procediento manova, el interes es del Referee basicamente. Lo >> que he hecho es hacer el manova: >> >> manova(cbind(A1, A2, A3) ~ Genotipo, data = a_Amil) >> >> Terms: >> Genotip Residuals >> resp 1 18.270 4512.025 >> resp 2 3.591 2117.572 >> resp 3 652.406 7632.530 >> Deg. of Freedom 2 665 >> >> Residual standard error: 2.604803 1.784466 3.387844 >> >> testar que hay dif.Significativas y con el residual estandard del error >> hacer el LSD. Mis dudas son, >> >> - Creo que no esta bien echo el procedimento, por que si hago el anova >> individual para cada variable los MSE son los mismos. >> >> - Hay algun procedimento para hacer comparacion multiple de medianas >> utilizando manova en R? >> >> - Hay alguna forma para hacer una representacion grafica de los >> resultados, he probado T2Hotteling, si mucho exito. >> >> - Puedo de alguna forma argumentar que el test anova es suficientemente >> valido para este proposito? >> >> Muchas gracias por adelantado. >> >> Marcal >> >> ______________________________**_________________ >> R-help-es mailing list >> R-help-es@r-project.org >> https://stat.ethz.ch/mailman/**listinfo/r-help-es<https://stat.ethz.ch/mailman/listinfo/r-help-es> >> > > > > -- > Saludos, > Carlos Ortega > www.qualityexcellence.es >-- Saludos, Carlos Ortega www.qualityexcellence.es [[alternative HTML version deleted]]