README.md

Ecological networks are variable both in time and space [@Poisot2015SpeWhy; @Trojelsgaard2016EcoNet] - this variability motivated the emergence of methodology to compare ecological networks, including in a way that meshes with the core concept for the comparison of ecological communities, namely $\beta$-diversity [@Poisot2012DisSpe]. The need to understand network variability through partitioning in components equivalent to $\alpha$, $\beta$, and $\gamma$ diversities is motivated by the prospect to further integrate the analysis of species interactions to the analysis of species compositions. Because species that make up the networks do not react to their environment in the same way, and because interactions are only expressed in subsets of the environments in which species co-occurr, the $\beta$-diversity of networks may behave in complex ways, and its quantification is likely to be ecologically informative.

@Poisot2012DisSpe and @Canard2014EmpEva have suggested an approach to $\beta$-diversity for ecological networks which is based on the comparison of the number of shared and unique links among species within a pair of networks. Their approach differentiates this sharing of links between those established between species occurring in both networks, and those established with at least one unique species. This framework is expressed as the decomposition $\beta_{wn} = \beta_{os} + \beta_{st}$, namely the fact that network dissimilarity ($\beta_{wn}$) has a component that can be calculated directly from the dissimilarity of interactions between shared species ($\beta_{os}$), and a component that cannot ($\beta_{st}$). The $\beta_{st}$ component differs slightly from the others, in that it is a quantification of the relative rewiring to overall dissimilarity, and not an absolute measure of interaction turnover. Presumably, the value of these components for a pair of networks can generate insights about the mechanisms involved in dissimilarity, when interpreted within the context of species turnover and differences in network connectance.

This approach has been widely adopted since its publication, with recent examples using it to understand the effect of fire on pollination systems [@Baronio2021NatFir]; the impact of rewiring on spatio-temporal network dynamics [@Campos-Moreno2021ImpInt]; the effects of farming on rural and urban landscapes on species interactions [@Olsson2021IntPla]; the impact of environment gradients on multi-trophic metacommunities [@Ohlmann2018MapImp]; and as a tool to estimate the sampling completeness of networks [@Souza2021PlaSam]. It has, similarly, received a number of extensions, including the ability to account for interaction strength [@Magrach2017PlaNet], the ability to handle probabilistic ecological networks [@Poisot2016StrPro], and the integration into the Local Contribution to Beta Diversity [@Legendre2013BetDiv] approach to understand how environment changes drive network dissimilarity [@Poisot2017HosPar].

{#fig:conceptual}

Yet, the precise meaning of $\beta_{st}$, namely the importance of species turnover in the overall dissimilarity, has been difficult to capture, and a source of confusion for some practitioners. This is not particularly surprising, as this component of the decomposition responds to unique species introducing their unique interactions both between themselves, and with species that are common to both networks (@fig:conceptual). For this reason, it is important to come up with guidelines for the interpretation of this measure, and how to use it to extract ecological insights.

Furthermore, much like the definition of $\beta$-diversity in all its forms is a contentious topic amongst community ecologists [see e.g. @Tuomisto2010DivBet], the $\beta$-diversity of networks has been submitted to methodological scrutiny over the years. A synthesis of some criticisms, related to the correct denominator to use to express the proportion of different links, has recently been published [@Frund2021DisSpe]. It argues that the calculation of network dissimilarity terms as originally outlined by @Poisot2012DisSpe is incorrect, as it can lead to over-estimating the role of interactions between shared species in a network ("rewiring"), and therefore underestimate the importance of species turnover across networks. As mist-understanding either of these quantities can lead to biased inferences about the mechanisms generating network dissimilarity, it is important to assess how the values (notably of $\beta_{os}$, and therefore of $\beta_{st}$) react to methodological choices.

Here, I present a mathematical analysis of the @Poisot2012DisSpe method, explain how information about species turnover and link rewiring can be extracted from its decomposition, and conduct numerical experiments to guide the interpretation of the $\beta$-diversity values thus obtained (with a specific focus on $\beta_{st}$). These numerical experiments establish three core facts. First, the decomposition adequately captures the relative roles of species turnover and interaction rewiring; second, the decomposition responds to differences in network structure (like connectance) as expected; finally, the decomposition more accurately captures rewiring than the proposed alternative using a different denominator put forth by @Frund2021DisSpe.

Partitioning network dissimilarity

The approach to quantifying the difference between pairs of networks established in @Poisot2012DisSpe is a simple extension of the overall method by @Koleff2003MeaBet for species dissimilarity based on presence-absence data. The objects to compare, $X_1$ and $X_2$, are partitioned into three values, $a = |X_1 \cup X_2|$, $b = |X_2 \setminus X_1|$, and $c = |X_1 \setminus X_2|$, where $|\cdot|$ is the cardinality of set $\cdot$ (the number of elements it contains), and $\setminus$ is the set substraction operation. In the perspective of species composition comparison, $X_1$ and $X_2$ are the sets of species in either community, so that if $X_1 = {x, y, z}$ and $X_2 = {v, w, x, y}$, we have $X_1 \cup X_2 = {v, w, x, y, z}$, $X_1 \cap X_2 = {x, y}$, $X_2 \setminus X_1 = {v, w}$, and $X_1 \setminus X_2 = {z}$. The core message of @Koleff2003MeaBet is that the overwheling majority of measures of $\beta$-diversity can be re-expressed as functions that operate on the cardinality of these sets -- this allows to focus on the number of unique and common elements, as outlined in @fig:conceptual.

Re-expressing networks as sets

Applying this framework to networks requires a few additional definitions. Although ecologists tend to think of networks as their adjacency matrix (as is presented in @fig:conceptual), this representation is not optimal to reach a robust understanding of which elements should be counted as part of which set when measuring network dissimilarity. For this reason, we need fall back on the definition of a graph as a pair of sets, wherein $\mathcal{G} = (V, E)$. These two components $V$ and $E$ represent vertices (nodes, species) and edges (interactions), where $V$ is specifically a set containing the vertices of $\mathcal{G}$, and $E$ is a set of ordered pairs, in which every pair is composed of two elements of $V$; an element ${i,j}$ in $E$ indicates that there is an interaction from species $i$ to species $j$ in the network $\mathcal{G}$. The adjancency matrix $\mathbf{A}$ of this network would therefore have a non-zero entry at $A_{ij}$.

In the context of networks comparison (assuming the networks to compare are $\mathcal{M}$ and $\mathcal{N}$), we can further decompose the contents of these sets as

$$\mathcal{M} = (V_c \cup V_m, E_c \cup E_{sm} \cup E_{um}) ,,$$

and

$$\mathcal{N} = (V_c \cup V_n, E_c \cup E_{sn} \cup E_{un}) ,,$$

where $V_c$ is the set of common species, $V_m$ and $V_n$ are the species belonging only to network $m$ and $n$ (respectively), $E_c$ are the common edges, and $E_{sm}$ and $E_{um}$ are the interactions unique to $k$ involving, respectively, only species in $V_c$, and at least one species from $V_m$ (the same notation applies for the subscript $_{n}$).

Defining the partitions from networks as sets

The metaweb [@Dunne2006NetStr], which is to say the entire regional species pool and their interaction, can be defined as $\mathcal{M} \cup \mathcal{N}$ (this operation is commutative), which is to say

$$\mathcal{M} \cup \mathcal{N} = (V_c \cup V_m \cup V_n, E_c \cup E_{sm} \cup E_{um} \cup E_{sn} \cup E_{un}) ,.$$

This operation gives us an equivalent to $\gamma$-diversity for networks, in that the set of vertices contains all species from the two networks, and the set of edges contains all the interactions between these species. If, further, we make the usual assumption that only species with at least one interaction are present in the set of vertices, then all elements of the set of vertices are present at least once in the set of edges, and the set of vertices can be entire reconstructed from the set of edges. Although measures of network $\beta$-diversity operate on interactions (not species), this property is maintained at every decomposition we will describe next.

We can similarly define the intersection (also commutative) of two networks:

$$\mathcal{M} \cap \mathcal{N} = (V_c, E_c),.$$

The decomposition of $\beta$-diversity from @Poisot2012DisSpe uses these components to measure $\beta_{os}$ ("rewiring"), and $\beta_{wn}$ (the overall dissimilarity including non-shared species). We can express the components $a$, $b$, and $c$ of @Koleff2003MeaBet as the cardinality of the following sets:

Component	$a$	$b$	$c$
$\beta_{os}$	$E_c$	$E_{sn}$	$E_{sm}$
$\beta_{wn}$	$E_c$	$E_{sn} \cup E_{un}$	$E_{sm} \cup E_{um}$

It is fundamental to note that these components can be measured entirely from the interactions, and that the number of species in either network are never directly involved.

In the following sections, I present a series of calculations aimed at expressing the values of $\beta_{os}$, $\beta_{wn}$, and therefore $\beta_{st}$ as a function of species sharing probability (as a proxy for mechanisms generating turnover), and link rewiring probability (as a proxy for mechanisms generating differences in interactions among shared species). These calculations are done using Symbolics.jl [@Gowda2021HigSym], and subsequently transformed in executable code for Julia [@Bezanson2017JulFre], used to produce the figures.

Quantifying the importance of species turnover

The difference between $\beta_{os}$ and $\beta_{wn}$ stems from the species dissimilarity between $\mathcal{M}$ and $\mathcal{N}$, and it is easier to understand the effect of turnover by picking a dissimilarity measure to work as an exemplar. We will use $\beta = (b+c)/(2a+b+c)$, which in the @Koleff2003MeaBet framework is [@Wilson1984MeaBet]. This measure returns values in $[0,1]$, with $0$ meaning complete similarity, and $1$ meaning complete dissimilarity.

Based on a partition between three sets of cardinality $a$, $b$, and $c$,

$$\beta_t = \frac{b+c}{2a+b+c},.$$

Note that this measure is written as $\beta_t$ for consistency with @Koleff2003MeaBet. So as to simplify the notation of the following section, I will introduce a series of new variables. Let $C = |E_c|$ be the number of links that are identical between networks (as a mnemonic, $C$ stands for "common"); $R = |E_{sn} \cup E_{sm}|$ be the number of links that are not shared, but only involve shared species (i.e. links from $\mathcal{M}\cup\mathcal{N}$ established between species from $\mathcal{M}\cap\mathcal{N}$; as a mnemonic, $R$ stands for "rewired"); and $T = |E_{un} \cup E_{um}|$ the number of links that are not shared, and involve at least one unique species (as a mnemonic, $T$ stands for "turnover").

There are two important points to note here. First, as mentionned earlier, the number or proportion of species that are shared is not involved in the calculation. Second, the connectance of either network is not involved in the calculation. That all links counted in e.g. $T$ come from $\mathcal{M}$, or that they are evenly distributed between $\mathcal{M}$ and $\mathcal{N}$, has no impact on the result. This is a desirable property of the approach: whatever quantitative value of the components of dissimilarity can be interpreted in the light of the connectance and species turnover without any risk of circularity; indeed, I present a numerical experiment where connectance varies independently later in this manuscript, reinforcing this point.

The final component of network dissimilarity in @Poisot2012DisSpe is $\beta_{st}$, i.e. the part of $\beta_{wn}$ that is not explained by changes in interactions between shared species ($\beta_{os}$), and therefore stems from species turnover. This fraction is defined as $\beta_{st} = \beta_{wn}-\beta_{os}$. The expression of $\beta_{st}$ does not involve a partition into sets that can be plugged into the framework of @Koleff2003MeaBet, because the part of $\mathcal{M}$ and $\mathcal{N}$ that are composed of their unique species cannot, by definition, share interactions. One could, theoretically, express these as $\mathcal{M} \setminus \mathcal{N} = (V_m, E_{um})$ and $\mathcal{N} \setminus \mathcal{M} = (V_v, E_{un})$ (note the non-commutativity here), but the dissimilarity between these networks is trivially maximal for the measures considered.

Using the $\beta_t$ measure of dissimilarity, we can re-write (using the notation with $R$, $C$, and $T$)

$$\beta_{os} = \frac{R}{2C+R},,$$

and

$$\beta_{wn} = \frac{R+T}{2C+R+T},.$$

Note that $\beta_{os}$ has the form $x/y$ with $x = S$ and $y = 2A+S$, and $\beta_{wn}$ has the form $(x+k)/(y+k)$, with $k = U$. As long as $k \ge 0$, it is guaranteed that $\beta_{wn} \ge \beta_{os}$, and therefore that $0 \ge \beta_{st} \ge 1$; as $C$, $T$, and $R$ are cardinalities of sets, they are necessarily satisfying this condition.

We can get an expression for $\beta_{st}$, by bringing $\beta_{os}$ and $\beta_{wn}$ to a common denominator and simplifying the numerator:

$$\beta_{st} = \frac{2CT}{(2C+R)(2C+R+T)},.$$

Note that this value varies in a non-monotonic way with regards to the number of interactions that are part of the common set of species -- this is obvious when developing the denominator into $4C^2 + R^2 + 4CR + 2CT + RT$. As such, we expect that the value of $\beta_{st}$ will vary in a hump-shaped way with the proportion of shared interactions. For this reason, @Poisot2012DisSpe suggest that $\beta_{st}/\beta{wn}$ (alt. $1-\beta_{os}/\beta_{wn}$) is a better indicator of the relative importance of turnover processes on network dissimilarity. This can be calculated as

$$\frac{\beta_{st}}{\beta_{wn}} = \frac{2CT}{(2C+S)(2C+R+T)}\times\frac{R+T}{2C+R+T},,$$

which reduces to

$$\frac{\beta_{st}}{\beta_{wn}} = \frac{2CT}{(2C+R)(R+T)},.$$

The roots of this expression are $C=0$ (the turnover of species has no contribution to the difference between $\beta_{wn}$ and $\beta_{os}$ if there are no shared species, and therefore no rewiring), and for $T = 0$ (the turnover of species has no contribution if all species are shared).

Quantifying the response of network beta-diversity to souces of variation

The relative effect of species turnover and link rewiring

As the decomposition of beta diversity into sets presented above reveals, the value of the components $\beta_{os}$ and $\beta_{st}$ will respond to two family of mechanisms: the probability of sharing a species between the two networks, noted $p$, which will impose bounds on the value of $T$; and the probability of an interactions between shared species not being rewired, noted $q$, which will impose bounds on the value of $C$. These two probabilities represent, respectively, mechanisms involved in species turnover and link turnover, as per @Poisot2015SpeWhy, and the aim of this numerical experiment is to describe how these families of processes drive network dissimilarity.

In order to simplify the calculations, I make the assumptions that the networks have equal species richness (noted $S$), so that $S_1 = S_2 = S$, and the same connectance (noted $\rho$), so that $\rho_1 = \rho_2 = \rho$. As a consequence, the two networks have the same number of links $L = \rho\times S_1^2 = \rho\times S_2^2$. The assumption of equal connectance will be relaxed in a subsequent numerical experiment. These simplifications allow to express the size of $C$, $R$, and $T$ only as functions of $p$ and $q$, as they would all be multiplied by $L$, which can therefore be dropped from the calculation.

{#fig:turnrew}

The value of $C$ is the proportion of shared species $p^2$, as per @fig:conceptual, times the proportion of shared links, $q$, giving $C = qp^2$. Each network has $r = p^2-(qp^2)$ rewired links, which leads to $R = 2r = 2p^2(1-q)$. Finally, we can get the number of unique links in each network $t$ by substracting $C+r$ from the total number of links (which, since we scale everything by $L$, is 1), yielding $t = 1 - qp^2 - p^2 + qp^2$, which is $t = 1-p^2$. The total number of unique links due to turnover is $T = 2t = 2(1-p^2)$. It is important to note that $C$ and $R$, namely the number of links that are kept or rewired, depends on species sharing ($p$), as the possible size of the overlap between the two networks does, but the quantity of links that are different due to turnover does not depends on rewiring.

With the values of $C$, $R$, and $T$, we can write

$$\beta_{os} = \frac{2p^2(1-q)}{2p^2q+2p^2(1-q)} = \frac{1-q}{q + 1 -q} = (1-q),.$$

This is a first noteworthy result: the value of $\beta_{os}$, in the ideal scenario of equal links and richness, is the probability of link re-wiring. Because this is true regardless of the value of $p$ (species turnover), this makes $\beta_{os}$ a strongly ecologically informative component.

Similarly, we can write

$$\beta_{wn} = \frac{2p^2(1-q)+2(1-p^2)}{2p^2q + 2p^2(1-q)+2(1-p^2)} = \frac{p^2(1-q)+(1-p^2)}{p^2q+p^2(1-q)+(1-p^2)} = 1-qp^2,.$$

The overall dissimilarity responds to $q$ (rewiring) linerarly, and to $p$ quadratically (which is expected assuming unipartite networks, in which species are present on both sides).

Expressing $\beta_{os}$ and $\beta_{wn}$ as functions of $p$ and $q$ trivializes the search for the expression of $\beta_{st}$, which is

$$\beta_{st} = 1 - p^2q - 1 + q = q\times(1-p^2),.$$

It is worth examining this solution in some detail. $\beta_{st}$ scales linearly with the probability that a link will not be rewired -- in other words, in a pair of networks for which rewiring is important ($q$ goes to 0), species turnover is going to be a relatively less important mechanism to dissimilarity. $\beta_{st}$ increases when turnover is important ($p$ goes to 0), and therefore $\beta_{st}$ represents a balance between species turnover and link rewiring. These three values, as well as $\beta_{st}/\beta_{wn}$, are represented in @fig:turnrew.

Sensibility of the decomposition to differences in connectance

The results presented in @fig:turnrew include the strong assumption that the two networks have equal connectance. Although the range of connectances in nature tends to be very strongly conserved within a system, we can relax this assumption, by letting one network have more interactions than the other. Note that for the sake of notation simplicity, I maintain the constraint that the two networks are equally species rich. Therefore, the sole variation in this numerical experiment is that one network has $L_1 = \rho\times a\times S^2$, and the other network has $L_2 = \rho\times S^2$; in other words, $L_1 = a\times L$ and $L_2 = L$. As one step of the components calculations involves a $\text{min}$ operation, I will add the constraint that $L_1 \le L_2$, which is to say $0 < a \le 1$. The value of $a$ is the ratio of connectances of the two networks, and the terms $S^2$ and $\rho$ being shared across all factors, they will be dropped from the calculations.

The maximal number of links that can be shared is $ap^2$ (i.e. $\text{min}(p^2, ap^2)$), as we cannot share more links than are in the sparsest of the two networks. Of these, $q$ are not rewired, leading to $C = aqp^2$. The number of links that are rewired in network 1 is the number of its links between shared species minus $C$, i.e. $r_1 = ap^2 - aqp^2 = ap^2(1-q)$, and similarly $r_2 = p^2 - aqp^2 = p^2(1-aq)$, leading to $R = r_1 + r_2 = p^2 \left[a(1-q)+1\right]$. Using the same approach, we can get $t_1 = a(1-p^2)$ and $t_2 = (1-p^2)$, leading to $T = t_1 + t_2 = (1-p^2)(1+a)$.

As in the previous section, we can use these values to write

$$\beta_{os} = 1 - 2\frac{aq}{1+a},,$$

$$\beta_{wn} = 1 - 2\frac{ap^2q}{1+a},,$$

and

$$\beta_{st} = 2aq\frac{(1-p^2)(1+a)}{a^2 + 2a + 1},.$$

{#fig:connectance}

The values of these components are visualized in @fig:connectance. The introduction of the connectance ratio makes these expressions marginally more complex than in the case without differences in connectance, but the noteworthy result remains that in the presence of differences of connectance, the value of $\beta_{os}$ is still independent from species turnover. In fact, there is an important conclusion to be drawn from this expression. The shared species component is by definition square, meaning that from an actual measurement of $\beta_{os}$ between two networks for which we know the connectance, noted $\mathbf{b}{os}$, we can get the probability of rewiring by reorganizing the terms of $\mathbf{b}{os} = 1 - 2aq/(1+a)$ as

$$q \approx \frac{(1-\mathbf{b}_{os})(a+1)}{2a},,$$

which gives the probability of rewiring as $1-q$; note that this is an approximation, as it assumes that the connectances of the entire network and the connectances of the shared components are the same.

Does the partition of network dissimilarity needs a new normalization?

One of the arguments put forth in a recent paper by @Frund2021DisSpe is that the decomposition outlined above will overestimate the effect of rewiring; I argue that this is based on a misunderstanding of what $\beta_{st}$ achieves. It is paramount to clarify that $\beta_{st}$ is not a direct measure of the importance of turnover: it is a quantification of the relative impact of rewiring to overall dissimilarity, which, all non-turnover mechanisms being accounted for in the decomposition, can be explained by turnover mechanisms. In this section, I present two numerical experiments showing (i) that the $\beta_{os}$ component is in fact an accurate measure of rewiring, and (ii) that $\beta_{st}$ captures the consequences of species turnover, and of the interactions brought by unique species.

Illustrations on arbitrarily small networks are biased

We can re-calculate the illustration of @Frund2021DisSpe, wherein a pair of networks with two shared interactions ($C = 2$) receive either an interaction in $T$, in $R$, or in both:

$C$	$T$	$R$	$\beta_{os}$	$\beta_{wn}$	$\beta_{st}$	$\beta_{st}/\beta_{wn}$
2	0	0	0	0	0
2	1	0	$1/5$	$1/5$	0	0
2	0	1	0	$1/5$	$1/5$	0
2	1	1	$1/5$	$1/3$	$2/15$	$2/5$

The over-estimation argument hinges on the fact that $\beta_{st} < \beta_{os}$ in the last situation (one interaction as rewiring, one as turnover). Reaching the conclusion of an overestimation from this is based on a mis-interpretation of what $\beta_{st}$ means. The correct interpretation is that, out of the entire network dissimilarity, only three-fifths are explained by re-wiring. The fact that this fraction is not exactly one-half comes from the fact that the @Wilson1984MeaBet measure counts shared interactions twice (i.e. it has a $2C$ term), which over-amplifies the effect of shared interactions as the network is really small. Running the same calculations with $C = 10$ gives a relative importance of the turnover processes of 47%, and $\beta_{st}$ goes to $1/2$ as $C/(T+R)$ increases. As an additional caveat, the value of $\beta_{st}$ will depend on the measure of beta-diversity used. Measures that do not count the shared interaction twice are not going to amplify the effect of rewiring.

Based on the arguments presented above, I do not think the suggestion of @Frund2021DisSpe to change the denominator of $\beta_{os}$ makes sense as a default; the strength of the original approach by @Poisot2012DisSpe is indeed that the effect of turnover is based on a rigorous definition of networks as graphs (as opposed to networks as matrices), in which the induction of vertices from the edgelist being compared gives rise to biologically meaningful denominators. The advantage of this approach is that at no time does the turnover of species itself (or indeed, as shown in many places in this manuscript, the network richness), or the connectance of the network, enter into the calculation of the beta-diversity components. As such, it is possible to use $\beta_{os}$ and $\beta_{wn}$ in relationship to these terms, calculated externally [as was recently done by e.g. @Higino2021BetPhy], without creating circularities.

Therefore the argument of @Frund2021DisSpe, whereby the $\beta_{os}$ component should decrease with turnover, and be invariant to connectance, does not hold: the very point of the approach is to provide measures that can be interpreted in the light of connectance and species turnover. Adopting the perspective developed in the previous section, wherein networks are sets and the measures of $\beta$-diversity operates on these sets, highlights the conceptual issue in the @Frund2021DisSpe alternative normalization: they are using components (namely, interactions) of the networks that are not directly part of the two networks being compared.

Using an alternative normalization trivializes the results

In this numerical experiment, we reproduce the results in @fig:turnrew, but using the alternative normalization described above. The results are presented in @fig:commden. Producing the analytical solutions for the various components, following the expressions for $C$, $T$, and $R$ given for @fig:turnrew, yields a similar value for $\beta_{wn}$ (i.e. the two approaches estimate the same value for total dissimiliarity), but different values for $\beta_{st}$ and $\beta_{os}$. Specifically, $\beta_{os}$ becomes $p^2(1-q)$, which becomes dependent on species turnover. This, from an ecological point of view, makes no sense: the quantification of how much shared species interact in a similar way should not depend on how much species actually overlap. The opposite problem arises for $\beta_{st}$, which becomes $1-p^2$. In short, the relative importance of species turnover is simply species turnover itself, and has no information on interaction dissimilarity. Therefore the core issue of the @Frund2021DisSpe alternative is that, by attempting to fix a non-issue (namely the over-estimate of the importance of re-wiring, which is only true in trivially small networks), it blurs the meaning of $\beta_{os}$, and renders $\beta_{st}$ useless as it is a re-expression of species beta-diversity.

{#fig:commden}

Measuring network beta-diversity: recommendations

Based on the numerical experiments and the derivations presented in this paper, we can establish a number of recommendations for the measurement and analysis of network dissimilarity. First, $\beta_{os}$ allows to estimate the rate of rewiring, which is an important ecological information to have; quantifying it properly can give insights as to how networks differ. Second, $\beta_{st}$ captures both turnover and rewiring mechanisms, but its interpretation is easier to accomplish in the context of total network dissimilarity, and therefore $\beta_{st}/\beta_{wn}$ should be interpreted more thoroughly. Finally, because the alternative denominator from @Frund2021DisSpe removes the interesting property of $\beta_{os}$ (independent estimate of rewiring rate), and trivializes the meaning of $\beta_{st}$ (by turning it into species dissimilarity), there seems to be no valid reason to use it.

Conflict of interest disclosure: the authors of this article declare that they have no financial conflict of interest with the content of this article; TP is one of the PCIEcology recommenders.

References