#
 #=====================================================================
 #                 Dados - Colírio - Crossover
 #======================================================================
 # 
 # O estudo em questão envolve pacientes com hipertensão ocular 
 # que nunca fizeram uso de colírios. O objetivo do estudo é comparar 
 # o efeito de dois colírios disponíveis no mercado: Betoptic e Timoptol.
 # 
 # O estudo foi conduzido de forma longitudinal. As seguintes variáveis 
 # foram medidas ao nível do  paciente: Sexo e Idade (em anos). 
 # E ao nível do episódio foram medidas as seguintes variáveis:
 #	           IOP: Pressão Intra Ocular (mm Hg);
 #
 # O estudo foi conduzido da seguinte forma para cada paciente: 
 # inicialmente a  variável IOP ao nível do episódio é medida, 
 # a seguir um colírio é escolhido de forma aleatória e submetido ao 
 # paciente pelo período de dois meses. Novamente a mesma variável 
 # ao nível do episódio é medida no paciente. Dois meses são 
 # concedidos ao paciente para limpar o efeito do colírio (período de 
 # “washout”) e em seguida é aplicado o outro colírio pelo mesmo período. 
 # Novamente ao final de dois meses a IOP é medida. 
 #  Trinta e dois pacientes foram utilizados no estudo. E a ordem foi
 #   0 - BETOP e a seguir TIMO
 #   1 - TIMO e a seguir BETOP
 #           
 #=======================================================================
 #           Leitura e Preliminares dos Dados 
 #========================================================================
 dados<-read.table("colirio_page.txt",dec=",",h=T) 
 names(dados)
 dim(dados)
 dados[1:15,]
 attach(dados) # o risco é todo seu....
#
#                  Descritiva dos dados
# 
 summary(IOP_sem)
 summary(IOP_betop)
 summary(IOP_timo)
 boxplot(IOP_sem, IOP_betop, IOP_timo,ylab="pressão ocular (mmHg)")
 axis(1, 1:3, c("sem","betop","timol"))
## Cuidado pois o Boxplot não considera o pareamento ####
#
#====================================================================
#       Suposição: Não Existe Efeito da Ordem
#       Teste t - Pareado
#====================================================================
 dif1<-IOP_sem - IOP_betop  
 dif2<-IOP_sem - IOP_timo
 dif3<-IOP_betop - IOP_timo
 dif<-dif3
 summary(dif)
 boxplot(dif)
# verificacao de normalidade           
  shapiro.test(dif)
  plot(qqnorm(dif),ylab="quantis")
# teste t para diferencas           
  t.test(dif)
#help.search("wilcoxon")
  wilcox.test(dif)
#=====================================================================