Struct rv::dist::Gaussian

pub struct Gaussian { /* fields omitted */ }

Gaussian / Normal distribution, N(μ, σ) over real values.

Examples

Compute the KL Divergence between two Gaussians.

use rv::prelude::*;

let gauss_1 = Gaussian::new(0.1, 2.3).unwrap();
let gauss_2 = Gaussian::standard();

// KL is not symmetric
let kl_12 = gauss_1.kl(&gauss_2);
let kl_21 = gauss_2.kl(&gauss_1);

// ... but kl_sym is because it's the sum of KL(P|Q) and KL(Q|P)
let kl_sym = gauss_1.kl_sym(&gauss_2);
assert!((kl_sym - (kl_12 + kl_21)).abs() < 1E-12);

Methods

impl Gaussian

pub fn set_mu(&mut self, val: f64) -> &mut Self

Mean

pub fn set_sigma(&mut self, val: f64) -> &mut Self

Standard deviation

impl Gaussian

pub fn new(mu: f64, sigma: f64) -> Result<Self>

Create a new Gaussian distribution

Aruments

• mu: mean
• sigma: standard deviation

pub fn new_unchecked(mu: f64, sigma: f64) -> Self

Creates a new Gaussian without checking whether the parameters are valid.

pub fn standard() -> Self

Standard normal

Example

let gauss = Gaussian::standard();

assert_eq!(gauss, Gaussian::new(0.0, 1.0).unwrap());

pub fn mu(&self) -> f64

Get mu parameter

Example

let gauss = Gaussian::new(2.0, 1.5).unwrap();

assert_eq!(gauss.mu(), 2.0);

pub fn sigma(&self) -> f64

Get sigma parameter

Example

let gauss = Gaussian::new(2.0, 1.5).unwrap();

assert_eq!(gauss.sigma(), 1.5);

Trait Implementations

impl Rv<f32> for Gaussian

fn ln_f(&self, x: &f32) -> f64

Probability function

fn draw<R: Rng>(&self, rng: &mut R) -> f32

Single draw from the Rv

fn sample<R: Rng>(&self, n: usize, rng: &mut R) -> Vec<f32>

Multiple draws of the Rv

fn f(&self, x: &X) -> f64

Probability function

impl Rv<f64> for Gaussian

fn ln_f(&self, x: &f64) -> f64

Probability function

fn draw<R: Rng>(&self, rng: &mut R) -> f64

Single draw from the Rv

fn sample<R: Rng>(&self, n: usize, rng: &mut R) -> Vec<f64>

Multiple draws of the Rv

fn f(&self, x: &X) -> f64

Probability function

impl Rv<Gaussian> for NormalGamma

fn ln_f(&self, x: &Gaussian) -> f64

Probability function

fn draw<R: Rng>(&self, rng: &mut R) -> Gaussian

Single draw from the Rv

fn f(&self, x: &X) -> f64

Probability function

fn sample<R: Rng>(&self, n: usize, rng: &mut R) -> Vec<X>

Multiple draws of the Rv

impl Support<f32> for Gaussian

fn supports(&self, x: &f32) -> bool

Returns true if x is in the support of the Rv

impl Support<f64> for Gaussian

fn supports(&self, x: &f64) -> bool

Returns true if x is in the support of the Rv

impl Support<Gaussian> for NormalGamma

fn supports(&self, x: &Gaussian) -> bool

Returns true if x is in the support of the Rv

impl ContinuousDistr<f32> for Gaussian

fn pdf(&self, x: &X) -> f64

The value of the Probability Density Function (PDF) at x

fn ln_pdf(&self, x: &X) -> f64

The value of the log Probability Density Function (PDF) at x

impl ContinuousDistr<f64> for Gaussian

fn pdf(&self, x: &X) -> f64

The value of the Probability Density Function (PDF) at x

fn ln_pdf(&self, x: &X) -> f64

The value of the log Probability Density Function (PDF) at x

impl ContinuousDistr<Gaussian> for NormalGamma

fn pdf(&self, x: &X) -> f64

The value of the Probability Density Function (PDF) at x

fn ln_pdf(&self, x: &X) -> f64

The value of the log Probability Density Function (PDF) at x

impl Cdf<f32> for Gaussian

fn cdf(&self, x: &f32) -> f64

The value of the Cumulative Density Function at x

fn sf(&self, x: &X) -> f64

Survival function, 1 - CDF(x)

impl Cdf<f64> for Gaussian

fn cdf(&self, x: &f64) -> f64

The value of the Cumulative Density Function at x

fn sf(&self, x: &X) -> f64

Survival function, 1 - CDF(x)

impl InverseCdf<f32> for Gaussian

fn invcdf(&self, p: f64) -> f32

The value of the x at the given probability in the CDF

fn quantile(&self, p: f64) -> X

Alias for invcdf

fn interval(&self, p: f64) -> (X, X)

Interval containing p proportion for the probability

impl InverseCdf<f64> for Gaussian

fn invcdf(&self, p: f64) -> f64

The value of the x at the given probability in the CDF

fn quantile(&self, p: f64) -> X

Alias for invcdf

fn interval(&self, p: f64) -> (X, X)

Interval containing p proportion for the probability

impl Mean<f32> for Gaussian

fn mean(&self) -> Option<f32>

Returns None if the mean is undefined

impl Mean<f64> for Gaussian

fn mean(&self) -> Option<f64>

Returns None if the mean is undefined

impl Median<f32> for Gaussian

fn median(&self) -> Option<f32>

Returns None if the median is undefined

impl Median<f64> for Gaussian

fn median(&self) -> Option<f64>

Returns None if the median is undefined

impl Mode<f32> for Gaussian

fn mode(&self) -> Option<f32>

Returns None if the mode is undefined or is not a single value

impl Mode<f64> for Gaussian

fn mode(&self) -> Option<f64>

Returns None if the mode is undefined or is not a single value

impl Variance<f64> for Gaussian

fn variance(&self) -> Option<f64>

Returns None if the variance is undefined

impl Entropy for Gaussian

fn entropy(&self) -> f64

The entropy, H(X)

impl Skewness for Gaussian

fn skewness(&self) -> Option<f64>

impl Kurtosis for Gaussian

fn kurtosis(&self) -> Option<f64>

impl KlDivergence for Gaussian

fn kl(&self, other: &Self) -> f64

The KL divergence, KL(P|Q) between this distribution, P, and another, Q

fn kl_sym(&self, other: &Self) -> f64

Symmetrised divergence, KL(P|Q) + KL(Q|P)

impl HasSuffStat<f32> for Gaussian

type Stat = GaussianSuffStat

fn empty_suffstat(&self) -> Self::Stat

impl HasSuffStat<f64> for Gaussian

type Stat = GaussianSuffStat

fn empty_suffstat(&self) -> Self::Stat

impl ConjugatePrior<f64, Gaussian> for NormalGamma

type Posterior = Self

fn posterior(&self, x: &DataOrSuffStat<f64, Gaussian>) -> Self

Computes the posterior distribution from the data

fn ln_m(&self, x: &DataOrSuffStat<f64, Gaussian>) -> f64

Log marginal likelihood

fn ln_pp(&self, y: &f64, x: &DataOrSuffStat<f64, Gaussian>) -> f64

Log posterior predictive of y given x

fn m(&self, x: &DataOrSuffStat<X, Fx>) -> f64

Marginal likelihood of x

fn pp(&self, y: &X, x: &DataOrSuffStat<X, Fx>) -> f64

Posterior Predictive distribution

impl<'_> From<&'_ Gaussian> for String

fn from(gauss: &Gaussian) -> String

Performs the conversion.

impl Clone for Gaussian

fn clone(&self) -> Gaussian

Returns a copy of the value.

fn clone_from(&mut self, source: &Self)

Performs copy-assignment from source.

impl Default for Gaussian

fn default() -> Self

Returns the "default value" for a type.

impl PartialEq<Gaussian> for Gaussian

fn eq(&self, other: &Gaussian) -> bool

This method tests for self and other values to be equal, and is used by ==.

fn ne(&self, other: &Gaussian) -> bool

This method tests for !=.

impl PartialOrd<Gaussian> for Gaussian

fn partial_cmp(&self, other: &Gaussian) -> Option<Ordering>

This method returns an ordering between self and other values if one exists.

fn lt(&self, other: &Gaussian) -> bool

This method tests less than (for self and other) and is used by the < operator.

fn le(&self, other: &Gaussian) -> bool

This method tests less than or equal to (for self and other) and is used by the <= operator.

fn gt(&self, other: &Gaussian) -> bool

This method tests greater than (for self and other) and is used by the > operator.

fn ge(&self, other: &Gaussian) -> bool

This method tests greater than or equal to (for self and other) and is used by the >= operator.

impl Display for Gaussian

fn fmt(&self, f: &mut Formatter) -> Result

Formats the value using the given formatter.

impl Debug for Gaussian

fn fmt(&self, f: &mut Formatter) -> Result

Formats the value using the given formatter.