cor_gcms: Compute and Plot Correlation Matrix for a Set of General Circulation Models

Description

This function computes and visualizes the correlation matrix for a set of General Circulation Models (GCMs) based on their variables.

Usage

cor_gcms(
  s,
  var_names = c("bio_1", "bio_12"),
  study_area = NULL,
  scale = TRUE,
  method = "pearson"
)

Value

A list containing two items: cor_matrix (the calculated correlations between GCMs) and cor_plot (a plot visualizing the correlation matrix).

Arguments

s: A list of stacks of General Circulation Models (GCMs).
var_names: Character. A vector with names of the variables to compare, or 'all' to include all variables.
study_area: An Extent object, or any object from which an Extent object can be extracted. Defines the study area for cropping and masking the rasters.
scale: Logical. Whether to apply centering and scaling to the data. Default is TRUE.
method: Character. The correlation method to use. Default is "pearson". Possible values are: "pearson", "kendall", or "spearman".

Author

Luíz Fernando Esser (luizesser@gmail.com) https://luizfesser.wordpress.com

Examples

Run this code

var_names <- c("bio_1", "bio_12")
s <- import_gcms(system.file("extdata", package = "chooseGCM"), var_names = var_names)
study_area <- terra::ext(c(-80, -30, -50, 10)) |> terra::vect(crs="epsg:4326")
cor_gcms(s, var_names, study_area, method = "pearson")