//predicate declarations
bloodtypeA(person)
bloodtypeB(person)
bloodtypeAB(person)
bloodtypeN(person)

pchromA(person)
pchromB(person)
pchromN(person)

mchromA(person)
mchromB(person)
mchromN(person)

father(person,person)
mother(person,person)
founder(person)

// hard formulas

pchromA(p1) => !(pchromB(p1) v pchromN(p1)).
pchromB(p1) => !(pchromA(p1) v pchromN(p1)).
pchromN(p1) => !(pchromA(p1) v pchromB(p1)).
pchromA(p1) v pchromB(p1) v pchromN(p1).

mchromA(p1) => !(mchromB(p1) v mchromN(p1)).
mchromB(p1) => !(mchromA(p1) v mchromN(p1)).
mchromN(p1) => !(mchromA(p1) v mchromB(p1)).
mchromA(p1) v mchromB(p1) v mchromN(p1).

bloodtypeA(p1) => !(bloodtypeB(p1) v bloodtypeAB(p1) v bloodtypeN(p1)).
bloodtypeB(p1) => !(bloodtypeA(p1) v bloodtypeAB(p1) v bloodtypeN(p1)).
bloodtypeAB(p1) => !(bloodtypeA(p1) v bloodtypeB(p1) v bloodtypeN(p1)).
bloodtypeN(p1) => !(bloodtypeA(p1) v bloodtypeB(p1) v bloodtypeAB(p1)).
bloodtypeA(p1) v bloodtypeB(p1) v bloodtypeAB(p1) v bloodtypeN(p1).

// weighted formulas


// P(bloodtype(p)=A|pchrom(p),mchrom(p))
// P(bloodtype(p)=A|pchrom(p)=A,mchrom(p)=A)
-0.105360516 pchromA(p) ^ mchromA(p) ^ bloodtypeA(p) //0.9
-2.302585093 pchromA(p) ^ mchromA(p) ^ !bloodtypeA(p) //0.1
// P(bloodtype(p)=A|pchrom(p)=A,mchrom(p)=B)
-3.506557897 pchromA(p) ^ mchromB(p) ^ bloodtypeA(p) //0.03
-0.030459207 pchromA(p) ^ mchromB(p) ^ !bloodtypeA(p) //0.97
// P(bloodtype(p)=A|pchrom(p)=A,mchrom(p)=N)
-0.105360516 pchromA(p) ^ mchromN(p) ^ bloodtypeA(p) //0.9
-2.302585093 pchromA(p) ^ mchromN(p) ^ !bloodtypeA(p) //0.1
// P(bloodtype(p)=A|pchrom(p)=B,mchrom(p)=A)
-3.506557897 pchromB(p) ^ mchromA(p) ^ bloodtypeA(p) //0.03
-0.030459207 pchromB(p) ^ mchromA(p) ^ !bloodtypeA(p) //0.97
// P(bloodtype(p)=A|pchrom(p)=B,mchrom(p)=B)
-3.218875825 pchromB(p) ^ mchromB(p) ^ bloodtypeA(p) //0.04
-0.040821995 pchromB(p) ^ mchromB(p) ^ !bloodtypeA(p) //0.96
// P(bloodtype(p)=A|pchrom(p)=B,mchrom(p)=N)
-3.912023005 pchromB(p) ^ mchromN(p) ^ bloodtypeA(p) //0.02
-0.020202707 pchromB(p) ^ mchromN(p) ^ !bloodtypeA(p) //0.98
// P(bloodtype(p)=A|pchrom(p)=N,mchrom(p)=A)
-0.105360516 pchromN(p) ^ mchromA(p) ^ bloodtypeA(p) //0.9
-2.302585093 pchromN(p) ^ mchromA(p) ^ !bloodtypeA(p) //0.1
// P(bloodtype(p)=A|pchrom(p)=N,mchrom(p)=B)
-3.912023005 pchromN(p) ^ mchromB(p) ^ bloodtypeA(p) //0.02
-0.020202707 pchromN(p) ^ mchromB(p) ^ !bloodtypeA(p) //0.98
// P(bloodtype(p)=A|pchrom(p)=N,mchrom(p)=N)
-3.912023005 pchromN(p) ^ mchromN(p) ^ bloodtypeA(p) //0.02
-0.020202707 pchromN(p) ^ mchromN(p) ^ !bloodtypeA(p) //0.98

// P(bloodtype(p)=B|pchrom(p),mchrom(p))
// P(bloodtype(p)=B|pchrom(p)=A,mchrom(p)=A)
-3.506557897 pchromA(p) ^ mchromA(p) ^ bloodtypeB(p) //0.03
-0.030459207 pchromA(p) ^ mchromA(p) ^ !bloodtypeB(p) //
// P(bloodtype(p)=B|pchrom(p)=A,mchrom(p)=B)
-3.506557897 pchromA(p) ^ mchromB(p) ^ bloodtypeB(p) //0.03
-0.030459207 pchromA(p) ^ mchromB(p) ^ !bloodtypeB(p) //
// P(bloodtype(p)=B|pchrom(p)=A,mchrom(p)=N)
-3.218875825 pchromA(p) ^ mchromN(p) ^ bloodtypeB(p) //0.04
-0.040821995 pchromA(p) ^ mchromN(p) ^ !bloodtypeB(p) //
// P(bloodtype(p)=B|pchrom(p)=B,mchrom(p)=A)
-3.506557897 pchromB(p) ^ mchromA(p) ^ bloodtypeB(p) //0.03
-0.030459207 pchromB(p) ^ mchromA(p) ^ !bloodtypeB(p) //
// P(bloodtype(p)=B|pchrom(p)=B,mchrom(p)=B)
-0.105360516 pchromB(p) ^ mchromB(p) ^ bloodtypeB(p) //0.9
-2.302585093 pchromB(p) ^ mchromB(p) ^ !bloodtypeB(p) //
// P(bloodtype(p)=B|pchrom(p)=B,mchrom(p)=N)
-0.105360516 pchromB(p) ^ mchromN(p) ^ bloodtypeB(p) //0.9
-2.302585093 pchromB(p) ^ mchromN(p) ^ !bloodtypeB(p) //
// P(bloodtype(p)=B|pchrom(p)=N,mchrom(p)=A)
-3.506557897 pchromN(p) ^ mchromA(p) ^ bloodtypeB(p) //0.03
-0.030459207 pchromN(p) ^ mchromA(p) ^ !bloodtypeB(p) //
// P(bloodtype(p)=B|pchrom(p)=N,mchrom(p)=B)
-0.105360516 pchromN(p) ^ mchromB(p) ^ bloodtypeB(p) //0.9
-2.302585093 pchromN(p) ^ mchromB(p) ^ !bloodtypeB(p) //
// P(bloodtype(p)=B|pchrom(p)=N,mchrom(p)=N)
-3.218875825 pchromN(p) ^ mchromN(p) ^ bloodtypeB(p) //0.04
-0.040821995 pchromN(p) ^ mchromN(p) ^ !bloodtypeB(p) //

// P(bloodtype(p)=AB|pchrom(p),mchrom(p))
// P(bloodtype(p)=AB|pchrom(p)=A,mchrom(p)=A)
-3.506557897 pchromA(p) ^ mchromA(p) ^ bloodtypeAB(p) //0.03
-0.030459207 pchromA(p) ^ mchromA(p) ^ !bloodtypeAB(p) //
// P(bloodtype(p)=AB|pchrom(p)=A,mchrom(p)=B)
-0.105360516 pchromA(p) ^ mchromB(p) ^ bloodtypeAB(p) //0.9
-2.302585093 pchromA(p) ^ mchromB(p) ^ !bloodtypeAB(p) //
// P(bloodtype(p)=AB|pchrom(p)=A,mchrom(p)=N)
-3.506557897 pchromA(p) ^ mchromN(p) ^ bloodtypeAB(p) //0.03
-0.030459207 pchromA(p) ^ mchromN(p) ^ !bloodtypeAB(p) //
// P(bloodtype(p)=AB|pchrom(p)=B,mchrom(p)=A)
-0.105360516 pchromB(p) ^ mchromA(p) ^ bloodtypeAB(p) //0.9
-2.302585093 pchromB(p) ^ mchromA(p) ^ !bloodtypeAB(p) //
// P(bloodtype(p)=AB|pchrom(p)=B,mchrom(p)=B)
-3.506557897 pchromB(p) ^ mchromB(p) ^ bloodtypeAB(p) //0.03
-0.030459207 pchromB(p) ^ mchromB(p) ^ !bloodtypeAB(p) //
// P(bloodtype(p)=AB|pchrom(p)=B,mchrom(p)=N)
-3.218875825 pchromB(p) ^ mchromN(p) ^ bloodtypeAB(p) //0.04
-0.040821995 pchromB(p) ^ mchromN(p) ^ !bloodtypeAB(p) //
// P(bloodtype(p)=AB|pchrom(p)=N,mchrom(p)=A)
-3.506557897 pchromN(p) ^ mchromA(p) ^ bloodtypeAB(p) //0.03
-0.030459207 pchromN(p) ^ mchromA(p) ^ !bloodtypeAB(p) //
// P(bloodtype(p)=AB|pchrom(p)=N,mchrom(p)=B)
-3.218875825 pchromN(p) ^ mchromB(p) ^ bloodtypeAB(p) //0.04
-0.040821995 pchromN(p) ^ mchromB(p) ^ !bloodtypeAB(p) //
// P(bloodtype(p)=AB|pchrom(p)=N,mchrom(p)=N)
-3.218875825 pchromN(p) ^ mchromN(p) ^ bloodtypeAB(p) //0.04
-0.040821995 pchromN(p) ^ mchromN(p) ^ !bloodtypeAB(p) //

// P(bloodtype(p)=N|pchrom(p),mchrom(p))
// P(bloodtype(p)=N|pchrom(p)=A,mchrom(p)=A)
-3.218875825 pchromA(p) ^ mchromA(p) ^ bloodtypeN(p) //0.04
-0.040821995 pchromA(p) ^ mchromA(p) ^ !bloodtypeN(p) //
// P(bloodtype(p)=N|pchrom(p)=A,mchrom(p)=B)
-3.218875825 pchromA(p) ^ mchromB(p) ^ bloodtypeN(p) //0.04
-0.040821995 pchromA(p) ^ mchromB(p) ^ !bloodtypeN(p) //
// P(bloodtype(p)=N|pchrom(p)=A,mchrom(p)=N)
-3.506557897 pchromA(p) ^ mchromN(p) ^ bloodtypeN(p) //0.03
-0.030459207 pchromA(p) ^ mchromN(p) ^ !bloodtypeN(p) //
// P(bloodtype(p)=N|pchrom(p)=B,mchrom(p)=A)
-3.218875825 pchromB(p) ^ mchromA(p) ^ bloodtypeN(p) //0.04
-0.040821995 pchromB(p) ^ mchromA(p) ^ !bloodtypeN(p) //
// P(bloodtype(p)=N|pchrom(p)=B,mchrom(p)=B)
-3.506557897 pchromB(p) ^ mchromB(p) ^ bloodtypeN(p) //0.03
-0.030459207 pchromB(p) ^ mchromB(p) ^ !bloodtypeN(p) //
// P(bloodtype(p)=N|pchrom(p)=B,mchrom(p)=N)
-3.218875825 pchromB(p) ^ mchromN(p) ^ bloodtypeN(p) //0.04
-0.040821995 pchromB(p) ^ mchromN(p) ^ !bloodtypeN(p) //
// P(bloodtype(p)=N|pchrom(p)=N,mchrom(p)=A)
-3.218875825 pchromN(p) ^ mchromA(p) ^ bloodtypeN(p) //0.04
-0.040821995 pchromN(p) ^ mchromA(p) ^ !bloodtypeN(p) //
// P(bloodtype(p)=N|pchrom(p)=N,mchrom(p)=B)
-3.218875825 pchromN(p) ^ mchromB(p) ^ bloodtypeN(p) //0.04
-0.040821995 pchromN(p) ^ mchromB(p) ^ !bloodtypeN(p) //
// P(bloodtype(p)=N|pchrom(p)=N,mchrom(p)=N)
-0.105360516 pchromN(p) ^ mchromN(p) ^ bloodtypeN(p) //0.9
-2.302585093 pchromN(p) ^ mchromN(p) ^ !bloodtypeN(p) //

// P(pchrom(p)=A|father(f,p),pchrom(f),mchrom(f))
-0.105360516 father(f,p) ^ pchromA(f) ^ mchromA(f) ^ pchromA(p) //0.9
-2.302585093 father(f,p) ^ pchromA(f) ^ mchromA(f) ^ !pchromA(p) //
-0.713349888 father(f,p) ^ pchromA(f) ^ mchromB(f) ^ pchromA(p) //0.49
-0.673344553 father(f,p) ^ pchromA(f) ^ mchromB(f) ^ !pchromA(p) //
-0.713349888 father(f,p) ^ pchromA(f) ^ mchromN(f) ^ pchromA(p) //0.49
-0.673344553 father(f,p) ^ pchromA(f) ^ mchromN(f) ^ !pchromA(p) //
-0.713349888 father(f,p) ^ pchromB(f) ^ mchromA(f) ^ pchromA(p) //0.49
-0.673344553 father(f,p) ^ pchromB(f) ^ mchromA(f) ^ !pchromA(p) //
-2.995732274 father(f,p) ^ pchromB(f) ^ mchromB(f) ^ pchromA(p) //0.05
-0.051293294 father(f,p) ^ pchromB(f) ^ mchromB(f) ^ !pchromA(p) //
-3.912023005 father(f,p) ^ pchromB(f) ^ mchromN(f) ^ pchromA(p) //0.02
-0.020202707 father(f,p) ^ pchromB(f) ^ mchromN(f) ^ !pchromA(p) //
-0.713349888 father(f,p) ^ pchromN(f) ^ mchromA(f) ^ pchromA(p) //0.49
-0.673344553 father(f,p) ^ pchromN(f) ^ mchromA(f) ^ !pchromA(p) //
-3.912023005 father(f,p) ^ pchromN(f) ^ mchromB(f) ^ pchromA(p) //0.02
-0.020202707 father(f,p) ^ pchromN(f) ^ mchromB(f) ^ !pchromA(p) //
-2.995732274 father(f,p) ^ pchromN(f) ^ mchromN(f) ^ pchromA(p) //0.05
-0.051293294 father(f,p) ^ pchromN(f) ^ mchromN(f) ^ !pchromA(p) //


// P(pchrom(p)=B|father(f,p),pchrom(f),mchrom(f))
-2.995732274 father(f,p) ^ pchromA(f) ^ mchromA(f) ^ pchromB(p) //0.05
-0.051293294 father(f,p) ^ pchromA(f) ^ mchromA(f) ^ !pchromB(p) //
-0.713349888 father(f,p) ^ pchromA(f) ^ mchromB(f) ^ pchromB(p) //0.49
-0.673344553 father(f,p) ^ pchromA(f) ^ mchromB(f) ^ !pchromB(p) //
-3.912023005 father(f,p) ^ pchromA(f) ^ mchromN(f) ^ pchromB(p) //0.02
-0.020202707 father(f,p) ^ pchromA(f) ^ mchromN(f) ^ !pchromB(p) //
-0.713349888 father(f,p) ^ pchromB(f) ^ mchromA(f) ^ pchromB(p) //0.49
-0.673344553 father(f,p) ^ pchromB(f) ^ mchromA(f) ^ !pchromB(p) //
-0.105360516 father(f,p) ^ pchromB(f) ^ mchromB(f) ^ pchromB(p) //0.9
-2.302585093 father(f,p) ^ pchromB(f) ^ mchromB(f) ^ !pchromB(p) //
-0.713349888 father(f,p) ^ pchromB(f) ^ mchromN(f) ^ pchromB(p) //0.49
-0.673344553 father(f,p) ^ pchromB(f) ^ mchromN(f) ^ !pchromB(p) //
-3.912023005 father(f,p) ^ pchromN(f) ^ mchromA(f) ^ pchromB(p) //0.02
-0.020202707 father(f,p) ^ pchromN(f) ^ mchromA(f) ^ !pchromB(p) //
-0.713349888 father(f,p) ^ pchromN(f) ^ mchromB(f) ^ pchromB(p) //0.49
-0.673344553 father(f,p) ^ pchromN(f) ^ mchromB(f) ^ !pchromB(p) //
-2.995732274 father(f,p) ^ pchromN(f) ^ mchromN(f) ^ pchromB(p) //0.05
-0.051293294 father(f,p) ^ pchromN(f) ^ mchromN(f) ^ !pchromB(p) //

// P(pchrom(p)=N|father(f,p),pchrom(f),mchrom(f))
-2.995732274 father(f,p) ^ pchromA(f) ^ mchromA(f) ^ pchromN(p) //0.05
-0.051293294 father(f,p) ^ pchromA(f) ^ mchromA(f) ^ !pchromN(p) //
-3.912023005 father(f,p) ^ pchromA(f) ^ mchromB(f) ^ pchromN(p) //0.02
-0.020202707 father(f,p) ^ pchromA(f) ^ mchromB(f) ^ !pchromN(p) //
-0.713349888 father(f,p) ^ pchromA(f) ^ mchromN(f) ^ pchromN(p) //0.49
-0.673344553 father(f,p) ^ pchromA(f) ^ mchromN(f) ^ !pchromN(p) //
-3.912023005 father(f,p) ^ pchromB(f) ^ mchromA(f) ^ pchromN(p) //0.02
-0.020202707 father(f,p) ^ pchromB(f) ^ mchromA(f) ^ !pchromN(p) //
-2.995732274 father(f,p) ^ pchromB(f) ^ mchromB(f) ^ pchromN(p) //0.05
-0.051293294 father(f,p) ^ pchromB(f) ^ mchromB(f) ^ !pchromN(p) //
-0.713349888 father(f,p) ^ pchromB(f) ^ mchromN(f) ^ pchromN(p) //0.49
-0.673344553 father(f,p) ^ pchromB(f) ^ mchromN(f) ^ !pchromN(p) //
-0.713349888 father(f,p) ^ pchromN(f) ^ mchromA(f) ^ pchromN(p) //0.49
-0.673344553 father(f,p) ^ pchromN(f) ^ mchromA(f) ^ !pchromN(p) //
-0.713349888 father(f,p) ^ pchromN(f) ^ mchromB(f) ^ pchromN(p) //0.49
-0.673344553 father(f,p) ^ pchromN(f) ^ mchromB(f) ^ !pchromN(p) //
-0.105360516 father(f,p) ^ pchromN(f) ^ mchromN(f) ^ pchromN(p) //0.9
-2.302585093 father(f,p) ^ pchromN(f) ^ mchromN(f) ^ !pchromN(p) //

// P(mchrom(p)=A|mother(f,p),pchrom(f),mchrom(f))
-0.105360516 mother(f,p) ^ pchromA(f) ^ mchromA(f) ^ mchromA(p) //0.9
-2.302585093 mother(f,p) ^ pchromA(f) ^ mchromA(f) ^ !mchromA(p) //
-0.713349888 mother(f,p) ^ pchromA(f) ^ mchromB(f) ^ mchromA(p) //0.49
-0.673344553 mother(f,p) ^ pchromA(f) ^ mchromB(f) ^ !mchromA(p) //
-0.713349888 mother(f,p) ^ pchromA(f) ^ mchromN(f) ^ mchromA(p) //0.49
-0.673344553 mother(f,p) ^ pchromA(f) ^ mchromN(f) ^ !mchromA(p) //
-0.713349888 mother(f,p) ^ pchromB(f) ^ mchromA(f) ^ mchromA(p) //0.49
-0.673344553 mother(f,p) ^ pchromB(f) ^ mchromA(f) ^ !mchromA(p) //
-2.995732274 mother(f,p) ^ pchromB(f) ^ mchromB(f) ^ mchromA(p) //0.05
-0.051293294 mother(f,p) ^ pchromB(f) ^ mchromB(f) ^ !mchromA(p) //
-3.912023005 mother(f,p) ^ pchromB(f) ^ mchromN(f) ^ mchromA(p) //0.02
-0.020202707 mother(f,p) ^ pchromB(f) ^ mchromN(f) ^ !mchromA(p) //
-0.713349888 mother(f,p) ^ pchromN(f) ^ mchromA(f) ^ mchromA(p) //0.49
-0.673344553 mother(f,p) ^ pchromN(f) ^ mchromA(f) ^ !mchromA(p) //
-3.912023005 mother(f,p) ^ pchromN(f) ^ mchromB(f) ^ mchromA(p) //0.02
-0.020202707 mother(f,p) ^ pchromN(f) ^ mchromB(f) ^ !mchromA(p) //
-2.995732274 mother(f,p) ^ pchromN(f) ^ mchromN(f) ^ mchromA(p) //0.05
-0.051293294 mother(f,p) ^ pchromN(f) ^ mchromN(f) ^ !mchromA(p) //

// P(mchrom(p)=B|mother(f,p),pchrom(f),mchrom(f))
-2.995732274 mother(f,p) ^ pchromA(f) ^ mchromA(f) ^ mchromB(p) //0.05
-0.051293294 mother(f,p) ^ pchromA(f) ^ mchromA(f) ^ !mchromB(p) //
-0.713349888 mother(f,p) ^ pchromA(f) ^ mchromB(f) ^ mchromB(p) //0.49
-0.673344553 mother(f,p) ^ pchromA(f) ^ mchromB(f) ^ !mchromB(p) //
-3.912023005 mother(f,p) ^ pchromA(f) ^ mchromN(f) ^ mchromB(p) //0.02
-0.020202707 mother(f,p) ^ pchromA(f) ^ mchromN(f) ^ !mchromB(p) //
-0.713349888 mother(f,p) ^ pchromB(f) ^ mchromA(f) ^ mchromB(p) //0.49
-0.673344553 mother(f,p) ^ pchromB(f) ^ mchromA(f) ^ !mchromB(p) //
-0.105360516 mother(f,p) ^ pchromB(f) ^ mchromB(f) ^ mchromB(p) //0.9
-2.302585093 mother(f,p) ^ pchromB(f) ^ mchromB(f) ^ !mchromB(p) //
-0.713349888 mother(f,p) ^ pchromB(f) ^ mchromN(f) ^ mchromB(p) //0.49
-0.673344553 mother(f,p) ^ pchromB(f) ^ mchromN(f) ^ !mchromB(p) //
-3.912023005 mother(f,p) ^ pchromN(f) ^ mchromA(f) ^ mchromB(p) //0.02
-0.020202707 mother(f,p) ^ pchromN(f) ^ mchromA(f) ^ !mchromB(p) //
-0.713349888 mother(f,p) ^ pchromN(f) ^ mchromB(f) ^ mchromB(p) //0.49
-0.673344553 mother(f,p) ^ pchromN(f) ^ mchromB(f) ^ !mchromB(p) //
-2.995732274 mother(f,p) ^ pchromN(f) ^ mchromN(f) ^ mchromB(p) //0.05
-0.051293294 mother(f,p) ^ pchromN(f) ^ mchromN(f) ^ !mchromB(p) //

// P(mchrom(p)=N|mother(f,p),pchrom(f),mchrom(f))
-2.995732274 mother(f,p) ^ pchromA(f) ^ mchromA(f) ^ mchromN(p) //0.05
-0.051293294 mother(f,p) ^ pchromA(f) ^ mchromA(f) ^ !mchromN(p) //
-3.912023005 mother(f,p) ^ pchromA(f) ^ mchromB(f) ^ mchromN(p) //0.02
-0.020202707 mother(f,p) ^ pchromA(f) ^ mchromB(f) ^ !mchromN(p) //
-0.713349888 mother(f,p) ^ pchromA(f) ^ mchromN(f) ^ mchromN(p) //0.49
-0.673344553 mother(f,p) ^ pchromA(f) ^ mchromN(f) ^ !mchromN(p) //
-3.912023005 mother(f,p) ^ pchromB(f) ^ mchromA(f) ^ mchromN(p) //0.02
-0.020202707 mother(f,p) ^ pchromB(f) ^ mchromA(f) ^ !mchromN(p) //
-2.995732274 mother(f,p) ^ pchromB(f) ^ mchromB(f) ^ mchromN(p) //0.05
-0.051293294 mother(f,p) ^ pchromB(f) ^ mchromB(f) ^ !mchromN(p) //
-0.713349888 mother(f,p) ^ pchromB(f) ^ mchromN(f) ^ mchromN(p) //0.49
-0.673344553 mother(f,p) ^ pchromB(f) ^ mchromN(f) ^ !mchromN(p) //
-0.713349888 mother(f,p) ^ pchromN(f) ^ mchromA(f) ^ mchromN(p) //0.49
-0.673344553 mother(f,p) ^ pchromN(f) ^ mchromA(f) ^ !mchromN(p) //
-0.713349888 mother(f,p) ^ pchromN(f) ^ mchromB(f) ^ mchromN(p) //0.49
-0.673344553 mother(f,p) ^ pchromN(f) ^ mchromB(f) ^ !mchromN(p) //
-0.105360516 mother(f,p) ^ pchromN(f) ^ mchromN(f) ^ mchromN(p) //0.9
-2.302585093 mother(f,p) ^ pchromN(f) ^ mchromN(f) ^ !mchromN(p) //

// P(pchrom(p)=A|founder(p))
-1.203972804 founder(p) ^ pchromA(p)
-0.356674944 founder(p) ^ !pchromA(p)
// P(pchrom(p)=B|founder(p))
-1.203972804 founder(p) ^ pchromB(p)
-0.356674944 founder(p) ^ !pchromB(p)
// P(pchrom(p)=N|founder(p))
-0.916290732 founder(p) ^ pchromN(p)
-0.510825624 founder(p) ^ !pchromN(p)

// P(mchrom(p)=A|founder(p))
-1.203972804 founder(p) ^ mchromA(p)
-0.356674944 founder(p) ^ !mchromA(p)
// P(mchrom(p)=B|founder(p))
-1.203972804 founder(p) ^ mchromB(p)
-0.356674944 founder(p) ^ !mchromB(p)
// P(mchrom(p)=N|founder(p))
-0.916290732 founder(p) ^ mchromN(p)
-0.510825624 founder(p) ^ !mchromN(p)