Object-Oriented Probabilistic Relational Model

import pyagrum.prm as gum
import pyagrum.lib.notebook as gnb

O3PRM is inspired by relational language, directly based on the BN and improved by object oriented paradigms, where the focus is set on classes of objects and by defining relations among these objects (see http://o3prm.gitlab.io/ for more details)

<img src=“/images/reference/o3prm.png”width=“350” alt=“PyAgrum inline image”>

O3PRM syntax

with open("res/Asia.o3prm", "r") as prm:
  for line in prm.readlines():
    print(line, end="")

class Asia {

    boolean visitToAsia {
        [ 0.99,  // False
          0.01 ] // True
    };

    boolean tuberculosis dependson visitToAsia {
        // False | True
        [    0.99, 0.95,  // False
             0.01, 0.05 ] // True
    };

    boolean smoking {
        [ 0.50,  // False
          0.50 ] // True
    };

    boolean lungCancer dependson smoking {
        // False | True => smoking
        [    0.99, 0.90,  // False
             0.01, 0.10 ] // True
    };

    boolean bronchitis dependson smoking {
        // False | True => smoking
        [    0.70, 0.40,  // False
             0.30, 0.60 ] // True
    };

    boolean tubOrCancer dependson tuberculosis, lungCancer {
        //      False    ||     True      => tuberculosis
        //  False | True || False | True  => lungCancer
        [     1.00, 0.00,     0.00, 0.00,  // False
              0.00, 1.00,     1.00, 1.00 ] // True
    };

    boolean positiveXRay dependson tubOrCancer {
        // False | True => tubOrCancer
        [ 0.95, 0.02,  // False
          0.05, 0.98 ] // True
    };

    boolean dyspnea dependson tubOrCancer, bronchitis {
        //      False    ||     True      => tubOrCancer
        //  False | True || False | True  => bronchitis
        [     0.90, 0.20,     0.30, 0.10,  // False
              0.10, 0.80,     0.70, 0.90 ] // True
    };
}

Using o3prm syntax for creating BayesNet

bn = gum.loadBN("res/Asia.o3prm", verbose=False)
bn

bn = gum.loadBN("res/aSys.o3prm")
bn

classpath = "res/ComplexPrinters"
filename = "res/ComplexPrinters/fr/lip6/printers/system.o3prm"

system = "Work"
bn = gum.loadBN(filename, system=system, classpath=classpath)

## the inference will take place in a rather large junction tree
gnb.showJunctionTreeMap(bn, scaleClique=0.1, scaleSep=0.05, lenEdge=1.2, size="8!")

gum.saveBN(bn, "out/bnprinters.bgum")

bn2 = gum.loadBN("out/bnprinters.bgum")
gnb.showInference(bn2, size="25!")

import pyagrum.lib.image as gimg

gimg.export(bn, filename="out/ComplexPrinters.pdf")

bn = gum.loadBN("res/aSys.o3prm")
gnb.sideBySide(gnb.getBN(bn, size="5"), gnb.getInference(bn, size="5"))

Exploring Probabilistic Relational Model

classpath = "res/ComplexPrinters"
filename = "res/ComplexPrinters/complexprinters_system.o3prm"
explor = gum.PRMexplorer()
explor.load(filename)

for cl in explor.classes():
  print("Class : " + cl)
  print("  - Super class : " + ("None" if explor.getSuperClass(cl) == None else explor.getSuperClass(cl)))
  print("  - Implemented interface : ")
  for inter in explor.classImplements(cl):
    print("      " + inter)
  print("  - Direct sub-types : ")
  for ext in explor.getDirectSubClass(cl):
    print("      " + ext)
  print("  - Attributes : ")
  for t, n, depensons in explor.classAttributes(cl):
    s = ""
    for depenson in depensons:
      s = s + depenson + " "
    print("      " + t + " " + n + " (" + s + ")")
  print("  - References : ")
  for t, n, isArray in explor.classReferences(cl):
    print("      " + t + ("[]" if isArray else "") + " " + n)
  print("  - Aggragates : ")
  for t, n, g, l, slots in explor.classAggregates(cl):
    s = ""
    for slot in slots:
      s = s + slot + " "
    print("      " + t + " " + n + " " + g + " " + ("NoLabel" if l == None else l) + " (" + s + ")")
  print("  - SlotChains : ")
  for t, n, isMultiple in explor.classSlotChains(cl):
    print("      " + t + " " + n + " " + ("[]" if isMultiple else ""))
  print("  - Parameters : ")
  for param in explor.classParameters(cl):
    print("      " + param)
  # print("  - Dag : ")
  # (dic, dotString) = explor.classDag(cl)
  # print(dic)
  # print(dotString)

  print()

Class : PowerSupply
  - Super class : None
  - Implemented interface :
  - Direct sub-types :
  - Attributes :
      t_state state ()
  - References :
  - Aggragates :
  - SlotChains :
  - Parameters :

Class : Room
  - Super class : None
  - Implemented interface :
  - Direct sub-types :
  - Attributes :
  - References :
      PowerSupply power
  - Aggragates :
  - SlotChains :
  - Parameters :

Class : ColorPrinter
  - Super class : None
  - Implemented interface :
      Printer
  - Direct sub-types :
  - Attributes :
      t_ink yellow ()
      t_degraded equipState (room.power.state hasPaper hasInk black )
      t_ink magenta ()
      t_paper hasPaper ()
      t_ink black ()
      t_ink cyan ()
  - References :
      Room room
  - Aggragates :
      boolean hasInk forall NotEmpty (yellow cyan black magenta )
  - SlotChains :
      t_state room.power.state
  - Parameters :

Class : SafeComputer
  - Super class : None
  - Implemented interface :
      Equipment
  - Direct sub-types :
  - Attributes :
      boolean can_print (working_printer equipState )
      t_degraded equipState (room.power.(t_state)state )
  - References :
      Printer[] printers
      Room room
  - Aggragates :
      boolean degraded_printer exists Degraded (printers.(t_degraded)equipState )
      boolean working_printer exists true (functional_printer degraded_printer )
      boolean functional_printer exists OK (printers.(t_state)equipState )
  - SlotChains :
      t_state room.power.(t_state)state
      t_state printers.(t_state)equipState []
      t_degraded printers.(t_degraded)equipState []
  - Parameters :

Class : ParamClass<lambda=0.001,t=4>
  - Super class : ParamClass
  - Implemented interface :
  - Direct sub-types :
  - Attributes :
      t_degraded equipState (room.power.state hasInk hasPaper )
      t_ink hasInk ()
      t_paper hasPaper ()
  - References :
      Room room
  - Aggragates :
  - SlotChains :
      t_state room.power.state
  - Parameters :
      lambda
      t

Class : BWPrinter
  - Super class : None
  - Implemented interface :
      Printer
  - Direct sub-types :
  - Attributes :
      t_paper hasPaper ()
      t_ink hasInk ()
      t_degraded equipState (room.power.state hasPaper hasInk )
  - References :
      Room room
  - Aggragates :
  - SlotChains :
      t_state room.power.state
  - Parameters :

Class : Computer
  - Super class : None
  - Implemented interface :
      Equipment
  - Direct sub-types :
  - Attributes :
      boolean can_print (equipState working_printer )
      t_degraded equipState (room.power.state )
  - References :
      Printer[] printers
      Room room
  - Aggragates :
      boolean working_printer exists true (degraded_printer functional_printer )
      mycount count_printers count true (degraded_printer functional_printer )
      boolean functional_printer exists OK (printers.equipState )
      boolean degraded_printer exists Degraded (printers.equipState )
  - SlotChains :
      t_state room.power.state
      t_degraded printers.equipState []
  - Parameters :

Class : ParamClass
  - Super class : None
  - Implemented interface :
  - Direct sub-types :
      ParamClass<lambda=0.001,t=4>
      ParamClass<lambda=0.4,t=4>
  - Attributes :
      t_degraded equipState (room.power.state hasInk hasPaper )
      t_paper hasPaper ()
      t_ink hasInk ()
  - References :
      Room room
  - Aggragates :
  - SlotChains :
      t_state room.power.state
  - Parameters :
      t
      lambda

Class : ParamClass<lambda=0.4,t=4>
  - Super class : ParamClass
  - Implemented interface :
  - Direct sub-types :
  - Attributes :
      t_ink hasInk ()
      t_paper hasPaper ()
      t_degraded equipState (room.power.state hasInk hasPaper )
  - References :
      Room room
  - Aggragates :
  - SlotChains :
      t_state room.power.state
  - Parameters :
      t
      lambda

print("The following lists the systems of the prm:\n")
systems = explor.getalltheSystems()
sys1 = systems[0]
print("Name of the system: " + sys1[0] + "\n")
print("Nodes : dict(id: [name,type])")
print(sys1[1])
print("\n")
print("Arcs : list[(tail, head),(tail, head)...]")
print(sys1[2])

The following lists the systems of the prm:

Name of the system: aSys

Nodes : dict(id: [name,type])
{0: ('pow', 'PowerSupply'), 1: ('r', 'Room'), 2: ('bw_printers[0]', 'BWPrinter'), 3: ('bw_printers[1]', 'BWPrinter'), 4: ('bw_printers[2]', 'BWPrinter'), 5: ('bw_printers[3]', 'BWPrinter'), 6: ('bw_printers[4]', 'BWPrinter'), 7: ('bw_printers[5]', 'BWPrinter'), 8: ('bw_printers[6]', 'BWPrinter'), 9: ('bw_printers[7]', 'BWPrinter'), 10: ('bw_printers[8]', 'BWPrinter'), 11: ('bw_printers[9]', 'BWPrinter'), 12: ('color_printers[0]', 'ColorPrinter'), 13: ('color_printers[1]', 'ColorPrinter'), 14: ('c1', 'Computer'), 15: ('c2', 'Computer'), 16: ('p', 'ParamClass<lambda=0.4,t=4>'), 17: ('paramBis', 'ParamClass<lambda=0.001,t=4>')}


Arcs : list[(tail, head),(tail, head)...]
[]

gnb.showTensor(explor.cpf("Computer", "equipState"))

	t_degraded-0x7835708500
room.power.state	OK	Dysfunctional	Degraded
OK	0.9000	0.0900	0.0100
NOK	0.0000	1.0000	0.0000

for cl in explor.types():
  print("Type : " + cl)
  print("  - Super type : " + ("None" if explor.getSuperType(cl) == None else explor.getSuperType(cl)))
  print("  - Direct sub-types : ")
  for name in explor.getDirectSubTypes(cl):
    print("      " + name)
  print("  - Labels : ")
  for t in explor.getLabels(cl):
    print("      " + t)
  print("  - Labels mapping : ")
  for key, val in dict().items() if explor.getLabelMap(cl) == None else explor.getLabelMap(cl).items():
    print("      " + key + " -> " + val)
  print()

Type : boolean
  - Super type : None
  - Direct sub-types :
      t_state
  - Labels :
      false
      true
  - Labels mapping :

Type : t_state
  - Super type : boolean-0x7835709d80
  - Direct sub-types :
      t_ink
      t_paper
      t_degraded
  - Labels :
      OK
      NOK
  - Labels mapping :
      OK -> true
      NOK -> false

Type : mycount
  - Super type : None
  - Direct sub-types :
  - Labels :
      0
      1
      2
      3
      4
      5
  - Labels mapping :

Type : t_ink
  - Super type : t_state-0x7839d00700
  - Direct sub-types :
  - Labels :
      NotEmpty
      Empty
  - Labels mapping :
      NotEmpty -> OK
      Empty -> NOK

Type : t_paper
  - Super type : t_state-0x7839d00700
  - Direct sub-types :
  - Labels :
      Ready
      Jammed
      Empty
  - Labels mapping :
      Ready -> OK
      Jammed -> NOK
      Empty -> NOK

Type : t_degraded
  - Super type : t_state-0x7839d00700
  - Direct sub-types :
  - Labels :
      OK
      Dysfunctional
      Degraded
  - Labels mapping :
      OK -> OK
      Dysfunctional -> NOK
      Degraded -> NOK

for cl in explor.interfaces():
  print("Interface : " + cl)
  print("  - Super interface : " + ("None" if explor.getSuperInterface(cl) == None else explor.getSuperInterface(cl)))
  print("  - Direct sub-interfaces : ")
  for name in explor.getDirectSubInterfaces(cl):
    print("      " + name)
  print("  - Implementations : ")
  for impl in explor.getImplementations(cl):
    print("      " + impl)
  print("  - Attributes : ")
  for t, n in explor.interAttributes(cl, allAttributes=True):
    print("      " + t + " " + n)
  print("  - References : ")
  for t, n, isArray in explor.interReferences(cl):
    print("      " + t + ("[]" if isArray else "") + " " + n)
  print()

Interface : Equipment
  - Super interface : None
  - Direct sub-interfaces :
      Printer
  - Implementations :
      SafeComputer
      Computer
  - Attributes :
      boolean equipState
      t_state equipState
      t_degraded equipState
  - References :
      Room room

Interface : Printer
  - Super interface : Equipment
  - Direct sub-interfaces :
  - Implementations :
      ColorPrinter
      BWPrinter
  - Attributes :
      boolean hasInk
      t_state equipState
      boolean equipState
      t_degraded equipState
      boolean hasPaper
  - References :
      Room room

print(explor.isType("fr.lip6.printers.base.Printer"))
print(explor.isClass("fr.lip6.printers.base.Printer"))
print(explor.isInterface("fr.lip6.printers.base.Printer"))

False
False
False