Get started

Design

The core of shar are functions to to simulate null model data by randomizing either the environmental data (i.e. raster data) or the locations of species (i.e. point pattern data). The null data is then used to analyse if significant species-habitat associations are present. Additionally, functions to visualize and analyse the results are available as well as some utility functions. The methods are mainly described in Harms et al. (2001), Plotkin et al. (2000) and Wiegand & Moloney (2014). The methods are not neessary complementary, but are rather different approaches for the same result.

Preprocessing of input data

All functions are designed for discrete habitat classes. Following, in case only continuous data is available, this has to be classified to discrete classes. classify_habitats provides several ways to classify the data such as the Fisher-Jenks algorithm (Fisher 1958, Jenks & Caspall 1971)

landscape_discrete <- classify_habitats(raster = landscape, classes = 5)

Randomize environmental data

There are two functions to randomize the environmental data: translate_raster() and randomize_raster(). The first function is a torus translation of the raster, shifting the habitat map in all four cardinal directions. This is only possible for rectangular observation areas and results in n_random <- (raster::nrow(landscape) + 1) * (raster::ncol(landscape) + 1) - 4 randomized rasters. The other function randomizes the environmental data using a random-walk algorithm. Here, the number of randomized rasters can be specified using the n_random argument.

landscape_torus <- translate_raster(raster = landscape_discrete)

landscape_random <- randomize_raster(raster = landscape_discrete, n_random = 19)

Randomize location data

To randomize the location data (i.e. the point pattern) either fit_point_process() or reconstruct_pattern() are available. The first fits either a Poisson process or a cluster process to the data. The difference to solutions from the spatstat package is that the number of points is always identical. The second functions reconstructs the spatial characteristics of the data using pattern reconstruction (Tscheschel & Stoyan 2006). This is advantageous for point patterns not describable by simple point process models. For both function, the number of patterns can be specified by the n_random argument.

fitted <- fit_point_process(pattern = species_a, n_random = 19, process = "cluster")

reconstructed <- reconstruct_pattern(pattern = species_a, max_runs = 500, n_random = 19) # takes some time

Analyse results

The most important function to analyse results is results_habitat_association(). This function compares the observed data to the null model data and by that is able to show significant species-habitat associations. The functions work for both, randomized environmental data or randomized location data.

results_habitat_association(pattern = species_a, raster = landscape_random)
#> > Input: randomized raster | Quantile thresholds: negative < 0.025 - positive > 0.975
#>   habitat count    lo    hi significance
#> 1       1    10  0.45  7.65     positive
#> 2       2    14 11.35 21.30         n.s.
#> 3       3    30 15.45 28.65     positive
#> 4       4     0 10.25 24.65     negative
#> 5       5    14  6.45 13.55     positive

results_habitat_association(pattern = reconstructed, raster = landscape_discrete)
#> > Input: randomized point pattern | Quantile thresholds: negative < 0.025 - positive > 0.975
#>   habitat count    lo    hi significance
#> 1       1    10  0.00  3.00     positive
#> 2       2    14  9.00 20.00         n.s.
#> 3       3    30 17.35 32.65         n.s.
#> 4       4     0  9.25 25.00     negative
#> 5       5    14  4.45 11.00     positive

There is also the possibility to visualize the randomized data using plot_randomized_raster() or plot_randomized_pattern().

plot_randomized_raster(landscape_random)

plot_randomized_pattern(reconstructed)

For the randomized point pattern data, it is also possible to show the “difference” in terms of the energy (Tscheschel & Stoyan 2006) between the patterns.

calculate_energy(pattern = fitted, return_mean = TRUE)
#> [1] 0.357844

calculate_energy(pattern = reconstructed, return_mean = TRUE)
#> [1] 0.07572195