[R-es] realizar ANOVAs en Loop

2023-07-18 Thread Yesica Pallavicini Fernandez 
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)


datosa-fict.xlsx
Description: MS-Excel 2007 spreadsheet
___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

[R-es] Supuestos de una ANOVA

2023-06-16 Thread Yesica Pallavicini Fernandez 
Buenos días y muchas gracias por adelantado
En cuando a probar gráficamente los supuestos de normalidad y
homocedasticidad de una ANOVA

Lo estoy haciendo de forma gráfica de la siguiente manera:
Primero hago el nova con: aov()
luego hago la comprobación de los supuestos con: plot(modelo)
y me salen 4 gráficos;
1 un scatter plot de los "residuals" vs " Fitted",
2"standardizez residuals" vs " Fitted",
3QQplot,
4 Residuals vs leverage
Pero ¿Qué es el leverage?
 1)¿teneis algun script que reemplace a plot(modelo) y que solo contenga el
QQplot y los valores residuales frente a los fitted?

2) Shapiro.test () no funciona para observaciones mayores a 500 y mis datos
tienen más de 3 observaciones. ¿sabeis de algun otro test para este
caso?

Muchísimas gracias y que acabéis fenomenal la semana

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

[R-es] loop para sumar una constante

2022-11-25 Thread Yesica Pallavicini Fernandez 
Hola.
En el mensaje anterior no he expresado correctamente lo que quiero.
A ver si ahora me explico mejor.

Tengo de partida un valor; por ejemplo 1
A ese 1 quiero sumarle 5 y luego al punto anterior quiero sumar otros 5 más
y así hasta llegar a 30.
El primer punto sería 6, el segundo 6+5=11, el siguiente 11+5= 16 y
así sucesivamente hasta llegar a 30.
Es decir que necesito que me devuelva como resultado
1, 6,11,16,21,26

Gracias por vuestra ayuda

Yésica

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

[R-es] Loop para sumar

2022-11-25 Thread Yesica Pallavicini Fernandez 
Buenos días.
Necesito dividir un mapa en celdas de 250m.
Tengo la UTM de partida y la de Llegada.

El comando que quisiera elaborar sería:
Desde este punto de partida, suma 250m y el posterior a 250m desde el punto
anterior y asi sucesivamente hasta llegar al punto de llegada. Solo voy a
utilizar la UTM X.
Y quisiera crear un vector que contenga todos esos puntos.
¿Me explico bien?
169631.4
568497.6
¿R puede aguantar mas de un millón quinientosmil datos que va a tener ese
vector?
Este es mi intento, pero me falta decirle que sume los 250 al punto
anterior.

datos=for(for(i in  169631.4: 568497.6  ){
 print(i^2)
}

Muchas gracias

Yésica

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

[R-es] Correlaciones con permutaciones.

2020-06-02 Thread Yesica Pallavicini Fernandez 
Hola Necesito hacer correlaciones con permutaciones de los scores de los
ejes de un PCA con ciertas variables.
¿Alguien sabe qué función utilizar?
Gracias
Yésica

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Graficos: como hacer que las etiquetas no estén sobrepuestas

2020-05-12 Thread Yesica Pallavicini Fernandez 
GRACIAS!
Un gran saludo
Yesica

El mar., 12 may. 2020 a las 13:15, Carlos Ortega ()
escribió:

> Hola,
> Sí, hay un paquete específicamente para tratar esto. "ggrepel":
> https://cloud.r-project.org/web/packages/ggrepel/index.html
>
> Claro tu gráfico tiene que ser ggplot.
>
> Si no es ggplot, te recomiendo entonces que lo generes con el paquete
> "factoextra" que sí que genera gráficos ggplot y puedes aplicar esta mejora
> de ggrepel.
>
> Gracias,
> Carlos Ortega
> www.qualityexcellence.es
>
> El mar., 12 may. 2020 a las 12:07, Yesica Pallavicini Fernandez (<
> yesipa...@gmail.com>) escribió:
>
>> Hola,
>> Estoy haciendo un PCA con el paquete ade4.
>> En el gráfico, las etiquetas de las especies suelen quedar
>> sobrepuestas unas sobre las otras y no se pueden distinguir
>> individualmente. ¿Hay alguna manera de generar un poco de espacio entre
>> ellas para poder visualizarlas todas?
>>
>> gracias
>> Yésica
>>
>> [[alternative HTML version deleted]]
>>
>> ___
>> R-help-es mailing list
>> R-help-es@r-project.org
>> https://stat.ethz.ch/mailman/listinfo/r-help-es
>>
>
>
> --
> Saludos,
> Carlos Ortega
> www.qualityexcellence.es
>

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

[R-es] Graficos: como hacer que las etiquetas no estén sobrepuestas

2020-05-12 Thread Yesica Pallavicini Fernandez 
Hola,
Estoy haciendo un PCA con el paquete ade4.
En el gráfico, las etiquetas de las especies suelen quedar
sobrepuestas unas sobre las otras y no se pueden distinguir
individualmente. ¿Hay alguna manera de generar un poco de espacio entre
ellas para poder visualizarlas todas?

gracias
Yésica

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

[R-es] ver los resutados cuando son muy extensos y normalmente se cortan

2018-03-08 Thread Yesica Pallavicini Fernandez 
Gracias Carlos y Xavy

Efectivamente, cuando la salida del resultado es muy grande y por defecto
se cortan y no los puedes ver, se puede usar el comando
options(max.print=9)


Saludos

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Correlacionar una variable con muchas a la vez

2018-03-08 Thread Yesica Pallavicini Fernandez 
 Hola,
Estoy intentando correlacionar una sola variable con 52 mas
 ¿Sabéis cómo se puede automatizar la función?
Gracias
Yésica

Ya lo he sacado
cor(data)[,1]
Con esta función se correlaciona la primer variable con todas las demás del
arreglo de datos.

¿Cómo me puedo contestar a mi misma para que todos los veáis?
Gracias



El 8 de marzo de 2018, 8:51, Yesica Pallavicini Fernandez <
yesipa...@gmail.com> escribió:

> Hola,
> Estoy intentando correlacionar una sola variable con 52 mas
>  ¿Sabéis cómo se puede automatizar la función?
> Gracias
> Yésica
>
>

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

[R-es] Correlacionar una variable con muchas a la vez

2018-03-07 Thread Yesica Pallavicini Fernandez 
Hola,
Estoy intentando correlacionar una sola variable con 52 mas
 ¿Sabéis cómo se puede automatizar la función?
Gracias
Yésica

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

[R-es] exportar resultados a CSV

2018-03-07 Thread Yesica Pallavicini Fernandez 
Hola,
tengo dos dudas
1) he corrido una matriz de correlación con demasiadas variables y R no me
muestra todos los resultados porque ha llegado a su máxima capacidad. ¿cómo
puedo ver tdos los resultados?

2) Quisiera exportar los resultados a un CVS para poder verlos mejor ¿Se
puede hacer?

Gracias
Yésica

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Resumen de R-help-es, Vol 108, Envío 24

2018-02-19 Thread Yesica Pallavicini Fernandez 
Hola Juan Bautista.
Yo instalé agricolae la semana pasada y sin problemas..
No se cómo ayudarte.
A no ser que desinstales R y lo vuelvas a instalar ( perderías todos los
paquetes instalados)

Saludos
Yésica.


El 17 de febrero de 2018, 13:21,  escribió:

> Envíe los mensajes para la lista R-help-es a
> r-help-es@r-project.org
>
> Para subscribirse o anular su subscripción a través de la WEB
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
> O por correo electrónico, enviando un mensaje con el texto "help" en
> el asunto (subject) o en el cuerpo a:
> r-help-es-requ...@r-project.org
>
> Puede contactar con el responsable de la lista escribiendo a:
> r-help-es-ow...@r-project.org
>
> Si responde a algún contenido de este mensaje, por favor, edite la
> linea del asunto (subject) para que el texto sea mas especifico que:
> "Re: Contents of R-help-es digest...". Además, por favor, incluya en
> la respuesta sólo aquellas partes del mensaje a las que está
> respondiendo.
>
>
> Asuntos del día:
>
>1. RV: no puedo cargar el paquete "agricolae" (Juan Bautista)
>2. Re: RV: no puedo cargar el paquete "agricolae" (Javier Marcuzzi)
>
> --
>
> Message: 1
> Date: Sat, 17 Feb 2018 11:43:47 +
> From: Juan Bautista 
> To: "'R Help Es  (r-help-es@r-project.org)'" 
> Subject: [R-es] RV: no puedo cargar el paquete "agricolae"
> Message-ID:
> 
> Content-Type: text/plain; charset="iso-8859-1"
>
>
> Cada vez que intento cargar el paquete "agricolae" me da el siguiente
> mensaje:
>
> Error: package or namespace load failed for 'agricolae' in
> loadNamespace(i, c(lib.loc, .libPaths()), versionCheck = vI[[i]]):
> there is no package called 'spData'
>
> He intentado todo actualizar los paquetes, volverlos a instalar pero no
> hay forma.
>
> Si sabeis la solución ??
>
> Un saludo.
>
> Juan Bautista Relloso Barrio
>
>
>  próxima parte 
> A non-text attachment was scrubbed...
> Name: winmail.dat
> Type: application/ms-tnef
> Size: 14747 bytes
> Desc: no disponible
> URL:  attachments/20180217/0724f79f/attachment-0001.bin>
>
>
> --
>
> Message: 2
> Date: Sat, 17 Feb 2018 09:20:58 -0300
> From: Javier Marcuzzi 
> To: Juan Bautista 
> Cc: "R Help Es (r-help-es@r-project.org)" 
> Subject: Re: [R-es] RV: no puedo cargar el paquete "agricolae"
> Message-ID:
>  3nn1...@mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> Estimado Juan Bautista
>
> Yo intenté instalarlo y cargarlo, en mi caso no sale ese error, copio y
> pego todo lo que informa mi ambiente de trabajo por su puedes duplicarlo.
>
> > R.version   _
> platform   x86_64-w64-mingw32
> arch   x86_64
> os mingw32
> system x86_64, mingw32
> status
> major  3
> minor  4.3
> year   2017
> month  11
> day30
> svn rev73796
> language   R
> version.string R version 3.4.3 (2017-11-30)
> nickname   Kite-Eating Tree>
> install.packages("agricolae")Installing package into
> ‘C:/Users/HP/Documents/R/win-library/3.4’
> (as ‘lib’ is unspecified)also installing the dependencies ‘gtools’,
> ‘gdata’, ‘combinat’, ‘sp’, ‘spData’, ‘LearnBayes’, ‘deldir’,
> ‘gmodels’, ‘expm’, ‘klaR’, ‘spdep’, ‘AlgDesign’
> trying URL 'https://mran.microsoft.com/snapshot/2018-01-01/bin/
> windows/contrib/3.4/gtools_3.5.0.zip'Content
> type 'application/zip' length 144747 bytes (141 KB)downloaded 141 KB
> trying URL 'https://mran.microsoft.com/snapshot/2018-01-01/bin/
> windows/contrib/3.4/gdata_2.18.0.zip'Content
> type 'application/zip' length 1186810 bytes (1.1 MB)downloaded 1.1 MB
> trying URL 'https://mran.microsoft.com/snapshot/2018-01-01/bin/
> windows/contrib/3.4/combinat_0.0-8.zip'Content
> type 'application/zip' length 29663 bytes (28 KB)downloaded 28 KB
> trying URL 'https://mran.microsoft.com/snapshot/2018-01-01/bin/
> windows/contrib/3.4/sp_1.2-5.zip'Content
> type 'application/zip' length 1539525 bytes (1.5 MB)downloaded 1.5 MB
> trying URL 'https://mran.microsoft.com/snapshot/2018-01-01/bin/
> windows/contrib/3.4/spData_0.2.6.8.zip'Content
> type 'application/zip' length 3986137 bytes (3.8 MB)downloaded 3.8 MB
> trying URL 'https://mran.microsoft.com/snapshot/2018-01-01/bin/
> windows/contrib/3.4/LearnBayes_2.15.zip'Content
> type 'application/zip' length 1131042 bytes (1.1 MB)downloaded 1.1 MB
> trying URL 'https://mran.microsoft.com/snapshot/2018-01-01/bin/
> windows/contrib/3.4/deldir_0.1-14.zip'Content
> type 'application/zip' length 237317 bytes (231 KB)downloaded 231 KB
> trying URL

[R-es] Tapply para medias entre mas de 2 factores

2018-02-14 Thread Yesica Pallavicini Fernandez 
Hola.
Tengo unos datos de rendimiento de distintas variedades en 2 localidades.
Quiero sacar la media de cada variedad en cada una de las localidades.
La función tapply solo me hace la media por variedad o por localidad, pero
no juntas.
¿Conoceis alguna función que calcule la media de cada variedad separada por
localidad?
Gracias
Yésica

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

[R-es] ANOVA factores anidados

2018-02-13 Thread Yesica Pallavicini Fernandez 
En cuanto meter un factor anidado en una ANOVA ¿puede quedar asi?
¿Será mejor usar modelos mixtos en vez de ANOVA?
Solo quiero testar la influencia de las repeticiones del ensayo.

A<-aov(rto~A*T*y*R*rep/A/T/R,dato))

Gracias

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

[R-es] ANOVA, otra pregunta

2018-02-13 Thread Yesica Pallavicini Fernandez 
 Si recordías, quiero hacer el siguiente modelo.
Freddy ya me dió una solución. PEro me surgió otra duda:
En el modelo completo, debo meter la repecición ( 3 repeticiones) en el
modelo. Estas repeticiones están anidadas
repeticion(ByA)
¿Cómo se pone en la función "aov" de R?

Mil gracias

Y~T,
Y~A
Y~B
Y~TxA
Y~AxB
Y~AxBxT
¿Está bien planteada la función? anova2<-lm(d$Y~d$T*d$A*d$B)

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

[R-es] ANOVA no da el mismo resultado

2018-02-13 Thread Yesica Pallavicini Fernandez 
Hola,
He realizado 2 anovas, una simple y la otra multifactorial

 anova1<-lm(d$Y~d$T)
anova2<-lm(d$Y~d$T*d$A*d$B)

Las preguntas son:
1) en anova2 quiero testar el modelo

Y~T,
Y~A
Y~B
Y~TxA
Y~AxB
Y~AxBxT
¿Está bien planteada la función? anova2<-lm(d$Y~d$T*d$A*d$B)

2) En la salida de ambas anovas, el resultado es diferente para la misma
variable (T). Es mas, en la anova 1 el resultado es no significativo y en
anova2 si lo es ¿A que se debe esto?

3) Estos análisis los he repetido en SAS ¿ Por que puede dar diferente
resultado en un programa y en otro?

Muchísimas gracias por vuestra valiosa colaboración.

Saludos
Yésica

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] incluir las replicas en anova

2018-01-05 Thread Yesica Pallavicini Fernandez 
Hola
1) Gracias José por tu comentario, tienes razón.
2) como prueba ad hoc estoy usando la función pairwise.t.test(data1,
data2,p.adjust="bonferroni")
¿Alguien sabe alguna otra función que me ponga letras (a ab b c etc.)
cuando haga la comparación?
gracias
Yesica

El 5 de enero de 2018, 12:36, Yesica Pallavicini Fernandez <
yesipa...@gmail.com> escribió:

> Hola.
> Estoy realizando una ANOVA con la función lm()
> El ensayo tiene 3 réplicas por tratamiento y no se cómo incluir estas
> réplicas en la función lm().
> Gracias por vuestra ayuda
> Yesica
>

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

[R-es] incluir las replicas en anova

2018-01-05 Thread Yesica Pallavicini Fernandez 
Hola.
Estoy realizando una ANOVA con la función lm()
El ensayo tiene 3 réplicas por tratamiento y no se cómo incluir estas
réplicas en la función lm().
Gracias por vuestra ayuda
Yesica

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es