[R-es] glm com etiquetas en las variables

[Jose Luis Cañadas Reche]

Hola.

Otra idea es utilizar el paquete memisc

library(MASS)
data(birthwt, package="MASS")
birthwt$low  <- factor(birthwt$low)
birthwt$race <- factor(birthwt$smoke)
REG_LOG <- glm (low ~ smoke, family = "binomial", data = birthwt)
summary (REG_LOG)

library(memisc)

mtable1 <- mtable(REG_LOG)

mtable1 <- relabel(mtable1,
                       "(Intercept)" = "Constante",
                       smoke = "Consumo tabaco embarazo"
)
mtable1

Calls:
REG_LOG: glm(formula = low ~ smoke, family = "binomial", data = birthwt)

==================================
Constante                -1.087***
                          (0.215)
Consumo tabaco embarazo   0.704*
                          (0.320)
----------------------------------
Aldrich-Nelson R-sq.        0.025
McFadden R-sq.              0.021
Cox-Snell R-sq.             0.025
Nagelkerke R-sq.            0.036
phi                         1.000
Likelihood-ratio            4.867
p                           0.027
Log-likelihood           -114.902
Deviance                  229.805
AIC                       233.805
BIC                       240.288
N                         189
==================================

# o un fichero separado por comas
write.mtable(mtable1, file = "resumen.csv", colsep=",")



El 21/07/15 a las 11:26, Carlos Ortega escribió:
> Hola,
>
> Esta es una alternativa sólo válida para el modelo de ejemplo, pero te da
> una idea de cómo se podría hacer...
>
> Se puede generalizar para que tengas un vector con las relaciones nombres
> de las variables y su descripción y al hacer la asignación (mySummary$call)
> las cambie programáticamente.
>
> #-----------------------
>
> mySummary <- summary (REG_LOG)
> mySummary$call <- c("glm(formula = low = 'Bajo peso...' ~ smoke = 'Consumo
> tabaco..', family = 'binomial', data = birthwt)")
> mySummary
> #-----------------------
>
>
> Saludos,
> Carlos Ortega
> www.qualityexcellence.es
>
> El 21 de julio de 2015, 10:57, Griera <griera en yandex.com> escribió:
>
>> Hola:
>>
>> Si aún hay alguien que no está de vacaciones, igual me pueden ayudar.
>>
>> Quiero ajustar unos modelos:
>>
>>      REG_LOG <- glm (low ~ X, family = "binomial", data = DATOS)
>>
>> Ejemplo:
>>    library(MASS)
>>    data(birthwt, package="MASS")
>>    birthwt$low  <- factor(birthwt$low)
>>    birthwt$race <- factor(birthwt$smoke)
>>    REG_LOG <- glm (low ~ smoke, family = "binomial", data = birthwt)
>>    summary (REG_LOG)
>>
>> Se pueden colocar etiquetas en las variable de la fórmula de manera que en
>> los resultados salgan las etiquetas y no el nombre de la variable? Seria
>> una cosa como:
>>
>> REG_LOG <- glm (low = "Bajo peso recién nacido" ~ smoke = "Consumo tabaco
>> embarazo", family = "binomial", data = birthwt)
>>
>> que, evidentemente no funciona.
>>
>> Muchas gracias y saludos
>>
>> _______________________________________________
>> R-help-es mailing list
>> R-help-es en r-project.org
>> https://stat.ethz.ch/mailman/listinfo/r-help-es
>>
>
>