Hi,
I am beginner in LR subjects especially in Ordinal LR. I would like to perform ordinal logistic regression with this data:
ordinal_data <- structure(list(ODIPain = structure(c(3L, 1L, 2L, 3L, 4L, 3L,
4L, 2L, 3L, 5L, 2L, 5L, 5L, 6L, 2L, 3L, 1L, 2L, 3L, 3L, 1L, 3L,
3L, 2L, 2L, 5L, 5L, 2L, 5L, 3L, 5L, 1L, 3L, 3L, 3L, 1L, 5L, 3L,
5L, 1L, 1L, 2L, 1L, 2L, 3L, 2L, 3L, 1L, 2L, 1L, 2L, 4L, 6L, 4L,
3L, 3L, 3L, 3L, 1L, 4L, 5L, 4L, 3L, 3L, 1L, 3L, 1L, 4L, 3L, 3L,
2L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 1L, 2L, 2L, 1L, 3L, 4L, 4L, 3L,
2L, 2L, 2L, 2L, 2L, 1L, 1L, 3L, 1L, 3L, 1L, 3L, 4L, 4L, 3L, 3L,
1L, 2L, 3L, 3L, 3L, 3L, 5L, 2L, 2L), levels = c("[0] No pain",
"[1] Very mild pain", "[2] Moderate pain", "[3] Fairly severe pain",
"[4] Very severe pain", "[5] Worst imaginable pain"), class = c("ordered",
"factor")), Arm = structure(c(2L, 2L, 1L, 2L, 1L, 2L, 1L, 2L,
2L, 2L, 1L, 1L, 1L, 2L, 2L, 1L, 1L, 2L, 1L, 2L, 2L, 1L, 2L, 1L,
1L, 2L, 1L, 1L, 2L, 2L, 2L, 1L, 1L, 1L, 2L, 1L, 1L, 2L, 1L, 2L,
2L, 2L, 1L, 1L, 2L, 1L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L,
1L, 2L, 2L, 2L, 2L, 1L, 2L, 2L, 1L, 1L, 1L, 2L, 2L, 2L, 1L, 1L,
1L, 1L, 1L, 2L, 2L, 2L, 1L, 2L, 2L, 1L, 2L, 1L, 1L, 2L, 1L, 2L,
1L, 1L, 2L, 1L, 2L, 1L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 1L, 2L, 2L,
1L, 1L, 2L, 1L, 2L, 2L, 2L), levels = c("A", "B"), class = "factor"),
Age_centered = c(-6.15315315315316, 12.8468468468468, -9.15315315315316,
14.8468468468468, 12.8468468468468, 2.84684684684684, -10.1531531531532,
-18.1531531531532, -1.15315315315316, 8.84684684684684, -17.1531531531532,
13.8468468468468, 9.84684684684684, 17.8468468468468, -19.1531531531532,
-7.15315315315316, -10.1531531531532, -19.1531531531532,
-7.15315315315316, 0.846846846846844, -17.1531531531532,
5.84684684684684, -25.1531531531532, -1.15315315315316, -15.1531531531532,
4.84684684684684, 1.84684684684684, 12.8468468468468, -11.1531531531532,
5.84684684684684, -6.15315315315316, -0.153153153153156,
20.8468468468468, 5.84684684684684, -0.153153153153156, 12.8468468468468,
-19.1531531531532, -11.1531531531532, 1.84684684684684, 0.846846846846844,
-21.1531531531532, 9.84684684684684, 15.8468468468468, 14.8468468468468,
-12.1531531531532, -11.1531531531532, -9.15315315315316,
5.84684684684684, -4.15315315315316, 12.8468468468468, 1.84684684684684,
-7.15315315315316, -3.15315315315316, 7.84684684684684, 0.846846846846844,
-4.15315315315316, 5.84684684684684, -0.153153153153156,
1.84684684684684, -7.15315315315316, 1.84684684684684, -9.15315315315316,
6.84684684684684, 9.84684684684684, 17.8468468468468, 5.84684684684684,
9.84684684684684, -10.1531531531532, -5.15315315315316, 18.8468468468468,
21.8468468468468, -0.153153153153156, 2.84684684684684, -8.15315315315316,
-5.15315315315316, 5.84684684684684, 2.84684684684684, -15.1531531531532,
2.84684684684684, 25.8468468468468, -11.1531531531532, 27.8468468468468,
2.84684684684684, 20.8468468468468, -0.153153153153156, -2.15315315315316,
12.8468468468468, -0.153153153153156, 0.846846846846844,
11.8468468468468, -8.15315315315316, 3.84684684684684, 22.8468468468468,
5.84684684684684, 12.8468468468468, 4.84684684684684, 11.8468468468468,
-5.15315315315316, -17.1531531531532, -7.15315315315316,
-16.1531531531532, 0.846846846846844, -13.1531531531532,
-13.1531531531532, -19.1531531531532, -15.1531531531532,
-11.1531531531532, -4.15315315315316, -0.153153153153156,
-4.15315315315316, -7.15315315315316)), row.names = c(NA,
-111L), class = "data.frame")
So far I usually work with SPSS to do it.
I would like to use rms package or ordinal package or any other suggested and perform ordinal logistic regression.
I would like to get an output in the form of data frame not a lists if possible.
I would be very grateful for any help here.
What I have tried so far:
library(rms)
ormys <- orm(ODIPain ~ Age_centered + Arm, data = ordinal_data)
summary_orm <- summary(ormys)
Difficult part in this field is how to properly interpret the results.
best regards,
Andrzej