#------------------------------------------------------------------------
#  Script exercício 6 - Capítulo 7
#------------------------------------------------------------------------
#  Exerc. 01
#
 t<-(3160-2829)/sqrt(((453^2)/19)+((708^2)/36))
 t
#
 gl<-(((453^2)/19)+((708^2)/36))^2/((((453^2)/19)^2/18)+
     (((708^2)/36)^2/35))
 gl
#
 2*(1-pt(t,gl))
#
# Intervalo de 95% de Confiança
#
 limsup<- (3160-2829)+(qt(0.975,gl)*sqrt(((453^2)/19)+((708^2)/36)))
 liminf<- (3160-2829)-(qt(0.975,gl)*sqrt(((453^2)/19)+((708^2)/36)))
 limsup # 647
 liminf # 15
#--------------------------------------------------------------------------
#       Exer. 03 
#   Para Faixa Etária de 40 a 50 anos 
#        letras a) e b)
#
 dados<-matrix(c(96,108,124,94),nc=2) # Mulheres 40-50 anos
 dados
 Qp<-chisq.test(dados,correct=F) # teste usual
 Qp
 OR<-(dados[1,1]*dados[2,2])/(dados[1,2]*dados[2,1])
 OR
 vf<-(1/dados[1,1])+(1/dados[1,2])+(1/dados[2,1]+(1/dados[2,2]))
 vf
 dpf<-sqrt(vf)
 dpf
 z<-qnorm(0.975)
 li<-exp(log(OR)-z*dpf)
 li # 0.46
 ls<-exp(log(OR)+z*dpf)
 ls # 0.99
#
#--------------------------------------------------------------------------
#       Exer. 03 
#   Para Faixa Etária superior a 50 anos
#        letra c)
#
 dados<-matrix(c(16,9,20,8),nc=2) # Mulheres >50 anos
 dados
 Qp<-chisq.test(dados,correct=F) # teste usual
 Qp
 OR<-(dados[1,1]*dados[2,2])/(dados[1,2]*dados[2,1])
 OR
 vf<-(1/dados[1,1])+(1/dados[1,2])+(1/dados[2,1]+(1/dados[2,2]))
 vf
 dpf<-sqrt(vf)
 dpf
 z<-qnorm(0.975)
 li<-exp(log(OR)-z*dpf)
 li # 0.22
 ls<-exp(log(OR)+z*dpf)
 ls # 2.26
#