ATIVIDADE PROPOSTA DO MÓDULO#2 - TURMA 2018.2

por Cleyton Tenório 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?

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

ATIVIDADE PROPOSTA DO MÓDULO#2 - TURMA 2018.2

ATIVIDADE PROPOSTA DO MÓDULO#2 - TURMA 2018.2

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