# How to perform ordinal LR in R ?

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

This topic was automatically closed 42 days after the last reply. New replies are no longer allowed.

If you have a query related to it or one of the replies, start a new topic and refer back with a link.