Struct statrs::distribution::StudentsT

pub struct StudentsT { /* fields omitted */ }

Implements the Student's T distribution

Examples

use statrs::distribution::{StudentsT, Continuous};
use statrs::statistics::Mean;
use statrs::prec;

let n = StudentsT::new(0.0, 1.0, 2.0).unwrap();
assert_eq!(n.mean(), 0.0);
assert!(prec::almost_eq(n.pdf(0.0), 0.353553390593274, 1e-15));

Methods

impl StudentsT

pub fn new(location: f64, scale: f64, freedom: f64) -> Result<StudentsT>

Constructs a new student's t-distribution with location location, scale scale, and freedom freedom.

Errors

Returns an error if any of location, scale, or freedom are NaN. Returns an error if scale <= 0.0 or freedom <= 0.0

Examples

use statrs::distribution::StudentsT;

let mut result = StudentsT::new(0.0, 1.0, 2.0);
assert!(result.is_ok());

result = StudentsT::new(0.0, 0.0, 0.0);
assert!(result.is_err());

pub fn location(&self) -> f64

Returns the location of the student's t-distribution

Examples

use statrs::distribution::StudentsT;

let n = StudentsT::new(0.0, 1.0, 2.0).unwrap();
assert_eq!(n.location(), 0.0);

pub fn scale(&self) -> f64

Returns the scale of the student's t-distribution

Examples

use statrs::distribution::StudentsT;

let n = StudentsT::new(0.0, 1.0, 2.0).unwrap();
assert_eq!(n.scale(), 1.0);

pub fn freedom(&self) -> f64

Returns the freedom of the student's t-distribution

Examples

use statrs::distribution::StudentsT;

let n = StudentsT::new(0.0, 1.0, 2.0).unwrap();
assert_eq!(n.freedom(), 2.0);

Trait Implementations

impl Univariate<f64, f64> for StudentsT

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

Calculates the cumulative distribution function for the student's t-distribution at x

Formula

if x < μ {
    (1 / 2) * I(t, v / 2, 1 / 2)
} else {
    1 - (1 / 2) * I(t, v / 2, 1 / 2)
}

where t = v / (v + k^2), k = (x - μ) / σ, μ is the location, σ is the scale, v is the freedom, and I is the regularized incomplete beta function

impl Continuous<f64, f64> for StudentsT

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

Calculates the probability density function for the student's t-distribution at x

Formula

Γ((v + 1) / 2) / (sqrt(vπ) * Γ(v / 2) * σ) * (1 + k^2 / v)^(-1 / 2 * (v
+ 1))

where k = (x - μ) / σ, μ is the location, σ is the scale, v is the freedom, and Γ is the gamma function

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

Calculates the log probability density function for the student's t-distribution at x

Formula

ln(Γ((v + 1) / 2) / (sqrt(vπ) * Γ(v / 2) * σ) * (1 + k^2 / v)^(-1 / 2 *
(v + 1)))

where k = (x - μ) / σ, μ is the location, σ is the scale, v is the freedom, and Γ is the gamma function

impl Min<f64> for StudentsT

fn min(&self) -> f64

Returns the minimum value in the domain of the student's t-distribution representable by a double precision float

Formula

-INF

impl Max<f64> for StudentsT

fn max(&self) -> f64

Returns the maximum value in the domain of the student's t-distribution representable by a double precision float

Formula

INF

impl Mean<f64> for StudentsT

fn mean(&self) -> f64

Returns the mean of the student's t-distribution

Panics

If freedom <= 1.0

Formula

μ

where μ is the location

impl CheckedMean<f64> for StudentsT

fn checked_mean(&self) -> Result<f64>

Returns the mean of the student's t-distribution

Errors

If freedom <= 1.0

Formula

μ

where μ is the location

impl Variance<f64> for StudentsT

fn variance(&self) -> f64

Returns the variance of the student's t-distribution

Panics

If freedom <= 1.0

Formula

if v == INF {
    σ^2
} else if freedom > 2.0 {
    v * σ^2 / (v - 2)
} else {
    INF
}

where σ is the scale and v is the freedom

fn std_dev(&self) -> f64

Returns the standard deviation of the student's t-distribution

Panics

If freedom <= 1.0

Formula

let variance = if v == INF {
    σ^2
} else if freedom > 2.0 {
    v * σ^2 / (v - 2)
} else {
    INF
}
sqrt(variance)

where σ is the scale and v is the freedom

impl CheckedVariance<f64> for StudentsT

fn checked_variance(&self) -> Result<f64>

Returns the variance of the student's t-distribution

Errors

If freedom <= 1.0

Formula

if v == INF {
    σ^2
} else if freedom > 2.0 {
    v * σ^2 / (v - 2)
} else {
    INF
}

where σ is the scale and v is the freedom

fn checked_std_dev(&self) -> Result<f64>

Returns the standard deviation of the student's t-distribution

Errors

If freedom <= 1.0

Formula

let variance = if v == INF {
    σ^2
} else if freedom > 2.0 {
    v * σ^2 / (v - 2)
} else {
    INF
}
sqrt(variance)

where σ is the scale and v is the freedom

impl Entropy<f64> for StudentsT

fn entropy(&self) -> f64

Returns the entropy for the student's t-distribution

Panics

If location != 0.0 && scale != 1.0

Formula

(v + 1) / 2 * (ψ((v + 1) / 2) - ψ(v / 2)) + ln(sqrt(v) * B(v / 2, 1 /
2))

where v is the freedom, ψ is the digamma function, and B is the beta function

impl CheckedEntropy<f64> for StudentsT

fn checked_entropy(&self) -> Result<f64>

Returns the entropy for the student's t-distribution

Errors

If location != 0.0 && scale != 1.0

Formula

(v + 1) / 2 * (ψ((v + 1) / 2) - ψ(v / 2)) + ln(sqrt(v) * B(v / 2, 1 /
2))

where v is the freedom, ψ is the digamma function, and B is the beta function

impl Skewness<f64> for StudentsT

fn skewness(&self) -> f64

Returns the skewness of the student's t-distribution

Panics

If x <= 3.0

Formula

0

impl CheckedSkewness<f64> for StudentsT

fn checked_skewness(&self) -> Result<f64>

Returns the skewness of the student's t-distribution

Errors

If x <= 3.0

Formula

0

impl Median<f64> for StudentsT

fn median(&self) -> f64

Returns the median of the student's t-distribution

Formula

μ

where μ is the location

impl Mode<f64> for StudentsT

fn mode(&self) -> f64

Returns the mode of the student's t-distribution

Formula

μ

where μ is the location

impl Clone for StudentsT

fn clone(&self) -> StudentsT

Returns a copy of the value.

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

Performs copy-assignment from source.

impl PartialEq<StudentsT> for StudentsT

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

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

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

This method tests for !=.

impl Copy for StudentsT

impl Debug for StudentsT

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

Formats the value using the given formatter.

impl Distribution<f64> for StudentsT

fn sample<R: Rng + ?Sized>(&self, r: &mut R) -> f64

Generate a random value of T, using rng as the source of randomness.

fn sample_iter<R>(self, rng: R) -> DistIter<Self, R, T> where
    R: Rng, 

Create an iterator that generates random values of T, using rng as the source of randomness.