metromc 0.2.0

Markov chain Monte Carlo sampling using the Independence Metropolis-Hastings algorithm.
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21

use statrs::distribution::Continuous;
use crate::Pdf;

pub struct Beta {
    dist: statrs::distribution::Beta,
}
impl Beta {
    #[inline(always)]
    pub fn new(alpha: f64, beta: f64) -> Self {
        Self {
            dist: statrs::distribution::Beta::new(alpha, beta).unwrap(),
        }
    }
}

impl Pdf for Beta {
    #[inline(always)]
    fn prob(&self, x: f64) -> f64 {
        self.dist.pdf(x)
    }
}