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

summary(caixeta)
names(caixeta)
unique(caixeta$local)

caixeta$dap <- caixeta$cap/pi

hist(caixeta$dap)


par(mfrow=c(2,2))
hist(caixeta$dap[caixeta$local == "chauas"])
hist(caixeta$dap[caixeta$local == "jureia"])
hist(caixeta$dap[caixeta$local == "retiro"])
dev.off()