#6.1
conj1 <- rnorm(10, mean=6, sd=3)
conj2 <- rnorm(10, mean=7.5, sd=3.2)
source("simula.r")
sim.dif=simula(conj1, conj2, nsim=2000)
sim.maior=simula(conj2, conj1, nsim=2000, teste="uni")
t.dif <- t.test(conj1,conj2)
t.maior <- t.test(conj2, conj1, alternative="greater")

#6.2
mudas <- read.table("altura-mudas.csv", header=T, sep=",", as.is=F)
mudas
str(mudas)
tamboril <- mudas[mudas$especie=="tamboril",]
str(tamboril)
tamboril$substrato <- as.factor(tamboril$substrato)
mean.tamboril <- mean(tamboril$altura)
mean.tamboril
dif.geral <- mean.tamboril - tamboril$altura
dif.geral
sum(dif.geral)
ss.total <- sum((dif.geral)^2)
ss.total
vetor.obs <- 1:60
vetor.dados <- c(tamboril$altura)
vetor.cor <- rep(1:60, each=6)

plot(vetor.obs,vetor.dados,ylim=c(0,80),pch=(rep(c(15,16,17,18,19,20,1,2,8,9),each=6)),col=vetor.cor,ylab="Variável Resposta", xlab="Observações")
	for(i in 1:60)
	{
	lines(c(i,i),c(vetor.dados[i],mean(vetor.dados)),col=vetor.cor[i])
	}
	abline(h=mean.tamboril)

vetor.medias <- rep(mean.sust, each=6)
mean.sust <- tapply(tamboril$altura, tamboril$substrato, mean)
mean.sust
ss.intra <- sum((tamboril$altura - rep(mean.sust, each=6))^2)
ss.entre <- sum((mean.tamboril-mean.sust)^2)*6
ss.entre
gl.total <- 60-1
gl.entre <- 10-1
gl.intra <- gl.total - gl.entre 
dm.intra <- ss.intra/gl.intra
dm.entre <- ss.entre/gl.entre
razao <- dm.entre / dm.intra
razao
prob <- pf(razao, 9, 50, lower.tail=F)
prob
porc <- (ss.entre/ss.total)* 100
porc

#Variacao Intra Grupos
plot(vetor.obs,vetor.dados,ylim=c(0,90),pch=(rep(c(15,16,17,18,19,20,1,2,8,9),each=6)),col=vetor.cor,main="Variação Intra Grupos",ylab="Variável Resposta", xlab="Observações")
	for(i in 1:60)
	{
	lines(c(i,i),c(vetor.medias[i],vetor.dados[i]),col=vetor.cor[i])
	}
	lines(c(1,6),c(mean.sust[1],mean.sust[1]),col=1)
	lines(c(7,12),c(mean.sust[2],mean.sust[2]),col=2)
	lines(c(13,18),c(mean.sust[3],mean.sust[3]),col=3)
	lines(c(19,24),c(mean.sust[4],mean.sust[4]),col=4)
	lines(c(25,30),c(mean.sust[5],mean.sust[5]),col=5)
	lines(c(31,36),c(mean.sust[6],mean.sust[6]),col=6)
	lines(c(37,42),c(mean.sust[7],mean.sust[7]),col=7)
	lines(c(43,48),c(mean.sust[8],mean.sust[8]),col=8)
	lines(c(47,54),c(mean.sust[9],mean.sust[9]),col=9)
	lines(c(55,60),c(mean.sust[10],mean.sust[10]),col=10)

#Variacao entre grupos
plot(vetor.obs,vetor.medias,ylim=c(10,90),pch=(rep(c(15,16,17,18,19,20,1,2,8,9),each=6)),col=vetor.cor,main="Variação Entre Grupos",ylab="Variável Resposta", xlab="Observações")
		for(i in 1:60)
		{
		lines(c(i,i),c(vetor.medias[i],mean(vetor.medias)),col=vetor.cor[i])
		}
	abline(h=mean.tamboril)
	points(vetor.obs,vetor.dados,ylim=c(0,90),pch=(rep(c(0,1,2),each=6)),col=vetor.cor,cex=0.5)

#6.3
altura.mudas <- read.table("altura-mudas.csv", header=T, sep=",", as.is=F)
altura.mudas
tamboril <- altura.mudas[altura.mudas$especie=="tamboril",]
tamboril$substrato <- as.factor(tamboril$substrato)
aov.tamboril <- aov(tamboril$altura~tamboril$substrato)
summary(aov.tamboril)