Jak obliczyć iloraz szans w r (z przykładem)

W statystyce iloraz szans informuje nas o stosunku szans na wystąpienie zdarzenia w grupie badanej do szans na wystąpienie zdarzenia w grupie kontrolnej.

Często obliczamy iloraz szans podczas przeprowadzania analizy tabeli 2 na 2, która ma następujący format:

Aby obliczyć iloraz szans w R, możemy użyć funkcji oddsratio() z pakietu epitools .

Poniższy przykład pokazuje, jak zastosować tę składnię w praktyce.

Przykład: Oblicz iloraz szans w R

Załóżmy, że 50 koszykarzy korzysta z nowego programu treningowego, a 50 zawodników korzysta ze starego programu treningowego. Na koniec programu testujemy każdego gracza, aby sprawdzić, czy zdał określony test umiejętności.

Poniższa tabela przedstawia liczbę graczy, którzy przeszli i nie, w oparciu o program, którego używali:

Załóżmy, że chcemy obliczyć iloraz szans, aby porównać szanse, że gracz przejdzie test umiejętności przy użyciu nowego programu w porównaniu ze starym programem.

Oto jak utworzyć tę macierz w 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

A oto jak obliczyć iloraz szans za pomocą funkcji oddsratio() z pakietu 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”

Okazuje się, że iloraz szans wynosi 0,6045506 .

Interpretujemy to w ten sposób, że szanse, że gracz zda test przy użyciu nowego programu, są tylko 0,6045506 razy większe od szans , że gracz zda test przy użyciu starego programu.

Innymi słowy, dzięki zastosowaniu nowego programu szanse na zdanie testu przez gracza faktycznie zmniejszają się o około 39,6%.

Możemy również użyć wartości w dolnej i górnej kolumnie wyniku, aby skonstruować następujący 95% przedział ufności dla ilorazu szans:

95% przedział ufności dla ilorazu szans: [0,24, 1,48] .

Jesteśmy w 95% pewni, że prawdziwy iloraz szans pomiędzy nowym i starym programem treningowym mieści się w tym przedziale.

Kolumna midp.exact w wynikach wyświetla również wartość p powiązaną z ilorazem szans.

Okazuje się, że ta wartość p wynosi 0,271899 . Ponieważ wartość ta jest nie mniejsza niż 0,05, można stwierdzić, że iloraz szans nie jest istotny statystycznie.

Innymi słowy, z ilorazu szans wiemy, że szansa gracza na sukces przy użyciu nowego programu jest mniejsza niż szansa na sukces przy użyciu starego programu, ale różnica między tymi szansami nie jest w rzeczywistości istotna statystycznie.

Dodatkowe zasoby

Poniższe samouczki dostarczają dodatkowych informacji na temat ilorazów szans:

Iloraz szans a ryzyko względne: jaka jest różnica?
Kompletny przewodnik: jak zgłaszać iloraz szans
Jak obliczyć przedział ufności dla ilorazu szans