Monise Terra Cerezini

Exercicios da aula 4

#Cervejas

Uma amostra de 30 estudantes foi indagada sobre seu tipo de cerveja preferida, com o seguinte resultado:

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")

1. Represente este resultado como um gráfico de barras e um dotplot (função dotchart). Qual tem maior razão dado/tinta? 

barplot(table(cervejas))
dotchart(table(cervejas))

#Rios

O conjunto de dado rivers está em um objeto do pacote datasets. e tem a extensão dos principais rios da América do Norte, em milhas (consulte a ajuda para detalhes). Para carregar um objeto de dados de um pacote em sua área de trabalho use a função data: 

data(rivers)

1. Qual a proporção do total de rios que tem extensão menor do que a média?

>mean(rivers<mean(rivers))
[1] 0.6666667

2. Qual é o quantil de 75%?

>quantile(rivers, prob=.75)
75% 
680 

3. Compare a média, média truncada a 25% e mediana. Há diferenças?

> mean(rivers,trim=25)
[1] 425
> mean(rivers)
[1] 591.1844
> median(rivers)
[1] 425


#Caixetais

caixeta <- read.csv("caixeta.csv", as.is=T)
str(caixeta)

1. Construa um histograma da altura dos fustes do caixetal.

hist(caixeta$h)

2. Construa histogramas da altura das árvores para os diferentes caixetais ('local').

hist(caixeta$h[caixeta$local=="chauas"])
hist(caixeta$h[caixeta$local=="retiro"])
hist(caixeta$h[caixeta$local=="jureia"])

3. Há diferenças entre as estruturas (distribuição de tamanhos) dos caixetais? 

par(mfrow=c(1,1))

histogram(~caixeta$h|caixeta$local) #Não parece haver diferença entre os caixetais


#Eucaliptos

1. Utilize o gráfico boxplot para analisar o DAP de árvores de E. grandis em função das variáveis região ('regiao') e rotação ('rot'). 

eucaliptos <- read.csv("egrandis.csv", as.is=T)
str(eucaliptos)

library(lattice)
bwplot(dap~regiao|rot,data=eucaliptos,ylab = "Dap")

2. Avalie a normalidade da altura do conjunto total de árvores com um gráfico quantil-quantil contra a distribuição normal.

qqnorm(eucaliptos$ht)
qqline(eucaliptos$ht)


#Mais Caixetais

1. Analise a relação dap-altura ('dap' e 'h') em função do caixetal, mas somente para os fustes de caixeta (Tabebuia cassinoides). 

caixeta.tabebuia=caixeta[caixeta$especie=="Tabebuia cassinoides",]
str(caixeta.tabebuia)

plot(caixeta.tabebuia$cap,caixeta.tabebuia$h,ylab="DAP", xlab="h")

2. Para a mesma relação do item anterior, verifique linearidade com a função scatter.smooth

scatter.smooth(caixeta.tabebuia$cap~caixeta.tabebuia$h,ylab="DAP", xlab="h")

3. Utilizando o pacote lattice, analise a relação dap-altura ('dap' e 'h') em função do caixetal, mas somente para os fustes de caixeta (Tabebuia cassinoides).

xyplot(cap~h|local,data=caixeta.tabebuia,ylab="DAP", xlab="h")