Re: [R-es] ANOVA no da el mismo resultado

2018-02-13 Por tema Freddy Omar López Quintero
Hola. El mar, 13-02-2018 a las 12:32 +0100, Yesica Pallavicini Fernandez escribió: > 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

[R-es] ANOVA, otra pregunta

2018-02-13 Por tema 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

[R-es] ANOVA no da el mismo resultado

2018-02-13 Por tema 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

Re: [R-es] problema de fechas

2018-02-13 Por tema Álvaro Hernández
Yo entiendo que lo que quiere es saber cómo hacer el 'ifelse' para que no le devuelva un vector numérico porque para transformarlo podría hacer simplemente 'as_date(fechas2)'. El problema yo creo que es por tema coerción con el 'NA' (sumado al comportamiento de 'ifelse'). Estaba buscando

Re: [R-es] problema de fechas

2018-02-13 Por tema Carlos Ortega
Hola, Sí, en el ejemplo de la ayuda de "ifelse" se advierte de este problema y ofrece dos alternativas: una la que he propuesto y otra. Con cualquiera de las dos el problema se resuelve. ## ifelse() strips attributes ## This is important when working with Dates and factors x <-

Re: [R-es] problema de fechas

2018-02-13 Por tema Patricio Suárez Gil
Muchas gracias a todos por vuestra ayuda. La solución que propone Carlos Ortega me funciona perfectamente. Gracias mil, Carlos. Patricio Suárez Gil Unidad de Investigación Área V-Gijón Planta 5ª Impar Hospital Universitario de Cabueñes C/Prado, 395 33394 Gijón (Asturias) Tfno: 985 185 000

[R-es] ANOVA factores anidados

2018-02-13 Por tema 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]]

Re: [R-es] Crear bucle

2018-02-13 Por tema Andrés Hirigoyen
Xavier el tema que no logro solucionar es pasar de subzona==1 a subzona==2. Probé con ifelse pero me cambia los valores de toda la columan cuando cambia de zona. Gracias por responder El 13 de febrero de 2018, 18:19, Xavier-Andoni Tibau Alberdi< xaviti...@gmail.com> escribió: > Creo que

[R-es] Crear bucle

2018-02-13 Por tema Andrés Hirigoyen
Buenas tardes para tod@s (de nuevo) Tengo el siguiente dataframe: margen<-c("A","B","A","B","A","B","A","B","A","B","A","B","A","B","A","B","A","B","A","B","A","B","A","B") margen<-sort(margen, decreasing=T) long<-c(9,12,22,36,36,98,58,47,54,87,26,36,35,91,10,20,36,57,30,69,32,33,36,37)

Re: [R-es] Crear bucle

2018-02-13 Por tema Xavier-Andoni Tibau Alberdi
Perdona, toda la razón, no miré bien el problema antes de responder. Seguro que hay una manera más eficiente. Pero yo lo haría con un bucle, luego usas una bariable dummy, digamos zonna_d, que incremente su valor en 1 cada vez que subzona es ==1. Entonces zona = zonna_d. Se entiende? Saludos,

Re: [R-es] Crear bucle

2018-02-13 Por tema Carlos Ortega
Hola, Una forma puede ser la siguiente... Para comparar, he creado una nueva columna "new_zona" que es la que se va rellenando con un bucle... #--- > cont <- 1 > new_zona <- 0 > for(i in 1:nrow(data)) { + new_zona[i] <- cont + if(data$Subzona[i] == 1 & i > 1) { + cont <-

Re: [R-es] Crear bucle

2018-02-13 Por tema Andrés Hirigoyen
Excelente Carlos, muchas gracias. Me maree con [i] en el bucle que hacia. Saludos El 13 de febrero de 2018, 18:38, Carlos Ortega escribió: > Hola, > > Una forma puede ser la siguiente... > Para comparar, he creado una nueva columna "new_zona" que es la que se va >

Re: [R-es] Crear bucle

2018-02-13 Por tema Xavier-Andoni Tibau Alberdi
Creo que necesitas usar la funcion: ifelse(). Saludos! El 13 feb. 2018 22:16, "Andrés Hirigoyen" escribió: > Buenas tardes para tod@s > (de nuevo) > > Tengo el siguiente dataframe: > margen<-c("A","B","A","B","A","B","A","B","A","B","A","B"," >

Re: [R-es] ANOVA factores anidados

2018-02-13 Por tema Javier Marcuzzi
Estimada Yesica Pallavicini Fernandez Su problema es de principiante de R, no importa, escriba lo siguiente (?) cada vez que necesita conocer como escribir algo, R la ayuda con algún ejemplo como ser aov(yield ~ block + N * P + K, npk) lo siguiente / tiene un significado ?aov?lm También

Re: [R-es] Crear bucle

2018-02-13 Por tema Proyecto R-UCA
Buenas: Una forma sin bucles: > data <- read.table("/tmp/a.csv", header=TRUE, sep="", na.strings="NA", dec=".", strip.white=TRUE) > new_zon <- as.numeric(diff(data$Subzona) < 0) > new_zon <- c(1, new_zon) > new_zon <- cumsum(new_zon) > data$new_zon <- new_zon > data Margen Zona Subzona Long