vglm en el paquete VGAM para variables ordinales:¿Qué significan las letras L y Q?
Pregunta
Estoy ajustando un modelo logit acumulativo usando un explicativo ordinal, los coeficientes que obtengo son los siguientes:
coefficients(cum_mod_podds_ord,matrix=TRUE)
logit(P[Y>=2]) logit(P[Y>=3])
(Intercept) -0.4106068 -1.2504301
depr_ord.L 1.3749694 1.3749694
depr_ord.Q -0.1039784 -0.1039784
Me gustaría saber qué significan la L y la Q que la función pone después del nombre de la variable.
gracias de antemano
Solución
Estos son los contrastes de un factor ordenado;contrastes basados en polinomios ortogonales.El L
es para el término polinómico "lineal", el Q
para el término cuadrático.Esto es un poco como x
y x^2
y así sucesivamente, pero sólo en términos de contrastes para el factor ordenado.
Si tiene tres niveles en el factor ordenado, los contrastes polinomiales para esos niveles son:
> zapsmall(contr.poly(3))
.L .Q
[1,] -0.7071068 0.4082483
[2,] 0.0000000 -0.8164966
[3,] 0.7071068 0.4082483
donde las filas de arriba son los 3 niveles del factor.Luego, cada nivel del factor se reemplaza por el nivel correspondiente. .L
y .Q
valores (por lo tanto, el factor se reemplaza por dos variables en la matriz del modelo).
Si necesita más información que esta, Stackoverflow probablemente no sea el lugar;intentar Validación cruzada en cambio.