library(devtools)
## Carregando pacotes exigidos: usethis
library(ecodados)
library(vegan)
## Carregando pacotes exigidos: permute
##
## Attaching package: 'permute'
## The following object is masked from 'package:devtools':
##
## check
## Carregando pacotes exigidos: lattice
## This is vegan 2.6-4
library(ggplot2)
library(BiodiversityR)
## Carregando pacotes exigidos: tcltk
## BiodiversityR 2.15-1: Use command BiodiversityRGUI() to launch the Graphical User Interface;
## to see changes use BiodiversityRGUI(changeLog=TRUE, backward.compatibility.messages=TRUE)
composicao_especies <- ecodados::composicao_anuros_div_taxonomica
rank_com6 <- rankabundance(composicao_especies[6, composicao_especies[6,] > 0])
## Warning in qt(0.975, df = n - 1): NaNs produzidos
rankabunplot(rank_com6, scale = "logabun", specnames = c(1),
pch = 19, col = "violet")
rank_com7 <- rankabundance(composicao_especies[7, composicao_especies[7,] > 0])
## Warning in qt(0.975, df = n - 1): NaNs produzidos
rankabunplot(rank_com7, scale = "logabun", specnames = c(1),
pch = 19, col = "Green")
rank_com8 <- rankabundance(composicao_especies[8, composicao_especies[8,] > 0])
## Warning in qt(0.975, df = n - 1): NaNs produzidos
rankabunplot(rank_com8, scale = "logabun", specnames = c(1),
pch = 19, col = "Blue")
rank_com9 <- rankabundance(composicao_especies[9, composicao_especies[9,] > 0])
## Warning in qt(0.975, df = n - 1): NaNs produzidos
rankabunplot(rank_com9, scale = "logabun", specnames = c(1),
pch = 19, col = "Red")
rank_com10 <- rankabundance(composicao_especies[10, composicao_especies[10,] > 0])
## Warning in qt(0.975, df = n - 1): NaNs produzidos
rankabunplot(rank_com10, scale = "logabun", specnames = c(1),
pch = 19, col = "Yellow")
riqueza_res <- specnumber(composicao_especies)
riqueza_res
## Com_1 Com_2 Com_3 Com_4 Com_5 Com_6 Com_7 Com_8 Com_9 Com_10
## 10 10 5 5 5 6 2 4 6 4
shannon_res <- diversity(composicao_especies, index = "shannon", MARGIN = 1)
shannon_res
## Com_1 Com_2 Com_3 Com_4 Com_5 Com_6 Com_7 Com_8
## 2.3025851 0.5002880 0.9580109 1.6068659 1.4861894 1.5607038 0.6931472 1.1058899
## Com_9 Com_10
## 1.7140875 1.2636544
simpson_res <- diversity(composicao_especies, index = "simpson", MARGIN = 1)
simpson_res
## Com_1 Com_2 Com_3 Com_4 Com_5 Com_6 Com_7 Com_8
## 0.9000000 0.1710000 0.4814815 0.7989636 0.7587500 0.7674858 0.5000000 0.5850000
## Com_9 Com_10
## 0.8088889 0.6942149
precipitacao <- ecodados::precipitacao_div_taxonomica
precipitacao
## prec
## Com_1 3200
## Com_2 3112
## Com_3 2800
## Com_4 1800
## Com_5 2906
## Com_6 3005
## Com_7 930
## Com_8 1000
## Com_9 1300
## Com_10 987
dados_div <- data.frame(precipitacao$prec, riqueza_res,shannon_res,
simpson_res)
colnames(dados_div) <- c("Precipitacao", "Riqueza", "Shannon", "Simpson")
anova_shan <- lm(Shannon ~ Precipitacao, data = dados_div)
anova(anova_shan)
## Analysis of Variance Table
##
## Response: Shannon
## Df Sum Sq Mean Sq F value Pr(>F)
## Precipitacao 1 0.10989 0.10989 0.3627 0.5637
## Residuals 8 2.42366 0.30296
plot(anova_shan)
anova_simpson <- lm(Simpson ~ Precipitacao, data = dados_div)
anova(anova_simpson)
## Analysis of Variance Table
##
## Response: Simpson
## Df Sum Sq Mean Sq F value Pr(>F)
## Precipitacao 1 0.00132 0.001325 0.0252 0.8778
## Residuals 8 0.42064 0.052580
plot(anova_simpson)
anova_riqueza <- lm(Riqueza ~ Precipitacao, data = dados_div)
anova(anova_riqueza)
## Analysis of Variance Table
##
## Response: Riqueza
## Df Sum Sq Mean Sq F value Pr(>F)
## Precipitacao 1 30.622 30.6224 8.9156 0.01744 *
## Residuals 8 27.478 3.4347
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
plot(anova_riqueza)
View(composicao_especies)