# *******************************************************************
# 1. Definição da trilha de dados
# *******************************************************************
setwd("C:\\dados\\")


# *******************************************************************
# 2. Leitura dos dados
# *******************************************************************
X<-read.table("C:\\dados\\plano13.dat",header=F)
# Listar dados
X


# *******************************************************************
# 3. Estatísticas descritivas
# *******************************************************************
summary(X)


# *******************************************************************
# 4.1 Teste t  - comparação entre médias
# *******************************************************************
variavel<-(X[,2])
npop<-max(X$V1)
fim<-npop-1
for(i in 1:fim)
{
inicio<-(i+1)
for(j in inicio:npop)
{
pop1 <-i
pop2 <-j
cat("populações : ",cbind(i,j))
print(t.test(variavel[X[,1]==pop1], variavel[X[,1]==pop2], alternative="two.sided"))
}
}


# *******************************************************************
# 4.2 Teste t  - comparação entre variâncias
# *******************************************************************
variavel<-(X[,2])
npop<-max(X$V1)
fim<-npop-1
for(i in 1:fim)
{
inicio<-(i+1)
for(j in inicio:npop)
{
pop1 <-i
pop2 <-j
cat("populações : ",cbind(i,j))
print(var.test(variavel[X[,1]==pop1], variavel[X[,1]==pop2], alternative="two.sided"))
}
}


# *******************************************************************
# 4.3 Teste t  - comparação entre variâncias e médias
# *******************************************************************
alpha <- 0.05
variavel<-(X[,2])
npop<-max(X$V1)
fim<-npop-1
for(i in 1:fim)
{
inicio<-(i+1)
for(j in inicio:npop)
{
pop1 <-i
pop2 <-j
cat("populações : ",cbind(i,j))
if (variancia.teste$p.value > alpha)
media.teste <- t.test(variavel[X[,1]==pop1], variavel[X[,1]==pop2], alternative="two.sided", var.equal = TRUE, conf.level = 1-alpha)
else
media.teste <-t.test(variavel[X[,1]==pop1], variavel[X[,1]==pop2], alternative="two.sided", var.equal = FALSE, conf.level = 1-alpha)
print(media.teste)
}
}


# *******************************************************************
# 5. Teste t  usando o pacote gregmisc
# *******************************************************************
library(gregmisc)
C<-combinations(max(X$V1),2)
for(i in 1:nrow(C))
{
}