Exercício módulo 2 resolvido

Olá a todos do fórum, gostaria de compartilhar meu exercício, é saber se faltou alguma coisa que foi solicitado pelo prof. Marcos
aceito sugestões de todos, forte abraço!

##############################################
####### M2_Atividades_propostas
####### Curso do R para biologia da conservação
####### Prof. Marcos Vital
####### Aluno: Julio Cesar Wojciechowski
####### Data: 9-11-2016
#############################################################################

# Para funcionar a pasta 'interative folder' é necessário rodar esta função inicial nos computadores macintosh

choose.dir <- function() {
   system("osascript -e 'tell app \"R\" to POSIX path of (choose folder with prompt \"Choose Folder:\")' > /tmp/R_folder",
         intern = FALSE, ignore.stderr = TRUE)
   p <- system("cat /tmp/R_folder && rm -f /tmp/R_folder", intern = TRUE)
   return(ifelse(length(p), p, NA))
}

setwd(choose.dir()) #Abre uma janela para se escolher a pasta 'interativy folder', ou

# 1 Setando a pasta de trabalho, conferindo, criando o objeto dados e lendo o arquivo "dados4_UTF_16.txt"

setwd("/Users/juliocesar/Dropbox/R/CursoR")
dir() # confere os arquivos da pasta
dados<-read.table("dados4_UTF_16.txt", h=T, fileEncoding="UTF-16") # Lê o arquivo de dados
summary(dados) # conferindo os dados com as estatísticas iniciais
str(dados) # conferindo a categoria dos dados
attach(dados) # Faz o R reconhecer os nomes das colunas

#####################################################################

# 2 Cálculo de média, desvio e mediana com quartis da variável Riqueza

mean(Riqueza)
sd(Riqueza)
source("f_variacao.R")
coef.var(Riqueza)
quantile(Riqueza, c(0, 0.1, 0.25, 0.5, 0.9, 1))
hist(Riqueza)

####################################################################

# 3 Histograma por tipo de área

hist(Riqueza [Área=="Desprotegido"])
dev.new() # abri uma nova janela, e evita que o R sobreponha os gráficos na mesma janela
hist(Riqueza [Área=="Sustentável"])
dev.new()
hist(Riqueza [Área=="Integral"])

# Histograma em uma janela só
par(mfrow=c(3,1))
hist(Riqueza [Área=="Desprotegido"], col="gray", las=1)
hist(Riqueza [Área=="Sustentável"], col="gray", las=1)
hist(Riqueza [Área=="Integral"], col="gray", las=1)

# 4 Boxplot geral e por tipo de área
par(mfrow=c(1,2))
boxplot(Riqueza, las=1, ylab="Riqueza")
boxplot(Riqueza ~ Área, las=1, ylab="Riqueza")

# 5 Utilizando os gráficos de médias com wiskas representando o desvio padrão

install.packages("sciplot")
library(sciplot)
source("f_meu_graf.R") # nunca esqueça de colocar a extenção .R do arquivo, se não, não vai funcionar
graf.media.desvio(Área, Riqueza, "Riqueza", "Área")

################################# fim