Next Article in Journal
Resistance of Lichens and Mosses of Regenerated Alpine Communities to Repeated Experimental Trampling in the Belianske Tatras, Northern Slovakia
Previous Article in Journal
Impact of Nesting Mode, Diet, and Taxonomy in Structuring the Associated Microbial Communities of Amazonian Ants
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Importance of Including Spatial Autocorrelation When Modelling Species Richness in Archipelagos: A Bayesian Approach

by
Diogo Duarte Barros
1,*,
Maria da Luz Mathias
1,
Paulo A. V. Borges
2 and
Luís Borda-de-Água
3,4,5,*
1
CESAM—Centro de Estudos do Ambiente e do Mar, Departamento de Biologia Animal, Faculdade de Ciências, Universidade de Lisboa, Edifício C2, 3º piso, Campo Grande, 1749-016 Lisboa, Portugal
2
cE3c—Centre for Ecology, Evolution and Environmental Changes, Azorean Biodiversity Group, CHANGE—Global Change and Sustainability Institute, Faculty of Agricultural Sciences and Environment, University of the Azores, 9700-042 Angra do Heroísmo, Portugal
3
CIBIO/InBio, Centro de Investigação em Biodiversidade e Recursos Genéticos, Laboratório Associado, Universidade do Porto, Campus Agrário de Vairão, 4485-661 Vairão, Portugal
4
CIBIO/InBio, Centro de Investigação em Biodiversidade e Recursos Genéticos, Laboratório Associado, Instituto Superior de Agronomia, Universidade de Lisboa, Tapada da Ajuda, 1349-017 Lisboa, Portugal
5
BIOPOLIS Program in Genomics, Biodiversity and Land Planning, CIBIO, Campus de Vairão, 4485-661 Vairão, Portugal
*
Authors to whom correspondence should be addressed.
Diversity 2023, 15(2), 127; https://doi.org/10.3390/d15020127
Submission received: 1 November 2022 / Revised: 6 January 2023 / Accepted: 13 January 2023 / Published: 17 January 2023
(This article belongs to the Topic Biodiversity in the Azores: A Whole Biota Assessment)

Abstract

:
One of the aims of island biogeography theory is to explain the number of species in an archipelago. Traditionally, the variables used to explain the species richness on an island are its area and distance to the mainland. However, increasing evidence suggests that accounting for other variables is essential for better estimates. In particular, the distance between islands should play a role in determining species richness. This work uses a Bayesian framework using Gaussian processes to assess whether distance to neighbouring islands (spatial autocorrelation) can better explain arthropod species richness patterns in the Azores Archipelago and in the Canary Islands. This method is flexible and allows the inclusion of other variables, such as maximum altitude above sea level (elevation). The results show that accounting for spatial autocorrelation provides the best results for both archipelagos, but overall, spatial autocorrelation seems to be more important in the Canary archipelago. Similarly, elevation plays a more important role in determining species richness in the Canary Islands. We recommend that spatial autocorrelation should always be considered when modelling an archipelago’s species richness.

1. Introduction

Islands and archipelagos have provided fertile ground to develop theories in ecology, evolution, and biogeography (e.g., [1,2,3,4,5]). Assessing species richness on islands has been one of the major topics addressed by island biogeographers and one that has had a major impact on the development of theories in ecology overall. One of the best-known theories to explain the patterns and processes determining species richness in islands is the Equilibrium Theory of Island Biogeography by MacArthur and Wilson [6,7], hereafter ETIB. ETIB explains the species richness in an island based on immigration (dispersal) and extinction. In turn, these two processes are a function of the area of the island and its distance to the mainland, the latter acting as a source pool. According to this theory, assuming all other factors (such as climate or topographic complexity) are equal, for two islands at the same distance from the mainland, the largest island has more species, because more migrants are able to reach it (target effect), and there are fewer extinctions (resource availability effect). Similarly, for islands of the same size, the one closer to the mainland has more species because the number of individuals, and thus of species, reaching it from the mainland is larger and there is also a higher probability of rescue effects.
The pioneering work of MacArthur and Wilson inspired multiple lines of research on biogeography and biodiversity, gaining both support and criticism. For instance, the original ETIB model was criticized for being incomplete because it did not consider the evolutionary time scale on oceanic islands. Such criticisms eventually led to the development of the General Dynamic Model [8], hereafter GDM. In the GDM, the biogeographical patterns and processes of oceanic islands are fundamentally related to the age of the island and its geological evolution. GDM helped to unravel the prominent role of single island endemics, hereafter SIE, as potential indicators of evolutionary dynamics in ocean island archipelagos [2,9,10], by suggesting that the number of SIE on ocean islands is strongly determined by the island’s geological age, as reviewed by Borregaard et al. [5]. Moreover, islands with volcanic origin have complex topographies, and maximum altitude above sea level (elevation) can be considered a surrogate of an island’s habitat diversity [11] and island speciation rates [12]. The geological and biological history of each island’s mountains influences its species richness, level of endemism, and biological dynamics [8]. Therefore, we should expect a high number of SIE in mountainous islands with complex topographies.
The island species–area relationship, hereafter ISAR, relates the number of species, S, on an island of area, A, and it is often expressed as a power law,
S = cAz,
where c and z are constants. This formula has been a favourite among theoreticians and has also received empirical support (e.g., [13,14]). Diamond [15] added an exponential term to describe the effect of the distance of the island to the mainland, D,
S = c exp(−dD)Az,
where d is a positive constant (for more recent developments see [1]). However, often islands do not occur in isolation but in archipelagos. It is, then, to be expected that the species richness on an island is not only influenced by its area and distance to the mainland but also by the distance to the other islands in the archipelago [16]. For instance, the distance among islands is likely to affect the ability of species to disperse within the archipelago, which may affect not only species richness but also species composition (and beta-diversity) [16,17,18]. The latter is not the object of this work; however, closeness may affect species richness because it is likely to be a proxy for similar environmental conditions and for the islands’ age [16]. Moreover, some islands may have been connected during glacial periods, which is more likely to have occurred for nearby islands, and their present communities may still reflect a previous species’ richer community [2]). Therefore, when trying to fit an ISAR, the distance among islands should be incorporated into the model.
Spatial autocorrelation is often disregarded when fitting the ISAR (e.g., [15]). However, failing to address the effects of spatial autocorrelation may lead, for example, to bias in the parameter estimates. There are several approaches to deal with spatial autocorrelation, such as the use of spatial linear models of the Conditional Autoregressive Model [19], frequently used to address first-order dependencies, or the Simultaneous Autoregressive Model [17], used when there are multiple order dependencies. Dormann et al. [20] extensively reviewed and proposed multiple approaches to deal with spatial autocorrelation, including spatial generalized linear mixed models, where the spatial autocorrelation is modelled by specifying the correlation structure for the residuals. For instance, Selmi and Boulinier [18] suggested the use of spherical, Gaussian, and exponential covariance models to account for the role of spatial covariance in regression models.
Here, we use a model developed by McElreath [21]. It uses a Gaussian process regression (or simply a Gaussian process) to explicitly incorporate spatial autocorrelation. A Gaussian process is a Bayesian method that adds a varying intercept that accounts for the non-independence of the data, and a multivariate prior for these intercepts [21,22]. This approach has seen multiple practical applications in ecology, namely, within species distribution models [23], inference of function-valued traits [24], and beta diversity predictions as functions of environmental distance [25].
In this work, we illustrate how to model an ISAR with a Gaussian process, thus considering spatial autocorrelation, using data on arthropod species richness from the Azores [26], an archipelago consisting of nine islands, and from the Canary Islands [27], an archipelago consisting of seven islands. In addition to the area of the islands, their distances to the mainland, and the distances between them, we will also consider more complicated models that include the elevation of the islands. Finally, because ISARs only deal with the number of species richness, we will only assess spatial autocorrelation in the context of species richness; for works on community (dis)similarity as a function of distance, see, e.g., [28,29].

2. Materials and Methods

2.1. Statistical Model

Bayesian methods require that we assign a likelihood for the data and choose priors for the parameters. To choose the likelihood, first notice that the number of species on island i, Si, is a discrete number. Two distributions are obvious candidates: the Poisson distribution, which has only one parameter, μ, the mean (which is equal to the variance),
S i Poisson ( μ i )
and the negative binomial, which has two parameters, μ, the mean, and ϕ that controls for over-dispersion,
S i Negative _ binomial ( μ i , ϕ i ) .          
Typically, the number of species is over-dispersed, that is, the variance is much larger than the mean. This seems to rule out the Poisson distribution; however, this concern is only valid when we do not include a Gaussian process. When we model the parameter μ of the Poisson distribution with a Gaussian process, the result is a mixture distribution, very much like the negative binomial; recall that a negative binomial is a mixture distribution where the mean of a Poisson distribution is distributed according to a gamma distribution (e.g., [22]). Therefore, for models without the Gaussian process, we use only the negative binomial distribution, but for models with a Gaussian process, we assess both the Poisson and the negative binomial distributions.
As usual, we model the mean, μi, as a function of predictor variables, xij. Here we assume power–law relationships,
μ i = β 0 · x i 1 β 1 · x i 2 β 2   x i j β j ,  
where β j are the parameters to be estimated. We chose power laws because this is the form often assumed to be the relationship between the area of an island and the number of species, Equation (1) (e.g., [30]). Nevertheless, as we will see later, we also assume an exponential decay for the distance to the mainland.
As mentioned previously, the ISAR is likely to be also determined by factors that are related to the distance between the islands, such as dispersal or similar characteristics shared due to proximity, implying that the data from different islands are not independent. Often, dependence is dealt with by considering multilevel/hierarchical models (e.g., [31,32]). In the case of an ISAR, this would consist of adding a varying intercept, γ[i], to Equation (5) (the square brackets emphasize the hierarchical nature of this parameter), and then assuming that the terms γ[i] come from a common distribution (e.g., normal), whose parameters would also need to be estimated. To preclude μi to become negative, γ[i] enters Equation (5) as a multiplicative term (e.g., [21]),
μ i = β 0 · e γ [ i ] x i 1 β 1 · x i 2 β 2   x i j β j .    
This approach allows the varying intercepts, γ[i], to “absorb” some of the variation associated with each island. However, it assumes that γ[i] are discrete categories; therefore, it ignores that spatial autocorrelation is a function of distance and, thus, should be a continuous variable. One way to handle cases where the varying intercepts, γ[i], exhibit a continuous dependence is to use a Gaussian Process (e.g., [22]). Specifically, we assume that the varying intercepts γ[i] come from a multivariate normal distribution (MVNormal)
γ ~ MVNormal ( 0 , K ) ,
where 0 stands for the vector of the means, all equal to zero, and K is a covariance matrix. The dependence of γ[i] on distance is established through the functional form of the covariance matrix K. Often, γ[i] is assumed to decay exponentially with the square of the distance (e.g., [21]), with the elements of K, kij, given by
k i j = η 2 exp ( ρ 2 D i j 2 ) + δ i j σ 2 ,  
where Dij stands for the distance between islands i and j, and η, ρ, and σ are parameters to be estimated. δij is the Kronecker delta, equal to 1 when i = j and zero otherwise, meaning that the term δ i j σ 2 corresponds to the variance within an island. However, because there is only one observation per island (its number of species), the term δ i j σ 2 is irrelevant (there is no variance among the values for an island); hence, following [21], we set this term equal to a constant (0.01). The above choice of the covariance matrix, Equation (7), is by no means the only possible one; see, for instance, [33]. We tested alternative formulations, but the results were similar.
Notice that the covariance matrix, K, plays a major role in the formulation of the model. It is through this matrix that the distances between the islands, Dij, are explicitly included in the model and its elements reflect the importance (or not) of the spatial autocorrelation among the islands. Although this is the matrix whose parameters are being estimated, Equation (7), it is easier to interpret, instead, the correlation matrix (e.g., see [21]). Given two islands, I and J, if ki,j is the corresponding element in the covariance and ki,i and kj,j are their variances, then the correlation is calculated using k i , j / k i , i k j , j ; naturally, the diagonal elements of the correlation matrix are equal to 1.
In summary, we assume two types of likelihood, Equations (3) and (4), with the means modelled by Equations (5) or (6) (the former only for the negative binomial). This implies that there is one parameter β0, j parameters βj (to be identified), the parameter ϕ if we used a negative binomial distribution, and ρ2 and η when we include the Gaussian process. To illustrate the above model, assume a Poisson distribution with the area of the island, Ai, as the only explanatory variable but including a Gaussian process. From Equations (3) and (6), we obtain
S i ~ Poisson ( μ i ) ,
with
μ i = β 0 e γ [ i ] A i z ,
and recover the classical power–law formulation of the ISAR, Equation (1), with   c = β 0 e γ [ i ] , but where c now includes an autocorrelation term.
Finally, a Bayesian approach also requires that each parameter has a prior distribution. However, the choice of priors depends on the specificities of the situation being analyzed, hence, we defer the discussion on the priors until after the case studies are introduced.

2.2. Case Study—The Azores and Canary Archipelagos

2.2.1. Study Area

The Azores and the Canary Islands are two of the archipelagos that make up Macaronesia. The Azores is an isolated northern Atlantic archipelago of nine volcanic islands, located roughly between 37° and 40° N and 24°and 31° W (WGS 84 coordinate system). It is situated across the mid-Atlantic Ridge, which separates the western group (Flores and Corvo) from the central (Faial, Pico, São Jorge, Terceira, and Graciosa) and the eastern (São Miguel and Santa Maria) groups, Figure 1a. The climate is temperate oceanic, and its relative humidity often reaches 95% in the high-elevation native semi-tropical evergreen laurel forest, which contributes to small temperature fluctuations throughout the year [34]. The Canary Islands comprise seven volcanic islands located between 27° and 30° N and 13° and 19° W (WGS 84 coordinate system), distancing roughly 1300 km from the Azores. The Canary Islands are divided into two groups: the western one (Tenerife, La Palma, El Hierro, Gomera, and Gran Canaria) and the eastern one (Fuerteventura and Lanzarote), Figure 1b. The climate is subtropical and desert with a great variation in relative humidity [35].
For each island, we compiled information on the area, geological age, maximum altitude above sea level (elevation) and distance to the mainland [36,37,38,39,40,41,42]; see also Tables S1 and S2 in Supplementary Material S1). The distances between the islands were estimated using the Haversine formula and the centroid of each island; see Tables S3 and S4 in Supplementary Material S1.

2.2.2. Arthropod Data

Although new species are added and described occasionally, we believe the current knowledge of the arthropod diversity of these two archipelagos is reasonably complete. We used data from the two most recently updated checklists. For the Azores, this is an updated version of [26] based on new additions from the literature, leading to a total of 2323 species; see Table S5. For the Canary archipelago, we used [27], which contains a total of 7402 species; see Table S6. Therefore, all calculated richness values for each island are for the sums of the lowest taxonomic rank (i.e., species or subspecies). Species were classified as natives (those that were deliberately, or not, introduced by humans) archipelago endemics, single island endemics (those endemic species that occur only on one island, SIE), and multiple island endemics (endemic species that occur on two or more islands, hereafter MIE).
Since the number of endemic species may be determined by processes different from those of the exotic and the native non-endemic species, we modelled separately the total number of species, the number of endemic species, the number of SIE, and the number of MIE species.

2.3. Selection of Explanatory Variables

We started by considering the following variables as determinants of species richness: area of the island, age of the island, distance to the mainland, and maximum altitude, i.e., elevation [2,8]. However, the age of an island and the distance to the mainland are correlated; see Figures S1 and S2 in Supplementary Material (a relationship that was stronger for the Canary Islands). This happens because islands are mainly formed along the mid-ocean Atlantic ridge and, as it widens, younger islands are farther away from the continents. Since we have more confidence in the value of the distance to the mainland than in the ages of the islands, we used the former. Thus, in summary, we used “area”, “elevation”, and “distance to the mainland” as explanatory variables in Equations (5) and (6).
Both archipelagos have an east–west orientation (see Figure 1) and the distance of the islands from the mainland varies greatly within the archipelagos. In the Azores, the island closest to the mainland (Europe) is São Miguel, at a distance of 1368 km, and the farthest is Flores, at 1864 km. In the Canary Islands, the island closest to the mainland (Africa) is Fuerteventura, at a distance of 95 km, and the farthest is La Palma, at 414 km. Some authors (e.g., [15]) have assumed that Si exhibits an exponential relationship with the distance to the mainland, expression (2). Therefore, in addition to expressions (5) and (6), we also considered:
μ i = β 0 · e γ [ i ] x i 1 β 1 · x i 2 β 2   x i j β j e d D i ,
where Di is the distance to the mainland of the island i and d is a parameter to be estimated.

2.4. Choice of Priors

When choosing the priors, we aimed at priors that were weakly regularizing, so that they were non-informative for the range of values of interest but allowed for efficient numerical estimation of the posterior distribution [21]. Furthermore, we checked whether the priors led to values that were within our expectations for the number of species in the archipelagos. We adopted the following priors for both archipelagos:
β 0 normal ( 100 , 20 ) ,
β j exponential ( 1 ) ,
d exponential ( 1 ) ,
ρ 2 exponential ( 0.5 ) ,
and
η exponential ( 1 ) .

2.5. Software Packages

All the analyses were performed in R version 4.2.0 [43] with the packages Rethinking [21] and RStan [44]. We checked the convergence of the models obtained by confirming that the R ^ (“rhat”, see Supplementary Information S2) was not greater than 1. R ^ measures variance within all MCMC chain samples, and if all chains are at equilibrium, they should have the same variance, thus, R ^ = 1 [44]. In addition, we also show the effective number of samples, “n_eff”, that is, the estimated number of independent samples [44].

3. Results

3.1. The Traditional ISAR with and without the GAUSSIAN Process

Because the traditional power–law form of the ISAR, S = cAz, has played an important role in ecological theory (e.g., [6]), our analyses start with this version with and without the Gaussian process. The main purpose of this analysis is to illustrate the advantage of including a Gaussian process. For ease of comparison, we use in both cases the negative binomial likelihood. Table 1 and Table 2 show the values of the parameters and Figure 2a,b show the fitting.
Interestingly, the means and standard deviations of the c parameters with and without the Gaussian Process are similar for both archipelagos, but the parameters z and ϕ have large coefficients of variation in the model with the Gaussian process. This reveals that when spatial autocorrelation is not included, we have undue certainty about some parameters.
We used the “Widely Applicable Information Criterion” (WAIC) to compare the fitting of the two models, and adopted the terminology suggested by [45] to compare the performance of the models. Despite the larger number of parameters, WAIC shows that the model with the Gaussian process performed best in both archipelagos; see Table 3. However, the improvement is greater for the Canary Islands than for the Azores. Indeed, as Figure 2a shows, the uncertainty of the fitting provided by the two models is similar for the Azores but very different for the Canary Islands.
Spatial autocorrelation is modelled through the covariance matrix, K; however, as we previously mentioned, the correlation matrix is easier to interpret [21]. To obtain the correlation matrix, because both parameters have skewed distributions, we use the median of the η2 and ρ2 posteriors. The correlation matrices are shown in Table 4 and Table 5 for the Azores and the Canary Islands, respectively. Figure 2c,d shows the archipelagos and the intensity of the color of the line between the islands reflects the value of the correlation; for some islands, the correlation is so weak that the lines are not easily discerned. Although the model does not consider the species’ identity, only species richness, visual inspection of Figure 2c,d reveals clusters of islands in the Azores and in the Canary Islands. These clusters are clearly based on proximity among the islands; thus, independently of species identity, species richness alone is highly influenced by proximity.

3.2. Considering All Explanatory Variables

We next considered models with each of the explanatory variables, “area”, “elevation”, and “distance to the mainland” separately and all together (we call the latter the “full” model). We fitted the models for different types of species: all species, native, endemic, SIE, and MIE species. Table 6, Table 7, Table 8, Table 9 and Table 10 show the results of the WAIC scores for these models; to simplify the presentation, we only show the results for models with WAIC “weights” larger than 0.1; Supplementary Materials S3 and S4 show the full results. These tables reveal that the best models are always the ones with the Gaussian process (despite the increased number of parameters) and Poisson likelihood. (Therefore, in the following, we identify the models solely by their explanatory variables and omit that they include the Gaussian process and that the likelihood is a Poisson distribution).
For all species (Table 6) and for the Azores, the model with the lowest WAIC was the full one with distance to the mainland modelled by an exponential. However, the ΔWAIC reveals that four other models also have considerable support. For the Canary Islands, the area model had the smallest WAIC, followed by the elevation model, and according to the ΔWAIC, the four other models also have considerable support.
For the native species (Table 7) and for the Azores, the model with the lowest WAIC was the elevation-alone model. However, all the other models with the Gaussian process have very similar WAIC, as revealed by the ΔWAIC and the respective “weights”. A similar situation occurs for the Canary Islands, but now the model with the smallest WAIC is the “fulle” one, immediately followed by the elevation-alone one.
For the endemic species (Table 8) and for the Azores, the area-alone model has the smallest WAIC. For the Canary Islands, the best is the full model, with distance to the mainland modelled by an exponential, but the area, elevation, and distance models also have substantial support. As before, the models with “elevation” have a smaller ΔWAIC for the Canaries than for the Azores.
For the SIE species (Table 9) and for the Azores, the full model, with distance to the mainland modelled by an exponential, has the smallest WAIC. The second-best model has much less support. As in the previous cases, for the Canary Islands, there are several models with similar WAICs. The model with lower WAIC is the distance to the mainland only, modelled with a power–law decay.
Finally, for the MIE species (Table 10) and for the Azores, the models with more support are the area-only and the “fulle”, and for the Canaries, the one with smallest WAIC is the elevation-only model, although all the others also show significant support. Again, the ΔWAIC among models is smaller for the Canaries than it is for the Azores. Continuing the previous trend, the model with “elevation” only tends to exhibit lower ΔWAIC for the Canary Islands than for the Azores.
Analysis of the correlation matrices of the best models for the MIE and SIE species further highlights the importance of considering spatial autocorrelation. Because MIE species occur in several islands, dispersal is likely to determine their geographical distributions and proximity between islands should facilitate dispersal. On the other hand, the number of SIE species in an island should be less autocorrelated with the number of SIE in surrounding islands, but it would be incorrect to assume that there is no correlation, because proximity among islands may imply more similar geological ages or environmental conditions that may influence SIE species richness. In any case, we expect MIE species to display higher levels of positive spatial autocorrelation than SIE species. If that is the case, the matrix of correlation of the MIE species should have larger values than that of the SIE species. This is indeed what is observed, a result that becomes more evident if we divide each element of the MIE correlation matrix by its corresponding element in the SIE correlation matrix; see Table 11 and Table 12 (see also Tables S7–S10).

4. Discussion

The biogeographical processes that shape the species richness of an archipelago are complex and disentangling the effects of the driving factors behind these processes is challenging [1,4]. While island area and distance to the mainland are well-studied driving factors, others, such as the distance between islands within an archipelago, are not. Although dispersal between islands is widely recognized as extremely important for determining their species richness [16,46,47], the original ETIB only considered the distance between islands and the mainland. The novelty of our work, in the context of species richness in island biogeography, consists of also including the distance between the islands, and hence considering spatial autocorrelation using a Gaussian process (see also [48]). Our analyses demonstrated that spatial autocorrelation should not be left unaddressed in island biogeography studies for at least the following four reasons: (i) it affects the estimation of the parameters (e.g., their standard errors), (ii) reveals characteristics of the archipelagos (e.g., clustering of islands), (iii) helps identify differences between different types of richness measures (SIE versus MIE), and (iv) reveals differences between archipelagos (those where spatial autocorrelation plays a more important role). Indeed, ignoring spatial autocorrelation can eventually lead to wrong conclusions and even opposite patterns [17,49].
The formulation adopted [21] allows for the inclusion of several variables. In fact, although important, island area and distances to the mainland or between islands are not the only determinants of species richness in an archipelago. Under the GDM [8,50], the age of an island determines its topographic complexity and area, that is, the latter are not constant over time. There are no data to assess the temporal dynamics of species richness in the Azores and Canary archipelagos at the time scales relevant for the GDM. However, an important insight from GDM theory is that topographical complexity determines species richness. In fact, there are a large number of theories arguing for a different range of determinants of species richness as a function of topographical isolation and elevation ([12,51,52,53]), with some arguing that “diversity begets diversity” [10] (but see [54]). Our work may add to this debate by providing a method to include spatial autocorrelation. Overall, elevation seems to have a more dominant role in the Canary Islands than it has in the Azores. This may happen because the average maximum altitude is higher in the Canary Islands (~1800 m) than in the Azores (~1000 m) and because, even for the same range of altitudes, the climatic conditions are more homogeneous in the Azores [13].
An additional advantage of including spatial autocorrelation using a Gaussian process was the identification of island clusters within archipelagos through the correlation matrix [21]. We suggest that this method can be used routinely in the future to identify clusters. The correlation matrices also allowed the identification of the relative importance of spatial autocorrelation for SIE and MIE species richness. As expected, the correlation among islands is higher for MIE species, likely reflecting the importance of dispersal. Nevertheless, the correlation was different from zero among some islands for SIE species, probably due to similar environmental conditions and histories among nearby islands.
Finally, note that our approach only dealt with the number of species, not with their identities. Future work assessing the degree of community similarity as a function of distance is also possible within a Bayesian framework [28,29] and we plan to address this in future work. In summary, we developed a methodology to estimate species richness in an archipelago using Bayesian methods, paying special attention to the spatial autocorrelation, i.e., the influence of the distance between islands, which we modelled using a Gaussian process. We applied our methodology to the Azores and Canary Islands arthropod diversity, in two archipelagos of Macaronesia, and our approach highlighted the differences between them. The existence of off-the-shelf methods to identify spatial autocorrelations, as our results showed, warrants the conclusion that spatial autocorrelation should always be initially considered when studying the species richness of archipelagos, even if this analysis ends up showing that it is not present, which in itself would be an important conclusion.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/d15020127/s1, Information on the geographical and geological characteristics, and species richness, of the islands of the Azores and Canary archipelagos: Supplementary_Material_S1.pdf; selection of the explanatory variables: Supplementary_Material_S2.pdf; results for the models with the lowest WAIC: Supplementary_Material_S3.pdf; results for all models: Supplementary_Material_S4.pdf; R and Stan codes: Supplementary_Material_S5.zip.

Author Contributions

D.D.B. and L.B.-d.-Á. developed the Bayesian modelling framework and analysed the model output. P.A.V.B. was responsible for data collection and curation. M.d.L.M. provided important expertise in biogeography and ecology. D.D.B. and L.B.-d.-Á. wrote the first draft of the manuscript, while all authors contributed substantially to the revisions. All authors have read and agreed to the published version of the manuscript.

Funding

D.B. was funded by an FCT grant as part of the Biology and Ecology of Global Changes doctoral program. L.B.A. was financed under the Norma Transitória—L57/2016/CP1440/CT0022. P.A.V.B. worked on the manuscript under the framework of the Project FCT—Fundação para a Ciência e a Tecnologia, I.P., under the project UIDP/05292/2020 and UIDB/05292/2020. M.L.M. was supported by FCT/MCTES (UIDP/50017/2020 + UIDB/50017/2020 + LA/P/0094/2020).

Data Availability Statement

All the data used in this work are listed in Supplementary Material S1.

Acknowledgments

We thank Michael Betancourt, Giovani Silva, and Filipe Dias for their discussions on Bayesian theory and its implementation.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Matthews, T.J.; Triantis, K.A.; Whittaker, R.J. The Species–Area Relationship: Theory and Application; Cambridge University Press: Cambridge, UK, 2021. [Google Scholar]
  2. Whittaker, R.J.; Fernández-Palacios, J.M. Island Biogeography: Ecology, Evolution, and Conservation; Oxford University Press: Oxford, UK, 2007. [Google Scholar]
  3. Warren, B.H.; Simberloff, D.; Ricklefs, R.E.; Aguilée, R.; Condamine, F.L.; Gravel, D.; Morlon, H.; Mouquet, N.; Rosindell, J.; Casquet, J.; et al. Islands as Model Systems in Ecology and Evolution: Prospects Fifty Years after MacArthur-Wilson. Ecol. Lett. 2015, 18, 200–217. [Google Scholar] [CrossRef] [PubMed]
  4. Santos, A.M.; Field, R.; Ricklefs, R.E. New Directions in Island Biogeography. Glob. Ecol. Biogeogr. 2016, 25, 751–768. [Google Scholar] [CrossRef]
  5. Borregaard, M.K.; Amorim, I.R.; Borges, P.A.; Cabral, J.S.; Fernández-Palacios, J.M.; Field, R.; Heaney, L.R.; Kreft, H.; Matthews, T.J.; Olesen, J.M.; et al. Oceanic Island Biogeography through the Lens of the General Dynamic Model: Assessment and Prospect. Biol. Rev. 2016, 92, 830–853. [Google Scholar] [CrossRef] [PubMed]
  6. MacArthur, R.H.; Wilson, E.O. The Theory of Island Biogeography: Monographs in Population Biology; Princeton University Press: Princeton, NJ, USA, 1967. [Google Scholar]
  7. MacArthur, R.H.; Wilson, E.O. Equilibrium-Theory of Insular Zoogeography. Evolution 1963, 17, 373. [Google Scholar] [CrossRef]
  8. Whittaker, R.J.; Triantis, K.A.; Ladle, R.J. A General Dynamic Theory of Oceanic Island Biogeography. J. Biogeogr. 2008, 35, 977–994. [Google Scholar] [CrossRef]
  9. Witt, C.C.; Maliakal-Witt, S. Why Are Diversity and Endemism Linked on Islands? Ecography 2007, 30, 331–333. [Google Scholar] [CrossRef]
  10. Emerson, B.C.; Kolm, N. Species Diversity Can Drive Speciation. Nature 2005, 434, 1015–1017. [Google Scholar] [CrossRef] [PubMed]
  11. Irl, S.D.; Harter, D.E.; Steinbauer, M.J.; Gallego Puyol, D.; Fernández-Palacios, J.M.; Jentsch, A.; Beierkuhnlein, C. Climate vs. Topography–Spatial Patterns of Plant Species Diversity and Endemism on a High-Elevation Island. J. Ecol. 2015, 103, 1621–1633. [Google Scholar] [CrossRef]
  12. Steinbauer, M.J.; Otto, R.; Naranjo-Cigala, A.; Beierkuhnlein, C.; Fernández-Palacios, J.-M. Increase of Island Endemism with Altitude–Speciation Processes on Oceanic Islands. Ecography 2012, 35, 23–32. [Google Scholar] [CrossRef] [Green Version]
  13. Triantis, K.A.; Hortal, J.; Amorim, I.; Cardoso, P.; Santos, A.M.; Gabriel, R.; Borges, P.A. Resolving the Azorean Knot: A Response to Carine & Schaefer (2010). J. Biogeogr. 2012, 39, 1179–1184. [Google Scholar]
  14. May, R.M. Patterns of Species Abundance and Diversity. In Ecology and Evolution of Communities; Cody, M.L., Diamond, J.M., Eds.; Harvard University Press: Cambridge, MA, USA, 1975; pp. 81–120. [Google Scholar]
  15. Diamond, J.M. Distributional Ecology of New Guinea Birds. Science 1973, 179, 759–769. [Google Scholar] [CrossRef]
  16. Carvalho, J.C.; Cardoso, P.; Rigal, F.; Triantis, K.A.; Borges, P.A. Modeling Directional Spatio-Temporal Processes in Island Biogeography. Ecol. Evol. 2015, 5, 4671–4682. [Google Scholar] [CrossRef] [PubMed]
  17. Kreft, H.; Jetz, W.; Mutke, J.; Kier, G.; Barthlott, W. Global Diversity of Island Floras from a Macroecological Perspective. Ecol. Lett. 2008, 11, 116–127. [Google Scholar] [CrossRef] [PubMed]
  18. Selmi, S.; Boulinier, T. Ecological Biogeography of Southern Ocean Islands: The Importance of Considering Spatial Issues. Am. Nat. 2001, 158, 426–437. [Google Scholar] [CrossRef]
  19. Keitt, T.H.; Bjørnstad, O.N.; Dixon, P.M.; Citron-Pousty, S. Accounting for Spatial Pattern When Modeling Organism-Environment Interactions. Ecography 2002, 25, 616–625. [Google Scholar] [CrossRef] [Green Version]
  20. Dormann, C.F.; McPherson, J.M.; Araújo, M.B.; Bivand, R.; Bolliger, J.; Carl, G.; Davies, R.G.; Hirzel, A.; Jetz, W.; Daniel Kissling, W.; et al. Methods to Account for Spatial Autocorrelation in the Analysis of Species Distributional Data: A Review. Ecography 2007, 30, 609–628. [Google Scholar] [CrossRef] [Green Version]
  21. McElreath, R. Statistical Rethinking. In Texts in Statistical Science; CRC Press: Boca Raton, FL, USA, 2020. [Google Scholar]
  22. Gelman, A.; Carlin, J.B.; Stern, H.S.; Dunson, D.B.; Vehtari, A.; Rubin, D.B. Bayesian Data Analysis, 3rd ed.; CRC Press: Boca Raton, FL, USA, 2014. [Google Scholar]
  23. Golding, N.; Purse, B.V. Fast and Flexible Bayesian Species Distribution Modelling Using Gaussian Processes. Methods Ecol. Evol. 2016, 7, 598–608. [Google Scholar] [CrossRef] [Green Version]
  24. Hadjipantelis, P.Z.; Jones, N.S.; Moriarty, J.; Springate, D.A.; Knight, C.G. Function-Valued Traits in Evolution. J. R. Soc. Interface 2013, 10, 20121032. [Google Scholar] [CrossRef] [Green Version]
  25. Talluto, M.V.; Mokany, K.; Pollock, L.J.; Thuiller, W. Multifaceted Biodiversity Modelling at Macroecological Scales Using Gaussian Processes. Divers. Distrib. 2018, 24, 1492–1502. [Google Scholar] [CrossRef] [Green Version]
  26. Borges, P.A.V.; Vieira, V.; Amorim, I.R.; Bicudo, N.; Fritzén, N.; Gaspar, C.; Heleno, R.; Hortal, J.; Lissner, J.; Logunov, D.; et al. List of arthropods (Arthropoda). In A List of the Terrestrial and Marine Biota from the Azores; Borges, P.A.V., Costa, A., Cunha, R., Gabriel, R., Gonçalves, V., Martins, A.F., Melo, I., Parente, M., Raposeiro, P., Rodrigues, P., Eds.; Princípia: Cascais, Portugal, 2010; pp. 179–246. [Google Scholar]
  27. Arechavaleta, M.; Rodríguez, S.; Zurita, N.; García, A. Lista de Especies Silvestres de Canarias. Hongos, Plantas y Animales Terrestres. 2009. Gobierno de Canararias 2010. Available online: https://mc-stan.org (accessed on 10 December 2021).
  28. Dias, F.S.; Betancourt, M.; Rodríguez-González, P.M.; Borda-de-Água, L. BetaBayes—A Bayesian Approach for Comparing Ecological Communities. Diversity 2022, 15, 858. [Google Scholar] [CrossRef]
  29. Dias, F.S.; Betancourt, M.; Rodríguez-González, P.M.; Borda-de-Água, L. Analysing the Distance Decay of Community Similarity in River Networks Using Bayesian Methods. Sci. Rep. 2021, 11, 21660. [Google Scholar] [CrossRef] [PubMed]
  30. Rosenzweig, M.L. Species Diversity in Space and Time, 1st ed.; Cambridge University Press: Cambridge, UK, 1995. [Google Scholar]
  31. Gelman, A.; Hill, J. Data Analysis Using Regression and Multilevel/Hierarchical Models; Cambridge University Press: Cambridge, UK, 2007. [Google Scholar]
  32. Edge, M.D. Statistical Thinking From Scratch: A Primer For Scientists; Oxford University Press: Oxford, UK, 2019; ISBN 0-19-882762-8. [Google Scholar]
  33. Rasmussen, C.E.; Williams, C.K.I. Gaussian Processes for Machine Learning; MIT Press: Cambridge, MA, USA, 2006. [Google Scholar]
  34. Borges, P.A.V.; Lobo, J.M.; de Azevedo, E.B.; Gaspar, C.S.; Melo, C.; Nunes, L.V. Invasibility and Species Richness of Island Endemic Arthropods: A General Model of Endemic vs. Exotic Species. J. Biogeogr. 2006, 33, 169–187. [Google Scholar] [CrossRef] [Green Version]
  35. Fernandopullé, D. Climatic Characteristics of the Canary Islands. In Biogeography and ecology in the Canary Islands; Springer: Berlin/Heidelberg, Germany, 1976; pp. 185–206. [Google Scholar]
  36. Ramalho, R.S.; Helffrich, G.; Madeira, J.; Cosca, M.; Thomas, C.; Quartau, R.; Hipólito, A.; Rovere, A.; Hearty, P.J.; Ávila, S.P. Emergence and Evolution of Santa Maria Island (Azores)—The Conundrum of Uplifted Islands Revisited. Bulletin 2017, 129, 372–390. [Google Scholar] [CrossRef] [Green Version]
  37. Azevedo, J.; Ferreira, M.P. The Volcanotectonic Evolution of Flores Island, Azores (Portugal). J. Volcanol. Geotherm. Res. 2006, 156, 90–102. [Google Scholar] [CrossRef] [Green Version]
  38. Prægel, N.-O.; Holm, P.M. Lithospheric Contributions to High-MgO Basanites from the Cumbre Vieja Volcano, La Palma, Canary Islands and Evidence for Temporal Variation in Plume Influence. J. Volcanol. Geotherm. Res. 2006, 149, 213–239. [Google Scholar] [CrossRef]
  39. van den Bogaard, P. The Origin of the Canary Island Seamount Province-New Ages of Old Seamounts. Sci. Rep. 2013, 3, 1–7. [Google Scholar] [CrossRef] [Green Version]
  40. Ancochea, E.; Huertas, M.J.; Hernán, F.; Brändle, J.L.; Alonso, M. Structure, Composition and Age of the Small Islands of Santa Luzia, Branco and Raso (Cape Verde Archipelago). J. Volcanol. Geotherm. Res. 2015, 302, 257–272. [Google Scholar] [CrossRef]
  41. Ramalho, R.S.; Brum da Silveira, A.; Fonseca, P.E.; Madeira, J.; Cosca, M.; Cachão, M.; Fonseca, M.M.; Prada, S.N. The Emergence of Volcanic Oceanic Islands on a Slow-Moving Plate: The Example of M Adeira I Sland, NE A Tlantic. Geochem. Geophys. Geosystems 2015, 16, 522–537. [Google Scholar] [CrossRef]
  42. Ávila, S.; Cachão, M.; Ramalho, R.; Botelho, A.; Madeira, P.; Rebelo, A.; Cordeiro, R.; Melo, C.; Hipólito, A.; Ventura, M.; et al. The Palaeontological Heritage of Santa Maria Island (Azores: NE Atlantic): A Re-Evaluation of Geosites in GeoPark Azores and Their Use in Geotourism. Geoheritage 2016, 8, 155–171. [Google Scholar] [CrossRef]
  43. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2015. [Google Scholar]
  44. Stan Development Team. Stan Modeling Language User’s Guide and Reference Manual; Version 2.17. 0 2017; Stan Development Team: Aberdeen, South Dakota, 2017. [Google Scholar]
  45. Burnham, K.P.; Anderson, D.R. Model Selection and Multimodel Inference. A Practical Information-Theoretic Approach Second Edition; Springer: Berlin/Heidelberg, Germany, 2002. [Google Scholar]
  46. Cowie, R.H.; Holland, B.S. Dispersal Is Fundamental to Biogeography and the Evolution of Biodiversity on Oceanic Islands. J. Biogeogr. 2006, 33, 193–198. [Google Scholar] [CrossRef]
  47. Williamson, M. Relationship of Species Number to Area, Distance and Other Variables. In Analytical Biogeography. An Integrated Approach to the Study of Animal and Plant Distributions; Myers, A.A., Giller, P.S., Eds.; Chapman and Hall: Boca Raton, FL, USA, 1988; pp. 91–115. [Google Scholar]
  48. Sanmartín, I.; Van Der Mark, P.; Ronquist, F. Inferring Dispersal: A Bayesian Approach to Phylogeny-Based Island Biogeography, with Special Reference to the Canary Islands. J. Biogeogr. 2008, 35, 428–449. [Google Scholar] [CrossRef]
  49. Kühn, I. Incorporating Spatial Autocorrelation May Invert Observed Patterns. Divers. Distrib. 2007, 13, 66–69. [Google Scholar] [CrossRef]
  50. Borregaard, M.K.; Matthews, T.J.; Whittaker, R.J. The General Dynamic Model: Towards a Unified Theory of Island Biogeography? Glob. Ecol. Biogeogr. 2016, 25, 805–816. [Google Scholar] [CrossRef]
  51. Borges, P.A.V.; Hortal, J. Time, Area and Isolation: Factors Driving the Diversification of Azorean Arthropods. J. Biogeogr. 2009, 36, 178–191. [Google Scholar] [CrossRef] [Green Version]
  52. Steinbauer, M.J.; Field, R.; Grytnes, J.-A.; Trigas, P.; Ah-Peng, C.; Attorre, F.; Birks, H.J.B.; Borges, P.A.V.; Cardoso, P.; Chou, C.-H.; et al. Topography-Driven Isolation, Speciation and a Global Increase of Endemism with Elevation. Glob. Ecol. Biogeogr. 2016, 25, 1097–1107. [Google Scholar] [CrossRef]
  53. Cutts, V.; Katal, N.; Löwer, C.; Algar, A.C.; Steinbauer, M.J.; Irl, S.D.; Beierkuhnlein, C.; Field, R. The Effect of Small-Scale Topography on Patterns of Endemism within Islands. Front. Biogeogr. 2019, 11, e43737. [Google Scholar] [CrossRef]
  54. Pereira, H.M.; Proenca, V.M.; Vicente, L. Does Species Diversity Really Drive Speciation? Ecography 2007, 30, 328–330. [Google Scholar] [CrossRef]
Figure 1. The archipelagos of the Azores (a) and of the Canary Islands (b); the total number of species in each island is shown in brackets.
Figure 1. The archipelagos of the Azores (a) and of the Canary Islands (b); the total number of species in each island is shown in brackets.
Diversity 15 00127 g001
Figure 2. Plots (a,b) show the fitting with the Gaussian process, GP, (blue line) and without (red line) for the island species–area relationship (ISAR) model. The shadow areas correspond to 95% credible intervals. The fitting lines were obtained with the mean values of the posterior of the parameters c and z of S = cAz. For the model with the Gaussian process, the fitting line corresponds only to S = cAz, i.e., without the Gaussian process term. Plots (c,d) are geographical locations of the archipelagos and the darkness of the lines among the islands represents the strength of the correlations (for some pairs of islands, the correlation value is so small that the line is not visible). The size of the dots is proportional to the area of the islands.
Figure 2. Plots (a,b) show the fitting with the Gaussian process, GP, (blue line) and without (red line) for the island species–area relationship (ISAR) model. The shadow areas correspond to 95% credible intervals. The fitting lines were obtained with the mean values of the posterior of the parameters c and z of S = cAz. For the model with the Gaussian process, the fitting line corresponds only to S = cAz, i.e., without the Gaussian process term. Plots (c,d) are geographical locations of the archipelagos and the darkness of the lines among the islands represents the strength of the correlations (for some pairs of islands, the correlation value is so small that the line is not visible). The size of the dots is proportional to the area of the islands.
Diversity 15 00127 g002aDiversity 15 00127 g002b
Table 1. Stan output for the ISAR model with and without the Gaussian process (GP) for the Azores. “SD” stands for standard deviation, “CV” for coefficient of variation, “2.50%” and “97.50%” correspond to the width of the credible interval, “n_eff” is the effective sample size, and “Rhat” stands for R ^ (see main text).
Table 1. Stan output for the ISAR model with and without the Gaussian process (GP) for the Azores. “SD” stands for standard deviation, “CV” for coefficient of variation, “2.50%” and “97.50%” correspond to the width of the credible interval, “n_eff” is the effective sample size, and “Rhat” stands for R ^ (see main text).
MeanSDCV2.50%97.50%n_effRhat
Without GPc100.6420.230.2063.33139.772231.00
z0.410.050.120.320.521721.00
ϕ4.531.990.441.669.361951.00
With GPc99.5319.160.1964.45136.282191.01
z0.420.060.140.310.561231.00
ϕ6.473.170.492.0514.542321.02
η20.230.482.090.001.393631.01
ρ21.561.831.170.026.224061.00
Table 2. Stan output for the ISAR model with and without the Gaussian process (GP) for the Canary Islands. “SD” stands for standard deviation, “CV” for coefficient of variation, “2.50%” and “97.50%” correspond to the width of the credible interval, “n_eff” is the effective sample size and “Rhat” stands for R ^ (see main text).
Table 2. Stan output for the ISAR model with and without the Gaussian process (GP) for the Canary Islands. “SD” stands for standard deviation, “CV” for coefficient of variation, “2.50%” and “97.50%” correspond to the width of the credible interval, “n_eff” is the effective sample size and “Rhat” stands for R ^ (see main text).
MeanSDCV2.50%97.50%n_effRhat
Without GPc101.9120.310.2061.69139.742301.00
z0.470.050.110.380.591721.01
ϕ2.931.470.500.896.002241.00
With GPc99.3120.40.2057.9136.312481.00
z0.460.100.220.180.63541.00
ϕ5.032.640.521.3011.492921.00
η20.580.801.340.022.691221.00
ρ21.371.961.430.027.792041.00
Table 3. Comparison of ISAR models using the “Widely Applicable Information Criterion” (WAIC). ΔWAIC for model i is calculated as WAICi—WAICmin, and “weight” is calculated as exp ( 0.5 Δ WAIC i ) / A l l   m o d e l s exp ( 0.5 Δ WAIC i ) . “GP” stands for Gaussian Process.
Table 3. Comparison of ISAR models using the “Widely Applicable Information Criterion” (WAIC). ΔWAIC for model i is calculated as WAICi—WAICmin, and “weight” is calculated as exp ( 0.5 Δ WAIC i ) / A l l   m o d e l s exp ( 0.5 Δ WAIC i ) . “GP” stands for Gaussian Process.
ISAR ModelWAICΔWAICWeight
AzoresWith GP127.000.000.56
Without GP127.500.500.44
Canary IslandsWith GP116.700.000.76
Without GP119.002.300.24
Table 4. Correlation matrix of the ISAR model with a Gaussian process for the Azores (median values). The grey shading areas highlight the three major groups of islands: the western, central, and eastern groups. The abbreviations mean: C—Corvo, Fl—Flores, Fa—Faial, P—Pico, G—Graciosa, SJ—São Jorge, T—Terceira, SM—São Miguel, SMa—Santa Maria.
Table 4. Correlation matrix of the ISAR model with a Gaussian process for the Azores (median values). The grey shading areas highlight the three major groups of islands: the western, central, and eastern groups. The abbreviations mean: C—Corvo, Fl—Flores, Fa—Faial, P—Pico, G—Graciosa, SJ—São Jorge, T—Terceira, SM—São Miguel, SMa—Santa Maria.
CFlFaPGSJTSMSMa
C1.0000.8310.0040.0010.0010.0000.0000.0000.000
Fl0.8311.0000.0050.0010.0010.0000.0000.0000.000
Fa0.0040.0051.0000.8100.4990.6540.1850.0000.000
P0.0010.0010.8101.0000.5670.8190.3460.0020.000
G0.0010.0010.4990.5671.0000.7400.5100.0020.000
SJ0.0000.0000.6540.8190.7401.0000.5600.0040.000
T0.0000.0000.1850.3460.5100.5601.0000.0400.001
SM0.0000.0000.0000.0020.0020.0040.0401.0000.371
SMa0.0000.0000.0000.0000.0000.0000.0010.3711.000
Table 5. Correlation matrix of the ISAR model with a Gaussian process for the Canary Islands (median values). The grey shading areas highlight the two major groups of islands: the western and eastern groups. The abbreviations mean: F—Fuerteventura, L—Lanzarote, GC—Gran Canaria, T—Tenerife, LG—La Gomera, EH—El Hierro, LP—La Palma.
Table 5. Correlation matrix of the ISAR model with a Gaussian process for the Canary Islands (median values). The grey shading areas highlight the two major groups of islands: the western and eastern groups. The abbreviations mean: F—Fuerteventura, L—Lanzarote, GC—Gran Canaria, T—Tenerife, LG—La Gomera, EH—El Hierro, LP—La Palma.
FLGCTLGEHLP
F1.0000.5060.0810.0030.0000.0000.000
L0.5061.0000.0080.0000.0000.0000.000
GC0.0810.0081.0000.3450.0780.0040.005
T0.0030.0000.3451.0000.5790.0970.166
LG0.0000.0000.0780.5791.0000.4490.436
EH0.0000.0000.0040.0970.4491.0000.317
LP0.0000.0000.0050.1660.4360.3171.000
Table 6. Comparison of the several models using the “Widely Applicable Information Criterion” (WAIC) criterion for all species. Only models with ‘weight’ greater than 0.01 are shown. See the caption of Table 3 for the explanation of the acronyms “ΔWAIC” and “weight”. “poisson” means that the model assumes a Poisson likelihood and “gp” means it included a Gaussian process (thus necessarily the distance among islands). “area”, “dist” and “elev” stand for the models that (in addition to the Gaussian process) also consider area, distance to the mainland, or elevation, respectively, while “full” stands for a model with all these three explanatory variables included. The “e” in “diste” and “fulle” means that the distance to the mainland entered the model through an exponential relationship, and “dist” and “full” through a power–law relationship.
Table 6. Comparison of the several models using the “Widely Applicable Information Criterion” (WAIC) criterion for all species. Only models with ‘weight’ greater than 0.01 are shown. See the caption of Table 3 for the explanation of the acronyms “ΔWAIC” and “weight”. “poisson” means that the model assumes a Poisson likelihood and “gp” means it included a Gaussian process (thus necessarily the distance among islands). “area”, “dist” and “elev” stand for the models that (in addition to the Gaussian process) also consider area, distance to the mainland, or elevation, respectively, while “full” stands for a model with all these three explanatory variables included. The “e” in “diste” and “fulle” means that the distance to the mainland entered the model through an exponential relationship, and “dist” and “full” through a power–law relationship.
ArchipelagoModelWAICΔWAICWeight
Azoresfulle.poisson.gp91.60.00.23
area.poisson.gp91.70.10.21
dist.poisson.gp91.90.30.20
diste.poisson.gp92.00.40.19
elev.poisson.gp92.20.50.17
Canary Islandsarea.poisson.gp78.30.00.20
elev.poisson.gp78.40.10.19
dist.poisson.gp78.50.30.17
fulle.poisson.gp78.70.40.16
full.poisson.gp78.90.70.14
diste.poisson.gp79.10.80.13
Table 7. Comparison of the several models using the “Widely Applicable Information Criterion” (WAIC) criterion for the native species. Only models with ‘weight’ greater than 0.01 are shown. See the captions of Table 3 and Table 6 for the explanation of the acronyms.
Table 7. Comparison of the several models using the “Widely Applicable Information Criterion” (WAIC) criterion for the native species. Only models with ‘weight’ greater than 0.01 are shown. See the captions of Table 3 and Table 6 for the explanation of the acronyms.
ArchipelagoModelWAICΔWAICWeight
Azoreselev.poisson.gp85.70.00.18
area.poisson.gp85.80.10.17
fulle.poisson.gp85.80.10.17
full.poisson.gp85.80.20.16
dist.poisson.gp85.90.20.16
diste.poisson.gp85.90.20.16
Canary Islandsfulle.poisson.gp62.10.00.20
elev.poisson.gp62.20.10.19
area.poisson.gp62.40.30.17
full.poisson.gp62.40.30.17
dist.poisson.gp62.80.80.14
diste.poisson.gp62.90.80.13
Table 8. Comparison of the several models using the “Widely Applicable Information Criterion” (WAIC) criterion for the endemic species. Only models with ‘weight’ greater than 0.01 are shown. See the captions of Table 3 and Table 6 for the explanation of the acronyms.
Table 8. Comparison of the several models using the “Widely Applicable Information Criterion” (WAIC) criterion for the endemic species. Only models with ‘weight’ greater than 0.01 are shown. See the captions of Table 3 and Table 6 for the explanation of the acronyms.
ArchipelagoModelWAICΔWAICWeight
Azoresarea.poisson.gp70.70.00.38
fulle.poisson.gp72.21.50.18
full.poisson.gp72.82.10.13
dist.poisson.gp73.22.50.11
diste.poisson.gp73.32.60.10
Canary Islandsfulle.poisson.gp70.80.00.20
area.poisson.gp70.90.10.18
elev.poisson.gp71.00.20.17
full.poisson.gp71.20.40.16
dist.poisson.gp71.30.50.15
diste.poisson.gp71.60.80.13
Table 9. Comparison of the several models using the “Widely Applicable Information Criterion” (WAIC) criterion for the single island endemic species (SIE). Only models with ‘weight’ greater than 0.01 are shown. See the captions of Table 3 and Table 6 for the explanation of the acronyms.
Table 9. Comparison of the several models using the “Widely Applicable Information Criterion” (WAIC) criterion for the single island endemic species (SIE). Only models with ‘weight’ greater than 0.01 are shown. See the captions of Table 3 and Table 6 for the explanation of the acronyms.
ArchipelagoModelWAICΔWAICWeight
Azoresfulle.poisson.gp57.80.00.45
diste.poisson.gp59.72.00.17
full.poisson.gp60.32.60.13
Canary Islandsdist.poisson.gp61.30.00.19
elev.poisson.gp61.40.10.18
diste.poisson.gp61.40.10.18
area.poisson.gp61.60.30.17
full.poisson.gp61.60.30.16
fulle.poisson.gp62.20.90.12
Table 10. Comparison of the several models using the “Widely Applicable Information Criterion” (WAIC) criterion for the multiple island endemic species (MIE). Only models with ‘weight’ greater than 0.01 are shown. See the captions of Table 3 and Table 6 for the explanation of the acronyms.
Table 10. Comparison of the several models using the “Widely Applicable Information Criterion” (WAIC) criterion for the multiple island endemic species (MIE). Only models with ‘weight’ greater than 0.01 are shown. See the captions of Table 3 and Table 6 for the explanation of the acronyms.
ArchipelagoModelWAICΔWAICWeight
Azoresarea.poisson.gp69.60.00.4.0
fulle.poisson.gp70.20.60.29
dist.poisson.gp72.12.60.11
elev.poisson.gp72.42.80.10
diste.poisson.gp72.42.90.10
Canary Islandselev.poisson.gp68.70.00.23
dist.poisson.gp68.80.10.22
fulle.poisson.gp68.90.20.20
area.poisson.gp68.90.30.20
diste.poisson.gp69.50.90.15
Table 11. Ratio of the correlation matrices for the best models for MIE and SIE species of the Azores. The abbreviations mean: C—Corvo, Fl—Flores, Fa—Faial, P—Pico, G—Graciosa, SJ—São Jorge, T—Terceira, SM—São Miguel, SMa—Santa Maria.
Table 11. Ratio of the correlation matrices for the best models for MIE and SIE species of the Azores. The abbreviations mean: C—Corvo, Fl—Flores, Fa—Faial, P—Pico, G—Graciosa, SJ—São Jorge, T—Terceira, SM—São Miguel, SMa—Santa Maria.
CFlFaPGSJTSMSMa
C1.001.140.000.000.000.000.000.000.00
Fl1.141.000.000.000.000.000.000.000.00
Fa0.000.001.001.202.741.7214.850.000.00
P0.000.001.201.002.201.175.130.000.00
G0.000.002.742.201.001.402.630.000.00
SJ0.000.001.721.171.401.002.250.000.00
T0.000.0014.855.132.632.251.00204.250.00
SM0.000.000.000.000.000.00204.251.004.54
SMa0.000.000.000.000.000.000.004.541.00
Table 12. Ratio of the correlation matrices for the best models for MIE and SIE species of the Canary Islands. The abbreviations mean: F—Fuerteventura, L—Lanzarote, GC—Gran Canaria, T—Tenerife, LG—La Gomera, EH—El Hierro, LP—La Palma.
Table 12. Ratio of the correlation matrices for the best models for MIE and SIE species of the Canary Islands. The abbreviations mean: F—Fuerteventura, L—Lanzarote, GC—Gran Canaria, T—Tenerife, LG—La Gomera, EH—El Hierro, LP—La Palma.
FLGCTLGEHLP
F1.001.445.0847.140.000.000.00
L1.441.0024.880.000.000.000.00
GC5.0824.881.001.885.2237.2335.4
T47.140.001.881.001.314.513.11
LG0.000.005.221.311.001.561.60
EH0.000.0037.234.511.561.001.99
LP0.000.0035.43.111.61.991.00
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Barros, D.D.; Mathias, M.d.L.; Borges, P.A.V.; Borda-de-Água, L. The Importance of Including Spatial Autocorrelation When Modelling Species Richness in Archipelagos: A Bayesian Approach. Diversity 2023, 15, 127. https://doi.org/10.3390/d15020127

AMA Style

Barros DD, Mathias MdL, Borges PAV, Borda-de-Água L. The Importance of Including Spatial Autocorrelation When Modelling Species Richness in Archipelagos: A Bayesian Approach. Diversity. 2023; 15(2):127. https://doi.org/10.3390/d15020127

Chicago/Turabian Style

Barros, Diogo Duarte, Maria da Luz Mathias, Paulo A. V. Borges, and Luís Borda-de-Água. 2023. "The Importance of Including Spatial Autocorrelation When Modelling Species Richness in Archipelagos: A Bayesian Approach" Diversity 15, no. 2: 127. https://doi.org/10.3390/d15020127

APA Style

Barros, D. D., Mathias, M. d. L., Borges, P. A. V., & Borda-de-Água, L. (2023). The Importance of Including Spatial Autocorrelation When Modelling Species Richness in Archipelagos: A Bayesian Approach. Diversity, 15(2), 127. https://doi.org/10.3390/d15020127

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop