Odds Ratio

# odds ratio ou razão de chances quantifica a probabilidade de um evento ocorrer em um grupo em relação a outro (ou outros) grupo (s)

# pode ser usada nas mesmas cituações que o teste qui-quadrado para independência

# Será considerado um exemplo de Sampaio (2010) página 118 onde será avaliado se a prenhez de éguas depende do diluente. Neste exemplo um diluente A foi avaliado em 30 éguas sendo os resultados de prenhez (17(+)/13(-)). Um diluente B foi avaliado também em 30 éguas sendo os resultados de prenhez (25(+)/5(-)). 

# entrada dos dados (forma de matriz)

dados=matrix(c(17,25,13,5), ncol=2)

dados

# pode-se inserir nomes na matriz (tabela)

colnames(dados)=c("Prenhez_Positiva","Prenhez_Negativa")

rownames(dados)=c("Diluente_A","Diluente_B")

dados

# para obtermos a odds vamos requerer o pacote "epitools" que deve estar instalado (ver Sugestões de pacotes)

require(epitools)

# odds ratio

oddsratio(dados)

# odds estimada foi de 0,27, ou seja, uma égua inseminada com o diluente B tem 0,27 menos chance de não prenhez que uma égua inseminada com o diluente A

# no entando é mais interessante interpretar a odds com valores maiores que 1 e para isso faremos

oddsratio(dados, rev="c")

# odds estimada foi de 3,68, ou seja, tem 3,68 mais chance de prenhez quando utilizado o diluente B em relação ao diluente A

# essa seria a melhor maneira de apresentar a odds, que foi significativa (p-value<0,05)

Referências

SAMPAIO, I. B. M. Estatistica aplicada a experimentacao animal. 3nd Edition. Belo Horizonte: Editora FEPMVZ, Fundacao de Ensino e Pesquisa em Medicina Veterinaria e Zootecnia, 2010. 264p.