Buenas tardes a todos. 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í. Os quería pedir que me dijerais el que está mal y porqué está mal. Los comandos que doy son estos:
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 <mailto:jbauti...@neiker.eus> | M. 688 62 98 14 www.neiker.eus<http://www.neiker.eus/> [cid:image001.png@01D9FB92.5EBA6090] <https://www.linkedin.com/company/neiker/> [cid:image002.png@01D9FB92.5EBA6090] <https://twitter.com/neiker_BRTA> [cid:image003.png@01D9FB92.5EBA6090] <https://www.facebook.com/neikerBRTA> [cid:image004.png@01D9FB92.5EBA6090] <https://www.youtube.com/user/neikertec> [cid:image005.png@01D9FB92.5EBA6090] [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-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 <mailto:jbauti...@neiker.eus> | M. 688 62 98 14 www.neiker.eus<http://www.neiker.eus/> [cid:image001.png@01D9FB92.5EBA6090] <https://www.linkedin.com/company/neiker/> [cid:image002.png@01D9FB92.5EBA6090] <https://twitter.com/neiker_BRTA> [cid:image003.png@01D9FB92.5EBA6090] <https://www.facebook.com/neikerBRTA> [cid:image004.png@01D9FB92.5EBA6090] <https://www.youtube.com/user/neikertec> [cid:image005.png@01D9FB92.5EBA6090] [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 <c...@qualityexcellence.es<mailto:c...@qualityexcellence.es>> Enviado el: lunes, 18 de septiembre de 2023 23:04 Para: Relloso Barrio, Juan Bautista <jbauti...@neiker.eus<mailto:jbauti...@neiker.eus>> CC: r-help-es@r-project.org<mailto: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<http://www.qualityexcellence.es> El lun, 18 sept 2023 a las 22:47, Relloso Barrio, Juan Bautista (<jbauti...@neiker.eus<mailto: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 <mailto:jbauti...@neiker.eus> | M. 688 62 98 14 www.neiker.eus<http://www.neiker.eus/> [cid:image001.png@01D9EA82.0D7BB740] <https://www.linkedin.com/company/neiker/> [cid:image002.png@01D9EA82.0D7BB740] <https://twitter.com/neiker_BRTA> [cid:image003.png@01D9EA82.0D7BB740] <https://www.facebook.com/neikerBRTA> [cid:image004.png@01D9EA82.0D7BB740] <https://www.youtube.com/user/neikertec> [cid:image005.png@01D9EA82.0D7BB740] [ej_ekonomia_garapen_lateral] [cid:image007.jpg@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@r-project.org<mailto:R-help-es@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@r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es
