########Exercicio 6.1#############
###Crie seus dados
conj1<-rnorm(10, mean=6, sd=3)
conj2<-rnorm(10, mean=7.5, sd=3.2)
source("simula.r")
dif<-mean(conj1)-mean(conj2)

#a diferenças absoluta entre as médias é 1.311, há uma diferença 
#entre as média, mas não sabemos se essa diferença é gerada pelo acaso.

#gerando uma simulaçao

sim.dif<-simula(conj1,conj2, nsim=1000)
dados.maior<-sum(sim.dif$diferencas>=dif)
dados.menor<-sum(sim.dif$diferencas<=dif)
pvalue<-(dados.maior/length(sim.dif$diferencas))+
  (dados.menor/length(sim.dif$diferencas))
#as amostras não são diferentes

#teste de médias
sim.maior<-simula(conj2,conj1,nsim=1000,"uni")
pvalue2<-(sum(sim.maior$diferencas>=dif))/length(sim.maior$diferencas)
#testt
t.dif<-t.test(conj1,conj2,alternative="two.sided")
t.maior<-t.test(conj2,conj1,alternative="greater")

#não existe diferenças entre as amostras

####Exercício 6.2#########

mudas=read.table("altura-mudas.csv", header=T, sep=",")
str(mudas)
head(mudas)
tail(mudas)
mudas$substrato
mudas$especie
mudas$bloco
mudas$altura
#1.
#a variável medida é a altura das especies. Como fator temos as especies
# e o substrato.
#temos duas réplicas do fator especie e dez réplicas do fator substratos
#transformando substrato em fator
substrato<-as.factor(mudas$substrato)
mudas$substrato<-substrato
str(mudas)

#ANOVA na unha
tamboril<-mudas[mudas$especie=="tamboril",]
tamboril
#somatório dos desvios quadraticos geral dos gréficos
media.altura=mean(tamboril$altura)
dif.geral<-tamboril$altura-media.altura
ss.total<-sum(dif.geral^2)
#media de altura de cada substrato
media.substrato<-c(tapply(tamboril$altura,INDEX=tamboril$substrato, 
                                FUN=mean))
#Somatório do desvio quadrático intragrupo
str(tamboril)
head(tamboril)
tamboril
vetor.media.substrato<-rep(media.substrato,each=6)
tamboril$media.grupo<-vetor.media.substrato
ss.intra=sum((tamboril$altura-tamboril$media.grupo)^2)

#somatório do desvio quadrático entregrupo
media.substrato
media.altura
ss.entre=6*sum((media.substrato-media.altura)^2)

#Calculo do desvio médio intra grupos dividindo o ss.intra 
#pelo gl intra, e o desio médio entre grupos dividindo o ss.entre 
#pelo gl entre.
gl.total=60-1
gl.entre=10-1
gl.intra=gl.total-gl.entre
dm.entre=ss.entre/gl.entre
dm.intra=ss.intra/gl.intra
#Calculo da razão de variancia F
razao=dm.entre/dm.intra
#Probabilidade da diferença
prob=pf(razao,9,50, lower.tail=FALSE)
prob
porc=100*(ss.entre/ss.total)

####Exercício 6.3#############

mudas<-read.table("altura-mudas.csv", head=T, sep=",")
head(mudas)
tamboril<-mudas[mudas$especie=="tamboril",]
tamboril$substrato<-as.factor(tamboril$substrato)
str(tamboril)
aov.tamboril<-aov(tamboril$altura~tamboril$substrato)
aov.tamboril
summary(aov.tamboril)