1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92
//! Provides random variables for probabilistic modeling. //! //! The [`dist`] module provides a number of probability distributions with //! various traits implemented in the [`traits`] module. You can do all the //! standard probability distribution stuff like evaluate the PDF/PMF and draw //! values of different types. //! //! The [`prelude`] module provides all the distributions, all the traits, and //! creates a few useful type aliases. //! //! # Design //! //! Random variables are designed to be flexible. For example, we don't just //! want a `Beta` distribution that works with `f64`; we want it to work with a //! bunch of things //! //! ``` //! extern crate rand; //! extern crate rv; //! //! use rv::prelude::*; //! //! // Beta(0.5, 0.5) //! let beta = Beta::jeffreys(); //! //! let mut rng = rand::thread_rng(); //! //! // 100 f64 weights in (0, 1) //! let f64s: Vec<f64> = beta.sample(100, &mut rng); //! let pdf_x = beta.ln_pdf(&f64s[42]); //! //! // 100 f32 weights in (0, 1) //! let f32s: Vec<f32> = beta.sample(100, &mut rng); //! let pdf_y = beta.ln_pdf(&f32s[42]); //! //! // 100 Bernoulli distributions -- Beta is a prior on the weight //! let berns: Vec<Bernoulli> = beta.sample(100, &mut rng); //! let pdf_bern = beta.ln_pdf(&berns[42]); //! ``` //! //! # Examples //! //! For more examples, check out the `examples` directory. //! //! ## Conjugate analysis of coin flips //! //! ```rust //! extern crate rand; //! extern crate rv; //! //! use rand::Rng; //! use rv::prelude::*; //! //! fn main() { //! let mut rng = rand::thread_rng(); //! //! // A sequence of observations //! let flips = vec![true, false, true, true, true, false, true]; //! //! // Construct the Jeffreys prior of Beta(0.5, 0.5) //! let prior = Beta::jeffreys(); //! //! // Packages the data in a wrapper that marks it as having come from //! // Bernoulli trials. //! let obs: BernoulliData<bool> = DataOrSuffStat::Data(&flips); //! //! // Generate the posterior distribution P(θ|x); the distribution of //! // probable coin weights //! let posterior: Beta = prior.posterior(&obs); //! //! // What is the probability that the next flip would come up heads //! // (true) given the observed flips (posterior predictive)? //! let p_heads = prior.pp(&true, &obs); //! } //! ``` #![feature(associated_type_defaults)] #![feature(test)] #[cfg(feature = "serde_support")] #[macro_use] extern crate serde_derive; pub mod consts; pub mod data; pub mod dist; pub mod misc; mod model; pub mod prelude; mod priors; pub mod traits; pub use model::ConjugateModel;