Como calcular odds ratio em r (com exemplo)
Nas estatísticas, uma razão de probabilidades diz-nos a razão entre as probabilidades de um evento ocorrer num grupo de tratamento e as probabilidades de um evento ocorrer num grupo de controlo.
Freqüentemente calculamos uma razão de chances ao realizar análises em uma tabela 2 por 2, que assume o seguinte formato:
Para calcular uma razão de chances em R, podemos usar a função oddsratio() do pacote epitools .
O exemplo a seguir mostra como usar essa sintaxe na prática.
Exemplo: Calcule uma razão de chances em R
Digamos que 50 jogadores de basquete estejam usando um novo programa de treinamento e 50 jogadores estejam usando um programa de treinamento antigo. No final do programa, testamos cada jogador para ver se ele passa em determinado teste de habilidade.
A tabela a seguir mostra o número de jogadores que passaram e falharam, com base no programa que utilizaram:
Digamos que queremos calcular uma razão de chances para comparar as chances de um jogador passar no teste de habilidades usando o novo programa versus usando o programa antigo.
Veja como criar esta matriz em 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
E aqui está como calcular a razão de chances usando a função oddsratio() do pacote 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”
A razão de chances acabou sendo 0,6045506 .
Interpretamos isso como significando que as chances de um jogador passar no teste usando o novo programa são apenas 0,6045506 vezes as chances de um jogador passar no teste usando o programa antigo.
Ou seja, as chances de um jogador passar no teste são reduzidas em aproximadamente 39,6% com o uso do novo programa.
Também podemos usar os valores nas colunas inferior e superior do resultado para construir o seguinte intervalo de confiança de 95% para a razão de chances:
Intervalo de confiança de 95% para a razão de chances: [0,24, 1,48] .
Temos 95% de certeza de que a verdadeira razão de probabilidade entre o novo e o antigo programa de treinamento está contida neste intervalo.
A coluna midp.exact na saída também exibe o valor p associado à razão de chances.
Esse valor p acaba sendo 0,271899 . Dado que este valor não é inferior a 0,05, concluiríamos que o odds ratio não é estatisticamente significativo.
Por outras palavras, sabemos pela razão de probabilidades que a probabilidade de sucesso de um jogador utilizando o novo programa é menor do que a probabilidade de sucesso utilizando o programa antigo, mas a diferença entre estas probabilidades não é estatisticamente significativa.
Recursos adicionais
Os tutoriais a seguir fornecem informações adicionais sobre odds ratio:
Razão de chances versus risco relativo: qual a diferença?
O guia completo: como relatar taxas de probabilidade
Como calcular um intervalo de confiança para uma razão de chances