Por si te fuera de utilidad, esta es la respuesta de chat GPT 4.0 Parece que tienes un par de problemas y preguntas aquí. Intentaré ayudarte paso a paso.
### Problema 1: Error en la Separación de Medias HSD por Tratamiento Sobre el código que estás utilizando para la separación de medias, el error específico que estás recibiendo `argument "trt" is missing, with no default` posiblemente se relaciona con la función `HSD.test`. Vamos a revisar el bloque de código que has compartido: ```r by(pppp, pppp$TRAT, function(x) { anova_result <- aov(RTOSE ~ CORTE, data = x) HSD_result <- HSD.test(anova_result) Tratamiento <- unique(x$TRAT) return(list(Tratamiento = TRAT, HSD_result = HSD_result)) }) ``` Si estás utilizando la función `HSD.test` del paquete `agricolae`, asegúrate de proporcionar todos los argumentos necesarios para esta función, que serían el objeto de anova, el nombre del tratamiento y el nivel de significancia. Por ejemplo: ```r HSD_result <- HSD.test(anova_result, "CORTE", group=TRUE, console=TRUE) ``` Además, asegúrate de que la variable `TRAT` está disponible en el scope de tu función; en tu código, pareces utilizar `Tratamiento` para almacenar el valor único de `x$TRAT`, pero intentas devolver una variable llamada `TRAT` que no está definida. Tal vez quisiste hacer esto: ```r return(list(Tratamiento = Tratamiento, HSD_result = HSD_result)) ``` ### Problema 2: Comparación de Dos Métodos para ANOVA de Split Plot Tienes dos bloques de código para ANOVA: 1. `model <- with(CEPAS, sp.plot(REP, CEPA, VARIEDAD, trAUDPC))` 2. `model2.aov <- aov(trAUDPC ~ CEPA*VARIEDAD + REP + Error(REP/CEPA))` Dado que no tengo toda la información sobre las funciones y los datos, es un poco difícil decir con certeza si ambos realizan un ANOVA de split plot de la misma manera. - La función `sp.plot` no es estándar en R y puede ser parte de un paquete específico o un script personalizado que no estoy familiarizado con. Si proporcionas más información sobre esta función, podría ayudarte mejor. - La segunda línea de código parece utilizar la función `aov` de una manera que podría ser coherente con un diseño de split plot, especificando un término de error estratificado `Error(REP/CEPA)`, pero sin tener acceso a los datos y más contexto, es un poco difícil de verificar. Si los resultados que obtienes son diferentes, aquí hay algunas cosas que querrás revisar: - **Datos**: Asegúrate de que los datos que se están utilizando son exactamente los mismos en ambos análisis. - **Modelo**: Verifica que los modelos que estás ajustando sean equivalentes y apropiados para tu diseño experimental. - **Hipótesis**: Confirma que las hipótesis que se están probando sean las mismas. Recuerda que proporcionar un ejemplo reproducible (datos de ejemplo y código) es la mejor manera de obtener ayuda específica y precisa. Si puedes compartir más detalles, estaré encantado de ayudarte más! El mié, 11 oct 2023 a las 10:32, Relloso Barrio, Juan Bautista (<jbauti...@neiker.eus>) escribió: > Buenos días. > > Perdón por insistir. Quiero hacer un ANOVA para unos datos que tienen un > diseño Split plot. > > Y estoy utilizando dos editores diferentes, que, si todo fuera correcto y > ya que toman los mismos datos, me debería dar el mismo resultado, … pero no > es así. > > Querría saber si básicamente los dos hacen un análisis de varianza SPLIT > PLOT ?. Y el mismo análisis? > > Los comandos que doy son estos: > > CEPAS, es el nombre del archivo de datos. > > model<-with(CEPAS, sp.plot(REP, CEPA, VARIEDAD, trAUDPC)) > > Y el otro es este > > model2.aov<- aov(trAUDPC~CEPA*VARIEDAD+REP+Error(REP/CEPA)) > > > > Un saludo. > > *Juan Bautista Relloso Barrio* > > Coordinador de Equipos e Infraestructuras | Técnico de Cultivos del > Departamento de Producción Vegetal > > Talde eta Azpiegitura Koordinatzailea | Laborantza teknikaria Landare > Ekoizpen Departamentua > > jbauti...@neiker.eus | M. 688 62 98 14 > > *www.neiker.eus* <http://www.neiker.eus/> > <https://www.linkedin.com/company/neiker/> > <https://twitter.com/neiker_BRTA> <https://www.facebook.com/neikerBRTA> > <https://www.youtube.com/user/neikertec> > > *[image: ej_ekonomia_garapen_lateral]* [image: > https://www.emakunde.euskadi.eus/contenidos/informacion/violencia_puntolila/eu_def/adjuntos/punto.lila.jpg] > > > > PRIBATUTASUN POLITIKA <https://neiker.eus/eu/pribatutasun-politika/> | > POLÍTICA > DE PRIVACIDAD <https://neiker.eus/es/politica-privacidad/> | LEGAL NOTICE > <https://neiker.eus/en/privacy-policy/> > > > > > > > > *De:* R-help-es <r-help-es-boun...@r-project.org> *En nombre de *Relloso > Barrio, Juan Bautista > *Enviado el:* viernes, 22 de septiembre de 2023 16:18 > *Para:* r-help-es@r-project.org > *Asunto:* [R-es] no encuentro el error > > > > Buenas tardes. > > Intento hacer la separación de medias “por tratamientos” > > Pero me da un error que no consigo corregir. No sé que hago mal. Os paso > los comandos y el error > > > > # Separación de medias HSD por Tratamiento No funciona > > by(pppp, pppp$TRAT, function(x) { > > anova_result <- aov(RTOSE ~ CORTE, data = x) > > HSD_result <- HSD.test(anova_result) > > Tratamiento <- unique(x$TRAT) > > return(list(Tratamiento = TRAT, HSD_result = HSD_result)) > > }) > > Y el error es este: > > Error in pmatch(trt, names(A)) : > > argument "trt" is missing, with no default > > > > Un saludo. > > *Juan Bautista Relloso Barrio* > > Coordinador de Equipos e Infraestructuras | Técnico de Cultivos del > Departamento de Producción Vegetal > > Talde eta Azpiegitura Koordinatzailea | Laborantza teknikaria Landare > Ekoizpen Departamentua > > jbauti...@neiker.eus | M. 688 62 98 14 > > *www.neiker.eus* <http://www.neiker.eus/> > <https://www.linkedin.com/company/neiker/> > <https://twitter.com/neiker_BRTA> <https://www.facebook.com/neikerBRTA> > <https://www.youtube.com/user/neikertec> > > *[image: ej_ekonomia_garapen_lateral]* [image: > https://www.emakunde.euskadi.eus/contenidos/informacion/violencia_puntolila/eu_def/adjuntos/punto.lila.jpg] > > > > PRIBATUTASUN POLITIKA <https://neiker.eus/eu/pribatutasun-politika/> | > POLÍTICA > DE PRIVACIDAD <https://neiker.eus/es/politica-privacidad/> | LEGAL NOTICE > <https://neiker.eus/en/privacy-policy/> > > > > > > > > *De:* Carlos Ortega <c...@qualityexcellence.es> > *Enviado el:* lunes, 18 de septiembre de 2023 23:04 > *Para:* Relloso Barrio, Juan Bautista <jbauti...@neiker.eus> > *CC:* r-help-es@r-project.org > *Asunto:* Re: [R-es] Boxplot ordenados... > > > > Hola Juan, > > > > Sí, te recomiendo que uses la librería "ggeasy" que simplifica mucho estas > modificaciones, mejor que trabajar directamente con "ggplot". > > Sería así: > > > > #------------- > > *library(ggeasy)* > > > > ggplot(data = AÑOS, aes(x = VARIEDAD, y = rAUDPCor, colour = ANO, > > group = ANO)) + > > stat_summary(fun = mean, geom = "point") + > > stat_summary(fun = mean, geom = "line") + > > labs(y = 'mean rAUDPCor') + > > labs(x = 'Variety') + > > theme_bw() + > > *easy_rotate_x_labels( angle = 90)* > > > > #------------ > > > > > > Gracias, > > Carlos Ortega > > www.qualityexcellence.es > > > > El lun, 18 sept 2023 a las 22:47, Relloso Barrio, Juan Bautista (< > jbauti...@neiker.eus>) escribió: > > Buenas noches. > > Quisiera saber si alguien me puede decir como en estos comandos para un > gráfico … > > > > > > ggplot(data = AÑOS, aes(x = VARIEDAD, y = rAUDPCor, colour = ANO, > > group = ANO)) + > > stat_summary(fun = mean, geom = "point") + > > stat_summary(fun = mean, geom = "line") + > > labs(y = 'mean rAUDPCor') + > > labs(x = 'Variety') + > > theme_bw() > > > > Qué debo de añadir para que la orientación de las los nombres de las > variedades del eje “X” me aparezcan en vertical en vez de en horizontal. > > Muchas gracias. > > Un saludo. > > *Juan Bautista Relloso Barrio* > > Coordinador de Equipos e Infraestructuras | Técnico de Cultivos del > Departamento de Producción Vegetal > > Talde eta Azpiegitura Koordinatzailea | Laborantza teknikaria Landare > Ekoizpen Departamentua > > jbauti...@neiker.eus | M. 688 62 98 14 > > *www.neiker.eus* <http://www.neiker.eus/> > <https://www.linkedin.com/company/neiker/> > <https://twitter.com/neiker_BRTA> <https://www.facebook.com/neikerBRTA> > <https://www.youtube.com/user/neikertec> > > *[image: ej_ekonomia_garapen_lateral]* > > > > PRIBATUTASUN POLITIKA <https://neiker.eus/eu/pribatutasun-politika/> | > POLÍTICA > DE PRIVACIDAD <https://neiker.eus/es/politica-privacidad/> | LEGAL NOTICE > <https://neiker.eus/en/privacy-policy/> > > > > _______________________________________________ > R-help-es mailing list > R-help-es@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es > > > > -- > > Saludos, > Carlos Ortega > www.qualityexcellence.es > _______________________________________________ > R-help-es mailing list > R-help-es@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es >
_______________________________________________ R-help-es mailing list R-help-es@r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es