#Exercícios 4 - Análises exploratórias

#4.2 - Cervejas

#1

cervejas <-c("chope","lata","garrafa","chope","garrafa", "garrafa","lata","lata","nenhuma","lata","garrafa","garrafa", "garrafa","lata","lata","lata","garrafa","lata","chope","nenhuma", "garrafa","garrafa","garrafa","chope","garrafa","garrafa","chope","garrafa","lata","lata")
cervejas
table(cervejas)
barplot(table(cervejas))
dotchart(as.matrix(table(cervejas))) #versão antiga do r


#2
# O dotchart tem a maior razão dado/tinta.



#4.3 - Caixetas


caixeta<-read.table("caixeta.csv", header=T, sep=",", as.is=T)
caixeta 

#1

hist(caixeta$cap/pi)

#2

hist(caixeta$h[caixeta$local=="retiro"])
hist(caixeta$h[caixeta$local=="chauas"])
hist(caixeta$h[caixeta$local=="jureia"])
par(mfrow=c(1,3))

#3

#Sim, há diferenças entre as estruturas dos caixetais. No local "jureia" a distribuição 
é normal, no "chauas" a distribuição é unicaldal e no "retiro" podemos considerar a
distribuição também normal. 


#4.4 - Eucaliptos

#1

egrandis$rotacao<-factor(egrandis$rotacao)
egrandis$regiao<-factor(egrandis$regiao)

par(mfrow=c(1,2))
boxplot(egrandis$dap[egrandis$rotacao==1]~egrandis$regiao[egrandis$rotacao==1])
boxplot(egrandis$dap[egrandis$rotacao==2]~egrandis$regiao[egrandis$rotacao==2])


#2

qqnorm(egrandis$h)
qqline(egrandis$h)


#4.5 - Mais caixetas


#1

caixeta$areabasal=pi*((caixeta$cap/2*pi)^2)
caixeta=caixeta[caixeta$especie=="Tabebuia cassinoides",]
caixeta2=aggregate(caixeta$areabasal, by=list(caixeta$arvore, caixeta$h, caixeta$local),sum)
colnames(caixeta2)=c("arvore", "h","local","area.basal")
caixeta2$dap=sqrt(caixeta2$area.basal/pi)*2
caixeta2

plot(caixeta2$dap[caixeta2$local=="chauas"]~caixeta2$h[caixeta2$local=="chauas"])
points(x=caixeta2$h[caixeta2$local=="jureia"],y=caixeta2$dap[caixeta2$local=="jureia"], col="blue")
points(x=caixeta2$h[caixeta2$local=="retiro"],y=caixeta2$dap[caixeta2$local=="retiro"], col="red")
#ou
par(mfrow=c(2,2))
plot(caixeta2$dap[caixeta2$local=="chauas"]~caixeta2$h[caixeta2$local=="chauas"])
plot(caixeta2$dap[caixeta2$local=="retiro"]~caixeta2$h[caixeta2$local=="retiro"])
plot(caixeta2$dap[caixeta2$local=="jureia"]~caixeta2$h[caixeta2$local=="jureia"])

#2

with(caixeta2[caixeta2$local=="retiro",],scatter.smooth(x=caixeta2$h[caixeta2$local=="retiro"],y=caixeta2$dap[caixeta2$local=="retiro"]))
with(caixeta2[caixeta2$local=="chauas",],scatter.smooth(x=caixeta2$h[caixeta2$local=="chauas"],y=caixeta2$dap[caixeta2$local=="chauas"]))
with(caixeta2[caixeta2$local=="jureia",],scatter.smooth(x=caixeta2$h[caixeta2$local=="jureia"],y=caixeta2$dap[caixeta2$local=="jureia"]))

#3

library(lattice)
xyplot(caixeta2$dap~caixeta2$h|caixeta2$local)