Avaliando os efeitos de mudanças climáticas no R

Modelagem no wallace

No seguinte vídeo vocês vão aprender a fazer a modelagem de distribuição de espécies no pacote wallace do R: https://youtu.be/EGjC0KEBzAY

Note que o vídeo usa uma versão anterior. Busque fazer o donlwoad da versão mais nova do wallace. Como no código de R abaixo:

install.packages("devtools")
devtools::install_github("wallaceEcoMod/wallace")

Uma vez que você fez todos os passos, podemos agora fazer o pós processamento dos dados e analisar numericamente o impacto das mudanças climáticas.

Dados espaciais

Baixe os dados em formato “raster” das previsões de distribuição binária do wallace para o presente e para o futuro. Escolha o modelo com o maior AUC ou menor AIC. Pode usar o botão download como na figura abaixo:

Presente:

Futuro:

Carregar os dados no R

Primeiro instale o pacote necessário no R conforme o código abaixo:

install.packages("raster")

Carregue o pacote no R:

library(raster)

Loading required package: sp

Indique o caminho completo para o seu arquivo, como no código abaixo:

presente <- raster("Panthera_onca_p10.asc")
futuro <- raster("Panthera_onca_xfer_p10.asc")

Olhe o número de células do presente e do futuro:

pres_cells <- sum(values(presente) == 1, na.rm = TRUE)
fut_cells <- sum(values(futuro) == 1, na.rm = TRUE)

Agora veja a proporção de perda ou ganho:

1 - (fut_cells / pres_cells)

[1] 0.1678436

Nesse caso, podemos dizer que, segundo a previsão, haverá uma perda de 16,78% da área de distribuição geográfica nesse cenário futuro.