ATIVIDADE PROPOSTA DO MÓDULO#2 - TURMA 2018.2

Ir em baixo

ATIVIDADE PROPOSTA DO MÓDULO#2 - TURMA 2018.2

Mensagem por Cleyton Tenório em Qui Out 18, 2018 1:22 pm

Boa tarde, pessoal. Estou criando esse tópico para discutirmos e compartilharmos nossas ideias! Se alguém conseguiu desenvolver algumas perguntas, ou outras atividades através da atividade, compartilhem por gentileza! confused
Código:
##############################
#Gráficos de média e desvio  #
#Cleyton Tenório Barros      #
#R para Bio da conservação   #
##############################

#Carregando dados:

setwd("D:/R/Rbio")
dir()
dados<-read.table("dados4.txt", header=T)
summary(dados)

#Ver os dados
View(dados)

######################################
#Calculando a media
mean(dados$Riqueza)
#Calculando o Desvio
sd(dados$Riqueza)
#Calculando a mediana
median(dados$Riqueza)
#Calculando os quartis
quantile(dados$Riqueza)
######################################

#Histograma da variavel Riqueza
hist(dados$Riqueza, col = "grey")

#histograma por tipo de area
#Area Desprotegido
hist(dados$Riqueza[dados$Area=="Desprotegido"])
#Area Integral
hist(dados$Riqueza[dados$Area=="Integral"])
#Area Sustentavel
hist(dados$Riqueza[dados$Area=="Sustentavel"])
#####################################

#Boxplot por tipo de Area:
#Area Desprotegido
boxplot(dados$Riqueza[dados$Area=="Desprotegido"])
#Area Integral
boxplot(dados$Riqueza[dados$Area=="Integral"])
#Area Sustentavel
boxplot(dados$Riqueza[dados$Area=="Sustentavel"])
######################################

#Grafico de medias e desvio
medias<-tapply(dados$Riqueza, dados$Area, mean)

desvio<-tapply(dados$Riqueza, dados$Area, sd)

plot(1:3, medias, ylim = range(c(medias-desvio, medias+desvio)), pch=16)
arrows(1:3, medias-desvio, 1:3, medias+desvio, angle = 90, code = 3 )

#Será que a Riqueza de especie varia em relação ao tipo de Area?
plot(dados$Riqueza~dados$Area)
Pessoal. Consegui fazer a atividade, entender como constrói cada gráfico, verifica às médias, mediana, quartis e tal, só que não entendo muito bem a leitura dos dados e sua interpretação! Existe algumas dicas ou técnicas para leitura?
avatar
Cleyton Tenório

Mensagens : 7
Data de inscrição : 18/09/2018
Idade : 35
Localização : Bom Conselho - PE

Ver perfil do usuário

Voltar ao Topo Ir em baixo

Re: ATIVIDADE PROPOSTA DO MÓDULO#2 - TURMA 2018.2

Mensagem por Giovana em Sex Out 19, 2018 3:07 am

Código:
################################
#Atividade proposta do módulo 2#
#Giovana                       #
#R para Bio da conservação     #
################################


#Determinar a pasta de trabalho:
setwd("C:/Users/Giovana/Documents/Mestrado/Marcos Vital - R para Biologia da conservação/Módulo 2/Módulo 2 - Introdução à Bioestatística no R")
getwd()     #Mostra qual pasta é a pasta de trabalho


#Conferir dados na pasta:
dir()


#Criar o objeto dados4:
dados4<-read.table("dados4.txt", header = T)
ls()    #Objeto salvo na memória do R


#Resumo do objeto:
summary(dados4)


#Chamar um conjunto de dados (janela):
View(dados4)


#Calcular a média da variável "Riqueza":
mean(dados4$Riqueza)    #Resposta: 13.84444
mean(dados4$Riqueza[dados4$Area=="Integral"])     #Resposta: 13.4
mean(dados4$Riqueza[dados4$Area=="Sustentavel"])  #Resposta: 14.93333
mean(dados4$Riqueza[dados4$Area=="Desprotegido"]) #Resposta: 13.2


#Calcular o desvio da variável "Riqueza":
sd(dados4$Riqueza)    #Resposta: 3.988367


#Calcular a mediana com quartis da variável "Riqueza":
quantile(dados4$Riqueza)    #Resposta: 0%  25%  50%  75% 100%
                            #Resposta: 6   11   14   17   20

#Criar um histograma da variável "Riqueza":
hist(dados4$Riqueza, col = "pink", main = "", xlab = "Riqueza", ylab = "Frequência", las=1, box())


#Criar um histograma por tipo de área:
hist(dados4$Riqueza[dados4$Area=="Desprotegido"], col = "pink", main = "Área desprotegida", xlab = "Riqueza", ylab = "Frequência", las=1)
hist(dados4$Riqueza[dados4$Area=="Integral"], col = "pink", main = "Área integral", xlab = "Riqueza", ylab = "Frequência", las=1)
hist(dados4$Riqueza[dados4$Area=="Sustentavel"], col = "pink", main = "Área sustentável", xlab = "Riqueza", ylab = "Frequência", las=1)
box()

#Criar um boxplot que mostre a riqueza por tipo de área:
boxplot(dados4$Riqueza ~ dados4$Area, col = "pink", main = "Riqueza de espécies", xlab = "Área", ylab = "Riqueza", las=1 )


#Criar um gráfico de média que mostre a riqueza por tipo de área:
medias<-tapply(dados4$Riqueza, dados4$Area, mean)
desvios<-tapply(dados4$Riqueza, dados4$Area, sd)
plot(1:3, medias, pch=8, ylim=range(c(medias-desvios, medias+desvios)))
arrows(1:3, medias-desvios, 1:3, medias+desvios, angle=90, code=3)

source("plotmÇdias.R")
plotmédias(dados4$Area, dados4$Riqueza, "Tipo de Área", "Riqueza")

Giovana

Mensagens : 5
Data de inscrição : 21/09/2018

Ver perfil do usuário

Voltar ao Topo Ir em baixo

Voltar ao Topo

- Tópicos similares

 
Permissão deste fórum:
Você não pode responder aos tópicos neste fórum