Estoy ejecutando una regresión logística con una variable pronóstica categórica con niveles 0-6. Por defecto, R considera el nivel 0 como el grupo de referencia.Cambio del grupo de referencia para la variable predictora categórica en la regresión logística
¿Cómo puedo saber R para utilizar, por ejemplo, el nivel 3 en lugar de nivel 0 como referencia sin cambiar el nombre de los niveles?
Utilice la función C a defina sus contrastes en el marco de datos.
Si su trama de datos es DF y la variable factor es fct, entonces
DF$fct <- C(DF$fct, contr.treatment, base=3)
(no probado).
Funcionó muy bien, gracias! –
Esto es muy fácil con los paquetes de Frank Harrell (que permite muchas funciones adicionales).
Por ejemplo, para una trama de datos llamado 'df'
library(Hmisc)
library(rms)
dd=datadist(df)
options(datadist='dd')
(m=lrm(y ~ catvar, data=df)) #uses the mode as the reference group
summary(m, catvar=3) #using level=3 as the reference group
summary(m, catvar=0) #using level=0 as the reference group
(Y, algunos datos de ejemplo habría estado bien - siempre se puede utilizar dput para dar salida a un pequeño conjunto de datos)
Solución flexible interesante, gracias. Tendrá que leer lo que datadist está representando. –
La función C ya se ha sugerido, también mira contrasts, relevel y reorder, entre otros.
También es muy útil. Gracias. –
Puede utilizar la función Vuelva a nivelar: trama de datos x1 $ < - Vuelva a nivelar (datafrmae $ x1 "escriba aquí la categoría de referencia")
¿Existe una solución que trabajar con las etiquetas? Por ejemplo, si los niveles son '" a "," b ", ...," g "', 'R' usan orden alfabético. ¿Cómo se usa '" d "' como referencia? He intentado 'C (F, contr.treatment, base =" a ")' pero no funciona. –
'C (F, contr.treatment, base = which (levels (F) ==" d "))' funciona pero no es elegante. –
@MathieuDubois, mi solución perezosa es engañar a la alfabetización dando nombres a las categorías como "aCategory1" "bCategory2", etc. –