Testes de correlação

Teste de correlação entre oito variáveis:

• Conservação do solo: proporçao/porcentagem da área do município com erosão muito baixa ou baixa. (P. ex. se o município tem 80% da area com e_b_mb tem conservação do solo alta)

• Estoque de carbono: porcentagem da área do município coberta por vegetação nativa

• Balanço hídrico: média década precipitação - média década evapotranspiração (mm)

• Biodiversidade: ICMS ecológico (R$ por mun)

• Café

• Feijão

• Mandioca

• Commodities (soja, trigo e milho)

library(tidyverse)

── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
✔ dplyr     1.1.2     ✔ readr     2.1.4
✔ forcats   1.0.0     ✔ stringr   1.5.0
✔ ggplot2   3.4.2     ✔ tibble    3.2.1
✔ lubridate 1.9.2     ✔ tidyr     1.3.0
✔ purrr     1.0.1     
── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
✖ dplyr::filter() masks stats::filter()
✖ dplyr::lag()    masks stats::lag()
ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors

library(tidyr)
library(vegan)

Warning: package 'vegan' was built under R version 4.3.1

Carregando pacotes exigidos: permute

Warning: package 'permute' was built under R version 4.3.1

Carregando pacotes exigidos: lattice
This is vegan 2.6-4

library(corrplot)

Warning: package 'corrplot' was built under R version 4.3.1

corrplot 0.92 loaded

dados_new<-readr::read_delim("todos8_cs_ec_bh_bio_gr_caf_mand.csv")

Rows: 399 Columns: 9
── Column specification ────────────────────────────────────────────────────────
Delimiter: ";"
chr (1): municipio
dbl (8): conservacao_do_solo, estoque_de_carbono, balanço_hidrico, café, fei...

ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.

dados1<-dados_new %>% select(-1)

new.std<-decostand(dados1, method = "stand")

cor_new<-cor(new.std)
cor_new

                       conservacao_do_solo estoque_de_carbono balanço_hidrico
conservacao_do_solo             1.00000000         0.45398733      0.18889583
estoque_de_carbono              0.45398733         1.00000000      0.50701346
balanço_hidrico                 0.18889583         0.50701346      1.00000000
café                           -0.04910527        -0.07496737     -0.03649571
feijao                         -0.02888701         0.21345968      0.17501100
mandioca                        0.09480600        -0.36437797     -0.19131038
biodiversidade                  0.24767580         0.37265813      0.18549585
commodities(mi_tri_so)         -0.25367733        -0.01706918      0.08387193
                              café      feijao    mandioca biodiversidade
conservacao_do_solo    -0.04910527 -0.02888701  0.09480600     0.24767580
estoque_de_carbono     -0.07496737  0.21345968 -0.36437797     0.37265813
balanço_hidrico        -0.03649571  0.17501100 -0.19131038     0.18549585
café                    1.00000000 -0.05657509 -0.02291156    -0.02662221
feijao                 -0.05657509  1.00000000 -0.19203725     0.08953718
mandioca               -0.02291156 -0.19203725  1.00000000    -0.15076292
biodiversidade         -0.02662221  0.08953718 -0.15076292     1.00000000
commodities(mi_tri_so) -0.04242487  0.37701523 -0.27788717     0.04090380
                       commodities(mi_tri_so)
conservacao_do_solo               -0.25367733
estoque_de_carbono                -0.01706918
balanço_hidrico                    0.08387193
café                              -0.04242487
feijao                             0.37701523
mandioca                          -0.27788717
biodiversidade                     0.04090380
commodities(mi_tri_so)             1.00000000

corrplot(cor_new, method ='color')

* p<0.05