#编程语言 #概率 #编程 #贝叶斯 #MCMC #生成 #随机

ferric

一个具有随机变量声明性语法的概率编程语言

2个版本

0.1.2 2022年6月25日
0.1.0 2022年1月9日

#601 in 数学

MIT/Apache

12KB
217

Github Actions Tests crates.io Coverage Status

Ferric

一个具有声明性语法的Rust概率编程语言。

用法

将此添加到您的 Cargo.toml

[dependencies]
ferric = "0.1"

示例

use std::time::Instant;
use ferric::make_model;

make_model! {
    mod grass;
    use ferric::distributions::Bernoulli;

    let rain : bool ~ Bernoulli::new( 0.2 );

    let sprinkler : bool ~
        if rain {
            Bernoulli::new( 0.01 )
        } else {
            Bernoulli::new( 0.4 )
        };

    let grass_wet : bool ~ Bernoulli::new(
        if sprinkler && rain { 0.99 }
        else if sprinkler && !rain { 0.9 }
        else if !sprinkler && rain { 0.8 }
        else { 0.0 }
    );

    observe grass_wet;
    query rain;
    query sprinkler;
}

fn main() {
    let model = grass::Model {grass_wet: true};
    let mut num_rain = 0;
    let mut num_sprinkler = 0;
    let num_samples = 100000;
    let start = Instant::now();
    for sample in model.sample_iter().take(num_samples) {
        if sample.rain {
            num_rain += 1;
        }
        if sample.sprinkler {
            num_sprinkler += 1;
        }
    }
    let num_samples = num_samples as f64;
    println!(
        "posterior: rain = {} sprinkler = {}. Elapsed {} millisec for {} samples",
        (num_rain as f64) / num_samples,
        (num_sprinkler as f64) / num_samples,
        start.elapsed().as_millis(),
        num_samples,
    );
}

许可证

根据您的选择,受以下任一许可证的许可:

贡献

除非您明确声明,否则您提交给作品以供包含在内的任何贡献,根据Apache-2.0许可证的定义,将如上所述双重许可,不附加任何额外的条款或条件。

依赖项

~2.5MB
~56K SLoC