Loading [MathJax]/jax/output/SVG/jax.js
 

Moments of von mises and fisher distributions and applications

  • Received: 09 May 2016 Accepted: 22 November 2016 Published: 01 June 2017
  • MSC : Primary: 35Q92; Secondary: 62E15, 92B05

  • The von Mises and Fisher distributions are spherical analogues to the Normal distribution on the unit circle and unit sphere, respectively. The computation of their moments, and in particular the second moment, usually involves solving tedious trigonometric integrals. Here we present a new method to compute the moments of spherical distributions, based on the divergence theorem. This method allows a clear derivation of the second moments and can be easily generalized to higher dimensions. In particular we note that, to our knowledge, the variance-covariance matrix of the three dimensional Fisher distribution has not previously been explicitly computed. While the emphasis of this paper lies in calculating the moments of spherical distributions, their usefulness is motivated by their relationship to population statistics in animal/cell movement models and demonstrated in applications to the modelling of sea turtle navigation, wolf movement and brain tumour growth.

    Citation: Thomas Hillen, Kevin J. Painter, Amanda C. Swan, Albert D. Murtha. Moments of von mises and fisher distributions and applications[J]. Mathematical Biosciences and Engineering, 2017, 14(3): 673-694. doi: 10.3934/mbe.2017038

    Related Papers:

    [1] Said G. Nassr, Amal S. Hassan, Rehab Alsultan, Ahmed R. El-Saeed . Acceptance sampling plans for the three-parameter inverted Topp–Leone model. Mathematical Biosciences and Engineering, 2022, 19(12): 13628-13659. doi: 10.3934/mbe.2022636
    [2] Rui Zhu, Yuhang Chen, Qianqian Yu, Siqing Liu, Jianjie Wang, Zhili Zeng, Liming Cheng . Effects of contusion load on cervical spinal cord: A finite element study. Mathematical Biosciences and Engineering, 2020, 17(3): 2272-2283. doi: 10.3934/mbe.2020120
    [3] Ibrahim Alkhairy . Classical and Bayesian inference for the discrete Poisson Ramos-Louzada distribution with application to COVID-19 data. Mathematical Biosciences and Engineering, 2023, 20(8): 14061-14080. doi: 10.3934/mbe.2023628
    [4] Jan Haskovec, Ioannis Markou . Exponential asymptotic flocking in the Cucker-Smale model with distributed reaction delays. Mathematical Biosciences and Engineering, 2020, 17(5): 5651-5671. doi: 10.3934/mbe.2020304
    [5] Amal S. Hassan, Najwan Alsadat, Christophe Chesneau, Ahmed W. Shawki . A novel weighted family of probability distributions with applications to world natural gas, oil, and gold reserves. Mathematical Biosciences and Engineering, 2023, 20(11): 19871-19911. doi: 10.3934/mbe.2023880
    [6] M. G. M. Ghazal, H. M. M. Radwan . A reduced distribution of the modified Weibull distribution and its applications to medical and engineering data. Mathematical Biosciences and Engineering, 2022, 19(12): 13193-13213. doi: 10.3934/mbe.2022617
    [7] Jiajun Zhang, Tianshou Zhou . Stationary moments, distribution conjugation and phenotypic regions in stochastic gene transcription. Mathematical Biosciences and Engineering, 2019, 16(5): 6134-6166. doi: 10.3934/mbe.2019307
    [8] Saralees Nadarajah . Remark on the Paper by Rao And Kakehashi (2005). Mathematical Biosciences and Engineering, 2006, 3(2): 385-387. doi: 10.3934/mbe.2006.3.385
    [9] Oualid Kafi, Nader El Khatib, Jorge Tiago, Adélia Sequeira . Numerical simulations of a 3D fluid-structure interaction model for blood flow in an atherosclerotic artery. Mathematical Biosciences and Engineering, 2017, 14(1): 179-193. doi: 10.3934/mbe.2017012
    [10] Ariel Cintrón-Arias, Carlos Castillo-Chávez, Luís M. A. Bettencourt, Alun L. Lloyd, H. T. Banks . The estimation of the effective reproductive number from disease outbreak data. Mathematical Biosciences and Engineering, 2009, 6(2): 261-282. doi: 10.3934/mbe.2009.6.261
  • The von Mises and Fisher distributions are spherical analogues to the Normal distribution on the unit circle and unit sphere, respectively. The computation of their moments, and in particular the second moment, usually involves solving tedious trigonometric integrals. Here we present a new method to compute the moments of spherical distributions, based on the divergence theorem. This method allows a clear derivation of the second moments and can be easily generalized to higher dimensions. In particular we note that, to our knowledge, the variance-covariance matrix of the three dimensional Fisher distribution has not previously been explicitly computed. While the emphasis of this paper lies in calculating the moments of spherical distributions, their usefulness is motivated by their relationship to population statistics in animal/cell movement models and demonstrated in applications to the modelling of sea turtle navigation, wolf movement and brain tumour growth.


    1. Introduction

    The motivation for the present paper stems from the modelling of oriented biological movement. In many studies of animal or cell movement, a primary point of concern is the cue (or cues) used to navigate, necessitating statistical correlation between a preferred movement direction and orientating information in the environment. Typical examples include the orientation of animals in response to the earth's magnetic field [31], or the migration of invasive tumour cells through their local extracellular environment [30]. In ecological applications a quasi two-dimensional scenario generally suffices -for example, the habitat of a land-bound population or the typical swimming depth of marine organisms -and in such cases an animal's bearing can be described by a single angular coordinate. Other situations, however, may require a fully three-dimensional description, such as the motion of parasitic nematode worms in the soil, the movement of plankton in the oceans, or the flight of insects and birds. For the migration of cell populations in tissues a three dimensional understanding of movement can be particularly crucial. For example, malignant cancers of the central nervous system (CNS) are characterized by their diffuse and heterogeneous growth, and discovering the pathways of invasion is necessary for predicting tumour extent. The organisation of the brain itself plays a significant part in invasion, with cells channelled along the aligned neural fiber tracts that characterise white matter. Consequently, a three-dimensional understanding of the brain's architecture and the movement of cells in response to it can potentially lead to targeted treatment protocols.

    To model these and other movement phenomena we will consider random walk models in the next section (Section 2). We will draw a connection between spherical distributions, also called directional distributions, and movement models for biological systems. We denote the spherical distribution by q(n), where n is a unit vector indicating an individuals' bearing. It turns out that the first moment of q(n) (i.e. the expectation E[q]) is proportional to the population drift, while the variance-covariance matrix Var[q] is proportional to the diffusion of the population. Hence the moments of q(n) provide immediate significance for interpreting population-level responses. In Section 3 we compute the first and second moments of the von Mises distribution in two dimensions. These moments are known, and we include our computation to introduce a method based on the divergence theorem. In Section 4 we apply this method to compute the first and second moment in three dimensions (in three dimensions the von Mises distribution is generally known as the Fisher or von Mises-Fisher distribution [20,1]1). To our knowledge, an explicit formula for the three dimensional variance is not known in the literature and here we provide such a formula. Moreover, this method generalises to higher space dimensions and other circular/spherical distributions. Finally, in Section 5 we illustrate the importance of these calculations to three applications of directional movement: the orientation of sea turtles in an ocean; the tendency of animals to follow or avoid linear landscape features such as roads or seismic lines; and, the three-dimensional anisotropic spread of glioma. In each of these cases we find that the von-Mises and/or Fisher distribution lies at the heart of the organizing centre for directed movement and that their moments provide key information on population-level distributions.

    1Not to be confused with the statistical F-distribution.


    2. The von Mises and Fisher distributions

    We consider distributions q(n) on the unit circle Sn1 in Rn. The von Mises, Fisher and n-dimensional von Mises-Fisher distributions relate to a given direction uSn1, which identifies a "preferred" direction. In two dimensions the von Mises distribution is given as

    q(n)=12πI0(k)eknu, (1)

    for nS1. Here k is the concentration parameter and I0(k), and later Ij(k), denote the modified Bessel functions of first kind of order j. A more familiar representation of the von-Mises distribution in two dimensions uses the polar angle of the vectors involved. If we represent n=(cosα,sinα)T and u=(cosμ,sinμ)T, then the two-dimensional von-Mises distribution is

    q(n)=12πI0(k)ekcos(αμ).

    We prefer the coordinate free vector notation in (1) as it generalizes easily and it is not coupled to a specific coordinate system. The same distribution in three dimensions is called the Fisher distribution and is given by

    q(n)=k4πsinh(k)eknu,nS2. (2)

    The extension to higher dimensions forms the n-dimensional von Mises-Fisher distribution, given as

    q(n)=kn/21(2π)n/2In/21(k)eknunSn1.

    The von Mises distribution occupies a predominant position in spherical statistics. As a form of 'wrapped' normal distribution it acts as the de facto standard for fitting oriented data and has received a great deal of attention (e.g. see [22]). Further, applying a linear combination (or mixture) of von Mises distributions can be used to fit more complicated datasets, such as skewed or multi-modal distributions, thereby extending its general utility. For a von Mises mixture with given unit vectors u1,,uN, we consider

    q(n)=Ni=1ci2πI0(ki)ekinui, (3)

    where the coefficients ci satisfy Ni=1ci=1. Of particular note is the symmetric, bimodal von Mises distribution

    qB(n)=14πI0(k)(eknu+eknu), (4)

    along with its three-dimensional bimodal Fisher version

    q(n)=k8πsinh(k)(eknu+eknu), (5)

    which define bidirectional distributions with two opposite and equal peaks. Figure 1 illustrates the standard von Mises distribution superimposed onto the unit circle, along with a bimodal von Mises distribution with peaks at (1,0)T and (1,0)T. Figure 2 illustrates the three dimensional distribution (2) for increasing values of k.

    Figure 1. Left: The unimodal von Mises distribution as a function of nS1 with a peak at (1,0)T. Right: The bimodal von Mises distribution as a function of nS1 with peaks at n=±(1,0)T.
    Figure 2. The Fisher distribution on the unit sphere. From left to right we plot the distribution for increasing values of concentration parameter k. In each case, the dominant direction is taken as u=(1/2,1/2,2/2).

    Given the prominent role of von Mises-Fisher distributions, we illustrate a novel method for computing their first and second moments. In the context of the application to animal/cell movement, this allows a direct line to be drawn between the statistical analysis of an experimentally-derived dataset and inputs for the predictive macroscopic-scale model (6). Typically, moments are derived via explicit trigonometric integration (e.g. [24,3]) whereas our method here utilises the divergence theorem: while the former calculations become increasingly cumbersome as the space dimension increases, the latter generalises relatively easily.

    To explain the use of the von Mises and Fisher distribution in population models, we consider a simple random walk model for oriented animal/cell movement paths in two and three dimensions. Random walks have an established history in the quantification and modelling of biological movement (e.g. see [27,37,7,6,38]). In the velocity-jump random walk, [28], a path is traced out as a series of 'runs' (continuous movement through space with constant velocity) punctuated with 'turns' (selecting a new velocity). The run and tumbling motion of flagellated bacteria such as E. coli provides the exemplar, yet more generally, it offers a natural description for movement, compatible with established methods of recording animal/cell tracks: datasets of turning angles, times etc. inform probability distributions for input into the velocity-jump model.

    In the formulation here (see [14,29,16,32] for further details and applications) we presume an individual (taken as a point) is characterized by its position (xRn for n=2,3) and velocity (vRn) at time t[0,). For simplicity we assume the individual moves with a constant speed s, but that it undergoes changes in its direction (nSn1) as a result of a Poisson process, where λ is the (assumed constant) turning rate parameter. At each turn, a new direction (nSn1) is selected according to a suitable directional distribution q(t,x,n). Then q(t,x,n) will satisfy

    Sn1q(t,x,n)dn=1.

    To obtain an efficient evolution equation for the macroscopic population density distribution c(x,t), scaling techniques can be employed: we refer to [16] for full details. Adopting a moment-closure approach yields a macroscopic drift-diffusion equation in the form

    c(t,x)t+(a(t,x)c(t,x))=:(D(t,x)c(t,x)),=ni,j=1xixj(Dijc), (6)

    where the colon symbol (:) is used to denote the contraction of two tensors. In the above D(t,x) defines a (symmetric) n×n-diffusion tensor and a(t,x) is the ndimensional drift (or advective) velocity. Formally, equation (6) describes the probability distribution for finding an individual at position x and time t; however, if the population is dispersed and interactions are negligible, it can be interpreted as a macroscopic population density distribution. Significantly, the macroscopic quantities a and D relate to the first two moments of the turning distribution q:

    • the drift a is proportional to the expectation of q,

    a(t,x)=sE[q]=sSn1nq(t,x,n)dn; (7)

    • the D is proportional to the variance-covariance matrix of q,

    D(t,x)=s2λVar[q]=s2λSn1(na(t,x)/s)(na(t,x)/s)Tq(t,x,n)dn. (8)

    Note that the product vvT denotes the dyadic product of two vectors (in tensor notation: vv) and therefore D forms a symmetric matrix.

    For the case in which q is unchanged across x and t, a and D will be given by a constant vector and diffusion tensor respectively. In this case we can calculate the fundamental solution of (6) as

    c(x,t)=1detD(4πt)nexp(14t(xat)TD1(xat)). (9)

    The above describes the probability density for finding a random walker (starting at location 0, time 0) at location x, time t.

    Clearly, the dynamics predicted by (6) pivot about the directional distribution q which, in turn, is dictated by data. Ideally, an explicit form for the moments of q can be calculated, allowing a direct translation between the inputs of the stochastic/individual model and the macroscopic model. Yet, while the statistical methodology for analysing oriented datasets has grown substantially in recent decades (e.g. see [3,10,22]), it remains substantially less advanced than that of classical (linear) statistics.


    3. Moments of the von Mises distribution

    Lemma 3.1. For a given unit vector uS1 consider the von-Mises distribution

    q(n)=12πI0(k)eknu,nS1.

    Then q has the following expectation and variance:

    E[q]=I1(k)I0(k)u; (10)
    Var[q]=12(1I2(k)I0(k))I2+(I2(k)I0(k)(I1(k)I0(k))2)uuT. (11)

    In the above, Ij denote modified Bessel functions of first kind and order j, I2 denotes the 2×2 identity matrix, and uuT denotes the dyadic product of two vectors (in tensor notation uu).

    Proof2. Before we compute the moments of q(n) we collect some properties of Bessel functions (see e.g. [15]). Let Jn(x) denote the Bessel functions of first kind, then

    2The computations in this section have already been published in an appendix to a book chapter [16]. We include them here for completeness and to introduce our method.

    In(x):=(i)nJn(ix)

    denote Bessel functions of first kind with purely imaginary argument, or the modified Bessel functions. The modified Bessel functions have an explicit representation as

    In(k)=12π2π0cos(nϕ)ekcosϕdϕ. (12)

    Bessel functions satisfy a differential relation

    ddx(xnJn(x))=xnJn1(x) (13)

    and, for n1, they satisfy a recurrence relation

    Jn+1(x)=2nxJn(x)Jn1(x). (14)

    The total mass of the unimodal von Mises distribution (1) can be directly computed from (12), where we denote the angle between n and u as ϕ:

    S1q(n)dn=12πI0(k)2π0ekcosϕdϕ=1.

    For the expectation,

    E[q]=S1nq(n)dn,

    we use a test vector bR2 and multiply it by E[q] (note that the expectation here is a vector quantity), using summation convention for repeated indices. Then we have

    bE[q]=bS1nq(n)dn,

    and

    b2πI0(k)E[q]=bS1neknudn,=S1ni(bieknu)dn,=B1(0)vi(biekvu)dv,=B1(0)bikuiekvudv,=bku102π0erkcosϕrdrdϕ,=bku102πrI0(rk)dr,=b2πku10rI0(rk)dr.

    In the first step we use the summation convention, which requires summation over repeated indices, i.e. aibi=ni=1aibi. In the third step we used the divergence theorem to transform a surface integral over S1 into a volume integral over the unit ball B1(0) of a divergence-term.

    The last integral can be solved via (13) and

    rI0(rk)=irkJ0(irk)ik=1ikddx(xJ1(x))|x=irk=1ikddr(rJ1(irk)).

    Then

    10rI0(rk)dr=1ikJ1(ik)=1ikiI1(k)=I1(k)k, (15)

    and due to the fact that b is an arbitrary test vector

    E[q]=I1(k)I0(k)u,

    which is (10).

    The variance-covariance matrix of q is given by

    Var[q]=S1(nE[q])(nE[q])Tq(n)dn=S1nnTq(n)dnE[q]E[q]T.

    To find the second moment of q we consider two test vectors a,bR2 and we use index notation for automatic summation over repeated indices:

    2πI0(k)aS1nnTq(n)dnb (16)
    =S1ainibjnjeknluldn,=S1ni(aibjnjeknlul)dn,=B1(0)vi(aibjvjekvlul)dv,=B1(0)aibiekvudv+B1(0)ai(vb)kuiekvudv,=abB1(0)ekvudv+kaubB1(0)vekvudv, (17)

    where, again, we use the summation convention and the divergence theorem. The first integral in (17) can be solved directly:

    B1(0)ekvudv=10S1erknurdrdn=102πrI0(rk)dr=2πI1(k)k,

    where we used (12) and (15) in the penultimate and ultimate step, respectively. Using (12) we can transform the second integral from (17) as follows:

    B1(0)vekvudv=10S1rnerknurdrdn=10r2S1nerknudndr,=2πu10r2I1(rk)dr, (18)

    where we used (10) in the last step.

    Now we use the differential recurrence relation (13) to write

    r2I1(rk)=1ik2(irk)2J1(irk)=1ik2ddx(x2J1(x))|x=irk=1kddr(r2J1(irk)).

    Continuing from (18) we find

    B1(0)vekvudv=2πu101kddr(r2J1(irk))dr=2πuJ2(ik)=2πuI2(k)k. (19)

    Substituting all the integrals back into equation (17):

    aS1nnTq(n)dnb=ab2πI1(k)k2πI0(k)+kau2πubI2(k)k2πI0(k),=a(1kI1(k)I0(k)I2+uuTI2(k)I0(k))b.

    Finally, we use the identity (14) for n=1 to replace

    1kI1(k)I0(k)=12(1I2(k)I0(k)),

    and the second moment is given by

    S1nnTq(n)dn=12I2+I2(k)I0(k)(uuT12I2). (20)

    Together with the formula for the expectation (10) we find

    Var[q]=S1nnTq(n)dnE[q]E[q]T,=12I2+I2(k)I0(k)(uuT12I2)(I1(k)I0(k))2uuT,=12(1I2(k)I0(k))I2+(I2(k)I0(k)(I1(k)I0(k))2)uuT,

    which is (11).

    In the two dimensional case, it is useful to formulate the result of Lemma 10 in terms of a polar angle. For u=(cosμ,sinμ)T we find

    E[q]=I1(k)I0(k)(cosμsinμ); (21)
     Var[q]=12(1I2(k)I0(k))(1001)+(I2(k)I0(k)(I1(k)I0(k))2)(cos2μcosμsinμcosμsinμsin2μ). (22)

    As outlined earlier, for many applications it is useful to compute the first and second moments for the bimodal von Mises distribution.

    Corollary 1. For given uS1 consider the bimodal von Mises distribution

    q(n)=14πI0(k)(eknu+eknu),

    then

    E[q]=0, (23)
    Var[q]=12(1I2(k)I0(k))I2+I2(k)I0(k)uuT. (24)

    Proof. Since the bimodal von Mises distribution is symmetric we have E[q]=0 and Var[q]=nnTq(n)dn. We apply formula (20) for each of the components eknu and eknu separately and sum.

    It is interesting to consider the limiting scenario as the parameter of concentration k becomes small (i.e. k0). Here, the directional bias becomes increasingly small and (1) or (4) limit to a uniform distribution. A straightforward exploration reveals that I1(k)/I0(k),I2(k)/I0(k)0. Consequently, E[q]0 while the variance-covariance Var[q]12I2 for both the unimodal and bimodal cases.


    4. The 3-D Fisher distribution

    The main result of this paper is as follows.

    Theorem 4.1. For given uS2 consider the Fisher distribution

    q(n)=k4πsinhkekun,nS2.

    Then

    E[q]=(cothk1k)u, (25)
    Var[q]=(cothkk1k2)I+(1cothkk+2k2coth2k)uuT. (26)

    Proof. To compute the mean direction

    E[q]=S2nq(n)dn,

    we use a test-vector bR3 and multiply it by E[q]. Using summation convention for repeated indices:

    bE[q]=bS2nq(n)dn,=S2ni(biq(n))dn,=B1(0)vi(bik4πsinhkekulvl)dv,=biB1(0)k2ui4πsinhkekuvdv,=bukB1(0)k4πsinhkekuvdv,=buk10S2k4πsinhkerkunr2drdn,=buk10r2k4πsinhk4πsinhrkrkdr,=buksinhk(1k2sinhk+1kcoshk),=bu(cothk1k),

    where in the 5th step we used the normalization condition for a Fisher distribution with concentration parameter rk, i.e.

    S2erkundn=4πsin(rk)rk.

    Since the test vector b was arbitrary, we proved (25). This formula coincides with formula (9.3.9) in Mardia [22].

    To compute the second moment, we take two test-vectors a,bR3 and compute

    aS2nnTq(n)dnb=S2aininjbjq(n)dn,=S2ni(ainjbjq(n))dn,=B1(0)vi(aivjbjk4πsinhkekulvl)dv,=B1(0)aibik4πsinhkekuvdv+B1(0)ai(vb)k2ui4πsinhkekuvdv,=abB1(0)k4πsinhkekuvdv(A)+aubB1(0)k24πsinhkvekuvdv(B)). (27)

    From the computation of the first moment, we find for the first integral (A) that

    B1(0)k4πsinhkekuvdv=1kcothk1k2. (28)

    The second integral (B) can be computed as

    B1(0)vekuvdv=10S2r3nerkundndr.

    From the first moment we obtain

    S2nerkundn=[4πsinhrkrk(cothrk1rk)]u.

    Hence

    B1(0)vekuvdv,=[10r34πsinhrkrk(cothrk1rk)dr]u,=4π[10r2coshrkkdr10rsinhrkk2dr]u,=4π[sinhkk23coshkk3+3sinhkk4]u.

    Then

    k24πsinhkB1(0)vekuvdv=(13kcothk+3k2)u. (29)

    Combining (28) and (29) with (27) we obtain

    annTq(n)dnb=ab(cothkk1k2)+auub(13cothkk+3k2),

    and the second moment is given by

    nnTq(n)dn=(cothkk1k2)I+(13cothkk+3k2)uuT. (30)

    To compute the full variance-covariance matrix of q we use the formula Var[q]=nnTq(n)dnE[q]E[q]T to obtain

    Var[q]=(cothkk1k2)I+(1cothkk+2k2coth2k)uuT,

    which is (26). We were not able to find this formula elsewhere in the literature.

    Corollary 2. The first two moments of the bimodal Fisher distribution

    q(n)=k8πsinhk(ekun+ekun)=k4πsinhkcosh(kun)

    are given by

    E[q]=0, (31)
    Var[q]=(cothkk1k2)I+(13cothkk+3k2)uuT. (32)

    Proof. As for the bimodal von Mises distribution, we exploit the symmetry properties by considering the components eknu and eknu separately and summing.

    It can be revealing to consider the coefficients arising in (25, 26, 32) and their dependence on the concentration parameter k. Consider the following four coefficients as functions of k:

    a(k)=cothk1k;b(k)=cothkk1k2;c(k)=1cothkk+2k2coth2k;d(k)=13cothkk+3k2. (33)

    Then for the unimodal case

    E[q]=a(k)u,Var[q]=b(k)I+c(k)uuT, (34)

    and for the bimodal case

    E[q]=0,Var[q]=b(k)I+d(k)uuT. (35)

    Thus, the variance-covariance matrix is split into an isotropic component, determined by b(k), and an anisotropic component, determined by c(k) or d(k), respectively. In two dimensions as the (unimodal and bimodal) von Mises distribution tends towards a uniform distribution (k0) we found that the expectation vanishes while the variance-covariance becomes an isotropic matrix. While substituting k=0 directly into (25, 26, 32) is not possible, the limit as k0 for the various coefficients above can be computed. Using the above notation,

    limk0+a(k)=limk0+cothk1k,

    which can be rewritten using l'Hȏpital's rule as

    limk0+a(k)=limk0+ktanhkktanhk,=limk0+1sech2ktanhk+ksech2k,=limk0+2sech2ktanhksech2k+(12sechk)sech2k,=0.

    Continuing in a similar fashion for b(k), we have

    limk0+b(k)=limk0+cothkk1k2.

    This can be rewritten with repeated use of l'Hȏpital's rule as

    limk0+b(k)=limk0+ktanhkk2tanhk,=limk0+1sech2k2ktanhk+k2sech2k,=limk0+2sech2ktanhk2ksech2k+2tanhk2k2sech2ktanhk+2ksech2k,=limk0+4sech2ktanh2k+2sech4k6sech2k12ksech2ktanhk2k2(sech4k2sech2ktanh2k),=13.

    Then we have

    limk0+c(k)=limk0+1cothkk+2k2coth2k,

    which we may write as

    limk0+1(cothkk1k2)b(k)+1k2coth2k=limk0+23+(1k+cothk)(1kcothk).

    In the above we have used the fact that the limit of b(k) is 1/3. Taking a factor of 1/k from the first factor to the second, we have

    limk0+c(k)=limk0+23+(1k+cothk)(1kcothk),=limk0+23(1+kcothk)(b(k)),=limk0+2313(1+ksinhk),=2323=0,

    where we used l'Hȏpital's rule in the last step. Finally, for d(k) we have

    limk0+d(k)=limk0+13cothkk+3k2,=limk0+13(cothkk1k2),=13limk0+b(k),=13(13)=0.

    Summarising, for either the unimodal or bimodal Fisher distribution, as k0 we find

    a=0andD=τs213I3

    as would be expected.

    These results are demonstrated numerically by plotting the coefficients as a function of k, see Figure 3: as k decreases to 0, the directional distribution tends to the uniform distribution and the variance-covariance matrix for either the unimodal or bimodal case limits to an isotropic tensor, as c(k) and d(k) tend to 0. Furthermore, we can clearly see from these plots that as anisotropy increases (k increases), the expectation of the uni-modal distribution becomes more strongly aligned with u.

    Figure 3. Coefficients a(k),b(k),c(k),d(k) from (33) as functions of k.

    As a final note of interest, it is remarked that for the unimodal case (34) the coefficient c(k) that describes the anisotropic diffusion component is negative. This effectively means reduced diffusion in the dominant direction u: the unimodal von-Mises distribution describes enhanced movement in direction u, with a correspondingly reduced movement in direction u, thereby reducing the variance in this direction.


    5. Numerical tests and applications

    We use this section to briefly highlight three interesting applications. In the interest of brevity we refer to the literature for biological and modelling details. We first compare our macroscopic PDE model (6) to an explicit random walk simulation in Section 5.1. The results show excellent overlap between the continuous model and the Monte Carlo simulation. In Section 5.2 we consider orientation of sea turtles as described by the unimodal two dimensional von Mises distribution. Section 5.3 considers the use of a bimodal two-dimensional von Mises distribution for animal movement patterns in disturbed forest environments. Finally, in Section 5.4, we use a bimodal three-dimensional Fisher distribution to model brain tumour growth. In all cases, simulations of the population-level animal/cell density distributions exploit the explicit moment calculations derived in Sections 3 and 4.


    5.1. Fixed, homogeneous orienteering cue

    To illustrate anisotropic diffusion we consider the movement of individuals released at x=0 at t=0 and subjected to a fixed, spatially uniform orientating cue. Consequently, q is fixed with respect to x and t and the macroscopic density distribution predicted by (6) will be determined by the fundamental solution (9), with a and D given by (7)-(8).

    We initially consider a 2D scenario, where q is given by a von Mises distribution with ϕ=π/4 and κ=10. In Figure 4 the density distribution for particle position is generated through repeated (100000) simulations of the individual model and compared against the fundamental solution (9) with a and D determined by equations (10) and (11). The correspondence between solutions validates both the derivation of (6) as well as the calculations of a and D. In Figure 5 we show a similar comparison for a 3-dimensional scenario, using the Fisher distribution (2) and correspondingly the moments (25-26) for a and D.

    Figure 4. Comparison between simulations of the stochastic velocity-jump model and the fundamental solution of (6), equation (9) in two dimensions. Results shown at times t=50 (left) and t=200 (right). Red lines indicate contours of the average distribution for the particle position, generated through repeated simulations of the stochastic velocity-jump (VJ) model, while blue lines indicate contours of the fundamental solution (FS), equation (9); contours are shown here for c(x,t)=104 and 102. For the stochastic simulations we initialize the individual at (0,0) and compute the random walk path using the von Mises distribution for q (with ϕ=π/4 and κ=10), s=10 and λ=100. a and D for (9) are subsequently determined using (10) and (11).
    Figure 5. Comparison between simulations of the stochastic velocity-jump model and the fundamental solution of (6) in three dimensions. All plots are generated at time t=50. In (a)-(b) we plot the three dimensional contour surface at the fixed density c(x,t)=105 for (a) the fundamental solution given by equation (9) and (b) the particle probability distribution generated from multiple simulations (200000) of the 3D velocity-jump model. In (c-d) we illustrate the quantitative match by plotting contour lines of fixed density c(x,t)=104 (solid line) and c(x,t)=106 (dashed line) for projections onto the (c) xz plane (fixing y=50) and (d) yz plane (fixing x=460). For the simulations we initialise the individual at (0,0,0) and compute the random walk path using the Fisher distribution (2) for q (with u=(1,0,0)) and κ=10) along with s=10 and λ=100. a and D for equation (9) are subsequently determined by equations (25) and (26).

    5.2. Sea turtle navigation

    The ability of animals to navigate has fascinated scientists for centuries. Certain green sea turtles (Chelonia mydas) are particularly renowned for swimming halfway across the Atlantic Ocean to nest at tiny Ascension Island [25], a feat noted already by Charles Darwin [8]. While the precise navigating cues remain uncertain, a combination of factors from geomagnetic information to wind or water-borne odours have been proposed to help turtles on their journey. In [31] we formulated an individual-based model for turtle navigation, where the turning direction for a simulated turtle from its current location was biased towards the island via a two dimensional von Mises distribution. The subsequent population-level macroscopic model is given by equation (6), although augmented by an additional advection term to account for ocean current drift, where the first two moments from equation (1) provide the forms for a and D. We refer to [31] for full details of the modelling.

    Figure 6 shows the results from two simulation scenarios. All simulations consider a start date of 1st of January 2014 (early in a typical nesting season) and a population of 1000 turtles that are initially distributed uniformly throughout a circular region centred on the island. Simulated turtles subsequently attempt to swim towards the island, with some fixed concentration parameter k that reflects the "navigational strength" of the cue being utilised. For these simulations, the additional drift due to ocean currents is taken from (post-validated) ocean surface flow data provided by an ocean forecasting model [5] over the studied time period. Figure 6 shows simulations for populations of (a) "weak navigators" (k=0.5) and (b) "strong navigators" (k=2.0). For each case we show results from (top panel) the underlying individual-based model and (bottom panel) the macroscopic model (6) for the population distribution: the close match between these distributions demonstrates that the macroscopic model (6) can accurately capture the individual-level data and stresses the usefulness of the moment calculations. Our simulations show that weak-navigators are generally unsuccessful in their attempts to reach the island: the orientating bias is not sufficiently strong to overcome the often powerful ocean currents encountered. Strong navigators, however, are able to plot a relatively straightforward course to the island. More details, in particular on the estimation of model parameters, can be found in [31].

    Figure 6. Simulations of green sea turtle navigation towards Ascension island: we refer to [31] for full details. We consider (a) weak navigators and (b) strong navigators, showing in each case (top row) the individual-based model and (bottom row) the macroscopic model. For the IBM turtle positions (white circles) are shown at the indicated times, superimposed on ocean currents as illustrated by its direction (arrows) and magnitude (arrow length/density map), while for the macroscopic model (6) the population density is indicated by the density map, where the scale indicates the number/km2 and white regions represent a density below 105. Simulations start on 1st of January 2014 and turtles swim with mean fixed swimming speed of 40 km/day and make turns once per hour. Ocean currents are obtained from a standard ocean forecasting model (HYCOM, [5]).

    5.3. Wolf movement on seismic lines

    A good example in the usage of the bimodal von Mises distribution comes from wolf movement in disturbed forest environments. Boreal landscapes in Northern Alberta are characterized by a network of narrow straight corridors of cleared forest: so-called "seismic lines" created by oil companies to allow seismic imaging of potential oil reservoirs. Approximately 5 metres wide, these lines criss-cross the habitats of caribou and wolf populations and have been found to impact on their spatial distribution, thereby potentially altering the predator-prey dynamics. Location data from radio collars attached to wolves and caribou found that (on average) the former tended to be closer to seismic lines and the latter further away [18], suggesting that wolves may utilise the seismic lines to facilitate their movement through the landscape (and hence expand their hunting range) while caribou avoid them, perhaps due to a greater chance of encounters with wolves. More recently, GPS collar data from wolves demonstrated that they tended to move along the seismic lines, and that a bimodal von Mises distribution provided a good fit for their movement properties [23]. Subsequent modelling in [23], based on various modes of biased movement in response to seismic lines, was used to predict how this would impact on encounter/predation rates.

    To show how biased movement along seismic lines can alter the population distribution, we consider simulations of the macroscopic model (6) in different landscapes, where the expectation and variance-covariance required for a and D are those of the bimodal von Mises distribution (23-24). Specifically, we consider (a) an artificial landscape of regular criss-crossing seismic lines and (b) a genuine landscape based on an aerial photograph. We assume distinct movement properties on and off the lines as follows.

    • Off-seismic lines the population is assumed to move in a (more or less) random manner, i.e. k=0 in (4).

    • On-seismic lines the population is biased into either following or avoiding the seismic lines:

    – for seismic line following we take k=2 in (4) and a direction of movement along the (local) axial direction of the seismic line;

    – for seismic line avoidance we take k=2 in (4) and a dominant direction of movement orthogonal to the (local) axial direction of the seismic line.

    Figure 7 plots the long-term distributions for populations following one of these scenarios. As to be expected, there is a heterogeneous pattern in the population distribution according to the underlying network of seismic lines, such that an individual is either more or less likely to be found on (or close to) a seismic line for following or avoidance strategies respectively. A more formal matching of this distribution data to the location data found in studies such as [18,23] would allow us to examine theoretically how different strategies and responses of both wolves and caribou to the seismic line network could impact on predation rates: we refer to [23] to a specific such study in this regard.

    Figure 7. Population density distributions describing the different responses of a population to linear landscape features such as seismic lines. In each case we solve (6), where a and D are calculated using the moments of the bimodal von Mises distribution, i.e. equations (23-24). In (a) we impose a regular array of criss-crossing seismic lines: the population density distribution is either (a1) raised or (a2) decreased along the seismic lines according to whether they tend to (a1) follow along the paths or (a2) move directly off of them. In (b) an identical set of simulations is conducted for a genuine landscape, based on an aerial photograph of a Northern Alberta landscape.

    5.4. An application to brain tumour modelling

    One application of the three-dimensional von Mises-Fisher distribution can be found in the modelling of brain tumours. Gliomas are brain tumours that develop from the glial cells of the central nervous system. The most aggressive variant is the Grade Ⅳ Astrocytoma (Glioblastoma multiforme). Current standard of care includes a combination of maximal safe surgical de-bulking, radiation therapy, and chemotherapy. The typical median survival is approximately 14 months from diagnosis. Despite intensive treatment, recurrence is expected in virtually all patients, and usually occurs within and adjacent of the original region of grossly visible disease seen on diagnostic imaging. Microscopic extension of tumor cells can be found in the adjacent normal appearing brain and are thought to contribute to recurrence and subsequent extension of the tumour.

    Brain tissue can be classified into two main components: gray matter and white matter. While the former is relatively isotropic in structure, a significant proportion of white matter stems from the myelinated nerve axons, bundled into long fibre tracts that generate a highly anisotropic tissue. Significant evidence suggests that invasive glioma cells spread along these white matter tracts [33,11,12] hence directing the heterogeneous spread of cancer. Consequently it is of key interest to understand how these fibres are connected in the brain. An MRI-based technique called diffusion tensor imaging (DTI) generates voxel-by-voxel diffusion tensors for the anisotropic diffusion of water molecules, which can subsequently be reinterpreted to provide a measurement of tissue alignment and hence information on the arrangement of white matter tracts [4].

    A number of groups have developed methods to include DTI measurements within mathematical glioma growth models [35,13,19,21,26,30,9] and in [30] we proposed to incorporate DTI information via its parametrisation of a bimodal von Mises-Fisher distribution. Consequently the macroscopic glioma growth follows an anisotropic diffusion model, augmented by a logistic growth term to describe proliferation of the cancer cells:

    c(x,t)t=:(D(x)c(x,t))+rc(x,t)(1c(x,t))

    where c(x,t) is the time and space dependent (normalised) cell density and r is the growth rate. The colon symbol (:) is used as contraction as in (6). D(x) is given from (8) using (32), where the dominant direction u and concentration parameter k is determined by the anisotropy revealed in the DTI data and will be spatially variant. The concentration parameter k is taken to be proportional to the fractional anisotropy (FA) at each location, with constant of proportionality κ, i.e. k(x)=κFA(x). The fractional anisotropy, or FA measures "how anisotropic" the DTI measured tensor is, with FA=0 corresponding to an isotropic tensor and FA=1 corresponding to a fully anisotropic tensor. The constant κ measures the degree to which the invasion process follows the anisotropic structure and can be tuned for a specific patient: choosing κ=0 will generate a uniform distribution, and hence an isotropic tensor, while a higher value of κ simulates a stronger tendency of the cancer cells to invade anisotropically along white matter tracts.

    A full description of the model and simulation procedure can be found in [34] (see also [30]) and a typical simulation is shown in Figure 8. In (a), we plot the FA map in a two-dimensional slice of brain for a particular patient, and augment it with a black dot representing the initial condition (i.e. proposed point of tumour initiation) for a simulation; note that yellow regions correspond to high FA, and hence imply highly aligned tissue in the form of bundled white matter tracts. (b) and (d) show two-dimensional slices through the three-dimensional cell density generated by model simulation, creating an artificial tumour based on patient-specific DTI data. While the (b) and (c) corresponds to a low κ value (κ=0.5), (d) and (e) use a high κ value (κ=15). Here the blue to yellow colormap indicates increasing cell density and the tumour profile is plotted 280 days post initiation. (c) and (e) show an isosurface for the three-dimensional cell density corresponding to c(x,t)=0.16, estimated to be the density observable on a T2 MRI image [36] and hence indicating the visible tumour boundary. We see that for higher values of κ the finger-like projections often seen in these types of tumours can be recreated.

    Figure 8. Simulations of brain tumour growth using real patient data: see [34] for details. These are test cases to show the effect of changing the concentration parameter. (a) shows the fractional anisotropy for a two-dimensional axial slice of a real patient brain. Yellow indicates high fractional anisotropy, and thus high alignment, while blue indicates isotropic tissue. The initial condition for the brain tumour simulation is indicated by a black dot. (b)-(e) show two artificial tumours generated using real patient DTI data for two different values of κ, and thus k. (b) and (c) use κ=0.5, and (d) and (e) use κ=15. (b) and (d) represent a two-dimensional slice of the cell density (where dark blue = low cancer cell density, yellow = high cancer cell density), while (c) and (e) represent an isosurface corresponding to c(x,t)=0.16. It has been estimated that this is the cell-denstity that shows up on T2-MRI images, thus the isosurfaces in (c) and (e) can be thought of as the visible tumour boundaries.

    6. Conclusion

    We have shown that the divergence theorem can be a powerful tool to compute the moments of spherical distributions. In this paper we focussed on the first two moments of two and three dimensional spherical analogues of the Normal distribution: the von Mises distribution (in 2-D) and the Fisher distribution (in 3-D). The method presented avoids tedious trigonometric integrals and we expect that it can also be employed for higher order moments, for higher dimensions, or for other spherical distributions such as the Kent-or Bingham distributions [20,22]. The key argument was the use of the divergence theorem, and this is also possible in higher dimensions. In those cases, higher spherical moment calculations are necessary, but there is hope that these can be solved. See for example [17], where general spherical moments were computed in arbitrary space dimensions for the uniform distribution.

    Our focus on the first two moments of the von Mises and Fisher distributions stems from their direct relevance to the modelling of oriented movement of animals and cells, highlighted here by applications to the modelling of sea turtle navigation, wolf movement and brain tumour invasion: such movement-based problems naturally reside in two or three dimensional space and one could apply similar methods to model a myriad of other such phenomena in biology, ecology or medicine. More generally, the n-dimensional von Mises-Fisher distribution has found applications in numerous other areas, such as describing the type of high dimensional data found in areas such as gene expression data [2]. Consequently, we hope that the method described here may prove to have wider utility for other applications requiring descriptions of directional data.


    Acknowledgments

    AS is grateful to support through a NSERC CGS D3 scholarship and an Alberta Innovates Technology Futures Graduate Student Scholarship. TH is supported through an NSERC discovery grant. The DTI images were acquired through an Alberta Cancer Foundation sponsored grant aimed at using DTI imaging to help predict the pattern of glioma growth.


    [1] [ wikipedia. com, https://en.wikipedia.org/wiki/VonMises%E2%80%93Fisher_distribution, last accessed on 11/22/2016.
    [2] [ A. Banerjee,I. S. Dhillon,J. Ghosh,S. Sra, Clustering on the unit hypersphere using von Mises-Fisher distributions, Journal of Machine Learning Research, 6 (2005): 1345-1382.
    [3] [ E. Batschelet, null, Circular Statistics in Biology, Academic Press, London, 1981.
    [4] [ C. Beaulieu, The basis of anisotropic water diffusion in the nervous system -a technical review, NMR Biomed., 15 (2002): 435-455.
    [5] [ R. Bleck, An oceanic general circulation model framed in hybrid isopycnic-cartesian coordinates, Ocean Mod., 4 (2002): 55-88.
    [6] [ E. A. Codling,N. Bearon,G. J. Thorn, Thorn, Diffusion about the mean drift location in a biased random walk, Ecology, 91 (2010): 3106-3113.
    [7] [ E. A. Codling,M. J. Plank,S. Benhamou, Random walk models in biology, J. Roy. Soc. Interface, 5 (2008): 813-834.
    [8] [ C. Darwin, Perception in the lower animals, Nature, 7 (1873): 360.
    [9] [ C. Engwer,T. Hillen,M. Knappitsch,C. Surulescu, Glioma follow white matter tracts: A multiscale DTI-based model, J. Math. Biol., 71 (2015): 551-582.
    [10] [ N. I. Fisher, null, Statistical Analysis of Circular Data, Cambridge University Press, Cambridge, 1993.
    [11] [ A. Giese,L. Kluwe,B. Laube,H. Meissner,M. E. Berens,M. Westphal, Migration of human glioma cells on myelin, Neurosurgery, 38 (1996): 755-764.
    [12] [ P. G. Gritsenko,O. Ilina,P. Friedl, Interstitial guidance of cancer invasion, Journal or Pathology, 226 (2012): 185-199.
    [13] [ H. Hatzikirou,A. Deutsch,C. Schaller,M. Simaon,K. Swanson, Mathematical modelling of glioblastoma tumour development: A review, Math. Models Meth. Appl. Sci., 15 (2005): 1779-1794.
    [14] [ T. Hillen, M5 mesoscopic and macroscopic models for mesenchymal motion, J. Math. Biol., 53 (2006): 585-616.
    [15] [ T. Hillen, E. Leonard and H. van Roessel, Partial Differential Equations; Theory and Completely Solved Problems, Wiley, Hoboken, NJ, 2012.
    [16] [ T. Hillen and K. Painter, Transport and anisotropic diffusion models for movement in oriented habitats, In: Lewis, M., Maini, P., Petrovskii, S. (Eds.), Dispersal, Individual Movement and Spatial Ecology: A Mathematical Perspective, Springer, Heidelberg, 2071 (2013), 177-222.
    [17] [ T. Hillen, On the L2-moment closure of transport equations: The general case, Discr.Cont.Dyn.Systems, Series B, 5 (2005): 299-318.
    [18] [ A. R. C. James,A. K. Stuart-Smith, Distribution of caribou and wolves in relation to linear corridors, The Journal of Wildlife Management, 64 (2000): 154-159.
    [19] [ A. Jbabdi,E. Mandonnet,H. Duffau,L. Capelle,K. Swanson,M. Pelegrini-Issac,R. Guillevin,H. Benali, Simulation of anisotropic growth of low-grade gliomas using diffusion tensor imaging, Mang. Res. Med., 54 (2005): 616-624.
    [20] [ J. Kent, The Fisher-Bingham Distribution on the Sphere, J. Royal. Stat. Soc., 1982.
    [21] [ E. Konukoglu,O. Clatz,P. Bondiau,H. Delignette,N. Ayache, Extrapolation glioma invasion margin in brain magnetic resonance images: Suggesting new irradiation margins, Medical Image Analysis, 14 (2010): 111-125.
    [22] [ K. V. Mardia and P. E. Jupp, Directional Statistics, Wiley, New York, 2000.
    [23] [ H. McKenzie,E. Merrill,R. Spiteri,M. Lewis, How linear features alter predator movement and the functional response, Royal Society Interface, 2 (2012): 205-216.
    [24] [ P. Moorcroft,M. A. Lewis, null, Mechanistic Home Range Analysis, Princeton University Press, USA, 2006.
    [25] [ J. A. Mortimer,A. Carr, Reproduction and migrations of the Ascension Island green turtle (chelonia mydas), Copeia, 1987 (1987): 103-113.
    [26] [ P. Mosayebi,D. Cobzas,A. Murtha,M. Jagersand, Tumor invasion margin on the Riemannian space of brain fibers, Medical Image Analysis, 16 (2011): 361-373.
    [27] [ A. Okubo and S. Levin, Diffusion and Ecological Problems: Modern Perspectives, Second edition. Interdisciplinary Applied Mathematics, 14. Springer-Verlag, New York, 2001.
    [28] [ H. Othmer,S. Dunbar,W. Alt, Models of dispersal in biological systems, Journal of Mathematical Biology, 26 (1988): 263-298.
    [29] [ K. Painter, Modelling migration strategies in the extracellular matrix, J. Math. Biol., 58 (2009): 511-543.
    [30] [ K. Painter,T. Hillen, Mathematical modelling of glioma growth: The use of diffusion tensor imaging DTI data to predict the anisotropic pathways of cancer invasion, J. Theor. Biol., 323 (2013): 25-39.
    [31] [ K. Painter,T. Hillen, Navigating the flow: Individual and continuum models for homing in flowing environments, Royal Society Interface, 12 (2015): 20150,647.
    [32] [ K. J. Painter, Multiscale models for movement in oriented environments and their application to hilltopping in butterflies, Theor. Ecol., 7 (2014): 53-75.
    [33] [ J. Rao, Molecular mechanisms of glioma invasiveness: The role of proteases, Nature Reviews Cancer, 3 (2003): 489-501.
    [34] [ A. Swan, T. Hillen, J. Bowman and A. Murtha, A patient-specific anisotropic diffusion model for brain tumor spread, Bull. Math. Biol., to appear (2017).
    [35] [ K. Swanson,E. A. Jr,J. Murray, A quantitative model for differential motility of gliomas in grey and white matter, Cell Proliferation, 33 (2000): 317-329.
    [36] [ K. Swanson,R. Rostomily,E. Alvord Jr, Predicting survival of patients with glioblastoma by combining a mathematical model and pre-operative MR imaging characteristics: A proof of principle, British Journal of Cancer, 98 (2008): 113-119.
    [37] [ P. Turchin, Quantitative Analysis of Movement, Sinauer Assoc., Sunderland, 1998.
    [38] [ C. Xue,H. Othmer, Multiscale models of taxis-driven patterning in bacterial populations, SIAM Journal for Applied Mathematics, 70 (2009): 133-167.
  • This article has been cited by:

    1. Karol Salazar-Ariza, Rafael Torres, Statistical theory of the polarization on the Poincaré sphere, 2019, 44, 0146-9592, 3318, 10.1364/OL.44.003318
    2. Amanda Swan, Thomas Hillen, John C. Bowman, Albert D. Murtha, A Patient-Specific Anisotropic Diffusion Model for Brain Tumour Spread, 2018, 80, 0092-8240, 1259, 10.1007/s11538-017-0271-8
    3. Juan D. Olarte-Plata, Fernando Bresme, Theoretical description of the thermomolecular orientation of anisotropic colloids, 2019, 21, 1463-9076, 1131, 10.1039/C8CP06780E
    4. Karl K. Sabelfeld, Application of the von Mises–Fisher distribution to Random Walk on Spheres method for solving high-dimensional diffusion–advection–reaction equations, 2018, 138, 01677152, 137, 10.1016/j.spl.2018.03.002
    5. K.J. Painter, A.Z. Plochocka, Efficiency of island homing by sea turtles under multimodal navigating strategies, 2019, 391, 03043800, 40, 10.1016/j.ecolmodel.2018.10.025
    6. Daria Stepanova, Helen M. Byrne, Philip K. Maini, Tomás Alarcón, Roeland M.H. Merks, A multiscale model of complex endothelial cell dynamics in early angiogenesis, 2021, 17, 1553-7358, e1008055, 10.1371/journal.pcbi.1008055
    7. Ion Bica, Thomas Hillen, Kevin J. Painter, Aggregation of biological particles under radial directional guidance, 2017, 427, 00225193, 77, 10.1016/j.jtbi.2017.05.039
    8. Paul C. Bressloff, Samuel R. Carroll, Stochastic neural fields as gradient dynamical systems, 2019, 100, 2470-0045, 10.1103/PhysRevE.100.012402
    9. Kevin J. Painter, Thomas Hillen, 2018, Chapter 5, 978-3-319-96841-4, 103, 10.1007/978-3-319-96842-1_5
    10. Rangaprasad Arun Srivatsan, Mengyun Xu, Nicolas Zevallos, Howie Choset, Probabilistic pose estimation using a Bingham distribution-based linear filter, 2018, 37, 0278-3649, 1610, 10.1177/0278364918778353
    11. S.T. Johnston, K.J. Painter, The impact of short- and long-range perception on population movements, 2019, 460, 00225193, 227, 10.1016/j.jtbi.2018.10.031
    12. Yaping LI, Xiaowei CHEN, Yongbo SHAO, Numerical Modeling of Natural Fracture Distributions in Shale, 2023, 1000-9515, 10.1111/1755-6724.15050
    13. N. LOY, T. HILLEN, K. J. PAINTER, Direction-dependent turning leads to anisotropic diffusion and persistence, 2022, 33, 0956-7925, 729, 10.1017/S0956792521000206
    14. James D. McLaren, Heiko Schmaljohann, Bernd Blasius, Predicting performance of naïve migratory animals, from many wrongs to self-correction, 2022, 5, 2399-3642, 10.1038/s42003-022-03995-5
    15. Ryan Thiessen, Thomas Hillen, Anisotropic Network Patterns in Kinetic and Diffusive Chemotaxis Models, 2021, 9, 2227-7390, 1561, 10.3390/math9131561
    16. Martina Conte, Nadia Loy, Multi-Cue Kinetic Model with Non-Local Sensing for Cell Migration on a Fiber Network with Chemotaxis, 2022, 84, 0092-8240, 10.1007/s11538-021-00978-1
    17. Romain Fayat, Viviana Delgado Betancourt, Thibault Goyallon, Mathieu Petremann, Pauline Liaudet, Vincent Descossy, Lionel Reveret, Guillaume P. Dugué, Inertial Measurement of Head Tilt in Rodents: Principles and Applications to Vestibular Research, 2021, 21, 1424-8220, 6318, 10.3390/s21186318
    18. Yaping Li, Xiaowei Chen, Yongbo Shao, 3D natural fracture model of shale reservoir based on petrophysical characterization, 2023, 166, 01918141, 104763, 10.1016/j.jsg.2022.104763
    19. Szabolcs Suveges, Kismet Hossain-Ibrahim, J. Douglas Steele, Raluca Eftimie, Dumitru Trucu, Mathematical Modelling of Glioblastomas Invasion within the Brain: A 3D Multi-Scale Moving-Boundary Approach, 2021, 9, 2227-7390, 2214, 10.3390/math9182214
    20. Gwangbin Bae, Ignas Budvytis, Roberto Cipolla, 2021, Estimating and Exploiting the Aleatoric Uncertainty in Surface Normal Estimation, 978-1-6654-2812-5, 13117, 10.1109/ICCV48922.2021.01289
    21. Dominik Friml, Pavel Vaclavek, 2022, Bayesian Inference of Total Least-Squares With Known Precision, 978-1-6654-6761-2, 203, 10.1109/CDC51059.2022.9992409
    22. Vitaly V Ganusov, Viktor S Zenkov, Barun Majumder, Correlation between speed and turning naturally arises for sparsely sampled cell movements, 2023, 20, 1478-3967, 025001, 10.1088/1478-3975/acb18c
    23. Qingchao Li, Mohammed El-Hajjar, Ibrahim Hemadeh, Deepa Jagyasi, Arman Shojaeifard, Lajos Hanzo, Performance Analysis of Active RIS-Aided Systems in the Face of Imperfect CSI and Phase Shift Noise, 2023, 72, 0018-9545, 8140, 10.1109/TVT.2023.3239398
    24. Mamoru Ota, Yuji Tanaka, Yoshiya Kasahara, Propagation Analysis Method Considering Angular Spread for Random Electromagnetic Waves in Magnetized Cold Plasma, 2023, 58, 0048-6604, 10.1029/2023RS007673
    25. Qingchao Li, Mohammed El-Hajjar, Yanshi Sun, Ibrahim Hemadeh, Arman Shojaeifard, Yuanwei Liu, Lajos Hanzo, Achievable Rate Analysis of the STAR-RIS-Aided NOMA Uplink in the Face of Imperfect CSI and Hardware Impairments, 2023, 71, 0090-6778, 6100, 10.1109/TCOMM.2023.3287995
    26. Shengxi Li, Danilo Mandic, Von Mises–Fisher Elliptical Distribution, 2023, 34, 2162-237X, 11006, 10.1109/TNNLS.2022.3160519
    27. Hanyu WANG, Qiang SHEN, Zilong DENG, Xinyi CAO, Xiaokang Wang, Absolute pose estimation of UAV based on large-scale satellite image, 2023, 10009361, 10.1016/j.cja.2023.12.028
    28. Jisan Yang, Jie Jiang, Jian Li, Yan Ma, Lingfeng Tian, Guangjun Zhang, Uncertainty Estimation and Multi-FOV Data Fusion for Star Sensors Based on Directional Statistics, 2024, 73, 0018-9456, 1, 10.1109/TIM.2023.3341142
    29. Qingchao Li, Mohammed El-Hajjar, Ibrahim Hemadeh, Arman Shojaeifard, Lajos Hanzo, Low-Overhead Channel Estimation for RIS-Aided Multi-Cell Networks in the Presence of Phase Quantization Errors, 2024, 73, 0018-9545, 6626, 10.1109/TVT.2023.3339968
    30. Sourav Chandra, Rajeev Singh, Rakesh Kumar Singh, Statistical insights of polarization speckle via von Mises–Fisher distribution on the Poincaré sphere, 2024, 41, 1084-7529, 1287, 10.1364/JOSAA.519685
    31. Linzhou Zeng, Xuewen Liao, Zhangfeng Ma, Hao Jiang, Zhen Chen, UAV-to-UAV MIMO Systems Under Multimodal Nonisotropic Scattering: Geometrical Channel Modeling and Outage Performance Analysis, 2024, 11, 2327-4662, 26266, 10.1109/JIOT.2024.3395524
    32. Linzhou Zeng, Xuewen Liao, Zhangfeng Ma, Wenwu Xie, Hao Jiang, Zhen Chen, Envelope Level Crossing Rate and Average Fade Duration of Low-Altitude UAV-to-UAV Channels, 2024, 13, 2162-2337, 2220, 10.1109/LWC.2024.3407865
    33. Qingchao Li, Mohammed El-Hajjar, Yanshi Sun, Lajos Hanzo, Performance Analysis of Reconfigurable Holographic Surfaces in the Near-Field Scenario of Cell-Free Networks Under Hardware Impairments, 2024, 23, 1536-1276, 11972, 10.1109/TWC.2024.3386850
    34. Andrei Ciprian Macarie, Szabolcs Suveges, Mohamed Okasha, Kismet Hossain-Ibrahim, J. Douglas Steele, Dumitru Trucu, Post–operative glioblastoma cancer cell distribution in the peritumoural oedema, 2024, 14, 2234-943X, 10.3389/fonc.2024.1447010
    35. Alexandra Shyntar, Thomas Hillen, Mathematical modeling of microtube-driven regrowth of gliomas after local resection, 2024, 22, 1551-0018, 52, 10.3934/mbe.2025003
    36. Thomas Hillen, Maria R. D’Orsogna, Jacob C. Mantooth, Alan E. Lindsay, Mean First Passage Times for Transport Equations, 2025, 85, 0036-1399, 78, 10.1137/24M1647667
  • Reader Comments
  • © 2017 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(8461) PDF downloads(1054) Cited by(36)

Article outline

Figures and Tables

Figures(8)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog