Calculo do Odds Ratio

<?xml version="1.0" encoding="ISO-8859-1"?>
<Rscript>
<titulo>Calculo do Odds Ratio</titulo>
<descricao>
Este script calcula o odds ratio e seu intervalo de confiançade 95% 
para uma tabela 2x2.

A funçao foi baseada em uma de autoria de Frank Harrell
http://biostat.mc.vanderbilt.edu/twiki/bin/view/Main/RS
</descricao>
<autor>Oswaldo G Cruz</autor>
<email>oswaldo@fiocruz.br</email>
<versao>1.0</versao>
<parametro></parametro>
<parametro2></parametro2>
<colsel></colsel>
<colsel2></colsel2>
<script>
library(stats)
odds.ratio <- function(x,  conf=0.95) {
	tmp <- list()
	
	if (dim(x)[1]!=2 || dim(x)[2] !=2) {
		stop(" x tem de ser Matrix 2x2 ")
		}	
	
	odds <- x[1,1] * x[2,2] / ( x[2,1] * x[1,2] )
	ASE <- sqrt(sum(1/x))
	IC <- exp(log(odds) + c(-1,1) * qnorm(0.5*(1+conf)) *ASE )
	
	tmp <- list(or=odds,ASE=ASE,intervalo=IC,conf=conf,tab=x)
		
	class(tmp) <- "or"
	tmp
}

print.or <-  function(x)
{
	cat (titulo,"\n")
	cat (subtitulo,"\n")
	cat ("\n\n  Odds Ratio  \n\n")
	print(x$tab) 
	cat ("\n valor do Odds Ratio =", x$or,"\n")
	cat ("\n Intervalo de Confiança (",x$conf,"%) \t",x$intervalo[1],"   a   ",x$intervalo[2],"\n\n",sep="" )
	
}


odds.ratio(tabela)
</script>
<exemplo>

</exemplo>
</Rscript>