Bivariate von Mises distribution

{{Short description|Probability distribution on a torus}}

File:Bivariate von Mises distribution cosine samples.svg

In probability theory and statistics, the bivariate von Mises distribution is a probability distribution describing values on a torus. It may be thought of as an analogue on the torus of the bivariate normal distribution. The distribution belongs to the field of directional statistics. The general bivariate von Mises distribution was first proposed by Kanti Mardia in 1975.{{Cite journal | jstor = 2984782|last= Mardia|first=Kanti |date=1975|title=Statistics of directional data|journal=J. R. Stat. Soc. B |volume=37| issue = 3|pages= 349–393|doi= 10.1111/j.2517-6161.1975.tb01550.x}}{{Cite book | doi = 10.1007/978-3-642-27225-7_6| chapter = Statistics of Bivariate von Mises Distributions| title = Bayesian Methods in Structural Bioinformatics| url = https://archive.org/details/bayesianmethodss00hame| url-access = limited| pages = [https://archive.org/details/bayesianmethodss00hame/page/n178 159]| series = Statistics for Biology and Health| year = 2012| last1 = Mardia | first1 = K. V. | last2 = Frellsen | first2 = J. | isbn = 978-3-642-27224-0}} One of its variants is today used in the field of bioinformatics to formulate a probabilistic model of protein structure in atomic detail, {{Cite journal | doi = 10.1073/pnas.0801715105| title = A generative, probabilistic model of local protein structure| journal = Proceedings of the National Academy of Sciences| volume = 105| issue = 26| pages = 8932–7| year = 2008| last1 = Boomsma | first1 = W. | last2 = Mardia | first2 = K. V. | last3 = Taylor | first3 = C. C. | last4 = Ferkinghoff-Borg | first4 = J. | last5 = Krogh | first5 = A. | last6 = Hamelryck | first6 = T. | pmid=18579771 | pmc=2440424| bibcode = 2008PNAS..105.8932B| doi-access = free}}{{cite journal |vauthors=Shapovalov MV, Dunbrack, RL|title=A smoothed backbone-dependent rotamer library for proteins derived from adaptive kernel density estimates and regressions|journal=Structure |volume=19|issue=6|pages=844–858|year=2011|doi=10.1016/j.str.2011.03.019 |pmid=21645855|pmc=3118414}} such as backbone-dependent rotamer libraries.

Definition

The bivariate von Mises distribution is a probability distribution defined on the torus, S^1 \times S^1 in \mathbb{R}^3.

The probability density function of the general bivariate von Mises distribution for the angles \phi, \psi \in [0, 2\pi] is given by

:

f(\phi, \psi) \propto \exp [ \kappa_1 \cos(\phi - \mu) + \kappa_2 \cos(\psi - \nu) + (\cos(\phi-\mu), \sin(\phi-\mu)) \mathbf{A} (\cos(\psi - \nu), \sin(\psi - \nu))^T ],

where \mu and \nu are the means for \phi and \psi, \kappa_1 and \kappa_2 their concentration and the matrix \mathbf{A} \in \mathbb{M}(2,2) is related to their correlation.

Two commonly used variants of the bivariate von Mises distribution are the sine and cosine variant.

The cosine variant of the bivariate von Mises distribution has the probability density function

:

f(\phi, \psi) = Z_c(\kappa_1, \kappa_2, \kappa_3) \ \exp [ \kappa_1 \cos(\phi - \mu) + \kappa_2 \cos(\psi - \nu) - \kappa_3 \cos(\phi - \mu - \psi + \nu) ],

where \mu and \nu are the means for \phi and \psi, \kappa_1 and \kappa_2 their concentration and \kappa_3 is related to their correlation. Z_c is the normalization constant. This distribution with \kappa_3=0 has been used for kernel density estimates of the distribution of the protein dihedral angles \phi and \psi.

The sine variant has the probability density function{{Cite journal | doi = 10.1093/biomet/89.3.719| title = Probabilistic model for two dependent circular variables| journal = Biometrika| volume = 89| issue = 3| pages = 719–723| year = 2002| last1 = Singh | first1 = H.}}

:

f(\phi, \psi) = Z_s(\kappa_1, \kappa_2, \kappa_3) \ \exp [ \kappa_1 \cos(\phi - \mu) + \kappa_2 \cos(\psi - \nu) + \kappa_3 \sin(\phi - \mu) \sin(\psi - \nu) ],

where the parameters have the same interpretation.

See also

References