ATIVIDADE PROPOSTA DO MÓDULO#3 - TURMA 2018.2

Ir em baixo

ATIVIDADE PROPOSTA DO MÓDULO#3 - TURMA 2018.2

Mensagem por Giovana em Ter Out 23, 2018 7:04 pm

Question pale Question
Código:
####################################
#Atividade proposta para o módulo 3#
#Giovana                           #
#R para Bio da Conservação         #
####################################


#Estabelecendo a pasta de trabalho:
setwd("C:/Users/Giovana/Documents/Mestrado/Marcos Vital - R para Biologia da conservação/Módulo 3/Scripts e dados M3/atividade")
getwd()
dir()

# ATIVIDADE 1:
#Criar o objeto sanhaco:
sanhaco<-read.table( "sanhacos.txt" , header = T)
ls()  
summary(sanhaco)
View(sanhaco)


#Tabela de contingência com a função table():
Tabela.sanhaco<-table(sanhaco$Sanhaco, sanhaco$Melastomatacea)  
Tabela.sanhaco2<-table(sanhaco$Melastomatacea, sanhaco$Sanhaco)


#Olhando a tabela:
Tabela.sanhaco  
Tabela.sanhaco2

#Obervação em campo: Árvores com pássaros.
#H1: Melastomatacea são mais frequentadas por sanhaços.
#Gráfico de barras mostrando a relação entre as varáveis (forma mais intuitiva de interpretar):
barplot(Tabela.sanhaco, beside = T)
barplot(Tabela.sanhaco2, beside = T)


#O gráfico escolhido em uma versão bacana:
cores<-c("light blue", "light gray")
barplot(Tabela.sanhaco, beside=T, ylim=c(0, 50), las=1, xlab="Espécie de Melastomatacea", ylab="Número de árvores", col=cores)
legend("topleft", legend=levels(sanhaco$Sanhaco), col=cores, pch=15, title="É um Sanhaço?")
box()


# ATIVIDADE 2:
#Criar o objeto agricultura:
agricultura<-read.table( "dados_atividade3.txt", header = T)
ls()  
summary(agricultura)
View(agricultura)


#Gráfico de dispersão entre as variáveis quantitativas:
#Como o adubo (X) afeta a produtividade (Y):
plot(agricultura$Produtividade ~ agricultura$Adubo)


#Gráfico diferenciando os pontos de acordo com a variável "Cultivo":
levels(agricultura$Cultivo)
cores<-c("light green", "brown")
plot(agricultura$Produtividade ~ agricultura$Adubo, las=1, xlab="Quantidade de adubo", ylab="Produtividade", col=cores[agricultura$Cultivo], pch=18)
legend("topleft", legend=levels(agricultura$Cultivo), col=cores, pch=18, title="Tipo de cultivo")


Giovana

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

Ver perfil do usuário

Voltar ao Topo Ir em baixo

Dúvida função Plot

Mensagem por Ricardo Matheus Pires em Sab Out 27, 2018 11:50 am

Bom dia, pessoal! Aproveitando essa postagem, estou um pouco atrasadinho e terminando o Mód. 3 ainda. Me surgiu uma dúvida besta:

Qual a diferença entre uma sintaxe e outra?
plot(dados$Produtividade ~ dados$Adubo)
plot(dados$Adubo, dados$Produtividade)

Ambos me dão o mesmo gráfico, vi nos Scripts do professor a forma 1, mas quando fiquei com dúvida, consultei o help da função e a sintaxe era a forma 2, não sei se perdi alguma explicação durante as aulas...

?help #hahhahaha

Obrigado

Ricardo Matheus Pires

Mensagens : 2
Data de inscrição : 17/09/2018

Ver perfil do usuário

Voltar ao Topo Ir em baixo

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

Mensagem por Prof. Marcos em Sab Nov 10, 2018 9:31 am

Olá, Ricardo, tudo bem?

Ricardo Matheus Pires escreveu:

Qual a diferença entre uma sintaxe e outra?
plot(dados$Produtividade ~ dados$Adubo)
plot(dados$Adubo, dados$Produtividade)


Então, em termos de funcionamento, não há diferença nenhuma, os gráficos serão iguais. Wink

A diferença está mais ligada à forma de se "anunciar" a relação esperada entre as varáveis, e a forma de construirmos modelos estatísticos no R - especialmente os modelos lineares, como ANOVAs e regressões de vários tipos.

Ao escrever "dados$Produtividade ~ dados$Adubo", estamos falando algo como "achamos que a variável produtividade responde à variável adubo". E ao dizer "dados$Adubo, dados$Produtividade", estamos falando algo como "queremos a variável adubo no eixo horizontal e a produtividade no vertical".

Como eu disse antes, o efeito no gráfico é idêntico. Eu gosto mais de usar o "~" sempre que posso (existem uma ou outra função, como é o caso no sciplot que não permitem isso, mas são exceções), pois assim nos habituamos ao uso nos modelos estatísticos.

Enfim, como muitas coisas no R, no fundo é apenas uma questão de preferência, pois são dois caminhos que geram exatamente o mesmo resultado, ok? Então siga com o que você preferir e seja feliz! Smile

Abraços
avatar
Prof. Marcos

Mensagens : 377
Data de inscrição : 23/12/2015
Idade : 37
Localização : Maceió, AL

Ver perfil do usuário https://cantinhodor.wordpress.com/

Voltar ao Topo Ir em baixo

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

Mensagem por Cleyton Tenório em Seg Nov 12, 2018 10:30 pm

Código:
##############################
#Atividade módulo 3          #
#Cleyton Tenório Barros      #
#R para Bio da conservação   #
##############################

#Atividade 1:

#Estabelendo a pasta de trabalho:

setwd("D:/R/Rbio")
getwd() # Verificando se estou na pasta certa.
dir() #Conferindos os arquivos da pasta.

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

#Importando dados:
dados<-read.table("sanhacos.txt", h=T)
summary(dados) # Verificando os dados.

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

#Criando tabelas de contigência com a função table():

tabelinha<-table(dados$Sanhaco, dados$Melastomatacea)

#Vendo a tabela:
tabelinha

#Criando a tabela em outra ordem para comparar:

tabelinha.2<-table(dados$Melastomatacea, dados$Sanhaco)

#Vendo a tabela:
tabelinha.2

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

#H1: É que os sanhaços tem um predileção por Melastomatacea.

#Criando gráfico de barras (Obeservando o modo mais, intuitivo de interpretação:
barplot(tabelinha, beside = T)
barplot(tabelinha.2, beside = T)

#Gráfico de barras organizado e com legenda:

cores<-c("green1", "gray")
barplot(tabelinha, beside=T, ylim=c(0, 70), las=1, ylab="Nº de Árvores", xlab="Melastomatacea X Outras", col=cores)

legend("topleft", legend = levels(dados$Sanhaco), col=cores, pch = 15, title="Predileção do Sanhaço")
box()

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

#Atividade 2:

setwd("D:/R/Rbio")
getwd() # Verificando se estou na pasta certa.
dir() #Conferindos os arquivos da pasta.

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

#Importando dados:
producao<-read.table("dados_atividade3.txt", h=T)
summary(producao) # Verificando os dados.

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

#H1: Com uma quatidade adequada de adubo eu tenho uma melhor produtividade.

#Criando o gráfico de dispersão entre as variáveis quantitativas:

plot(producao$Produtividade~producao$Adubo)

#Criando gráfico, diferenciando os pontos de acordo com a variável "Cultivo":

levels(producao$Cultivo)

cores<-c("pink", "purple")

plot(producao$Produtividade~producao$Adubo, las=1, xlab="Quantidade de Adubo", ylab="Produtividade", col=cores[producao$Cultivo], pch=20, ylim=c(12, 26))

legend("topleft", legend=levels(producao$Cultivo), pch=20, col=cores, title="Tipos de Cultivo")
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#3 - TURMA 2018.2

Mensagem por Conteúdo patrocinado


Conteúdo patrocinado


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