Ayúdame a entender el odds ratio ajustado en regresión logística


20

Me ha costado mucho tratar de comprender el uso de la regresión logística en un artículo. El documento disponible aquí utiliza la regresión logística para predecir la probabilidad de complicaciones durante la cirugía de cataratas.

Lo que me confunde es que el documento presenta un modelo que asigna una razón de probabilidad de 1 a la línea de base que se describe a continuación:

Un paciente cuyo perfil de riesgo estaba en el grupo de referencia para todos los indicadores de riesgo (es decir, OR ajustado = 1.00 para todos en la Tabla 1) puede considerarse que tiene un "perfil de riesgo de referencia", y el modelo de regresión logística indica una "probabilidad predicha de referencia" para PCR o VL o ambos = 0.736%.

Por lo tanto, la probabilidad de 0.00736 se presenta con una razón de probabilidades de 1. Basado en la transformación de probabilidades a razones de probabilidades: , esto no puede ser igual a 1: .o=pag1-pag0.00741=0.007361-0.00736

Se vuelve aún más confuso. Las razones de probabilidades compuestas que representan múltiples covariables que tienen valores diferentes a la línea de base se usan para calcular el riesgo predicho.

... el OR compuesto de la Tabla 1 sería 1.28 X 1.58 X 2.99 X 2.46 X 1.45 X 1.60 = 34.5, y del gráfico de la Figura 1, vemos que este OR corresponde con una probabilidad predicha de PCR o VL o ambas alrededor del 20%

La única forma de llegar a los valores que el documento está dando como ejemplos es multiplicar la probabilidad de referencia con probabilidades compuestas como esta: .0.2025=(34,50 × 0.00736)1 + (34,50 × 0.00736)

Entonces, ¿qué está pasando aquí? ¿Cuál es la lógica para asignar odds ratio 1 a una probabilidad de referencia que no sea 0.5? La fórmula de actualización que se me ocurrió anteriormente presenta las probabilidades correctas para los ejemplos en el documento, pero esta no es la multiplicación directa de la razón de probabilidades que esperaría. ¿Entonces que es?


8
Es posible que tenga una simple confusión sobre la terminología: es una probabilidad , no una razón de probabilidades. Un odds ratio es la división de una de esas expresiones por otra. pag/ /(1-pag)
whuber

Respuestas:


35

Las probabilidades son una forma de expresar oportunidades. Las razones de probabilidades son solo eso: una cuota dividida por otra. Eso significa que una razón de probabilidades es lo que multiplica una probabilidad para producir otra. Veamos cómo funcionan en esta situación común.

Conversión entre probabilidades y probabilidad

Y10 Pr ( Y = 0 )Pr(Y=1)0 0Pr(Y=0 0)

Posibilidades(Y)=Pr(Y=1)Pr(Y=0 0)=Pr(Y=1)1-Pr(Y=1).

La expresión equivalente a la derecha muestra que es suficiente modelar para encontrar las probabilidades. Por el contrario, tenga en cuenta que podemos resolverPr(Y=1)

Pr(Y=1)=Posibilidades(Y)1+Posibilidades(Y)=1-11+Posibilidades(Y).

Regresión logística

La regresión logística modela el logaritmo de las probabilidades de como una función lineal de variables explicativas. En general, escribiendo estas variables como , e incluyendo un posible término constante en la función lineal, podemos nombrar los coeficientes (que se estimarán a partir de los datos) como y . Formalmente esto produce el modelox 1 , ... , x p β 1 , ... , β p β 0YX1,...,Xpagβ1,...,βpagβ0 0

Iniciar sesión(Posibilidades(Y))=β0 0+β1X1++βpagXpag.

Las probabilidades mismas pueden recuperarse deshaciendo el logaritmo:

Posibilidades(Y)=Exp(β0 0+β1X1++βpagXpag).

Usando variables categóricas

Las variables categóricas, como el grupo de edad, el sexo, la presencia de glaucoma, etc. , se incorporan mediante "codificación ficticia". Para mostrar que no importa cómo se codifica la variable, proporcionaré un ejemplo simple de un grupo pequeño; su generalización a múltiples grupos debería ser obvia. En este estudio, una variable es "tamaño de la pupila", con tres categorías, "Grande", "Mediano" y "Pequeño". (El estudio los trata como puramente categóricos, aparentemente sin prestar atención a su orden inherente). Intuitivamente, cada categoría tiene sus propias probabilidades, digamos para "Grande", para "Medio" y para "Pequeño" . Esto significa que, todas las demás cosas iguales,α M α SαLαMETROαS

Posibilidades(Y)=Exp(αL+β0 0+β1X1++βpagXpag)

para cualquiera en la categoría "Grande",

Posibilidades(Y)=Exp(αMETRO+β0 0+β1X1++βpagXpag)

para cualquiera en la categoría "Medio", y

Posibilidades(Y)=Exp(αS+β0 0+β1X1++βpagXpag)

para aquellos en la categoría "Pequeño".

Crear coeficientes identificables

He coloreado los dos primeros coeficientes para resaltarlos, porque quiero que note que permiten que ocurra un cambio simple: podríamos elegir cualquier número y, al sumarlo a y restarlo de cada uno de , y , no cambiaríamos ninguna probabilidad pronosticada. Esto se debe a las equivalencias obvias de la forma.β 0 α L α M α Sγβ0 0αLαMETROαS

αL+β0 0=(αL-γ)+(γ+β0 0),

etc. Aunque esto no presenta problemas para el modelo, todavía predice exactamente las mismas cosas, muestra que los parámetros no son en sí mismos interpretables. Lo que permanece igual cuando hacemos esta maniobra de suma y resta son las diferencias entre los coeficientes. Convencionalmente, para abordar esta falta de identificabilidad, las personas (y, por defecto, el software) eligen una de las categorías en cada variable como "base" o "referencia" y simplemente estipulan que su coeficiente será cero. Esto elimina la ambigüedad.

El artículo enumera primero las categorías de referencia; "Grande" en este caso. Por lo tanto, se resta de cada uno de y , y se agrega a para compensar.α L , α M , α S β 0αLαL,αMETRO,αSβ0 0

Las probabilidades de registro para un individuo hipotético que cae en todas las categorías base, por lo tanto, es igual a más un grupo de términos asociados con todas las demás "covariables", las variables no categóricas:β0 0

Probabilidades (categoría base)=Exp(β0 0+β1X1++βpagXpag).

No hay términos asociados con las variables categóricas se presentan aquí. (He cambiado ligeramente la notación en este punto: las betas ahora son los coeficientes solo de las covariables , mientras que el modelo completo incluye los para las diversas categorías).α jβyoαj

Comparación de probabilidades

Comparemos las probabilidades. Supongamos que un individuo hipotético es un

paciente masculino de 80 a 89 años con catarata blanca, sin vista fundamental y una pequeña pupila operada por un registrador especializado, ...

Asociados con este paciente (llamémoslo Charlie) se encuentran los coeficientes estimados para cada categoría: para su grupo de edad, para ser hombre, y así sucesivamente. Dondequiera que su atributo sea la base de su categoría, el coeficiente es cero por convención , como hemos visto. Debido a que este es un modelo lineal, los coeficientes se suman. Por lo tanto, para las probabilidades de registro base dadas anteriormente, las probabilidades de registro para este paciente se obtienen sumandoα machoα80-89αmasculino

α80-89+αmasculino+αsin glaucoma++αregistrador especialista.

Esta es precisamente la cantidad por la cual las probabilidades de registro de este paciente varían de la base. Para convertir de las probabilidades de registro, deshaga el logaritmo y recuerde que esto convierte la suma en multiplicación. Por lo tanto, las probabilidades básicas deben multiplicarse por

Exp(α80-89)Exp(αmasculino)Exp(αsin glaucoma)Exp(αregistrador especialista).

Estos son los números dados en la tabla bajo "OR ajustado" (odds ratio ajustado). (Se llama "ajustado" porque las covariables se incluyeron en el modelo. No juegan ningún papel en ninguno de nuestros cálculos, como verá. Se llama "relación" porque es precisamente la cantidad por cuyas probabilidades básicas se deben multiplicar para producir las probabilidades pronosticadas del paciente: consulte el primer párrafo de esta publicación.) En orden en la tabla, son , , , y así sucesivamente. Según el artículo, su producto funciona a . Por lo tantoX1,...,XpagExp(α80-89)=1,58Exp(αmasculino)=1,28Exp(αsin glaucoma)=1.0034,5

Cuotas (Charlie)=34,5×Probabilidades (Base).

(Observe que todas las categorías base tienen razones de probabilidad de , porque incluir en el producto lo deja sin cambios. Así es como puede detectar las categorías base en la tabla). 1.00=Exp(0 0)1

Reexpresar los resultados como probabilidades

Finalmente, convirtamos este resultado en probabilidades. Nos dijeron que la probabilidad prevista de referencia es . Por lo tanto, utilizando las fórmulas que relacionan las probabilidades y las probabilidades derivadas al principio, podemos calcular0,736%=0.00736

Probabilidades (Base)=0.007361-0.00736=0.00741.

En consecuencia, las probabilidades de Charlie son

Cuotas (Charlie)=34,5×0.00741=0.256.

Finalmente, convertir esto de nuevo a probabilidades da

Pr(Y(Charlie)=1)=1-11+0.256=0.204.

3
whuber: ponerse delante de mi computadora después de un día anterior muy agotador y encontrar esta respuesta extraordinaria de usted es simplemente genial. Me has ayudado mucho en una situación muy difícil. Muchas gracias. (de alguna manera @ whuber no aparecerá ...)
mahonya
Al usar nuestro sitio, usted reconoce que ha leído y comprende nuestra Política de Cookies y Política de Privacidad.
Licensed under cc by-sa 3.0 with attribution required.