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 en neiker.eus <mailto:jbautista en neiker.eus> | M. 688 62 98 14
www.neiker.eus<http://www.neiker.eus/> [cid:image001.png en
01D9FC2D.E9B85580] <https://www.linkedin.com/company/neiker/>
[cid:image002.png en 01D9FC2D.E9B85580] <https://twitter.com/neiker_BRTA>
[cid:image003.png en 01D9FC2D.E9B85580]
<https://www.facebook.com/neikerBRTA> [cid:image004.png en
01D9FC2D.E9B85580] <https://www.youtube.com/user/neikertec>
[cid:image005.png en 01D9FC2D.E9B85580] [ej_ekonomia_garapen_lateral]
[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 en r-project.org<mailto:r-help-es-bounces
en r-project.org>> En nombre de Relloso Barrio, Juan Bautista
Enviado el: viernes, 22 de septiembre de 2023 16:18
Para: r-help-es en r-project.org<mailto:r-help-es en 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 en neiker.eus <mailto:jbautista en neiker.eus> | M. 688 62 98 14
www.neiker.eus<http://www.neiker.eus/> [cid:image001.png en
01D9FC2D.E9B85580] <https://www.linkedin.com/company/neiker/>
[cid:image002.png en 01D9FC2D.E9B85580] <https://twitter.com/neiker_BRTA>
[cid:image003.png en 01D9FC2D.E9B85580]
<https://www.facebook.com/neikerBRTA> [cid:image004.png en
01D9FC2D.E9B85580] <https://www.youtube.com/user/neikertec>
[cid:image005.png en 01D9FC2D.E9B85580] [ej_ekonomia_garapen_lateral]
[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 en qualityexcellence.es<mailto:cof en
qualityexcellence.es>>
Enviado el: lunes, 18 de septiembre de 2023 23:04
Para: Relloso Barrio, Juan Bautista <jbautista en
neiker.eus<mailto:jbautista en neiker.eus>>
CC: r-help-es en r-project.org<mailto:r-help-es en 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
en neiker.eus<mailto:jbautista en 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 en neiker.eus <mailto:jbautista en neiker.eus> | M. 688 62 98 14
www.neiker.eus<http://www.neiker.eus/> [cid:image001.png en
01D9EA82.0D7BB740] <https://www.linkedin.com/company/neiker/>
[cid:image002.png en 01D9EA82.0D7BB740] <https://twitter.com/neiker_BRTA>
[cid:image003.png en 01D9EA82.0D7BB740]
<https://www.facebook.com/neikerBRTA> [cid:image004.png en
01D9EA82.0D7BB740] <https://www.youtube.com/user/neikertec>
[cid:image005.png en 01D9EA82.0D7BB740] [ej_ekonomia_garapen_lateral]
[cid:image007.jpg en 01D9EA82.0D7BB740]
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 en r-project.org<mailto:R-help-es en r-project.org>
https://stat.ethz.ch/mailman/listinfo/r-help-es
--
Saludos,
Carlos Ortega
www.qualityexcellence.es<http://www.qualityexcellence.es>
------------ próxima parte ------------
Se ha borrado un adjunto en formato HTML...
URL:
<https://stat.ethz.ch/pipermail/r-help-es/attachments/20231011/c4a1f08c/attachment-0001.html>
------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: image001.png
Type: image/png
Size: 1298 bytes
Desc: image001.png
URL:
<https://stat.ethz.ch/pipermail/r-help-es/attachments/20231011/c4a1f08c/attachment-0005.png>
------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: image002.png
Type: image/png
Size: 1318 bytes
Desc: image002.png
URL:
<https://stat.ethz.ch/pipermail/r-help-es/attachments/20231011/c4a1f08c/attachment-0006.png>
------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: image003.png
Type: image/png
Size: 1193 bytes
Desc: image003.png
URL:
<https://stat.ethz.ch/pipermail/r-help-es/attachments/20231011/c4a1f08c/attachment-0007.png>
------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: image004.png
Type: image/png
Size: 1305 bytes
Desc: image004.png
URL:
<https://stat.ethz.ch/pipermail/r-help-es/attachments/20231011/c4a1f08c/attachment-0008.png>
------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: image005.png
Type: image/png
Size: 21076 bytes
Desc: image005.png
URL:
<https://stat.ethz.ch/pipermail/r-help-es/attachments/20231011/c4a1f08c/attachment-0009.png>
------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: image006.jpg
Type: image/jpeg
Size: 14174 bytes
Desc: image006.jpg
URL:
<https://stat.ethz.ch/pipermail/r-help-es/attachments/20231011/c4a1f08c/attachment-0002.jpg>
------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: image007.jpg
Type: image/jpeg
Size: 5590 bytes
Desc: image007.jpg
URL:
<https://stat.ethz.ch/pipermail/r-help-es/attachments/20231011/c4a1f08c/attachment-0003.jpg>
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 en 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 en 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 en r-project.org> *En nombre de
*Relloso
> Barrio, Juan Bautista
> *Enviado el:* viernes, 22 de septiembre de 2023 16:18
> *Para:* r-help-es en 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 en 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 en qualityexcellence.es>
> *Enviado el:* lunes, 18 de septiembre de 2023 23:04
> *Para:* Relloso Barrio, Juan Bautista <jbautista en neiker.eus>
> *CC:* r-help-es en 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 en 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 en 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 en 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 en 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>
Relloso Barrio, Juan Bautista
2023-Nov-21 07:46 UTC
[R-es] Cambiar el intervalo de confianza en un anova
Buenos días.
Normalmente por defecto, todos los ANOVAS (análisis de varianza) que he
realizado han sido con un intervalo de confianza del 95%.
¿Cuál sería la orden que debería introducir si deseo modificar ese 95% y, por
ejemplo, bajarlo a un 90%?
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 en neiker.eus <mailto:jbautista en neiker.eus> | M. 688 62 98 14
www.neiker.eus<http://www.neiker.eus/> [cid:image001.png en
01DA1C57.393B21C0] <https://www.linkedin.com/company/neiker/>
[cid:image002.png en 01DA1C57.393B21C0] <https://twitter.com/neiker_BRTA>
[cid:image003.png en 01DA1C57.393B21C0]
<https://www.facebook.com/neikerBRTA> [cid:image004.png en
01DA1C57.393B21C0] <https://www.youtube.com/user/neikertec>
[cid:image005.png en 01DA1C57.393B21C0] [ej_ekonomia_garapen_lateral]
[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 en r-project.org> En nombre de Relloso
Barrio, Juan Bautista
Enviado el: miércoles, 11 de octubre de 2023 10:32
Para: r-help-es en r-project.org
Asunto: Re: [R-es] Hacer un Split plot
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 en neiker.eus <mailto:jbautista en neiker.eus> | M. 688 62 98 14
www.neiker.eus<http://www.neiker.eus/> [cid:image001.png en
01DA1C57.393B21C0] <https://www.linkedin.com/company/neiker/>
[cid:image002.png en 01DA1C57.393B21C0] <https://twitter.com/neiker_BRTA>
[cid:image003.png en 01DA1C57.393B21C0]
<https://www.facebook.com/neikerBRTA> [cid:image004.png en
01DA1C57.393B21C0] <https://www.youtube.com/user/neikertec>
[cid:image005.png en 01DA1C57.393B21C0] [ej_ekonomia_garapen_lateral]
[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 en r-project.org<mailto:r-help-es-bounces
en r-project.org>> En nombre de Relloso Barrio, Juan Bautista
Enviado el: viernes, 22 de septiembre de 2023 16:18
Para: r-help-es en r-project.org<mailto:r-help-es en 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 en neiker.eus <mailto:jbautista en neiker.eus> | M. 688 62 98 14
www.neiker.eus<http://www.neiker.eus/> [cid:image001.png en
01DA1C57.393B21C0] <https://www.linkedin.com/company/neiker/>
[cid:image002.png en 01DA1C57.393B21C0] <https://twitter.com/neiker_BRTA>
[cid:image003.png en 01DA1C57.393B21C0]
<https://www.facebook.com/neikerBRTA> [cid:image004.png en
01DA1C57.393B21C0] <https://www.youtube.com/user/neikertec>
[cid:image005.png en 01DA1C57.393B21C0] [ej_ekonomia_garapen_lateral]
[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 en qualityexcellence.es<mailto:cof en
qualityexcellence.es>>
Enviado el: lunes, 18 de septiembre de 2023 23:04
Para: Relloso Barrio, Juan Bautista <jbautista en
neiker.eus<mailto:jbautista en neiker.eus>>
CC: r-help-es en r-project.org<mailto:r-help-es en 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
en neiker.eus<mailto:jbautista en 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 en neiker.eus <mailto:jbautista en neiker.eus> | M. 688 62 98 14
www.neiker.eus<http://www.neiker.eus/> [cid:image001.png en
01D9EA82.0D7BB740] <https://www.linkedin.com/company/neiker/>
[cid:image002.png en 01D9EA82.0D7BB740] <https://twitter.com/neiker_BRTA>
[cid:image003.png en 01D9EA82.0D7BB740]
<https://www.facebook.com/neikerBRTA> [cid:image004.png en
01D9EA82.0D7BB740] <https://www.youtube.com/user/neikertec>
[cid:image005.png en 01D9EA82.0D7BB740] [ej_ekonomia_garapen_lateral]
[cid:image007.jpg en 01D9EA82.0D7BB740]
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 en r-project.org<mailto:R-help-es en r-project.org>
https://stat.ethz.ch/mailman/listinfo/r-help-es
--
Saludos,
Carlos Ortega
www.qualityexcellence.es<http://www.qualityexcellence.es>
------------ próxima parte ------------
Se ha borrado un adjunto en formato HTML...
URL:
<https://stat.ethz.ch/pipermail/r-help-es/attachments/20231121/614f089d/attachment-0001.html>
------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: image001.png
Type: image/png
Size: 1298 bytes
Desc: image001.png
URL:
<https://stat.ethz.ch/pipermail/r-help-es/attachments/20231121/614f089d/attachment-0005.png>
------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: image002.png
Type: image/png
Size: 1318 bytes
Desc: image002.png
URL:
<https://stat.ethz.ch/pipermail/r-help-es/attachments/20231121/614f089d/attachment-0006.png>
------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: image003.png
Type: image/png
Size: 1193 bytes
Desc: image003.png
URL:
<https://stat.ethz.ch/pipermail/r-help-es/attachments/20231121/614f089d/attachment-0007.png>
------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: image004.png
Type: image/png
Size: 1305 bytes
Desc: image004.png
URL:
<https://stat.ethz.ch/pipermail/r-help-es/attachments/20231121/614f089d/attachment-0008.png>
------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: image005.png
Type: image/png
Size: 21076 bytes
Desc: image005.png
URL:
<https://stat.ethz.ch/pipermail/r-help-es/attachments/20231121/614f089d/attachment-0009.png>
------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: image006.jpg
Type: image/jpeg
Size: 14174 bytes
Desc: image006.jpg
URL:
<https://stat.ethz.ch/pipermail/r-help-es/attachments/20231121/614f089d/attachment-0002.jpg>
------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: image007.jpg
Type: image/jpeg
Size: 5590 bytes
Desc: image007.jpg
URL:
<https://stat.ethz.ch/pipermail/r-help-es/attachments/20231121/614f089d/attachment-0003.jpg>