Estoy haciendo un modelo de regresión logística ordinal en R, pero no sé por qué me cambia los nombres de las variables si pongo una variable ordinal en las varibles independientes de la fórmula al usar polr del paquete MASS. (Supongo tiene sentido que cambien los valores cuando la variable es ordinal):
Por ejemplo con este código, al poner la variable clarity como categórica:
data(diamonds)
muestra_diamantes <- diamonds[1:100,]
muestra_diamantes$cut <- factor(muestra_diamantes$cut, levels=c('Fair', 'Good', 'Very Good', 'Premium', 'Ideal'), ordered=TRUE)
muestra_diamantes$clarity <- factor(muestra_diamantes$clarity, levels=c('I1', 'SI2', 'SI1', 'VS2', 'VS1', 'VVS2', 'VVS1', 'IF'), ordered=FALSE)
muestra_lro <- MASS::polr(formula=cut ~ x + y + z + clarity, data=muestra_diamantes, Hess=TRUE, method="logistic")
exp(coef(summary(muestra_lro)))
Coeficientes exponenciados:
| | Value | Std Error | t |
|-------------------|-------------|------------|------------|
| x | 0.127456184 | 124.139126 | 0.65229528 |
| y | 12.95810574 | 114.487089 | 1.71668465 |
| z | 0.196285788 | 14.70743 | 0.54572311 |
| claritySI2 | 0.63312405 | 4.733888 | 0.7452797 |
| claritySI1 | 0.512574993 | 4.579121 | 0.64452544 |
| clarityVS2 | 0.28127597 | 4.76252 | 0.44366549 |
| clarityVS1 | 0.20029496 | 4.718363 | 0.35472276 |
| clarityVVS2 | 0.454695091 | 4.88625 | 0.60847767 |
| clarityVVS1 | 0.432025719 | 4.807342 | 0.58595175 |
| Fair|Good | 0.001374046 | 11.163314 | 0.06512454 |
| Good|Very Good | 0.012629681 | 10.195421 | 0.15216611 |
| Very Good|Premium | 0.072818199 | 9.860981 | 0.3183131 |
| Premium|Ideal | 0.221585212 | 9.815635 | 0.51695916 |
Al poner la variable clarity como ordinal:
# muestra_diamantes$clarity <- factor(muestra_diamantes$clarity, levels=c('I1', 'SI2', 'SI1', 'VS2', 'VS1', 'VVS2', 'VVS1', 'IF'), ordered=TRUE)
| | Value | Std Error | t |
|-------------------|-------------|-------------|-----------|
| x | 0.123991625 | 123.701862 | 0.6483718 |
| y | 13.31223779 | 113.894553 | 1.7275099 |
| z | 0.196626228 | 14.693605 | 0.5459594 |
| clarity.L | 0.001855322 | 3427798.281 | 0.6583678 |
| clarity.Q | 0.008909861 | 3816358.193 | 0.7323537 |
| clarity.C | 0.009782939 | 158278.5153 | 0.6794367 |
| clarity^4 | 0.034829154 | 2709.986301 | 0.6539515 |
| clarity^5 | 0.191235962 | 66.482591 | 0.6742492 |
| clarity^6 | 0.807139755 | 5.494582 | 0.8818299 |
| Fair|Good | 0.011340425 | 43.987691 | 0.3061142 |
| Good|Very Good | 0.104355146 | 41.774658 | 0.5457936 |
| Very Good|Premium | 0.601682289 | 41.716317 | 0.8726966 |
| Premium|Ideal | 1.830864643 | 41.823995 | 1.1758496 |
¿Alguien sabe por qué pasa? Gracias.