Estoy llevando a cabo un estudio de una serie de pacientes con una enfermedad, y utilizando una escala ordinal de evaluación del estado funcional en 3 momentos diferentes.
En este ejemplo podemos decir que tengo 100 pacientes, y cada uno de ellos se mide en los tres puntos temporales (sin necesidad de imputación).
Véase la representación gráfica del cambio en las frecuencias de cada uno de estos 6 valores ordinales a lo largo de los tres puntos temporales (el grupo superior no tiene pacientes con puntuación ordinal 3,5,6):
Quiero comparar la diferencia global entre el punto temporal 1 y el punto temporal 2, el punto temporal 1 y el punto temporal 3, y el punto temporal 2 y el punto temporal 3 para ver si los cambios observados son significativos.
Además, si hay una forma de determinar los cambios entre los niveles ordinales específicos, también puede ser útil.
Estoy planeando hacer el análisis en R, pero me gustaría asegurarme de que estoy utilizando la prueba adecuada. ¿Sería apropiada alguna forma de regresión ordinal? ¿O es mejor alguna versión de la suma de rangos de Wilcoxin para este paradigma de estudio?
Gracias por su ayuda.
* ACTUALIZACIÓN CON CONJUNTO DE DATOS
Todavía he estado trabajando para analizar esto de la mejor manera.
El código de mi conjunto de datos figura a continuación como referencia:
mrsstat <-tibble(
Patient = c(1:21),
pMRS = c(0, 0, 0, 0, 0, 0, 4, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 2, 0),
dMRS = c(4, 2, 2, 4, 3, 6, 4, 1, 5, 0, 3, 0, 6, 1, 5, 3, 3, 1, 3, 4, 3),
fMRS = c(1, 2, 2, 3, 2, 6, 3, 1, 2, 0, 3, 0, 6, 1, 5, 2, 0, 1, 3, 4, 0))
mrsstat<- mrsstat %>%
pivot_longer(cols = c("pMRS", "dMRS", "fMRS"), names_to = "timepoint") %>%
mutate(timepoint = factor(timepoint,
levels= c("fMRS",
"dMRS",
"pMRS")),
value= factor(value, levels = c("0", "1","2","3","4","5","6")),
Patient= factor(Patient, levels = c(1:21)))
> mrsstat
# A tibble: 63 × 3
Patient timepoint value
<fct> <fct> <fct>
1 1 pMRS 0
2 1 dMRS 4
3 1 fMRS 1
4 2 pMRS 0
5 2 dMRS 2
6 2 fMRS 2
7 3 pMRS 0
8 3 dMRS 2
9 3 fMRS 2
10 4 pMRS 0
# … with 53 more rowsIntenté hacer un análisis de regresión ordinal con métodos utilizando el paquete ordinal et al. en R como se describe aquí:
Regresión ordinal bidireccional con CLM
No estoy seguro de que esto sea factible con mi conjunto de datos, ya que cada paciente tiene sólo 1 registro de la puntuación del resultado funcional (variable etiquetada "valor") en cada punto de tiempo, por lo que cuando se intenta llegar a un análisis ANOVA de los datos como se ve en el enlace anterior para la regresión ordinal de dos vías, no funcionará ya que obviamente no hay varianza con una medida por paciente por punto de tiempo.
Sigo pensando que debería haber alguna prueba para evaluar los cambios de significación entre las proporciones de "valor" entre los 3 puntos temporales y para las comparaciones múltiples entre los grupos ordinales de "valor".
