Sí llamamos DF al data.frame que tienes y DF$ID es la variable con los ID de
los lotes puede hacer con el código de que dispones lo siguiente (es
equivalente a un bucle for):
L <- unique(DF$ID)
Hay que tener cuidado que ID no sea un factor, si lo es
L <- unique(as.character(DF$ID))
Resultados <- lapply(L, function(x){
X <- subset(DF, ID == x)
Aquí vendrían las operaciones que quieras hacer con X (el data.frame
para un solo ID o, como lo llamas, lote)
hasta que obtienes el objeto que persigues, al que llamaré Z.
Y, por último
return(Z)
})
names(Resultados) <- L
Esta última instrucción te permite conocer a que lote corresponden los
resultados.
Resultados es una lista y puedes ver uno de ellos con Resultados[["C001"]] (por
ejemplo).
Espero haber sido de ayuda.
-----Mensaje original-----
De: R-help-es [mailto:[email protected]] En nombre de Andres
Hirigoyen
Enviado el: lunes, 30 de mayo de 2016 16:20
Para: Carlos Ortega
CC: Lista R
Asunto: Re: [R-es] Codigo
Carlos, capaz no me explique bien, va de nuevo...
Tengo un dataFrame con varias variables que describen Lotes (127), he creado un
código que calcula nuevas variables, re codifica otras y ordenas otras tantas.
Es muy laborioso correr ese código 127 veces por lo cual quería hacer un
sentencia (del tipo if() o repeat()), para obtener al final todos los lotes
procesados.
Saludos y gracias por tan rápida respuesta
El 30 de mayo de 2016, 11:06, Carlos Ortega <[email protected]>
escribió:
> 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
> <[email protected]>
> 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 <[email protected]> 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" <[email protected]>
>> > Para: [email protected]
>> > 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
>> los
>> > cuales debo efectuares un GLM y con los parámetros obtenidos
>> > calcular
>> otras
>> > variables para luego hacer gráficas, razón por la cual no puedo
>> procesarlos
>> > 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]]
>> >
>> > _______________________________________________
>> > R-help-es mailing list
>> > [email protected]
>> > https://stat.ethz.ch/mailman/listinfo/r-help-es
>> >
>>
>> [[alternative HTML version deleted]]
>>
>> _______________________________________________
>> R-help-es mailing list
>> [email protected]
>> https://stat.ethz.ch/mailman/listinfo/r-help-es
>>
>
>
>
> --
> Saludos,
> Carlos Ortega
> www.qualityexcellence.es
>
[[alternative HTML version deleted]]
_______________________________________________
R-help-es mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-help-es
_______________________________________________
R-help-es mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-help-es
