Buenos días y gracias de antemano por vuestra ayuda. Necesito realizar una serie de ANOVAS en loop. Os adjunto unos datos ficticios en este email. Dichos datos tienen 3 variables: 1)Valor: corresponde a la variable dependiente y es numérica 2) Grupo: Corresponde a la variable independiente y es u factor 3) Iteracción: Corresponde a la variable sobre la cual hay que repetir los ANOVAs con las variables anteriores y es un factor. Abajo os pego el código con el que he estado trabajando pero que no lo tengo bien, porque no puedo acceder a los resultados. Os agradeceria mucho si: -Me podéis ayudar a mejorar este código para que funcione -Si me podéis sugerir alguna fuente que explique bien cómo hacer un loop porque no he dado con los blogs adecuados. -Que me digáis cuál es vuestro libro/web de estadística favorito para profundizar en la matemática detrás de los análisis. Mil gracias compañeros Yésica library(agricolae) library(readxl) datosa_fict <- read_excel("datosa-fict.xlsx") #Cambiar nombre a la base de datos datos=datosa_fict #Copiar la variable sobre la que hacer el loop iter=datos$iteraccion #Crear el loop for(i in iter) { res=aov(valor~Grupo,data=datos) } #Salvar los resultados ANOVA(res) ------------ próxima parte ------------ Se ha borrado un adjunto en formato HTML... URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20230718/304fcf60/attachment.html> ------------ próxima parte ------------ A non-text attachment was scrubbed... Name: datosa-fict.xlsx Type: application/vnd.openxmlformats-officedocument.spreadsheetml.sheet Size: 10992 bytes Desc: no disponible URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20230718/304fcf60/attachment.xlsx>
Buenas, Yesica: La variable iteraccion parece una variable factor, en cuyo caso ¿no se debería incluir en el anova? en vez de repetir muchas veces el anova. Si se trata de repeticiones también hay modelos anova con repeticiones. En tu código estás haciendo el mismo anova muchas veces pues no cambias los datos en el bucle, además estás haciendo un análisis por cada fila sin eliminar las repeticiones. Por otra parte no imprimes ni guardas nada dentro del bucle con lo que cuando sales lo que calculas ya no está. Prueba lo siguiente. ## Extraer los valores de iteraccion sin repeticiones iter <- levels(factor(data$iteraccion)) for(i in iter) { ## Hacemos el análisis quedándonos con las filas correspondientes res = aov(valor~Grupo, data=datos[datos$iteraccion == iter,]) print(summary(res)) } Un saludo. P.D.: No lo he probado, puede que tenga alguna errata. El mar, 18-07-2023 a las 08:25 +0200, Yesica Pallavicini Fernandez escribió:> Buenos días y gracias de antemano por vuestra ayuda. > > Necesito realizar una serie de ANOVAS en loop. > Os adjunto unos datos ficticios en este email. > Dichos datos tienen 3 variables: > 1)Valor: corresponde a la variable dependiente y es numérica > 2) Grupo: Corresponde a la variable independiente y es u factor > 3) Iteracción: Corresponde a la variable sobre la cual hay que repetir los ANOVAs con las variables anteriores y es un factor. > > Abajo os pego el código con el que he estado trabajando pero que no lo tengo bien, porque no puedo acceder a los resultados. > Os agradeceria mucho si: > -Me podéis ayudar a mejorar este código para que funcione > -Si me podéis sugerir alguna fuente que explique bien cómo hacer un loop porque no he dado con los blogs adecuados. > -Que me digáis cuál es vuestro libro/web de estadística favorito para profundizar en la matemática detrás de los análisis. > > Mil gracias compañeros > Yésica > > library(agricolae) > library(readxl) > datosa_fict <- read_excel("datosa-fict.xlsx") > #Cambiar nombre a la base de datos > datos=datosa_fict > #Copiar la variable sobre la que hacer el loop > > iter=datos$iteraccion > #Crear el loop > for(i in iter) { > res=aov(valor~Grupo,data=datos) > > } > #Salvar los resultados > ANOVA(res) > _______________________________________________ > R-help-es mailing list > R-help-es en r-project.org > https://urldefense.com/v3/__https://stat.ethz.ch/mailman/listinfo/r-help-es__;!!D9dNQwwGXtA!VjLSaDKJejCZoZy1TgKKYQML4mA7DfvcDtUBYj5iRpdFNSsGUU5lEYtfkt0EYWTPD_-5HztE7U6zv48$ >
Hola Yésica, A lo mejor hay otras formas, pero esta funciona bien: R> u <- unique(datos$iteraccion) R> resultado <- lapply(u, function(i){ aov(valor ~ Grupo, data = subset(datos, iteraccion == i)) }) R> names(resultado) <- u Para acceder a la tabla ANOVA para iteraccion T simplemente haces> resultado[['T']]Call: aov(formula = valor ~ Grupo, data = subset(datos, iteraccion = i)) Terms: Grupo Residuals Sum of Squares 0.70326 62.83939 Deg. of Freedom 2 13 Residual standard error: 2.19859 Estimated effects may be unbalanced Si quieres obtener los estadísticos correspondientes basta con hacer R> summary(resultado[['T']]) Df Sum Sq Mean Sq F value Pr(>F) Grupo 2 0.70 0.352 0.073 0.93 Residuals 13 62.84 4.834 Espero sea de utilidad. Saludos, Jorge.- On Tue, Jul 18, 2023 at 1:25?AM Yesica Pallavicini Fernandez < yesipalla en gmail.com> wrote:> Buenos días y gracias de antemano por vuestra ayuda. > > Necesito realizar una serie de ANOVAS en loop. > Os adjunto unos datos ficticios en este email. > Dichos datos tienen 3 variables: > 1)Valor: corresponde a la variable dependiente y es numérica > 2) Grupo: Corresponde a la variable independiente y es u factor > 3) Iteracción: Corresponde a la variable sobre la cual hay que repetir los > ANOVAs con las variables anteriores y es un factor. > > Abajo os pego el código con el que he estado trabajando pero que no lo > tengo bien, porque no puedo acceder a los resultados. > Os agradeceria mucho si: > -Me podéis ayudar a mejorar este código para que funcione > -Si me podéis sugerir alguna fuente que explique bien cómo hacer un loop > porque no he dado con los blogs adecuados. > -Que me digáis cuál es vuestro libro/web de estadística favorito para > profundizar en la matemática detrás de los análisis. > > Mil gracias compañeros > Yésica > > library(agricolae) > library(readxl) > datosa_fict <- read_excel("datosa-fict.xlsx") > #Cambiar nombre a la base de datos > datos=datosa_fict > #Copiar la variable sobre la que hacer el loop > > iter=datos$iteraccion > #Crear el loop > for(i in iter) { > res=aov(valor~Grupo,data=datos) > > } > #Salvar los resultados > ANOVA(res) > _______________________________________________ > R-help-es mailing list > R-help-es en r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es >[[alternative HTML version deleted]]
Hola, library(dplyr) library(purrr) library(tidyr) datos_anova <- datosa_fict |> nest(data = c(valor, Grupo)) |> mutate(anova = map(.x = data, .f = ~aov(valor ~ Grupo, data = .x)), summary = map(.x = anova, .f = ~summary(.x))) # Acceder a los resultados datos_anova$anova[datos_anova$iteraccion == "Z"] datos_anova$summary[datos_anova$iteraccion == "Z"] El mar, 18 jul 2023 a las 1:25, Yesica Pallavicini Fernandez (< yesipalla en gmail.com>) escribió:> Buenos días y gracias de antemano por vuestra ayuda. > > Necesito realizar una serie de ANOVAS en loop. > Os adjunto unos datos ficticios en este email. > Dichos datos tienen 3 variables: > 1)Valor: corresponde a la variable dependiente y es numérica > 2) Grupo: Corresponde a la variable independiente y es u factor > 3) Iteracción: Corresponde a la variable sobre la cual hay que repetir los > ANOVAs con las variables anteriores y es un factor. > > Abajo os pego el código con el que he estado trabajando pero que no lo > tengo bien, porque no puedo acceder a los resultados. > Os agradeceria mucho si: > -Me podéis ayudar a mejorar este código para que funcione > -Si me podéis sugerir alguna fuente que explique bien cómo hacer un loop > porque no he dado con los blogs adecuados. > -Que me digáis cuál es vuestro libro/web de estadística favorito para > profundizar en la matemática detrás de los análisis. > > Mil gracias compañeros > Yésica > > library(agricolae) > library(readxl) > datosa_fict <- read_excel("datosa-fict.xlsx") > #Cambiar nombre a la base de datos > datos=datosa_fict > #Copiar la variable sobre la que hacer el loop > > iter=datos$iteraccion > #Crear el loop > for(i in iter) { > res=aov(valor~Grupo,data=datos) > > } > #Salvar los resultados > ANOVA(res) > _______________________________________________ > R-help-es mailing list > R-help-es en r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es >-- *Aviso legal:* El contenido de este mensaje y los archivos adjuntos son confidenciales y de uso exclusivo de la Universidad Nacional de Colombia. Se encuentran dirigidos sólo para el uso del destinatario al cual van enviados. La reproducción, lectura y/o copia se encuentran prohibidas a cualquier persona diferente a este y puede ser ilegal. Si usted lo ha recibido por error, infórmenos y elimínelo de su correo. Los Datos Personales serán tratados conforme a la Ley 1581 de 2012 y a nuestra Política de Datos Personales que podrá consultar en la página web www.unal.edu.co <http://www.unal.edu.co/>.* *Las opiniones, informaciones, conclusiones y cualquier otro tipo de dato contenido en este correo electrónico, no relacionados con la actividad de la Universidad Nacional de Colombia, se entenderá como personales y de ninguna manera son avaladas por la Universidad. [[alternative HTML version deleted]]