In [1]:

    

import numpy as np


    

In [2]:

    

import abtools as ab


    

In [3]:

    

m1 = ab.Lognormal(mu=1, std=1.83, n=5000)


    

In [4]:

    

b1 = ab.Bernoulli(alpha=5000, beta=90000)
b2 = ab.Bernoulli(alpha=6000, beta=90000)


    

In [5]:

    

(b1 * 20 >> b2).mean()


    

    
Out[5]:





0.062497944507175045

In [29]:

    

bl = ab.BLModel(mu=1.29, std=1.86, alpha=3000, n=10000)


    

In [30]:

    

bl2 = ab.BLModel(mu=1.24, std=1.81, alpha=1000, n=10000)


    

In [31]:

    

bl.mean()


    

    
Out[31]:





6.1532570916786105

In [32]:

    

bl2.mean()


    

    
Out[32]:





1.7838088202520972

In [33]:

    

(bl2 >> bl).mean()


    

    
Out[33]:





3.9594894351879102