Skip to content
Prev 12346 / 15379 Next

[R-es] Incongruencias entre ANOVA y comparación por pares

Antes de hablar de incongruencia habria que especificar si es de la interpretacion o del algoritmo dr R. Para descartar lo ?ltimo compara si puedes con otro softaware, o intenta, si no es muy complicado montar el algoritmo co formulas excel



Enviado desde mi smartphone Samsung Galaxy.


-------- Mensaje original --------
De: Yesica Pallavicini Fernandez <yesipalla en gmail.com>
Fecha: 20/2/18 13:48 (GMT+01:00)
Para: r-help-es en r-project.org
Asunto: [R-es] Incongruencias entre ANOVA y comparaci?n por pares

Hola, he hecho una anova con su posterior comparaci?n de Tukey pero me da
un resultado incongruente.
Mientras la ANOVA me dice que SI hay diferencias significativas entre
variedades ( resaltado en amarillos), la comparaci?n me da todo en la misma
letra ( lo que significa que no hay diferencias).
Yo esperaba al menos que  1 variedad sea diferente.

Ya he revisado los datos y est?n bien.
He realizado los mismos an?lisis con otros datos y las variedades me dan
diferentes entre si.

?A que se puede deber esta incongruencia?
?Se puede confiar en el P-valor de este an?lisis?
 ?Me recomiendan alguna otra funcion/librer?a para este an?lisis)

Gracias

summary(Za<-aov(yield~variety*treat*Y+rep,z))

                 Df Sum Sq Mean Sq  F value   Pr(>F)
variety          23   30.7     1.3    2.396 0.000680 ***
treat             1  753.8   753.8 1354.597  < 2e-16 ***
Y                 1  392.6   392.6  705.634  < 2e-16 ***
rep               1    6.1     6.1   10.894 0.001151 **
variety:treat    23   31.1     1.4    2.432 0.000553 ***
variety:Y        23   26.9     1.2    2.101 0.003570 **
treat:Y           1   14.8    14.8   26.608 6.21e-07 ***
variety:treat:Y  23   13.2     0.6    1.030 0.429382
Residuals       191  106.3     0.6
---
Signif. codes:
0 ?***? 0.001 ?**? 0.01 ?*? 0.05 ?.? 0.1 ? ? 1

out<-HSD.test(Za,"variety", group=TRUE,console=TRUE)

               yield groups
Avispa      7.447222      a
Pelayo      7.446759      a
Regallo     7.423611      a
Olivadur    7.395833      a
Martinur    7.345833      a
Iride       7.315278      a
Amilcar     7.303241      a
Gallareta   7.223611      a
Ramirez     7.127315      a
Sula        6.983796      a
Sculptur    6.983449      a
Claudio     6.971759      a
Simeto      6.910648      a
Bolo        6.907407      a
Vitron      6.896759      a
Bolido      6.878704      a
Mexa        6.804630      a
D.Pedro     6.779630      a
Dorondon    6.724537      a
D.Ricardo   6.656019      a
Burgos      6.582870      a
Tussur      6.529514      a
D.Sebastian 6.412500      a
Kiko.Nick   6.374074      a


_______________________________________________
R-help-es mailing list
R-help-es en r-project.org
https://nam01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fstat.ethz.ch%2Fmailman%2Flistinfo%2Fr-help-es&data=02%7C01%7C%7Cf82210b3ff3844fce81808d578604957%7C84df9e7fe9f640afb435aaaaaaaaaaaa%7C1%7C0%7C636547277298744079&sdata=AQrzzbSl1AQ42lBMJrpw74W59CDXEr%2FOkdgCqFV6V5k%3D&reserved=0
Message-ID: <VI1PR05MB1213E80ED36A54C55570A8D3CDCF0@VI1PR05MB1213.eurprd05.prod.outlook.com>
In-Reply-To: <CAD11z9xyhhjPh9PJMUQB2NOxKZ34B1G0YOrKA105TiXEOsHYgg@mail.gmail.com>