compiled-knowledge 4.0.0a20__cp313-cp313-macosx_10_13_x86_64.whl

ck/example/alarm.py

@@ -0,0 +1,366 @@
+from ck.pgm import PGM
+
+
+class Alarm(PGM):
+    """
+    This PGM is the well known, pedagogical 'Alarm' Bayesian network.
+    """
+
+    def __init__(self):
+        super().__init__(self.__class__.__name__)
+
+        MinVol = self.new_rv('MinVol', ('Zero', 'Low', 'Normal', 'High'))
+        Press = self.new_rv('Press', ('Zero', 'Low', 'Normal', 'High'))
+        PAP = self.new_rv('PAP', ('Low', 'Normal', 'High'))
+        ExpCO2 = self.new_rv('ExpCO2', ('Zero', 'Low', 'Normal', 'High'))
+        HRBP = self.new_rv('HRBP', ('Low', 'Normal', 'High'))
+        ErrLowOutput = self.new_rv('ErrLowOutput', ('True', 'False'))
+        HRSat = self.new_rv('HRSat', ('Low', 'Normal', 'High'))
+        HREKG = self.new_rv('HREKG', ('Low', 'Normal', 'High'))
+        ErrCauter = self.new_rv('ErrCauter', ('True', 'False'))
+        BP = self.new_rv('BP', ('Low', 'Normal', 'High'))
+        History = self.new_rv('History', ('True', 'False'))
+        CO = self.new_rv('CO', ('Low', 'Normal', 'High'))
+        HR = self.new_rv('HR', ('Low', 'Normal', 'High'))
+        Catechol = self.new_rv('Catechol', ('Normal', 'High'))
+        ArtCO2 = self.new_rv('ArtCO2', ('Low', 'Normal', 'High'))
+        TPR = self.new_rv('TPR', ('Low', 'Normal', 'High'))
+        Anaphylaxis = self.new_rv('Anaphylaxis', ('True', 'False'))
+        SaO2 = self.new_rv('SaO2', ('Low', 'Normal', 'High'))
+        PVSat = self.new_rv('PVSat', ('Low', 'Normal', 'High'))
+        FiO2 = self.new_rv('FiO2', ('Low', 'Normal'))
+        VentAlv = self.new_rv('VentAlv', ('Zero', 'Low', 'Normal', 'High'))
+        VentLung = self.new_rv('VentLung', ('Zero', 'Low', 'Normal', 'High'))
+        VentTube = self.new_rv('VentTube', ('Zero', 'Low', 'Normal', 'High'))
+        Disconnect = self.new_rv('Disconnect', ('True', 'False'))
+        VentMach = self.new_rv('VentMach', ('Zero', 'Low', 'Normal', 'High'))
+        MinVolSet = self.new_rv('MinVolSet', ('Low', 'Normal', 'High'))
+        KinkedTube = self.new_rv('KinkedTube', ('True', 'False'))
+        Shunt = self.new_rv('Shunt', ('Normal', 'High'))
+        Intubation = self.new_rv('Intubation', ('Normal', 'Esophageal', 'OneSided'))
+        PulmEmbolus = self.new_rv('PulmEmbolus', ('True', 'False'))
+        InsuffAnesth = self.new_rv('InsuffAnesth', ('True', 'False'))
+        PCWP = self.new_rv('PCWP', ('Low', 'Normal', 'High'))
+        CVP = self.new_rv('CVP', ('Low', 'Normal', 'High'))
+        StrokeVolume = self.new_rv('StrokeVolume', ('Low', 'Normal', 'High'))
+        LVEDVolume = self.new_rv('LVEDVolume', ('Low', 'Normal', 'High'))
+        LVFailure = self.new_rv('LVFailure', ('True', 'False'))
+        Hypovolemia = self.new_rv('Hypovolemia', ('True', 'False'))
+
+        pgm_factor0 = self.new_factor(MinVol, VentLung, Intubation)
+        pgm_function_2759353723488 = pgm_factor0.set_cpt()
+        pgm_function_2759353723488.set_cpd((0, 0), (0.97, 0.01, 0.01, 0.01))
+        pgm_function_2759353723488.set_cpd((1, 0), (0.01, 0.97, 0.01, 0.01))
+        pgm_function_2759353723488.set_cpd((2, 0), (0.01, 0.01, 0.97, 0.01))
+        pgm_function_2759353723488.set_cpd((3, 0), (0.01, 0.01, 0.01, 0.97))
+        pgm_function_2759353723488.set_cpd((0, 1), (0.97, 0.01, 0.01, 0.01))
+        pgm_function_2759353723488.set_cpd((1, 1), (0.6, 0.38, 0.01, 0.01))
+        pgm_function_2759353723488.set_cpd((2, 1), (0.5, 0.48, 0.01, 0.01))
+        pgm_function_2759353723488.set_cpd((3, 1), (0.5, 0.48, 0.01, 0.01))
+        pgm_function_2759353723488.set_cpd((0, 2), (0.97, 0.01, 0.01, 0.01))
+        pgm_function_2759353723488.set_cpd((1, 2), (0.01, 0.97, 0.01, 0.01))
+        pgm_function_2759353723488.set_cpd((2, 2), (0.01, 0.01, 0.97, 0.01))
+        pgm_function_2759353723488.set_cpd((3, 2), (0.01, 0.01, 0.01, 0.97))
+        pgm_factor1 = self.new_factor(Press, KinkedTube, Intubation, VentTube)
+        pgm_function_2759353724528 = pgm_factor1.set_cpt()
+        pgm_function_2759353724528.set_cpd((0, 0, 0), (0.97, 0.01, 0.01, 0.01))
+        pgm_function_2759353724528.set_cpd((1, 0, 0), (0.97, 0.01, 0.01, 0.01))
+        pgm_function_2759353724528.set_cpd((0, 1, 0), (0.97, 0.01, 0.01, 0.01))
+        pgm_function_2759353724528.set_cpd((1, 1, 0), (0.97, 0.01, 0.01, 0.01))
+        pgm_function_2759353724528.set_cpd((0, 2, 0), (0.97, 0.01, 0.01, 0.01))
+        pgm_function_2759353724528.set_cpd((1, 2, 0), (0.97, 0.01, 0.01, 0.01))
+        pgm_function_2759353724528.set_cpd((0, 0, 1), (0.01, 0.49, 0.3, 0.2))
+        pgm_function_2759353724528.set_cpd((1, 0, 1), (0.01, 0.97, 0.01, 0.01))
+        pgm_function_2759353724528.set_cpd((0, 1, 1), (0.1, 0.84, 0.05, 0.01))
+        pgm_function_2759353724528.set_cpd((1, 1, 1), (0.4, 0.58, 0.01, 0.01))
+        pgm_function_2759353724528.set_cpd((0, 2, 1), (0.01, 0.29, 0.3, 0.4))
+        pgm_function_2759353724528.set_cpd((1, 2, 1), (0.01, 0.9, 0.08, 0.01))
+        pgm_function_2759353724528.set_cpd((0, 0, 2), (0.01, 0.01, 0.08, 0.9))
+        pgm_function_2759353724528.set_cpd((1, 0, 2), (0.01, 0.01, 0.97, 0.01))
+        pgm_function_2759353724528.set_cpd((0, 1, 2), (0.05, 0.25, 0.25, 0.45))
+        pgm_function_2759353724528.set_cpd((1, 1, 2), (0.2, 0.75, 0.04, 0.01))
+        pgm_function_2759353724528.set_cpd((0, 2, 2), (0.01, 0.01, 0.08, 0.9))
+        pgm_function_2759353724528.set_cpd((1, 2, 2), (0.01, 0.01, 0.38, 0.6))
+        pgm_function_2759353724528.set_cpd((0, 0, 3), (0.01, 0.01, 0.01, 0.97))
+        pgm_function_2759353724528.set_cpd((1, 0, 3), (0.01, 0.01, 0.01, 0.97))
+        pgm_function_2759353724528.set_cpd((0, 1, 3), (0.01, 0.15, 0.25, 0.59))
+        pgm_function_2759353724528.set_cpd((1, 1, 3), (0.2, 0.7, 0.09, 0.01))
+        pgm_function_2759353724528.set_cpd((0, 2, 3), (0.01, 0.01, 0.01, 0.97))
+        pgm_function_2759353724528.set_cpd((1, 2, 3), (0.01, 0.01, 0.01, 0.97))
+        pgm_factor2 = self.new_factor(PAP, PulmEmbolus)
+        pgm_function_2759353723904 = pgm_factor2.set_cpt()
+        pgm_function_2759353723904.set_cpd((0,), (0.01, 0.19, 0.8))
+        pgm_function_2759353723904.set_cpd((1,), (0.05, 0.9, 0.05))
+        pgm_factor3 = self.new_factor(ExpCO2, ArtCO2, VentLung)
+        pgm_function_2759353724112 = pgm_factor3.set_cpt()
+        pgm_function_2759353724112.set_cpd((0, 0), (0.97, 0.01, 0.01, 0.01))
+        pgm_function_2759353724112.set_cpd((1, 0), (0.97, 0.01, 0.01, 0.01))
+        pgm_function_2759353724112.set_cpd((2, 0), (0.97, 0.01, 0.01, 0.01))
+        pgm_function_2759353724112.set_cpd((0, 1), (0.01, 0.97, 0.01, 0.01))
+        pgm_function_2759353724112.set_cpd((1, 1), (0.01, 0.01, 0.97, 0.01))
+        pgm_function_2759353724112.set_cpd((2, 1), (0.01, 0.01, 0.01, 0.97))
+        pgm_function_2759353724112.set_cpd((0, 2), (0.01, 0.97, 0.01, 0.01))
+        pgm_function_2759353724112.set_cpd((1, 2), (0.01, 0.01, 0.97, 0.01))
+        pgm_function_2759353724112.set_cpd((2, 2), (0.01, 0.01, 0.01, 0.97))
+        pgm_function_2759353724112.set_cpd((0, 3), (0.01, 0.97, 0.01, 0.01))
+        pgm_function_2759353724112.set_cpd((1, 3), (0.01, 0.01, 0.97, 0.01))
+        pgm_function_2759353724112.set_cpd((2, 3), (0.01, 0.01, 0.01, 0.97))
+        pgm_factor4 = self.new_factor(HRBP, ErrLowOutput, HR)
+        pgm_function_2759353721200 = pgm_factor4.set_cpt()
+        pgm_function_2759353721200.set_cpd((0, 0), (0.98, 0.01, 0.01))
+        pgm_function_2759353721200.set_cpd((1, 0), (0.98, 0.01, 0.01))
+        pgm_function_2759353721200.set_cpd((0, 1), (0.4, 0.59, 0.01))
+        pgm_function_2759353721200.set_cpd((1, 1), (0.01, 0.98, 0.01))
+        pgm_function_2759353721200.set_cpd((0, 2), (0.3, 0.4, 0.3))
+        pgm_function_2759353721200.set_cpd((1, 2), (0.01, 0.01, 0.98))
+        pgm_factor5 = self.new_factor(ErrLowOutput)
+        pgm_function_2759353722448 = pgm_factor5.set_cpt()
+        pgm_function_2759353722448.set_cpd((), (0.05, 0.95))
+        pgm_factor6 = self.new_factor(HRSat, HR, ErrCauter)
+        pgm_function_2759353722656 = pgm_factor6.set_cpt()
+        pgm_function_2759353722656.set_cpd((0, 0), (0.33333333, 0.33333333, 0.33333333))
+        pgm_function_2759353722656.set_cpd((1, 0), (0.33333333, 0.33333333, 0.33333333))
+        pgm_function_2759353722656.set_cpd((2, 0), (0.33333333, 0.33333333, 0.33333333))
+        pgm_function_2759353722656.set_cpd((0, 1), (0.98, 0.01, 0.01))
+        pgm_function_2759353722656.set_cpd((1, 1), (0.01, 0.98, 0.01))
+        pgm_function_2759353722656.set_cpd((2, 1), (0.01, 0.01, 0.98))
+        pgm_factor7 = self.new_factor(HREKG, HR, ErrCauter)
+        pgm_function_2759353722032 = pgm_factor7.set_cpt()
+        pgm_function_2759353722032.set_cpd((0, 0), (0.33333333, 0.33333333, 0.33333333))
+        pgm_function_2759353722032.set_cpd((1, 0), (0.33333333, 0.33333333, 0.33333333))
+        pgm_function_2759353722032.set_cpd((2, 0), (0.33333333, 0.33333333, 0.33333333))
+        pgm_function_2759353722032.set_cpd((0, 1), (0.98, 0.01, 0.01))
+        pgm_function_2759353722032.set_cpd((1, 1), (0.01, 0.98, 0.01))
+        pgm_function_2759353722032.set_cpd((2, 1), (0.01, 0.01, 0.98))
+        pgm_factor8 = self.new_factor(ErrCauter)
+        pgm_function_2759353722864 = pgm_factor8.set_cpt()
+        pgm_function_2759353722864.set_cpd((), (0.1, 0.9))
+        pgm_factor9 = self.new_factor(BP, CO, TPR)
+        pgm_function_2759353721408 = pgm_factor9.set_cpt()
+        pgm_function_2759353721408.set_cpd((0, 0), (0.98, 0.01, 0.01))
+        pgm_function_2759353721408.set_cpd((1, 0), (0.98, 0.01, 0.01))
+        pgm_function_2759353721408.set_cpd((2, 0), (0.9, 0.09, 0.01))
+        pgm_function_2759353721408.set_cpd((0, 1), (0.98, 0.01, 0.01))
+        pgm_function_2759353721408.set_cpd((1, 1), (0.1, 0.85, 0.05))
+        pgm_function_2759353721408.set_cpd((2, 1), (0.05, 0.2, 0.75))
+        pgm_function_2759353721408.set_cpd((0, 2), (0.3, 0.6, 0.1))
+        pgm_function_2759353721408.set_cpd((1, 2), (0.05, 0.4, 0.55))
+        pgm_function_2759353721408.set_cpd((2, 2), (0.01, 0.09, 0.9))
+        pgm_factor10 = self.new_factor(History, LVFailure)
+        pgm_function_2759353720992 = pgm_factor10.set_cpt()
+        pgm_function_2759353720992.set_cpd((0,), (0.9, 0.1))
+        pgm_function_2759353720992.set_cpd((1,), (0.01, 0.99))
+        pgm_factor11 = self.new_factor(CO, HR, StrokeVolume)
+        pgm_function_2759353723280 = pgm_factor11.set_cpt()
+        pgm_function_2759353723280.set_cpd((0, 0), (0.98, 0.01, 0.01))
+        pgm_function_2759353723280.set_cpd((1, 0), (0.95, 0.04, 0.01))
+        pgm_function_2759353723280.set_cpd((2, 0), (0.8, 0.19, 0.01))
+        pgm_function_2759353723280.set_cpd((0, 1), (0.95, 0.04, 0.01))
+        pgm_function_2759353723280.set_cpd((1, 1), (0.04, 0.95, 0.01))
+        pgm_function_2759353723280.set_cpd((2, 1), (0.01, 0.04, 0.95))
+        pgm_function_2759353723280.set_cpd((0, 2), (0.3, 0.69, 0.01))
+        pgm_function_2759353723280.set_cpd((1, 2), (0.01, 0.3, 0.69))
+        pgm_function_2759353723280.set_cpd((2, 2), (0.01, 0.01, 0.98))
+        pgm_factor12 = self.new_factor(HR, Catechol)
+        pgm_function_2759353723072 = pgm_factor12.set_cpt()
+        pgm_function_2759353723072.set_cpd((0,), (0.1, 0.89, 0.01))
+        pgm_function_2759353723072.set_cpd((1,), (0.01, 0.09, 0.9))
+        pgm_factor13 = self.new_factor(Catechol, InsuffAnesth, SaO2, TPR, ArtCO2)
+        pgm_function_2759353724320 = pgm_factor13.set_cpt()
+        pgm_function_2759353724320.set_cpd((0, 0, 0, 0), (0.01, 0.99))
+        pgm_function_2759353724320.set_cpd((1, 0, 0, 0), (0.05, 0.95))
+        pgm_function_2759353724320.set_cpd((0, 1, 0, 0), (0.01, 0.99))
+        pgm_function_2759353724320.set_cpd((1, 1, 0, 0), (0.1, 0.9))
+        pgm_function_2759353724320.set_cpd((0, 2, 0, 0), (0.01, 0.99))
+        pgm_function_2759353724320.set_cpd((1, 2, 0, 0), (0.95, 0.05))
+        pgm_function_2759353724320.set_cpd((0, 0, 1, 0), (0.01, 0.99))
+        pgm_function_2759353724320.set_cpd((1, 0, 1, 0), (0.05, 0.95))
+        pgm_function_2759353724320.set_cpd((0, 1, 1, 0), (0.01, 0.99))
+        pgm_function_2759353724320.set_cpd((1, 1, 1, 0), (0.95, 0.05))
+        pgm_function_2759353724320.set_cpd((0, 2, 1, 0), (0.05, 0.95))
+        pgm_function_2759353724320.set_cpd((1, 2, 1, 0), (0.99, 0.01))
+        pgm_function_2759353724320.set_cpd((0, 0, 2, 0), (0.01, 0.99))
+        pgm_function_2759353724320.set_cpd((1, 0, 2, 0), (0.05, 0.95))
+        pgm_function_2759353724320.set_cpd((0, 1, 2, 0), (0.05, 0.95))
+        pgm_function_2759353724320.set_cpd((1, 1, 2, 0), (0.95, 0.05))
+        pgm_function_2759353724320.set_cpd((0, 2, 2, 0), (0.05, 0.95))
+        pgm_function_2759353724320.set_cpd((1, 2, 2, 0), (0.95, 0.05))
+        pgm_function_2759353724320.set_cpd((0, 0, 0, 1), (0.01, 0.99))
+        pgm_function_2759353724320.set_cpd((1, 0, 0, 1), (0.05, 0.95))
+        pgm_function_2759353724320.set_cpd((0, 1, 0, 1), (0.01, 0.99))
+        pgm_function_2759353724320.set_cpd((1, 1, 0, 1), (0.1, 0.9))
+        pgm_function_2759353724320.set_cpd((0, 2, 0, 1), (0.01, 0.99))
+        pgm_function_2759353724320.set_cpd((1, 2, 0, 1), (0.95, 0.05))
+        pgm_function_2759353724320.set_cpd((0, 0, 1, 1), (0.01, 0.99))
+        pgm_function_2759353724320.set_cpd((1, 0, 1, 1), (0.05, 0.95))
+        pgm_function_2759353724320.set_cpd((0, 1, 1, 1), (0.01, 0.99))
+        pgm_function_2759353724320.set_cpd((1, 1, 1, 1), (0.95, 0.05))
+        pgm_function_2759353724320.set_cpd((0, 2, 1, 1), (0.05, 0.95))
+        pgm_function_2759353724320.set_cpd((1, 2, 1, 1), (0.99, 0.01))
+        pgm_function_2759353724320.set_cpd((0, 0, 2, 1), (0.01, 0.99))
+        pgm_function_2759353724320.set_cpd((1, 0, 2, 1), (0.05, 0.95))
+        pgm_function_2759353724320.set_cpd((0, 1, 2, 1), (0.05, 0.95))
+        pgm_function_2759353724320.set_cpd((1, 1, 2, 1), (0.95, 0.05))
+        pgm_function_2759353724320.set_cpd((0, 2, 2, 1), (0.05, 0.95))
+        pgm_function_2759353724320.set_cpd((1, 2, 2, 1), (0.99, 0.01))
+        pgm_function_2759353724320.set_cpd((0, 0, 0, 2), (0.01, 0.99))
+        pgm_function_2759353724320.set_cpd((1, 0, 0, 2), (0.01, 0.99))
+        pgm_function_2759353724320.set_cpd((0, 1, 0, 2), (0.01, 0.99))
+        pgm_function_2759353724320.set_cpd((1, 1, 0, 2), (0.1, 0.9))
+        pgm_function_2759353724320.set_cpd((0, 2, 0, 2), (0.01, 0.99))
+        pgm_function_2759353724320.set_cpd((1, 2, 0, 2), (0.3, 0.7))
+        pgm_function_2759353724320.set_cpd((0, 0, 1, 2), (0.01, 0.99))
+        pgm_function_2759353724320.set_cpd((1, 0, 1, 2), (0.01, 0.99))
+        pgm_function_2759353724320.set_cpd((0, 1, 1, 2), (0.01, 0.99))
+        pgm_function_2759353724320.set_cpd((1, 1, 1, 2), (0.3, 0.7))
+        pgm_function_2759353724320.set_cpd((0, 2, 1, 2), (0.01, 0.99))
+        pgm_function_2759353724320.set_cpd((1, 2, 1, 2), (0.99, 0.01))
+        pgm_function_2759353724320.set_cpd((0, 0, 2, 2), (0.01, 0.99))
+        pgm_function_2759353724320.set_cpd((1, 0, 2, 2), (0.01, 0.99))
+        pgm_function_2759353724320.set_cpd((0, 1, 2, 2), (0.01, 0.99))
+        pgm_function_2759353724320.set_cpd((1, 1, 2, 2), (0.3, 0.7))
+        pgm_function_2759353724320.set_cpd((0, 2, 2, 2), (0.01, 0.99))
+        pgm_function_2759353724320.set_cpd((1, 2, 2, 2), (0.3, 0.7))
+        pgm_factor14 = self.new_factor(ArtCO2, VentAlv)
+        pgm_function_2759353720784 = pgm_factor14.set_cpt()
+        pgm_function_2759353720784.set_cpd((0,), (0.01, 0.01, 0.98))
+        pgm_function_2759353720784.set_cpd((1,), (0.01, 0.01, 0.98))
+        pgm_function_2759353720784.set_cpd((2,), (0.04, 0.92, 0.04))
+        pgm_function_2759353720784.set_cpd((3,), (0.9, 0.09, 0.01))
+        pgm_factor15 = self.new_factor(TPR, Anaphylaxis)
+        pgm_function_2759353721616 = pgm_factor15.set_cpt()
+        pgm_function_2759353721616.set_cpd((0,), (0.98, 0.01, 0.01))
+        pgm_function_2759353721616.set_cpd((1,), (0.3, 0.4, 0.3))
+        pgm_factor16 = self.new_factor(Anaphylaxis)
+        pgm_function_2759353721824 = pgm_factor16.set_cpt()
+        pgm_function_2759353721824.set_cpd((), (0.01, 0.99))
+        pgm_factor17 = self.new_factor(SaO2, Shunt, PVSat)
+        pgm_function_2759353722240 = pgm_factor17.set_cpt()
+        pgm_function_2759353722240.set_cpd((0, 0), (0.98, 0.01, 0.01))
+        pgm_function_2759353722240.set_cpd((1, 0), (0.98, 0.01, 0.01))
+        pgm_function_2759353722240.set_cpd((0, 1), (0.01, 0.98, 0.01))
+        pgm_function_2759353722240.set_cpd((1, 1), (0.98, 0.01, 0.01))
+        pgm_function_2759353722240.set_cpd((0, 2), (0.01, 0.01, 0.98))
+        pgm_function_2759353722240.set_cpd((1, 2), (0.69, 0.3, 0.01))
+        pgm_factor18 = self.new_factor(PVSat, VentAlv, FiO2)
+        pgm_function_2759355657712 = pgm_factor18.set_cpt()
+        pgm_function_2759355657712.set_cpd((0, 0), (0.98, 0.01, 0.01))
+        pgm_function_2759355657712.set_cpd((1, 0), (0.98, 0.01, 0.01))
+        pgm_function_2759355657712.set_cpd((2, 0), (0.95, 0.04, 0.01))
+        pgm_function_2759355657712.set_cpd((3, 0), (0.95, 0.04, 0.01))
+        pgm_function_2759355657712.set_cpd((0, 1), (0.98, 0.01, 0.01))
+        pgm_function_2759355657712.set_cpd((1, 1), (0.98, 0.01, 0.01))
+        pgm_function_2759355657712.set_cpd((2, 1), (0.01, 0.95, 0.04))
+        pgm_function_2759355657712.set_cpd((3, 1), (0.01, 0.01, 0.98))
+        pgm_factor19 = self.new_factor(FiO2)
+        pgm_function_2759355653760 = pgm_factor19.set_cpt()
+        pgm_function_2759355653760.set_cpd((), (0.01, 0.99))
+        pgm_factor20 = self.new_factor(VentAlv, Intubation, VentLung)
+        pgm_function_2759355654384 = pgm_factor20.set_cpt()
+        pgm_function_2759355654384.set_cpd((0, 0), (0.97, 0.01, 0.01, 0.01))
+        pgm_function_2759355654384.set_cpd((1, 0), (0.97, 0.01, 0.01, 0.01))
+        pgm_function_2759355654384.set_cpd((2, 0), (0.97, 0.01, 0.01, 0.01))
+        pgm_function_2759355654384.set_cpd((0, 1), (0.01, 0.97, 0.01, 0.01))
+        pgm_function_2759355654384.set_cpd((1, 1), (0.01, 0.97, 0.01, 0.01))
+        pgm_function_2759355654384.set_cpd((2, 1), (0.03, 0.95, 0.01, 0.01))
+        pgm_function_2759355654384.set_cpd((0, 2), (0.01, 0.01, 0.97, 0.01))
+        pgm_function_2759355654384.set_cpd((1, 2), (0.01, 0.01, 0.97, 0.01))
+        pgm_function_2759355654384.set_cpd((2, 2), (0.01, 0.94, 0.04, 0.01))
+        pgm_function_2759355654384.set_cpd((0, 3), (0.01, 0.01, 0.01, 0.97))
+        pgm_function_2759355654384.set_cpd((1, 3), (0.01, 0.01, 0.01, 0.97))
+        pgm_function_2759355654384.set_cpd((2, 3), (0.01, 0.88, 0.1, 0.01))
+        pgm_factor21 = self.new_factor(VentLung, KinkedTube, VentTube, Intubation)
+        pgm_function_2759355650848 = pgm_factor21.set_cpt()
+        pgm_function_2759355650848.set_cpd((0, 0, 0), (0.97, 0.01, 0.01, 0.01))
+        pgm_function_2759355650848.set_cpd((1, 0, 0), (0.97, 0.01, 0.01, 0.01))
+        pgm_function_2759355650848.set_cpd((0, 1, 0), (0.95, 0.03, 0.01, 0.01))
+        pgm_function_2759355650848.set_cpd((1, 1, 0), (0.01, 0.97, 0.01, 0.01))
+        pgm_function_2759355650848.set_cpd((0, 2, 0), (0.4, 0.58, 0.01, 0.01))
+        pgm_function_2759355650848.set_cpd((1, 2, 0), (0.01, 0.01, 0.97, 0.01))
+        pgm_function_2759355650848.set_cpd((0, 3, 0), (0.3, 0.68, 0.01, 0.01))
+        pgm_function_2759355650848.set_cpd((1, 3, 0), (0.01, 0.01, 0.01, 0.97))
+        pgm_function_2759355650848.set_cpd((0, 0, 1), (0.97, 0.01, 0.01, 0.01))
+        pgm_function_2759355650848.set_cpd((1, 0, 1), (0.97, 0.01, 0.01, 0.01))
+        pgm_function_2759355650848.set_cpd((0, 1, 1), (0.97, 0.01, 0.01, 0.01))
+        pgm_function_2759355650848.set_cpd((1, 1, 1), (0.97, 0.01, 0.01, 0.01))
+        pgm_function_2759355650848.set_cpd((0, 2, 1), (0.97, 0.01, 0.01, 0.01))
+        pgm_function_2759355650848.set_cpd((1, 2, 1), (0.97, 0.01, 0.01, 0.01))
+        pgm_function_2759355650848.set_cpd((0, 3, 1), (0.97, 0.01, 0.01, 0.01))
+        pgm_function_2759355650848.set_cpd((1, 3, 1), (0.97, 0.01, 0.01, 0.01))
+        pgm_function_2759355650848.set_cpd((0, 0, 2), (0.97, 0.01, 0.01, 0.01))
+        pgm_function_2759355650848.set_cpd((1, 0, 2), (0.97, 0.01, 0.01, 0.01))
+        pgm_function_2759355650848.set_cpd((0, 1, 2), (0.95, 0.03, 0.01, 0.01))
+        pgm_function_2759355650848.set_cpd((1, 1, 2), (0.01, 0.97, 0.01, 0.01))
+        pgm_function_2759355650848.set_cpd((0, 2, 2), (0.5, 0.48, 0.01, 0.01))
+        pgm_function_2759355650848.set_cpd((1, 2, 2), (0.01, 0.01, 0.97, 0.01))
+        pgm_function_2759355650848.set_cpd((0, 3, 2), (0.3, 0.68, 0.01, 0.01))
+        pgm_function_2759355650848.set_cpd((1, 3, 2), (0.01, 0.01, 0.01, 0.97))
+        pgm_factor22 = self.new_factor(VentTube, VentMach, Disconnect)
+        pgm_function_2759355653968 = pgm_factor22.set_cpt()
+        pgm_function_2759355653968.set_cpd((0, 0), (0.97, 0.01, 0.01, 0.01))
+        pgm_function_2759355653968.set_cpd((1, 0), (0.97, 0.01, 0.01, 0.01))
+        pgm_function_2759355653968.set_cpd((2, 0), (0.97, 0.01, 0.01, 0.01))
+        pgm_function_2759355653968.set_cpd((3, 0), (0.97, 0.01, 0.01, 0.01))
+        pgm_function_2759355653968.set_cpd((0, 1), (0.97, 0.01, 0.01, 0.01))
+        pgm_function_2759355653968.set_cpd((1, 1), (0.01, 0.97, 0.01, 0.01))
+        pgm_function_2759355653968.set_cpd((2, 1), (0.01, 0.01, 0.97, 0.01))
+        pgm_function_2759355653968.set_cpd((3, 1), (0.01, 0.01, 0.01, 0.97))
+        pgm_factor23 = self.new_factor(Disconnect)
+        pgm_function_2759355654592 = pgm_factor23.set_cpt()
+        pgm_function_2759355654592.set_cpd((), (0.05, 0.95))
+        pgm_factor24 = self.new_factor(VentMach, MinVolSet)
+        pgm_function_2759355652720 = pgm_factor24.set_cpt()
+        pgm_function_2759355652720.set_cpd((0,), (0.01, 0.97, 0.01, 0.01))
+        pgm_function_2759355652720.set_cpd((1,), (0.01, 0.01, 0.97, 0.01))
+        pgm_function_2759355652720.set_cpd((2,), (0.01, 0.01, 0.01, 0.97))
+        pgm_factor25 = self.new_factor(MinVolSet)
+        pgm_function_2759355654800 = pgm_factor25.set_cpt()
+        pgm_function_2759355654800.set_cpd((), (0.01, 0.98, 0.01))
+        pgm_factor26 = self.new_factor(KinkedTube)
+        pgm_function_2759355651056 = pgm_factor26.set_cpt()
+        pgm_function_2759355651056.set_cpd((), (0.04, 0.96))
+        pgm_factor27 = self.new_factor(Shunt, PulmEmbolus, Intubation)
+        pgm_function_2759355656464 = pgm_factor27.set_cpt()
+        pgm_function_2759355656464.set_cpd((0, 0), (0.1, 0.9))
+        pgm_function_2759355656464.set_cpd((1, 0), (0.95, 0.05))
+        pgm_function_2759355656464.set_cpd((0, 1), (0.1, 0.9))
+        pgm_function_2759355656464.set_cpd((1, 1), (0.95, 0.05))
+        pgm_function_2759355656464.set_cpd((0, 2), (0.01, 0.99))
+        pgm_function_2759355656464.set_cpd((1, 2), (0.05, 0.95))
+        pgm_factor28 = self.new_factor(Intubation)
+        pgm_function_2759355657920 = pgm_factor28.set_cpt()
+        pgm_function_2759355657920.set_cpd((), (0.92, 0.03, 0.05))
+        pgm_factor29 = self.new_factor(PulmEmbolus)
+        pgm_function_2759355652304 = pgm_factor29.set_cpt()
+        pgm_function_2759355652304.set_cpd((), (0.01, 0.99))
+        pgm_factor30 = self.new_factor(InsuffAnesth)
+        pgm_function_2759355650224 = pgm_factor30.set_cpt()
+        pgm_function_2759355650224.set_cpd((), (0.2, 0.8))
+        pgm_factor31 = self.new_factor(PCWP, LVEDVolume)
+        pgm_function_2759355658128 = pgm_factor31.set_cpt()
+        pgm_function_2759355658128.set_cpd((0,), (0.95, 0.04, 0.01))
+        pgm_function_2759355658128.set_cpd((1,), (0.04, 0.95, 0.01))
+        pgm_function_2759355658128.set_cpd((2,), (0.01, 0.04, 0.95))
+        pgm_factor32 = self.new_factor(CVP, LVEDVolume)
+        pgm_function_2759355656256 = pgm_factor32.set_cpt()
+        pgm_function_2759355656256.set_cpd((0,), (0.95, 0.04, 0.01))
+        pgm_function_2759355656256.set_cpd((1,), (0.04, 0.95, 0.01))
+        pgm_function_2759355656256.set_cpd((2,), (0.01, 0.29, 0.7))
+        pgm_factor33 = self.new_factor(StrokeVolume, LVFailure, Hypovolemia)
+        pgm_function_2759355650432 = pgm_factor33.set_cpt()
+        pgm_function_2759355650432.set_cpd((0, 0), (0.98, 0.01, 0.01))
+        pgm_function_2759355650432.set_cpd((1, 0), (0.95, 0.04, 0.01))
+        pgm_function_2759355650432.set_cpd((0, 1), (0.5, 0.49, 0.01))
+        pgm_function_2759355650432.set_cpd((1, 1), (0.05, 0.9, 0.05))
+        pgm_factor34 = self.new_factor(LVEDVolume, Hypovolemia, LVFailure)
+        pgm_function_2759355652512 = pgm_factor34.set_cpt()
+        pgm_function_2759355652512.set_cpd((0, 0), (0.95, 0.04, 0.01))
+        pgm_function_2759355652512.set_cpd((1, 0), (0.01, 0.09, 0.9))
+        pgm_function_2759355652512.set_cpd((0, 1), (0.98, 0.01, 0.01))
+        pgm_function_2759355652512.set_cpd((1, 1), (0.05, 0.9, 0.05))
+        pgm_factor35 = self.new_factor(LVFailure)
+        pgm_function_2759355655216 = pgm_factor35.set_cpt()
+        pgm_function_2759355655216.set_cpd((), (0.05, 0.95))
+        pgm_factor36 = self.new_factor(Hypovolemia)
+        pgm_function_2759355656672 = pgm_factor36.set_cpt()
+        pgm_function_2759355656672.set_cpd((), (0.2, 0.8))

ck/example/asia.py

@@ -0,0 +1,28 @@
+from ck.pgm import PGM
+
+
+class Asia(PGM):
+    """
+    This PGM is the well known, pedagogical 'Asia' Bayesian network.
+    """
+
+    def __init__(self):
+        super().__init__(self.__class__.__name__)
+
+        asia = self.new_rv('asia', ('yes', 'no'))
+        tub = self.new_rv('tub', ('yes', 'no'))
+        smoke = self.new_rv('smoke', ('yes', 'no'))
+        lung = self.new_rv('lung', ('yes', 'no'))
+        bronc = self.new_rv('bronc', ('yes', 'no'))
+        either = self.new_rv('either', ('yes', 'no'))
+        xray = self.new_rv('xray', ('yes', 'no'))
+        dysp = self.new_rv('dysp', ('yes', 'no'))
+
+        self.new_factor(asia).set_dense().set_flat(0.01, 0.99)
+        self.new_factor(tub, asia).set_dense().set_flat(0.05, 0.01, 0.95, 0.99)
+        self.new_factor(smoke).set_dense().set_flat(0.5, 0.5)
+        self.new_factor(lung, smoke).set_dense().set_flat(0.1, 0.01, 0.9, 0.99)
+        self.new_factor(bronc, smoke).set_dense().set_flat(0.6, 0.3, 0.4, 0.7)
+        self.new_factor(either, lung, tub).set_dense().set_flat(1.0, 1.0, 1.0, 0.0, 0.0, 0.0, 0.0, 1.0)
+        self.new_factor(xray, either).set_dense().set_flat(0.98, 0.05, 0.02, 0.95)
+        self.new_factor(dysp, bronc, either).set_dense().set_flat(0.9, 0.8, 0.7, 0.1, 0.1, 0.2, 0.3, 0.9)

ck/example/binary_clique.py

@@ -0,0 +1,32 @@
+import random as _random
+
+from ck.pgm import PGM
+from ck.utils.iter_extras import pairs as _pairs
+
+
+class BinaryClique(PGM):
+    """
+    This PGM is a factor graph with binary factors and a fully connected set of variables
+    """
+
+    def __init__(
+            self,
+            vars_per_clique: int = 5,
+            states_per_var: int = 2,
+            include_unaries: bool = False,
+            random_seed: int = 123456,
+    ):
+        params = (vars_per_clique, states_per_var, include_unaries)
+        super().__init__(f'{self.__class__.__name__}({",".join(str(param) for param in params)})')
+
+        random_stream = _random.Random(random_seed).random
+
+        rvs = [self.new_rv(f'x{i}', states_per_var) for i in range(vars_per_clique)]
+
+        for x_i, x_j in _pairs(rvs):
+            self.new_factor(x_i, x_j).set_dense().set_stream(random_stream)
+
+        if include_unaries:
+            if include_unaries:
+                for rv in rvs:
+                    self.new_factor(rv).set_dense().set_stream(random_stream)

ck/example/bow_tie.py

@@ -0,0 +1,33 @@
+import random as _random
+
+from ck.pgm import PGM
+
+
+class BowTie(PGM):
+    """
+    This PGM is the 'BowTie' factor graph.
+
+    A BowTie is a factor graph with five random variables (x1, ..., x5).
+    One factor connects: x1, x2, x3.
+    Another factor connects: x1, x4, x5.
+    """
+
+    def __init__(
+            self,
+            states_per_var: int = 2,
+            random_seed: int = 123456,
+    ):
+        params = (states_per_var,)
+        super().__init__(f'{self.__class__.__name__}({",".join(str(param) for param in params)})')
+
+        random_stream = _random.Random(random_seed).random
+        rand_iter = iter(random_stream, None)
+
+        x1 = self.new_rv('x1', states_per_var)
+        x2 = self.new_rv('x2', states_per_var)
+        x3 = self.new_rv('x3', states_per_var)
+        x4 = self.new_rv('x4', states_per_var)
+        x5 = self.new_rv('x5', states_per_var)
+
+        self.new_factor(x1, x2, x3).set_dense().set_iter(rand_iter)
+        self.new_factor(x1, x4, x5).set_dense().set_iter(rand_iter)

ck/example/cancer.py

@@ -0,0 +1,37 @@
+from ck.pgm import PGM
+
+
+class Cancer(PGM):
+    """
+    This PGM is the well known, pedagogical 'Cancer' Bayesian network.
+    """
+
+    def __init__(self):
+        super().__init__(self.__class__.__name__)
+
+        pollution = self.new_rv('pollution', ('low', 'high'))
+        smoker = self.new_rv('smoker', ('True', 'False'))
+        cancer = self.new_rv('cancer', ('True', 'False'))
+        xray = self.new_rv('xray', ('positive', 'negative'))
+        dyspnoea = self.new_rv('dyspnoea', ('True', 'False'))
+
+        pgm_factor0 = self.new_factor(pollution)
+        pgm_function_2511325233376 = pgm_factor0.set_cpt()
+        pgm_function_2511325233376.set_cpd((), (0.9, 0.1))
+        pgm_factor1 = self.new_factor(smoker)
+        pgm_function_2511322114176 = pgm_factor1.set_cpt()
+        pgm_function_2511322114176.set_cpd((), (0.3, 0.7))
+        pgm_factor2 = self.new_factor(cancer, pollution, smoker)
+        pgm_function_2511324995136 = pgm_factor2.set_cpt()
+        pgm_function_2511324995136.set_cpd((0, 0), (0.03, 0.97))
+        pgm_function_2511324995136.set_cpd((1, 0), (0.05, 0.95))
+        pgm_function_2511324995136.set_cpd((0, 1), (0.001, 0.999))
+        pgm_function_2511324995136.set_cpd((1, 1), (0.02, 0.98))
+        pgm_factor3 = self.new_factor(xray, cancer)
+        pgm_function_2511331552432 = pgm_factor3.set_cpt()
+        pgm_function_2511331552432.set_cpd((0,), (0.9, 0.1))
+        pgm_function_2511331552432.set_cpd((1,), (0.2, 0.8))
+        pgm_factor4 = self.new_factor(dyspnoea, cancer)
+        pgm_function_2511325139344 = pgm_factor4.set_cpt()
+        pgm_function_2511325139344.set_cpd((0,), (0.65, 0.35))
+        pgm_function_2511325139344.set_cpd((1,), (0.3, 0.7))

ck/example/chain.py

@@ -0,0 +1,38 @@
+import math
+import random as _random
+
+from ck.pgm import PGM
+
+
+class Chain(PGM):
+    """
+    This PGM is the 'Chain' factor graph.
+
+    The Chain factor graph consists of a chain of random variables, x0, x1, x2, ...
+    where adjacent random variables in the chain are connected by a binary factor.
+    If include_unaries then, also includes one unary factor per random variable.
+    """
+
+    def __init__(
+            self,
+            vars_per_chain: int = 240,
+            states_per_var: int = 35,
+            include_unaries: bool = True,
+            random_seed: int = 123456,
+    ):
+        params = (vars_per_chain, states_per_var, include_unaries)
+        super().__init__(f'{self.__class__.__name__}({",".join(str(param) for param in params)})')
+
+        scale = 1 + math.log2(vars_per_chain)
+        random_stream = _random.Random(random_seed).random
+        binary_iter = map(lambda x: x / scale, iter(random_stream, None))
+        unary_iter = iter(random_stream, None)
+
+        rvs = [self.new_rv(f'x{i}', states_per_var) for i in range(vars_per_chain)]
+
+        for i in range(1, len(rvs)):
+            self.new_factor(rvs[i], rvs[i - 1]).set_dense().set_iter(binary_iter)
+
+        if include_unaries:
+            for rv in rvs:
+                self.new_factor(rv).set_dense().set_iter(unary_iter)