<html><head></head><body><div class="ydpdfcb6fbcyahoo-style-wrap" style="font-family:courier new, courier, monaco, monospace, sans-serif;font-size:16px;"><div></div>
<div dir="ltr" data-setdir="false">Hola</div><div dir="ltr" data-setdir="false">otra opción, a ver si te sirve.</div><div dir="ltr" data-setdir="false">Acabarás teniendo 6 archivos txt y 6 archivos xlsx cada uno con el resultado de cada una de las iteraciones.</div><div dir="ltr" data-setdir="false"><br></div><div dir="ltr" data-setdir="false"><div><div>library(purrr)</div><div>library(broom)</div></div><br></div><div dir="ltr" data-setdir="false"><br></div><div dir="ltr" data-setdir="false"><br></div><div dir="ltr" data-setdir="false"> <div><div>iter= unique(datos$iteraccion)</div><div>#Crear el loop</div><div>for(i in iter) {</div><div><br></div><div>a<-tidy(aov(valor~Grupo,data=subset(datos, iteraccion==i)))</div><div><br></div><div><br></div><div>capture.output(a,i, file =paste0(i, ".txt"))</div><div>print(a)</div><div>openxlsx:: write.xlsx(a,file = paste0(i,'.xlsx'),</div><div> sheetName = i, append = TRUE)</div><div><br></div><div>}</div></div><br></div><div dir="ltr" data-setdir="false"><img title="Imagen en el mensaje" alt="Imagen en el mensaje" src="cid:297e2bee-6ec2-637f-8647-4aee5d5f7089@yahoo.com" class="yahoo-inline-image" draggable="false" style="max-width: 643px; width: 100%;" data-inlineimagemanipulating="true" data-id="<297e2bee-6ec2-637f-8647-4aee5d5f7089@yahoo.com>"><br><span></span><br></div><div dir="ltr" data-setdir="false"><br></div><div><br></div>
</div><div id="yahoo_quoted_0309463671" class="yahoo_quoted">
<div style="font-family:'Helvetica Neue', Helvetica, Arial, sans-serif;font-size:13px;color:#26282a;">
<div>
En martes, 18 de julio de 2023, 16:53:30 CEST, Jorge I Velez <jorgeivanvelez@gmail.com> escribió:
</div>
<div><br></div>
<div><br></div>
<div>Hola Yésica,<br clear="none"><br clear="none">A lo mejor hay otras formas, pero esta funciona bien:<br clear="none"><br clear="none">R> u <- unique(datos$iteraccion)<br clear="none">R> resultado <- lapply(u, function(i){<br clear="none"> aov(valor ~ Grupo, data = subset(datos, iteraccion == i))<br clear="none">})<br clear="none">R> names(resultado) <- u<br clear="none"><br clear="none">Para acceder a la tabla ANOVA para iteraccion T simplemente haces<br clear="none"><br clear="none">> resultado[['T']]<br clear="none">Call:<br clear="none"> aov(formula = valor ~ Grupo, data = subset(datos, iteraccion ==<br clear="none"> i))<br clear="none"><br clear="none">Terms:<br clear="none"> Grupo Residuals<br clear="none">Sum of Squares 0.70326 62.83939<br clear="none">Deg. of Freedom 2 13<br clear="none"><br clear="none">Residual standard error: 2.19859<br clear="none">Estimated effects may be unbalanced<br clear="none"><br clear="none">Si quieres obtener los estadísticos correspondientes basta con hacer<br clear="none"><br clear="none">R> summary(resultado[['T']])<br clear="none">Df Sum Sq Mean Sq F value Pr(>F)<br clear="none">Grupo 2 0.70 0.352 0.073 0.93<br clear="none">Residuals 13 62.84 4.834<br clear="none"><br clear="none">Espero sea de utilidad.<br clear="none"><br clear="none">Saludos,<br clear="none">Jorge.-<br clear="none"><br clear="none"><br clear="none">On Tue, Jul 18, 2023 at 1:25 AM Yesica Pallavicini Fernandez <<br clear="none"><a shape="rect" ymailto="mailto:yesipalla@gmail.com" href="mailto:yesipalla@gmail.com">yesipalla@gmail.com</a>> wrote:<br clear="none"><br clear="none">> Buenos días y gracias de antemano por vuestra ayuda.<br clear="none">><br clear="none">> Necesito realizar una serie de ANOVAS en loop.<br clear="none">> Os adjunto unos datos ficticios en este email.<br clear="none">> Dichos datos tienen 3 variables:<br clear="none">> 1)Valor: corresponde a la variable dependiente y es numérica<br clear="none">> 2) Grupo: Corresponde a la variable independiente y es u factor<br clear="none">> 3) Iteracción: Corresponde a la variable sobre la cual hay que repetir los<br clear="none">> ANOVAs con las variables anteriores y es un factor.<br clear="none">><br clear="none">> Abajo os pego el código con el que he estado trabajando pero que no lo<br clear="none">> tengo bien, porque no puedo acceder a los resultados.<br clear="none">> Os agradeceria mucho si:<br clear="none">> -Me podéis ayudar a mejorar este código para que funcione<br clear="none">> -Si me podéis sugerir alguna fuente que explique bien cómo hacer un loop<br clear="none">> porque no he dado con los blogs adecuados.<br clear="none">> -Que me digáis cuál es vuestro libro/web de estadística favorito para<br clear="none">> profundizar en la matemática detrás de los análisis.<br clear="none">><br clear="none">> Mil gracias compañeros<br clear="none">> Yésica<br clear="none">><br clear="none">> library(agricolae)<br clear="none">> library(readxl)<br clear="none">> datosa_fict <- read_excel("datosa-fict.xlsx")<br clear="none">> #Cambiar nombre a la base de datos<br clear="none">> datos=datosa_fict<br clear="none">> #Copiar la variable sobre la que hacer el loop<br clear="none">><br clear="none">> iter=datos$iteraccion<br clear="none">> #Crear el loop<br clear="none">> for(i in iter) {<br clear="none">> res=aov(valor~Grupo,data=datos)<br clear="none">><br clear="none">> }<br clear="none">> #Salvar los resultados<br clear="none">> ANOVA(res)<br clear="none">> _______________________________________________<br clear="none">> R-help-es mailing list<br clear="none">> <a shape="rect" ymailto="mailto:R-help-es@r-project.org" href="mailto:R-help-es@r-project.org">R-help-es@r-project.org</a><br clear="none">> <a shape="rect" href="https://stat.ethz.ch/mailman/listinfo/r-help-es" target="_blank">https://stat.ethz.ch/mailman/listinfo/r-help-es</a><br clear="none">><br clear="none"><br clear="none"> [[alternative HTML version deleted]]<div class="yqt1985650995" id="yqtfd95031"><br clear="none"><br clear="none">_______________________________________________<br clear="none">R-help-es mailing list<br clear="none"><a shape="rect" ymailto="mailto:R-help-es@r-project.org" href="mailto:R-help-es@r-project.org">R-help-es@r-project.org</a><br clear="none"><a shape="rect" href="https://stat.ethz.ch/mailman/listinfo/r-help-es" target="_blank">https://stat.ethz.ch/mailman/listinfo/r-help-es</a><br clear="none"></div></div>
</div>
</div></body></html>