[R-es] Hacer un Split plot

Mie Oct 11 17:24:18 CEST 2023

Estimado Juan Bautista

¿A que se refiere con dos editores diferentes? ¿Editor de texto, de código R, o en código R escribir dos funciones diferentes?

Observando usted escribe model y luego model2.aov, antes de continuar, ¿que posibilidad hay que el algoritmo de cálculo sea diferente? Porque si internamente hay referencia a dos algoritmos diferentes, librerías o paquetes, el resultado puede variar. 

En lo personal me pasó recrear el ejemplo de un libro, y los resultados no me daban, consulté al autor, y la respuesta fue que en el libro había unos 5 datos, este tenía impreso todos los pasos con las matrices, numerito por numerito, en cambio el algoritmo estaba preparado para una gran cantidad de números, en otras palabras, al aumentar los cálculos los valores se aproximaban unos a otros, pero en una gran cantidad de cálculos la memoria o poder de procesamiento del procesador informático quedaría bloqueado, de ahí la diferencia entre el algoritmo del libro y el optimizado. ¿Que posibilidades hay que usted tenga alguna semejanza a esto?

Javier Rubén Marcuzzi

> El 11 oct 2023, a las 08:36, Manuel Mendoza <mmendoza using fulbrightmail.org> escribió:
> 
> 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  <mailto: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>
>>                   
>> 
>>  
>> 
>> 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 <mailto: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 <mailto: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  <mailto: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>
>>                   
>> 
>>  
>> 
>> 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 <mailto:cof using qualityexcellence.es>> 
>> Enviado el: lunes, 18 de septiembre de 2023 23:04
>> Para: Relloso Barrio, Juan Bautista <jbautista using neiker.eus <mailto:jbautista using neiker.eus>>
>> CC: r-help-es using r-project.org <mailto: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 <http://www.qualityexcellence.es/>
>>  
>> 
>> El lun, 18 sept 2023 a las 22:47, Relloso Barrio, Juan Bautista (<jbautista using neiker.eus <mailto: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  <mailto: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>
>>                   
>> 
>>  
>> 
>> 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 <mailto:R-help-es using r-project.org>
>> https://stat.ethz.ch/mailman/listinfo/r-help-es
>> 
>> 
>> 
>> --
>> 
>> Saludos,
>> Carlos Ortega
>> www.qualityexcellence.es <http://www.qualityexcellence.es/>_______________________________________________
>> R-help-es mailing list
>> R-help-es using r-project.org <mailto:R-help-es using r-project.org>
>> https://stat.ethz.ch/mailman/listinfo/r-help-es
> <image005.png><image003.png><image005.png><image001.png><image006.jpg><image002.png><image004.png><image007.jpg>_______________________________________________
> R-help-es mailing list
> R-help-es using r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es

	[[alternative HTML version deleted]]