#Exercício 4 - Uso da Linguagem R
#Liedson Tavares de Sousa Carneiro

#4.2
##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")
factor(cervejas)
barplot(table(cervejas))
dotchart(table(cervejas))

##2 Dotchart. No entanto, ambos exibem a mesma informação.

#4.3
##1
head(caixeta)
tail(caixeta)
caixeta$dap<- caixeta$cap/pi
hist(caixeta$dap)

##2
altura<-data.frame(caixeta$h,caixeta$local)
head(altura)
chauas<-altura[altura$caixeta.local==unique(altura$caixeta.local)[1],]
jureia<-altura[altura$caixeta.local==unique(altura$caixeta.local)[2],]
retiro<-altura[altura$caixeta.local==unique(altura$caixeta.local)[3],]
par(mfrow=c(2,2))
hist(chauas$caixeta.h)
hist(jureia$caixeta.h)
hist(retiro$caixeta.h)
###Extra
boxplot(chauas$caixeta.h,jureia$caixeta.h,retiro$caixeta.h)

##3
A diferença é pouco acentuada. As três localidades possuem dados com valores extremos, produzindo outliers superiores em chauas e jureia (como visto nos boxsplot).

#4.4
egrandis<-read.table("egrandis.csv",header=T,sep=";",as.is=T)
head(egrandis)
unique(egrandis$regiao)
factor(egrandis$regiao)
regiaobox<-data.frame(egrandis$dap,egrandis$regiao)
head(regiaobox)
Salto<-regiaobox[regiaobox$egrandis.regiao==unique(regiaobox$egrandis.regiao)[1],]
Bofete<-regiaobox[regiaobox$egrandis.regiao==unique(regiaobox$egrandis.regiao)[2],]
Itatinga<-regiaobox[regiaobox$egrandis.regiao==unique(regiaobox$egrandis.regiao)[3],]
Botucatu<-regiaobox[regiaobox$egrandis.regiao==unique(regiaobox$egrandis.regiao)[4],]
par=(mfrow=c(1,1))
boxplot(Salto$egrandis.dap,Bofete$egrandis.dap,Itatinga$egrandis.dap,Botucatu$egrandis.dap)
unique(egrandis$rotacao)
rotacaobox<-data.frame(egrandis$dap,egrandis$rotacao)
um<-rotacaobox[rotacaobox$egrandis.rotacao==unique(rotacaobox$egrandis.rotacao)[1],]
dois<-rotacaobox[rotacaobox$egrandis.rotacao==unique(rotacaobox$egrandis.rotacao)[2],]
boxplot(um$egrandis.dap,dois$egrandis.dap)
hist(egrandis$h, prob=T)
curve(expr=dnorm(x,mean=mean(egrandis$h),sd=sd(egrandis$h)),add=T, col="blue")
qqnorm(egrandis$h)
qqline(egrandis$h)

#4.5
caixeta<-read.table("caixeta.csv",head=T,sep=",",as.is=T)
head(caixeta)
caixeta$especie<-as.factor(caixeta$especie)
caixeta$local<-as.factor(caixeta$local)
unique(caixeta$especie)
table(caixeta$local)
tabebuia<- data.frame(caixeta$especie, caixeta$local, caixeta$dap, caixeta$h)
head(tabebuia)
tabapenas<-tabebuia[caixeta$especie==unique(caixeta$especie)[3],]
head(tabapenas)
par(mfrow=c(2,2))
plot(caixeta.dap~caixeta.h, data=tabapenas, subset=caixeta.local=="chauas",xlab="Altura",ylab="DAP", main="Chauas")
plot(caixeta.dap~caixeta.h, data=tabapenas, subset=caixeta.local=="jureia",xlab="Altura",ylab="DAP",main="Jureia")
plot(caixeta.dap~caixeta.h, data=tabapenas, subset=caixeta.local=="retiro",xlab="Altura",ylab="DAP", main="Retiro")