Hola,
Puedes utilizar el código de Oliver utilizando la librerÃa "broom" que
permite obtener los parámetros (coeficientes) y caracterÃsticas del modelo
(null.deviance, AIC, BIC) en un cómodo data.frame
El código serÃa algo asà como esto.
library(broom)
L <- list(L1, L2..., L120)
ajuste <- function(L) {
g <- glm(y ~., data = L)
g_df <- tidy(g) #tidy es una función del paquete "broom" que
devuelve los coeficientes y sus errores.
return(g_df)
}
fit <- apply(L, ajuste)
Saludos,
Carlos Ortega
www.qualityexcellence.es
El 30 de mayo de 2016, 15:15, Andres Hirigoyen <andreshirigoyen en gmail.com>
escribió:
Oliver gracias por le pista!
Es un poco mas complicado, el código tiene varios pasos incluidos calculo
de nuevas variables con los betas del glm y otras vueltas. NecesitarÃa una
sentencia que tome cada lote haga los cálculos por separado y luego los
une
en el dataframe de salida
El 30 de mayo de 2016, 10:02, Olivier Nuñez <onunez en unex.es> escribió:
Mira lapply ....
Si L= list(L1,L2,...,L120) es una lista de tus lotes
ajuste <- function(L) glm(y~x,data=L)
fit=lapply(L1,ajuste)
donde "fit" es la lista de 120 ajustes.
Un saludo. Olivier
----- Mensaje original -----
De: "Andres Hirigoyen" <andreshirigoyen en gmail.com>
Para: r-help-es en r-project.org
Enviados: Lunes, 30 de Mayo 2016 14:55:19
Asunto: [R-es] Codigo
Buenos dÃas tengo una consulta sobre un código rutinario que quiero
ejecutar.
Tengo 120 Lotes con un ID diferente (del C001 al C120) a cada uno de
cuales debo efectuares un GLM y con los parámetros obtenidos calcular
variables para luego hacer gráficas, razón por la cual no puedo
todos juntos.
Intente hacer una sentencia repeat() o con for() para que me ejecute de
forma individual cada lote y no tener que hacerlo a mano, pero no me
sale...
¿¿Alguna idea??
[[alternative HTML version deleted]]