Atividade Módulo 3

Eiei, bom dia.
Na Atividade proposta para o módulo 3 há o exercício de construir um gráfico de dispersão onde há diferenciação entre o Cultivo Agroflorestal e Comum. É possível fazer essa diferenciação assim:

# Como o aduno afeta a produtividade?
## Adubo explicativa e produtividade resposta

levels(dados2$Cultivo)
cores2 <- c("red3", "blue3")

plot( dados2 $ Produtividade ~ dados2 $ Adubo,
las = 1,
pch = 16,
col = cores2)

Fala um detalhe bem importante aí: associar as cores aos níveis da variável. Da forma como está construído, o R vai entender que deve usar as duas cores de forma alternada - aí o gráfico não vai fazer sentido...

No seu caso, ficaria assim:

plot( dados2 $ Produtividade ~ dados2 $ Adubo,
las = 1,
pch = 16,
col = cores2[dados2 $ Cultivo])

Aqui está o script da minha tarefa 1 do módulo 3. Não sei se está correta, mas fiz dessa forma.

# Atividade módulo 3 - Tarefa 1

# carregando os dados:

setwd("//Users/gilmarajunqueira/Documents/Programa R/Programa R para Biologia da Conservação/Modulo_3/Scripts e dados M3/atividade")
dir()
dados<-read.table("sanhacos.txt", h=T)
summary(dados)

#Criando uma tabela de contingência (linha, colunas)
tabela1<-table(dados$Melastomatacea, dados$Sanhaco)
tabela2<-table(dados$Sanhaco, dados$Melastomatacea)

tabela1
tabela2

#criando um gráfico de barras, com cores e nome dos eixos:
cores<-c("green4", "red4")
barplot(tabela1, beside=T, ylim=c(0, 60), las=1, col= cores, xlab="Presente em uma Melatomatacea", ylab="Número de sanhaços")

#adicionando a legenda:
legend("topleft", legend=levels(dados$Sanhaco), col=cores, pch=15)
box()

Aqui está o script da minha tarefa 2 do módulo 3. Essa achei mais fácil de fazer do que o gráfico de barras.

# Atividade módulo 3 - Tarefa 2

# carregando os dados:

setwd("//Users/gilmarajunqueira/Documents/Programa R/Programa R para Biologia da Conservação/Modulo_3/Scripts e dados M3/atividade")
dir()
dados<-read.table("dados_atividade3.txt", h=T)
summary(dados)

#criando o gráfico:
plot(dados$Produtividade ~ dados$Adubo, xlab="Adubo", ylab="Produtividade", las=1, pch=16)

#diferenciando os pontos de acordo com o tipo de cultivo:
levels(dados$Cultivo)
cores<-c("green4", "red4")

plot(dados$Produtividade ~ dados$Adubo, xlab="Adubo", ylab="Produtividade", las=1, pch=16, col=cores[dados$Cultivo])

# adicionando a legenda:
legend("topleft", pch=16, col=cores, legend=c("Agroflorestal" ,"Comum" ))