Как рассчитать коэффициенты шансов в r (с примером)

В статистике отношение шансов показывает нам отношение шансов события, произошедшего в экспериментальной группе, к шансам события, произошедшего в контрольной группе.

Мы часто рассчитываем отношение шансов при выполнении анализа таблицы 2 на 2, которая имеет следующий формат:

Чтобы вычислить отношение шансов в R, мы можем использовать функцию Oddsratio() из пакета epitools .

В следующем примере показано, как использовать этот синтаксис на практике.

Пример: вычислить отношение шансов в R

Допустим, 50 баскетболистов используют новую программу тренировок, а 50 игроков — старую программу тренировок. В конце программы мы проверяем каждого игрока, чтобы увидеть, прошел ли он определенный тест на навыки.

В следующей таблице показано количество игроков, сдавших и не сдавших экзамен, в зависимости от используемой ими программы:

Допустим, мы хотим рассчитать отношение шансов, чтобы сравнить шансы игрока пройти тест на навыки с использованием новой программы и с использованием старой программы.

Вот как создать эту матрицу в R:

 #create matrix
program <- c(' New Program ', ' Old Program ')
outcome <- c(' Pass ', ' Fail ')
data <- matrix(c(34, 16, 39, 11), nrow= 2 , ncol= 2 , byrow= TRUE )
dimnames(data) <- list(' Program '=program, ' Outcome '=outcome)

#view matrix
data

             Outcome
Program Pass Fail
  New Program 34 16
  Old Program 39 11

А вот как вычислить отношение шансов с помощью функции Oddsratio() из пакета epitools :

 install. packages (' epitools ')

library (epitools)

#calculate odds ratio
oddsratio(data)

$measure
             odds ratio with 95% CI
Program estimate lower upper
  New Program 1.0000000 NA NA
  Old Program 0.6045506 0.2395879 1.480143

$p.value
             two-sided
Program midp.exact fisher.exact chi.square
  New Program NA NA NA
  Old Program 0.271899 0.3678219 0.2600686

$correction
[1] FALSE

attr(,"method")
[1] “median-unbiased estimate & mid-p exact CI”

Отношение шансов оказывается равным 0,6045506 .

Мы интерпретируем это так, что шансы игрока пройти тест с использованием новой программы всего в 0,6045506 раза превышают шансы игрока пройти тест с использованием старой программы.

Другими словами, благодаря использованию новой программы шансы игрока пройти тест фактически снижаются примерно на 39,6%.

Мы также можем использовать значения в нижнем и верхнем столбцах результата, чтобы построить следующий 95% доверительный интервал для отношения шансов:

95% доверительный интервал для отношения шансов: [0,24, 1,48] .

Мы на 95% уверены, что истинное соотношение шансов между новой и старой программой обучения находится в пределах этого интервала.

Столбец Midp.exact в выходных данных также отображает значение p, связанное с отношением шансов.

Это значение p оказывается равным 0,271899 . Поскольку это значение не менее 0,05, мы можем заключить, что отношение шансов не является статистически значимым.

Другими словами, из отношения шансов мы знаем, что шансы игрока на успех при использовании новой программы ниже, чем шансы на успех при использовании старой программы, но разница между этими шансами на самом деле не является статистически значимой.

Дополнительные ресурсы

Следующие руководства предоставляют дополнительную информацию об отношениях шансов:

Отношение шансов и относительный риск: в чем разница?
Полное руководство: как сообщить о коэффициентах шансов
Как рассчитать доверительный интервал для отношения шансов