Estatísticas para tabelas 2 x 2

<?xml version="1.0" encoding="ISO-8859-1"?>
<Rscript>
<titulo>Estatísticas para tabelas 2 x 2</titulo>
<descricao>

</descricao>
<autor>Ernani Bandarra</autor>
<email></email>
<versao>1.0</versao>
<parametro></parametro>
<parametro2></parametro2>
<colsel></colsel>
<colsel2></colsel2>
<script>
tabela = tabela[1:2]

if (dim(tabela)[1]!=2 || dim(tabela)[2] !=2) {
  cat(" ATENÇÃO: Este esquema de análise somente se aplica para tabelas com duas linhas e duas colunas ");q()
}

cat(titulo,"\n\n")
tabela

#Risco Relativo
r1 = tabela[1,1] / (tabela[1,1] + tabela [1,2])
r2 = tabela[2,1] / (tabela[2,1] + tabela [2,2])
rr = r1/r2
cat("\nRisco relativo     = ",rr,"\n")
dr = (r1-r2)*100
cat("Diferença de risco = ",dr,"%\n\n")

#Odds ratio
odds = (tabela[1,1] * tabela[2,2]) / (tabela[1,2] * tabela[2,1])
cat("Razão dos produtos = ",odds,"\n")

# Cálculo do Qui-QUadrado
library(ctest)
a=chisq.test(tabela)
cat("\n")
cat("Qui Quadrado       = ",a$statistic,"\n")
cat("Valor p            = ",a$p.value,"\n")
cat("Graus de liberdade = ",a$parameter,"\n")
cat("Método             = ",a$method,"\n")

#Gráfico de áreas proporcionais a freqüência
grafico.tabwin()
a=numeric(4)
a[1] = tabela[1,1]
a[2] = tabela[1,2]
a[3] = tabela[2,2]
a[4] = tabela[2,1]
b = sqrt(a)
s = sum(b)
old.par = par(no.readonly = TRUE)
par(pty="s")
plot(c(1,s),c(1,s), type = "n", xlab=names(tabela)[1], ylab=rownames(tabela)[1])
z = s/2
polygon(c(z-b[1],z,z,z-b[1]), c(z+b[1],z+b[1],z,z), col="red",    border = "red")
polygon(c(z+b[2],z,z,z+b[2]), c(z+b[2],z+b[2],z,z), col="orange", border = "orange")
polygon(c(z+b[3],z,z,z+b[3]), c(z-b[3],z-b[3],z,z), col="green",  border = "green")
polygon(c(z-b[4],z,z,z-b[4]), c(z-b[4],z-b[4],z,z), col="yellow", border = "yellow")
title(titulo)
par(old.par)

#Prova exata de Fisher
if (length(a[a < 6]) > 0 ) fisher.test(tabela)
</script>
<exemplo>

</exemplo>
</Rscript>