[R-es] Hacer un Split plot

Mie Oct 11 13:36:29 CEST 2023

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
(<jbautista using 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
>
> jbautista using 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-bounces using r-project.org> *En nombre de *Relloso
> Barrio, Juan Bautista
> *Enviado el:* viernes, 22 de septiembre de 2023 16:18
> *Para:* r-help-es using 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
>
> jbautista using 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 <cof using qualityexcellence.es>
> *Enviado el:* lunes, 18 de septiembre de 2023 23:04
> *Para:* Relloso Barrio, Juan Bautista <jbautista using neiker.eus>
> *CC:* r-help-es using 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 (<
> jbautista using 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
>
> jbautista using 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 using 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 using r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>

------------ próxima parte ------------
Se ha borrado un adjunto en formato HTML...
URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20231011/5f6a482b/attachment-0001.html>

------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: image001.png
Type: image/png
Size: 1298 bytes
Desc: no disponible
URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20231011/5f6a482b/attachment-0006.png>

------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: image002.png
Type: image/png
Size: 1318 bytes
Desc: no disponible
URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20231011/5f6a482b/attachment-0007.png>

------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: image003.png
Type: image/png
Size: 1193 bytes
Desc: no disponible
URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20231011/5f6a482b/attachment-0008.png>

------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: image004.png
Type: image/png
Size: 1305 bytes
Desc: no disponible
URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20231011/5f6a482b/attachment-0009.png>

------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: image005.png
Type: image/png
Size: 21076 bytes
Desc: no disponible
URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20231011/5f6a482b/attachment-0010.png>

------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: image006.jpg
Type: image/jpeg
Size: 14174 bytes
Desc: no disponible
URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20231011/5f6a482b/attachment-0002.jpg>

------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: image007.jpg
Type: image/jpeg
Size: 5590 bytes
Desc: no disponible
URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20231011/5f6a482b/attachment-0003.jpg>

------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: image005.png
Type: image/png
Size: 21076 bytes
Desc: no disponible
URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20231011/5f6a482b/attachment-0011.png>