Calculate correlation
Arguments
- metrics
Tibble with results of as returned by the landscapemetrics package.
- method
Type of correlation. See
link{cor}
for details.- diag
If FALSE, values on the diagonal will be NA.
- simplify
If TRUE and only one level is present, only a tibble is returned.
Details
The functions calculates the correlation between all metrics. In order to calculate correlations, for the landscape level more than one landscape needs to be present. All input must be structured as returned by the landscapemetrics package.
Examples
landscape <- terra::rast(landscapemetrics::landscape)
metrics <- calculate_lsm(landscape, what = c("patch", "class"))
#> Warning: Please use 'check_landscape()' to ensure the input data is valid.
#> Warning: Class 1: PAFRAC = NA for class with < 10 patches
#> Warning: Class 3: PAFRAC = NA for class with < 10 patches
calculate_correlation(metrics, method = "pearson")
#> $patch
#> # A tibble: 78 × 3
#> metric_1 metric_2 value
#> <chr> <chr> <dbl>
#> 1 area area 1
#> 2 cai area 0.886
#> 3 circle area 0.511
#> 4 contig area 0.742
#> 5 core area 0.987
#> 6 enn area -0.443
#> 7 frac area 0.318
#> 8 gyrate area 0.952
#> 9 ncore area 0.702
#> 10 para area -0.709
#> # ℹ 68 more rows
#>
#> $class
#> # A tibble: 1,540 × 3
#> metric_1 metric_2 value
#> <chr> <chr> <dbl>
#> 1 ai ai 1
#> 2 area_cv ai -0.926
#> 3 area_mn ai 1.00
#> 4 area_sd ai 0.986
#> 5 ca ai 0.986
#> 6 cai_cv ai -0.973
#> 7 cai_mn ai 0.997
#> 8 cai_sd ai 0.984
#> 9 circle_cv ai -0.983
#> 10 circle_mn ai 0.990
#> # ℹ 1,530 more rows
#>