Citation: Thomas M. H. Hope. The Bilingual Cognitive Advantage: No Smoke without Fire?[J]. AIMS Neuroscience, 2015, 2(2): 58-65. doi: 10.3934/Neuroscience.2015.2.58
[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 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
1Not to be confused with the statistical F-distribution.
We consider distributions
$ q({\bf n}) = \frac{1}{2\pi I_0(k)} e^{k{\bf n}\cdot{\bf u}}, $ | (1) |
for
$ q({\bf n}) = \frac{1}{2\pi I_0(k)} e^{k\cos(\alpha-\mu)}. $ |
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({\bf n}) = \frac{k}{4\pi \sinh(k)} e^{k{\bf n}\cdot{\bf u}},{\bf n} \in \mathbb{S}^2. $ | (2) |
The extension to higher dimensions forms the
$ q({\bf n}) = \frac{k^{n/2-1}}{(2\pi)^{n/2} I_{n/2-1}(k)} e^{k {\bf n} \cdot {\bf u}}\, {\bf n}\in \mathbb{S}^{n-1}. $ |
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
$ q({{\bf n}}) = \sum\limits_{i=1}^N \dfrac{c_i}{2\pi I_0(k_i)} e^{k_i {\bf n}\cdot{\bf u}_i}\,, $ | (3) |
where the coefficients
$ q_B({{\bf n}}) = \dfrac{1}{4\pi I_0(k)} \Bigl(e^{k {\bf n}\cdot{\bf u}} + e^{-k {\bf n}\cdot{\bf u}}\Bigr)\,, $ | (4) |
along with its three-dimensional bimodal Fisher version
$ q({\bf n}) = \frac{k}{8\pi \sinh(k)} \Bigl( e^{k{\bf n}\cdot{\bf u}} + e^{-k{\bf n}\cdot{\bf u}} \Bigr)\,, $ | (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
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 (
$ \int_{\mathbb{S}^{n-1}} q(t,{\bf x},{\bf n}) d{\bf n} = 1\,. $ |
To obtain an efficient evolution equation for the macroscopic population density distribution
$ c(t,{\bf x})_t +\nabla \cdot \Bigl({\bf a} (t,{\bf x}) c(t,{\bf x})\Bigr) = \nabla\nabla : \Bigl(D(t,{\bf x}) c(t,{\bf x})\Bigr)\,, \label{driftdiffusion} \\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;= \sum\limits_{i,j =1 }^n \frac{\partial}{\partial x_i}\frac{\partial}{\partial x_j} (D^{ij} c) \,, \nonumber $ | (6) |
where the colon symbol (:) is used to denote the contraction of two tensors. In the above
• the drift
$ {\bf a} (t,{\bf x}) = s {\bf E}[q] = s \int_{\mathbb{S}^{n-1}} {\bf n} q(t,{\bf x},{\bf n}) d{\bf n}\,; \label{exp} $ | (7) |
• the
$ D (t,{\bf x}) = \frac{s^2}{\lambda} \mbox{Var}[q] = \frac{s^2}{\lambda} \int_{\mathbb{S}^{n-1}}({\bf n}- {\bf a} (t,{\bf x})/s)({\bf n}- {\bf a} (t,{\bf x})/s)^T q(t,{\bf x},{\bf n}) d{\bf n}\,. \label{var} $ | (8) |
Note that the product
For the case in which
$ \label{fundamental} c({\bf x},t) = \frac{1}{\sqrt{\mbox{det}D ( 4\pi t)^n}}\exp\left(-\frac{1}{4t} ({\bf x}-{\bf a} t)^T D^{-1} ({\bf x}-{\bf a} t)\right). $ | (9) |
The above describes the probability density for finding a random walker (starting at location 0, time
Clearly, the dynamics predicted by (6) pivot about the directional distribution
Lemma 3.1. For a given unit vector
$ q({{\bf n}}) = \frac{1}{2\pi I_0(k)} e^{k{{\bf n}}\cdot{\bf u}}, {\bf n}\in\mathbb{S}^1. $ |
Then
$ {\bf E}[q] = \frac{I_1(k)}{I_0(k)} \; {\bf u}\,;\\ $ | (10) |
$ {\rm Var}[q] = \frac{1}{2}\left(1 -\frac{I_2(k)}{I_0(k)} \right){\mathbb{I}}_2 + \left(\frac{I_2(k)}{I_0(k)} -\left(\frac{I_1(k)}{I_0(k)}\right)^2 \right) {\bf u}{\bf u}^T\,. $ | (11) |
In the above,
Proof2. Before we compute the moments of
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.
$I_n(x):= (-i)^{-n}J_n (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
$ I_n(k) = \frac{1}{2\pi}\int_0^{2\pi} \cos(n\phi) e^{k\cos\phi} d\phi \,. $ | (12) |
Bessel functions satisfy a differential relation
$ \frac{d}{dx}\left(x^n J_n(x) \right) = x^n J_{n-1} (x) $ | (13) |
and, for
$ J_{n+1}(x) = \frac{2n}{x} J_n(x) - J_{n-1}(x). $ | (14) |
The total mass of the unimodal von Mises distribution (1) can be directly computed from (12), where we denote the angle between
$ \int_{\mathbb{S}^1} q({{\bf n}}) d{{\bf n}} = \frac{1}{2\pi I_0(k) } \int_0^{2\pi} e^{k\cos \phi}d\phi =1.$ |
For the expectation,
$ {\bf E}[q]=\int_{\mathbb{S}^1} {\bf n} q({\bf n}) d{\bf n}\,, $ |
we use a test vector
$ {\bf b} \cdot {\bf E}[q]={\bf b} \cdot \int_{\mathbb{S}^1} {\bf n} q({\bf n}) d{\bf n}, $ |
and
$ {\bf b} \cdot 2\pi I_0(k) E[q] = {\bf b} \cdot \int_{\mathbb{S}^1} {{\bf n}} e^{k{{\bf n}}\cdot {\bf u}} d{\bf n}\,,\\ = \int_{\mathbb{S}^1} {{\bf n}^i}({\bf b}_i e^{k{{\bf n}}\cdot {\bf u}}) d{\bf n}\,, \\ = \int_{\mathbb{B}_1(0)} \frac{\partial }{\partial {\bf v}^i} ({\bf b}_i e^{k{\bf v}\cdot{\bf u}})d{\bf v}\,,\\ = \int_{\mathbb{B}_1(0)} {\bf b}_i k {\bf u}_i e^{k{\bf v}\cdot{\bf u}} d{\bf v}\,,\\ = {\bf b} \cdot k{\bf u} \int_0^1 \int_0^{2\pi} e^{rk cos\phi} r dr d\phi\,,\\ = {\bf b} \cdot k {\bf u} \int_0^1 2\pi r I_0(rk) dr\,,\\ = {\bf b} \cdot 2\pi k {\bf u}\int_0^1 r I_0(rk) dr \,. $ |
In the first step we use the summation convention, which requires summation over repeated indices, i.e.
The last integral can be solved via (13) and
$ r I_0(rk) = \frac{irk J_0(irk)}{ik} = \frac{1}{ik}\frac{d}{dx} (x J_1(x)) |_{x=irk} = \frac{1}{ik} \frac{d}{dr} (r J_1(irk)).$ |
Then
$ \int_0^1 r I_0(rk) dr = \frac{1}{ik} J_1(ik) = \frac{1}{ik} i I_1(k) = \frac{I_1(k)}{k}, $ | (15) |
and due to the fact that
$ {\bf E}[q] = \frac{I_1(k)}{I_0(k)} \; {\bf u}, $ |
which is (10).
The variance-covariance matrix of
$ \mbox{Var}[q] = \int_{\mathbb{S}^1}({{\bf n}}-{\bf E}[q])({\bf n}-{\bf E}[q])^T q({{\bf n}}) d{{\bf n}} = \int_{\mathbb{S}^1}{{\bf n}} {{\bf n}}^T q({{\bf n}}) d{{\bf n}} - {\bf E}[q] {\bf E}[q] ^T.$ |
To find the second moment of
$ 2\pi I_0(k) {\bf a}\int_{\mathbb{S}^1} {{\bf n}}{{\bf n}}^T q({{\bf n}})d{{\bf n}} \; {\bf b} $ | (16) |
$ = \int_{\mathbb{S}^1} {\bf a}_i{{\bf n}}^i {\bf b}_j{{\bf n}}^j e^{k{{\bf n}}^l{\bf u}_l} d{{\bf n}}\nonumber,\\ = \int_{\mathbb{S}^1} {{\bf n}}^i({\bf a}_i {\bf b}_j{{\bf n}}^j e^{k{{\bf n}}^l{\bf u}_l}) d{{\bf n}}\nonumber,\\ = \int_{\mathbb{B}_1(0)} \frac{\partial}{\partial {\bf v}^i} ({\bf a}_i {\bf b}_j{\bf v}^j e^{k{\bf v}^l{\bf u}_l}) d{\bf v}\nonumber,\\ = \int_{\mathbb{B}_1(0)} {\bf a}_i {\bf b}_i e^{k {\bf v}\cdot{\bf u}} d{\bf v} + \int_{\mathbb{B}_1(0)} {\bf a}_i({\bf v}\cdot {\bf b}) k {\bf u}_i e^{k{\bf v}\cdot {\bf u}} d{\bf v}\nonumber, \\ = {\bf a}\cdot {\bf b} \int_{\mathbb{B}_1(0)} e^{k{\bf v}\cdot{\bf u}} d{\bf v} + k {\bf a}\cdot{\bf u} \; {\bf b}\cdot \int_{\mathbb{B}_1(0)} {\bf v} e^{k{\bf v}\cdot{\bf u}} d{\bf v}, $ | (17) |
where, again, we use the summation convention and the divergence theorem. The first integral in (17) can be solved directly:
$ \int_{\mathbb{B}_1(0)}e^{k{\bf v}\cdot{\bf u}} d{\bf v} = \int_0^1\int_{\mathbb{S}^1} e^{rk{{\bf n}}\cdot{\bf u}} r dr d{{\bf n}} = \int_0^1 2 \pi r I_0(rk) dr = 2\pi \frac{I_1(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:
$ \int_{\mathbb{B}_1(0)} {\bf v} e^{k{\bf v}\cdot {\bf u}} d{\bf v} = \int_0^1\int_{\mathbb{S}^1} r{{\bf n}} e^{r k {{\bf n}}\cdot{\bf u}}r dr d{{\bf n}} = \int_0^1 r^2 \int_{\mathbb{S}^1} {{\bf n}} e^{rk{{\bf n}}\cdot{\bf u}} d{{\bf n}} dr \nonumber, \\ = 2\pi{\bf u} \int_0^1 r^2 I_1(rk) dr,\label{3rdBessel} $ | (18) |
where we used (10) in the last step.
Now we use the differential recurrence relation (13) to write
$ r^2 I_1(rk) = -\frac{1}{ik^2}(irk)^2 J_1(irk) = -\frac{1}{ik^2} \frac{d}{dx} (x^2 J_1(x))|_{x=irk} = -\frac{1}{k} \frac{d}{dr}(r^2 J_1(irk)).$ |
Continuing from (18) we find
$ \int_{\mathbb{B}_1(0)} {\bf v} e^{k{\bf v}\cdot {\bf u}} d{\bf v} = -2\pi{\bf u} \int_0^1 \frac{1}{k}\frac{d}{dr}(r^2 J_1(irk)) dr = -2\pi{\bf u} J_2(ik) = 2\pi {\bf u} \frac{I_2(k)}{k}. $ | (19) |
Substituting all the integrals back into equation (17):
$ {\bf a}\int_{\mathbb{S}^1} {{\bf n}}{{\bf n}}^T q({{\bf n}})d{{\bf n}}\; {\bf b} = {\bf a}\cdot {\bf b} \frac{2\pi \frac{I_1(k)}{k}}{2\pi I_0(k)} + k {\bf a}\cdot {\bf u} \frac{2\pi {\bf u}\cdot {\bf b} \frac{I_2(k)}{k}}{2\pi I_0(k)}\nonumber, \\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;= {\bf a}\left( \frac{1}{k} \frac{I_1(k)}{I_0(k)} {\mathbb{I}}_2 +{\bf u}{\bf u}^T \frac{I_2(k)}{I_0(k)} \right) {\bf b}\,.\nonumber $ |
Finally, we use the identity (14) for
$ \frac{1}{k}\frac{I_1(k)}{I_0(k)} = \frac{1}{2}\left(1-\frac{I_2(k)}{I_0(k)}\right),$ |
and the second moment is given by
$ \int_{\mathbb{S}^1} {{\bf n}}{{\bf n}}^T q({{\bf n}})d{{\bf n}} = \frac{1}{2}{\mathbb{I}}_2 + \frac{I_2(k)}{I_0(k)}\left({\bf u}{\bf u}^T-\frac{1}{2}{\mathbb{I}}_2\right). $ | (20) |
Together with the formula for the expectation (10) we find
$ \mbox{Var}[q] = \int_{\mathbb{S}^1}{{\bf n}} {{\bf n}}^T q({{\bf n}}) d{{\bf n}} - {\bf E}[q] {\bf E}[q] ^T \,,\\ \;\;\;\;\;\;\;\;\;\; = \frac{1}{2}{\mathbb{I}}_2 + \frac{I_2(k)}{I_0(k)}\left({\bf u}{\bf u}^T-\frac{1}{2}{\mathbb{I}}_2\right)- \left(\frac{I_1(k)}{I_0(k)}\right)^2 {\bf u}{\bf u}^T \,, \\ \;\;\;\;\;\;\;\;\;\;= \frac{1}{2}\left(1 -\frac{I_2(k)}{I_0(k)} \right){\mathbb{I}}_2 + \left(\frac{I_2(k)}{I_0(k)} -\left(\frac{I_1(k)}{I_0(k)}\right)^2 \right) {\bf u}{\bf u}^T\,, $ |
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
$
{\bf E}[q] = \frac{I_1(k)}{I_0(k)} \left(cosμsinμ \right);
$
|
(21) |
$
\mbox{ Var} [q]=
\frac{1}{2}\left(1 -\frac{I_2(k)}{I_0(k)} \right)\left(1001 \right) \nonumber \\
+ \left(\frac{I_2(k)}{I_0(k)} -\left(\frac{I_1(k)}{I_0(k)}\right)^2 \right) \left(cos2μcosμsinμcosμsinμsin2μ \right)\,. \label{Besselvariancepolar}
$
|
(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
$q({{\bf n}}) = \frac{1}{4\pi I_0(k)}\left(e^{k{{\bf n}}\cdot{\bf u}} + e^{-k{{\bf n}}\cdot{\bf u}}\right),$ |
then
$ {\bf E}[q] = 0\,, $ | (23) |
$ \mbox{Var}[q] = \frac{1}{2}\left(1- \frac{I_2(k)}{I_0(k)}\right) {\mathbb{I}}_2 + \frac{I_2(k)}{I_0(k)} {\bf u}{\bf u}^T. $ | (24) |
Proof. Since the bimodal von Mises distribution is symmetric we have
It is interesting to consider the limiting scenario as the parameter of concentration
The main result of this paper is as follows.
Theorem 4.1. For given
$ q({{\bf n}}) = \frac{k}{4\pi\sinh k} e^{k{\bf u}\cdot {{\bf n}}}, {{\bf n}} \in \mathbb{S}^2.$ |
Then
$ {\bf E}[q] = \left(\coth k - \frac{1}{k}\right) \; {\bf u}\,,\label{Exp3}\\ $ | (25) |
$ Var[q] = \left(\frac{\coth k }{k} - \frac{1}{k^2} \right){\mathbb{I}} + \left(1-\frac{ \coth k}{k} + \frac{2}{k^2}-\coth^2 k\right)\;{\bf u}{\bf u}^T. $ | (26) |
Proof. To compute the mean direction
$ {\bf E}[q] = \int_{\mathbb{S}^2} {{\bf n}} q({{\bf n}}) d{{\bf n}}, $ |
we use a test-vector
${\bf b}\cdot {\bf E}[q] = {\bf b}\cdot \int_{\mathbb{S}^2} {{\bf n}} q({{\bf n}}) d{{\bf n}}\,,\\ = \int_{\mathbb{S}^2} {{\bf n}}^i({\bf b}_i q({{\bf n}}))d{{\bf n}}\,, \\ = \int_{\mathbb{B}_1(0)} \frac{\partial}{\partial {\bf v}^i}\left({\bf b}_i \frac{k}{4\pi \sinh k } e^{k{\bf u}_l {\bf v}^l }\right) d{\bf v}\,, \\ = {\bf b}^i \int_{\mathbb{B}_1(0)} \frac{k^2{\bf u}_i}{4\pi\sinh k} e^{k {\bf u}\cdot {\bf v}} d{\bf v}\,, \\ = {\bf b}\cdot{\bf u} k \int_{\mathbb{B}_1(0)} \frac{k}{4\pi\sinh k} e^{k{\bf u}\cdot {\bf v}} d{\bf v}\,,\\ = {\bf b}\cdot{\bf u} k \int_0^1 \int_{\mathbb{S}^2} \frac{k}{4\pi \sinh k} e^{rk{\bf u}\cdot {\bf n}} r^2 dr d{{\bf n}}\,, \\ = {\bf b}\cdot{\bf u} k \int_0^1 \frac{r^2 k}{4\pi \sinh k}\frac{4\pi \sinh rk}{rk} dr\,,\\ = \frac{{\bf b}\cdot{\bf u} k}{\sinh k} \left(-\frac{1}{k^2}\sinh k + \frac{1}{k} \cosh k\right)\,,\\ = {\bf b}\cdot{\bf u} \left(\coth k - \frac{1}{k}\right), $ |
where in the 5th step we used the normalization condition for a Fisher distribution with concentration parameter
$ \int_{\mathbb{S}^2} e^{rk{\bf u}\cdot{\bf n}} d{\bf n} = \frac{4\pi\sin(rk)}{rk}.$ |
Since the test vector
To compute the second moment, we take two test-vectors
$ {\bf a}\cdot\int_{\mathbb{S}^2} {{\bf n}}{{\bf n}}^T q({{\bf n}}) d{{\bf n}} \cdot {\bf b} \nonumber \\ = \int_{\mathbb{S}^2} {\bf a}_i{{\bf n}}^i{{\bf n}}^j {\bf b}_j q({{\bf n}}) d{{\bf n}}\,, \nonumber\\ = \int_{\mathbb{S}^2} {{\bf n}}^i({\bf a}_i {{\bf n}}^j {\bf b}_j q({{\bf n}})) d{{\bf n}}\,,\nonumber\\ = \int_{\mathbb{B}_1(0)} \frac{\partial}{\partial {\bf v}^i}\left( {\bf a}_i {\bf v}^j {\bf b}_j \frac{k}{4\pi \sinh k} e^{k{\bf u}_l {\bf v}^l}\right) d{\bf v}\,,\nonumber \\ = \int_{\mathbb{B}_1(0)} {\bf a}_i {\bf b}_i \frac{k}{4\pi\sinh k} e^{k{\bf u}\cdot {\bf v}} d{\bf v} + \int_{\mathbb{B}_1(0)} {\bf a}_i ({\bf v}\cdot {\bf b}) \frac{k^2{\bf u}^i}{4\pi \sinh k} e^{k{\bf u}\cdot {\bf v}} d{\bf v}\,,\nonumber \\ = \underbrace{{\bf a}\cdot {\bf b} \int_{\mathbb{B}_1(0)} \frac{k}{4\pi \sinh k} e^{k{\bf u}\cdot {\bf v}} d{\bf v}}_{(A)} +\underbrace{{\bf a}\cdot {\bf u} \; {\bf b}\cdot \int_{\mathbb{B}_1(0)} \frac{k^2}{4\pi \sinh k} {\bf v} e^{k{\bf u}\cdot {\bf v}} d{\bf v}}_{(B)}\label{4.1}). $ | (27) |
From the computation of the first moment, we find for the first integral (A) that
$ \int_{\mathbb{B}_1(0)}\frac{k}{4\pi\sinh k} e^{k{\bf u}\cdot {\bf v}} d{\bf v} = \frac{1}{k}\coth k -\frac{1}{k^2}. $ | (28) |
The second integral (B) can be computed as
$ \int_{\mathbb{B}_1(0)} {\bf v} e^{k {\bf u}\cdot {\bf v}} d{\bf v} = \int_0^1 \int_{\mathbb{S}^2} r^3 {{\bf n}} e^{rk{\bf u}\cdot {\bf n}} d{{\bf n}} dr.$ |
From the first moment we obtain
$ \int_{\mathbb{S}^2} {{\bf n}} e^{rk{\bf u}\cdot {{\bf n}}}d{{\bf n}} = \left[ \frac{4\pi \sinh rk}{rk}\left(\coth rk -\frac{1}{rk}\right) \right] {\bf u}.$ |
Hence
$ \int_{\mathbb{B}_1(0)} {\bf v} e^{k{\bf u}\cdot {\bf v}} d{\bf v}\,, = \left[ \int_0^1 r^3\frac{4\pi\sinh rk}{rk}\left(\coth rk - \frac{1}{rk}\right)dr \right] {\bf u}\,,\\ = 4\pi\left[\int_0^1\frac{ r^2 \cosh rk}{k} dr - \int_0^1\frac{ r \sinh rk}{k^2}dr\right] {\bf u}\,,\\ = 4\pi\left[ \frac{\sinh k}{k^2}-\frac{3 \cosh k}{k^3} + \frac{3\sinh k}{k^4} \right] {\bf u}. $ |
Then
$ \frac{k^2}{4\pi \sinh k} \int_{\mathbb{B}_1(0)} {\bf v} e^{k{\bf u}\cdot {\bf v}}d{\bf v} = \left(1 - \frac{3}{k}\coth k+ \frac{3}{k^2}\right){\bf u}. $ | (29) |
Combining (28) and (29) with (27) we obtain
$ {\bf a}\cdot\int{{\bf n}}{{\bf n}}^Tq({{\bf n}})d{{\bf n}}\cdot {\bf b} = {\bf a}\cdot {\bf b}\left(\frac{\coth k }{k} - \frac{1}{k^2} \right) + {\bf a}\cdot {\bf u}\; {\bf u}\cdot {\bf b}\left(1-\frac{3 \coth k}{k} + \frac{3}{k^2}\right),$ |
and the second moment is given by
$ \int{{\bf n}} {{\bf n}}^T q({{\bf n}}) d{{\bf n}} = \left(\frac{\coth k }{k} - \frac{1}{k^2} \right){\mathbb{I}} + \left(1-\frac{3 \coth k}{k} + \frac{3}{k^2}\right)\;{\bf u}{\bf u}^T. $ | (30) |
To compute the full variance-covariance matrix of
$ \mbox{Var}[q] = \left(\frac{\coth k }{k} - \frac{1}{k^2} \right){\mathbb{I}} + \left(1-\frac{ \coth k}{k} + \frac{2}{k^2}-\coth^2 k\right)\;{\bf u}{\bf u}^T, $ |
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({{\bf n}}) = \frac{k}{8\pi\sinh k}\left( e^{k{\bf u}\cdot {{\bf n}}} + e^{-k{\bf u}\cdot{{\bf n}}}\right) = \frac{k}{4\pi\sinh k}\cosh(k{\bf u}\cdot {{\bf n}})$ |
are given by
$ {\bf E}[q] = 0, $ | (31) |
$ \mbox{Var}[q] = \left(\frac{\coth k }{k} - \frac{1}{k^2} \right){\mathbb{I}} + \left(1-\frac{3 \coth k}{k} + \frac{3}{k^2}\right)\;{\bf u}{\bf u}^T. $ | (32) |
Proof. As for the bimodal von Mises distribution, we exploit the symmetry properties by considering the components
It can be revealing to consider the coefficients arising in (25, 26, 32) and their dependence on the concentration parameter
$
a(k)=cothk−1k;b(k)=cothkk−1k2;c(k)=1−cothkk+2k2−coth2k;d(k)=1−3cothkk+3k2.
$
|
(33) |
Then for the unimodal case
$ {\bf E}[q] = a(k) {\bf u}, \mbox{Var}[q] = b(k) {\mathbb{I}} + c(k) {\bf u} {\bf u}^T, $ | (34) |
and for the bimodal case
$ {\bf E}[q] = 0, \mbox{Var}[q] = b(k) {\mathbb{I}} + d(k) {\bf u} {\bf u}^T. $ | (35) |
Thus, the variance-covariance matrix is split into an isotropic component, determined by
$\lim\limits_{k \to 0^+} a(k) = \lim\limits_{k \to 0^+} \coth k -\frac{1}{k}\,,$ |
which can be rewritten using l'Hȏpital's rule as
$limk→0+a(k)=limk→0+k−tanhkktanhk,=limk→0+1−sech2ktanhk+ksech2k,=limk→0+2sech2ktanhksech2k+(1−2sechk)sech2k,=0. $
|
Continuing in a similar fashion for
$\lim\limits_{k \to 0^+} b(k) = \lim\limits_{k \to 0^+} \frac{\coth k}{k} - \frac{1}{k^2}\,.$ |
This can be rewritten with repeated use of l'Hȏpital's rule as
$limk→0+b(k)=limk→0+k−tanhkk2tanhk,=limk→0+1−sech2k2ktanhk+k2sech2k,=limk→0+2sech2ktanhk2ksech2k+2tanhk−2k2sech2ktanhk+2ksech2k,=limk→0+−4sech2ktanh2k+2sech4k6sech2k−12ksech2ktanhk−2k2(sech4k−2sech2ktanh2k),=13. $
|
Then we have
$\lim\limits_{k \to 0^+} c(k) = \lim\limits_{k \to 0^+} 1-\frac{\coth k}{k}+\frac{2}{k^2}-\coth^2 k\,, $ |
which we may write as
$\lim\limits_{k \to 0^+} 1-\underbrace{\left( \frac{\coth k}{k}-\frac{1}{k^2}\right)}_{b(k)}+\frac{1}{k^2}-\coth^2 k = \lim\limits_{k \to 0^+} \frac{2}{3} +\left( \frac{1}{k}+\coth k \right) \left( \frac{1}{k}-\coth k\right)\,.$ |
In the above we have used the fact that the limit of
$limk→0+c(k)=limk→0+23+(1k+cothk)(1k−cothk),=limk→0+23−(1+kcothk)(b(k)),=limk→0+23−13(1+ksinhk),=23−23=0, $
|
where we used l'Hȏpital's rule in the last step. Finally, for
$limk→0+d(k)=limk→0+1−3cothkk+3k2,=limk→0+1−3(cothkk−1k2),=1−3limk→0+b(k),=1−3(13)=0. $
|
Summarising, for either the unimodal or bimodal Fisher distribution, as
$ {\bf a} = {\bf{0}} \;\;\;\mbox{and}\;\;\; D=\tau s^2 \frac{1}{3} \mathbb{I}_3 $ |
as would be expected.
These results are demonstrated numerically by plotting the coefficients as a function of
As a final note of interest, it is remarked that for the unimodal case (34) the coefficient
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.
To illustrate anisotropic diffusion we consider the movement of individuals released at
We initially consider a 2D scenario, where
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
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
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
• Off-seismic lines the population is assumed to move in a (more or less) random manner, i.e.
• On-seismic lines the population is biased into either following or avoiding the seismic lines:
– for seismic line following we take
– for seismic line avoidance we take
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.
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({\bf x},t)_t=\nabla \nabla : (D({\bf x}) c({\bf x},t))+rc({\bf x},t)(1-c({\bf x},t)) $ |
where
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
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
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] |
Paap KR, Sawi OM, Dalibar C, et al. (2014) The Brain Mechanisms Underlying the Cognitive Benefits of Bilingualism may be Extraordinarily Difficult to Discover. AIMS neurosci 1: 245-256. doi: 10.3934/Neuroscience.2014.3.245
![]() |
[2] | Saidi LG, Ansaldo AI (2015) Can a Second Language Help You in More Ways Than One? AIMS neurosci 1: 52-57. |
[3] |
Duyck W, Diependaele K, Drieghe D, et al. (2004) The Size of the Cross-Lingua Masked Phonological Priming Effect Does Not Depend on Second Language Proficiency. Exp Psychol 51: 116-124. doi: 10.1027/1618-3169.51.2.116
![]() |
[4] |
Lalor E, Kirsner K (2000) Cross-lingual transfer effects between English and Italian cognates and noncognates. Int J Biling 4: 385-398. doi: 10.1177/13670069000040030501
![]() |
[5] |
Lemhöfer K, Dijkstra T, Michel M (2004) Three languages, one ECHO: Cognate effects in trilingual word recognition. Lang Cogn Process 19: 585-611. doi: 10.1080/01690960444000007
![]() |
[6] | Kroll JF (2008) Juggling two languages in one mind. Psychological Science Agenda, American Psychological Association 22. |
[7] | Paap KR, Johnson HA, Sawi O (2014) Are bilingual advantages dependent upon specific tasks or specific bilingual experiences? J CognPsychol 26: 615-639. |
[8] |
Green DW (2008) Bilingual Aphasia: Adapted Language Networks and their Control. Annu Rev Appl Linguist 28: 25-48. doi: 10.1017/S0267190508080057
![]() |
[9] |
Green DW, Abutalebi J (2013) Language control in bilinguals: The adaptive control hypothesis. J Cogn Psychol 25: 515-530. doi: 10.1080/20445911.2013.796377
![]() |
[10] | Bialystok E, Majumder S, Martin MM (2003) Developing phonological awareness: Is there a bilingual advantage? Appl Psycholinguistics 24: 27-44. |
[11] |
Bialystok E, Martin MM (2004) Attention and inhibition in bilingual children: Evidence from the dimensional change card sort task. Developmental science 7: 325-339. doi: 10.1111/j.1467-7687.2004.00351.x
![]() |
[12] |
Bialystok E, Craik FIM, Freedman M (2007) Bilingualism as a protection against the onset of symptoms of dementia. Neuropsychologia 45: 459-464. doi: 10.1016/j.neuropsychologia.2006.10.009
![]() |
[13] |
Fernandes MA, Craik F, Bialystok E, et al. (2007) Effects of bilingualism, aging, and semantic relatedness on memory under divided attention. Can J Exp Psychol 61: 128-141. doi: 10.1037/cjep2007014
![]() |
[14] |
Bialystok E, Craik F, Luk G (2008) Cognitive control and lexical access in younger and older bilinguals. J Exp Psychol Learn Mem Cogn 34: 859-873. doi: 10.1037/0278-7393.34.4.859
![]() |
[15] |
Sterne JA, Davey Smith G (2001) Sifting the evidence-what's wrong with significance tests? Bmj 322: 226-231. doi: 10.1136/bmj.322.7280.226
![]() |
[16] |
Button KS, Ioannidis JPA, Mokrysz C, et al. (2013) Power failure: why small sample size undermines the reliability of neuroscience. Nat Rev Neurosci 14: 365-376. doi: 10.1038/nrn3475
![]() |
[17] | Hilchey MD, Klein RM (2011) Are there bilingual advantages on nonlinguistic interference tasks? Implications for the plasticity of executive control processes. Psychon Bull Rev 18: 625-658. |
[18] | Anton E, Dunabeitia JA, Estevez A, et al. (2014) Is there a bilingual advantage in the ANT task? Evidence from children. Front Psychol 5: 398. |
[19] |
Dunabeitia JA, Hernandez JA, Anton E, et al. (2014) The inhibitory advantage in bilingual children revisited: myth or reality? Exp Psychol 61: 234-251. doi: 10.1027/1618-3169/a000243
![]() |
[20] |
Paap KR, Greenberg ZI (2013) There is no coherent evidence for a bilingual advantage in executive processing. Cogn Psychol 66: 232-258. doi: 10.1016/j.cogpsych.2012.12.002
![]() |
[21] | Paap KR, Sawi O (2014) Bilingual advantages in executive functioning: problems in convergent validity, discriminant validity, and the identification of the theoretical constructs. Front Psychol 5: 962. |
[22] | Morton JB, Harper SN (2007) What did Simon say? Revisiting the bilingual advantage. Dev Sci 10: 719-726. |
[23] |
Kirk NW, Fiala L, Scott-Brown KC, et al. (2014) No evidence for reduced Simon cost in elderly bilinguals and bidialectals. J Cogn Psychol (Hove) 26: 640-648. doi: 10.1080/20445911.2014.929580
![]() |
[24] |
Calvo A, Bialystok E (2014) Independent effects of bilingualism and socioeconomic status on language ability and executive functioning. Cognition 130: 278-288. doi: 10.1016/j.cognition.2013.11.015
![]() |
[25] |
Ljungberg JK, Hansson P, Andrés P, et al. (2013) A Longitudinal Study of Memory Advantages in Bilinguals. PLoS ONE 8: e73029. doi: 10.1371/journal.pone.0073029
![]() |
[26] |
Sheng L, Lu Y, Gollan TH (2014) Assessing language dominance in Mandarin-English bilinguals: Convergence and divergence between subjective and objective measures. Biling (Camb Engl) 17: 364-383. doi: 10.1017/S1366728913000424
![]() |
[27] |
Ratcliff R, McKoon G (2008) The Diffusion Decision Model: Theory and Data for Two-Choice Decision Tasks. Neural comput 20: 873-922. doi: 10.1162/neco.2008.12-06-420
![]() |
[28] | Luk G, Anderson JAE, Craik FIM, et al. (2011) Distinct neural correlates for two types of inhibition in bilinguals: Response inhibition versus interference suppression. Brain Cogn 74: 347-357. |
[29] |
Hunt LT, Kolling N, Soltani A, et al. (2012) Mechanisms underlying cortical activity during value-guided choice. Nat Neurosci 15: 470-476. doi: 10.1038/nn.3017
![]() |
[30] |
Abutalebi J, Della Rosa PA, Green DW, et al. (2012) Bilingualism tunes the anterior cingulate cortex for conflict monitoring. Cereb Cortex 22: 2076-2086. doi: 10.1093/cercor/bhr287
![]() |
[31] |
Shenhav A, Botvinick MM, Cohen JD (2013) The expected value of control: an integrative theory of anterior cingulate cortex function. Neuron 79: 217-240. doi: 10.1016/j.neuron.2013.07.007
![]() |
[32] |
di Pellegrino G, Ciaramelli E, Ladavas E (2007) The regulation of cognitive control following rostral anterior cingulate cortex lesion in humans. J Cogn Neurosci 19: 275-286. doi: 10.1162/jocn.2007.19.2.275
![]() |
[33] |
Glascher J, Adolphs R, Damasio H, et al. (2012) Lesion mapping of cognitive control and value-based decision making in the prefrontal cortex. Proc Natl Acad Sci U S A 109: 14681-14686. doi: 10.1073/pnas.1206608109
![]() |
[34] |
Fellows LK, Farah MJ (2005) Is anterior cingulate cortex necessary for cognitive control? Brain 128: 788-796. doi: 10.1093/brain/awh405
![]() |
[35] |
Abutalebi J, Della Rosa PA, Green DW, et al. (2012) Bilingualism Tunes the Anterior Cingulate Cortex for Conflict Monitoring. Cereb Cortex 22: 2076-2086. doi: 10.1093/cercor/bhr287
![]() |
[36] |
Mechelli A, Crinion JT, Noppeney U, et al. (2004) Neurolinguistics: Structural plasticity in the bilingual brain. Nature 431: 757-757. doi: 10.1038/431757a
![]() |
[37] |
Mårtensson J, Eriksson J, Bodammer NC, et al. (2012) Growth of language-related brain areas after foreign language learning. NeuroImage 63: 240-244. doi: 10.1016/j.neuroimage.2012.06.043
![]() |
[38] |
Kovelman I, Baker SA, Petitto L-A (2008) Bilingual and Monolingual Brains Compared: A Functional Magnetic Resonance Imaging Investigation of Syntactic Processing and a Possible “Neural Signature” of Bilingualism. J cogn neurosci20: 153-169. doi: 10.1162/jocn.2008.20011
![]() |
[39] | Olsen RK, Pangelinan MM, Bogulski C, et al. (2015) The effect of lifelong bilingualism on regional grey and white matter volume. Brain Res. |
[40] |
Moreno S, Bialystok E, Wodniecka Z, et al. (2010) Conflict Resolution in Sentence Processing by Bilinguals. J Neurolinguistics 23: 564-579. doi: 10.1016/j.jneuroling.2010.05.002
![]() |
[41] | Calabria M, Hernandez M, Branzi FM, et al. (2012) Qualitative differences between bilingual language control and executive control: evidence from task switching. Front Psychol 2: 399. |
[42] |
Sagi D (2011) Perceptual learning in Vision Research. Vision Research 51: 1552-1566. doi: 10.1016/j.visres.2010.10.019
![]() |
[43] | Webster J, Whitworth A, Morris J (2015) Is it time to stop “fishing”? A review of generalisation following aphasia intervention. Aphasiology: 1-25. |
[44] |
Welch DB, Seitz AR (2013) Processing Irrelevant Location Information: Practice and Transfer Effects in a Simon Task. PLoS ONE 8: e64993. doi: 10.1371/journal.pone.0064993
![]() |
[45] |
Macnamara BN, Conway AR (2014) Novel evidence in support of the bilingual advantage: influences of task demands and experience on cognitive control and working memory. Psychon Bull Rev 21: 520-525. doi: 10.3758/s13423-013-0524-y
![]() |
[46] |
Prior A, Gollan TH (2013) The elusive link between language control and executive control: A case of limited transfer. J Cogn Psychol (Hove) 25: 622-645. doi: 10.1080/20445911.2013.821993
![]() |
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 |