djc.net

数据挖掘中的一算法 ines算法 c下实现的。适合初学习数据挖掘者借鉴

/*----------------------------------------------------------------------
  File    : djc.net
  Contents: Danish Jersey Cattle Blood Group Determination Network
            (reduced version, i.e., without parental_error node)
  Source  : Lene Kolind Rasmussen
----------------------------------------------------------------------*/

/*----------------------------------------------------------------------
  domains
----------------------------------------------------------------------*/
dom(dam_correct)        = { y, n };
dom(sire_correct)       = { y, n };
dom(true_dam_pg1)       = { f1, v1, v2 };
dom(true_dam_pg2)       = { f1, v1, v2 };
dom(true_sire_pg1)      = { f1, v1, v2 };
dom(true_sire_pg2)      = { f1, v1, v2 };
dom(stated_dam_pg1)     = { f1, v1, v2 };
dom(stated_dam_pg2)     = { f1, v1, v2 };
dom(stated_sire_pg1)    = { f1, v1, v2 };
dom(stated_sire_pg2)    = { f1, v1, v2 };
dom(offspring_pg1)      = { f1, v1, v2 };
dom(offspring_pg2)      = { f1, v1, v2 };
dom(offspring_genotype) = { f1f1, f1v1, f1v2, v1v1, v1v2, v2v2 };
dom(factor_40)          = { y, n };
dom(factor_41)          = { y, n };
dom(factor_42)          = { y, n };
dom(factor_43)          = { y, n };
dom(lysis_40)           = { 0, 1, 2, 3, 4, 5, 6, 7 };
dom(lysis_41)           = { 0, 1, 2, 3, 4, 5, 6, 7 };
dom(lysis_42)           = { 0, 1, 2, 3, 4, 5, 6, 7 };
dom(lysis_43)           = { 0, 1, 2, 3, 4, 5, 6, 7 };

/*----------------------------------------------------------------------
  graphical model
----------------------------------------------------------------------*/
probnet = {
  prob(dam_correct) =
    { y: 99, n: 1 };

  prob(sire_correct) =
    { y: 97, n: 3 };

  prob(true_dam_pg1) =
    { f1: 57.4, v1: 10.2, v2: 32.3 };

  prob(true_dam_pg2) =
    { f1: 57.4, v1: 10.2, v2: 32.3 };

  prob(true_sire_pg1) =
    { f1: 57.4, v1: 10.2, v2: 32.3 };

  prob(true_sire_pg2) =
    { f1: 57.4, v1: 10.2, v2: 32.3 };

  prob(stated_dam_pg1|dam_correct,true_dam_pg1) =
    { y: { f1: { f1:  1,   v1:  0,   v2:  0   },
           v1: { f1:  0,   v1:  1,   v2:  0   },
           v2: { f1:  0,   v1:  0,   v2:  1   } },
      n: { f1: { f1: 57.4, v1: 10.2, v2: 32.3 },
           v1: { f1: 57.4, v1: 10.2, v2: 32.3 },
           v2: { f1: 57.4, v1: 10.2, v2: 32.3 } } };

  prob(stated_dam_pg2|dam_correct,true_dam_pg2) =
    { y: { f1: { f1:  1,   v1:  0,   v2:  0   },
           v1: { f1:  0,   v1:  1,   v2:  0   },
           v2: { f1:  0,   v1:  0,   v2:  1   } },
      n: { f1: { f1: 57.4, v1: 10.2, v2: 32.3 },
           v1: { f1: 57.4, v1: 10.2, v2: 32.3 },
           v2: { f1: 57.4, v1: 10.2, v2: 32.3 } } };

  prob(stated_sire_pg1|sire_correct,true_sire_pg1) =
    { y: { f1: { f1:  1,   v1:  0,   v2:  0   },
           v1: { f1:  0,   v1:  1,   v2:  0   },
           v2: { f1:  0,   v1:  0,   v2:  1   } },
      n: { f1: { f1: 57.4, v1: 10.2, v2: 32.3 },
           v1: { f1: 57.4, v1: 10.2, v2: 32.3 },
           v2: { f1: 57.4, v1: 10.2, v2: 32.3 } } };

  prob(stated_sire_pg2|sire_correct,true_sire_pg2) =
    { y: { f1: { f1:  1,   v1:  0,   v2:  0   },
           v1: { f1:  0,   v1:  1,   v2:  0   },
           v2: { f1:  0,   v1:  0,   v2:  1   } },
      n: { f1: { f1: 57.4, v1: 10.2, v2: 32.3 },
           v1: { f1: 57.4, v1: 10.2, v2: 32.3 },
           v2: { f1: 57.4, v1: 10.2, v2: 32.3 } } };

  prob(offspring_pg1|true_dam_pg1,true_dam_pg2) =
    { f1: { f1: { f1: 1, v1: 0, v2: 0 },
            v1: { f1: 1, v1: 1, v2: 0 },
            v2: { f1: 1, v1: 0, v2: 1 } },
      v1: { f1: { f1: 1, v1: 1, v2: 0 },
            v1: { f1: 0, v1: 1, v2: 0 },
            v2: { f1: 0, v1: 1, v2: 1 } },
      v2: { f1: { f1: 1, v1: 0, v2: 1 },
            v1: { f1: 0, v1: 1, v2: 1 },
            v2: { f1: 0, v1: 0, v2: 1 } } };

  prob(offspring_pg2|true_sire_pg1,true_sire_pg2) =
    { f1: { f1: { f1: 1, v1: 0, v2: 0 },
            v1: { f1: 1, v1: 1, v2: 0 },
            v2: { f1: 1, v1: 0, v2: 1 } },
      v1: { f1: { f1: 1, v1: 1, v2: 0 },
            v1: { f1: 0, v1: 1, v2: 0 },
            v2: { f1: 0, v1: 1, v2: 1 } },
      v2: { f1: { f1: 1, v1: 0, v2: 1 },
            v1: { f1: 0, v1: 1, v2: 1 },
            v2: { f1: 0, v1: 0, v2: 1 } } };

  prob(offspring_genotype|offspring_pg1,offspring_pg2) =
    { f1: { f1: { f1f1: 1, f1v1: 0, f1v2: 0, v1v1: 0, v1v2: 0, v2v2: 0 },
            v1: { f1f1: 0, f1v1: 1, f1v2: 0, v1v1: 0, v1v2: 0, v2v2: 0 },
            v2: { f1f1: 0, f1v1: 0, f1v2: 1, v1v1: 0, v1v2: 0, v2v2: 0 } },
      v1: { f1: { f1f1: 0, f1v1: 1, f1v2: 0, v1v1: 0, v1v2: 0, v2v2: 0 },
            v1: { f1f1: 0, f1v1: 0, f1v2: 0, v1v1: 1, v1v2: 0, v2v2: 0 },
            v2: { f1f1: 0, f1v1: 0, f1v2: 0, v1v1: 0, v1v2: 1, v2v2: 0 } },
      v2: { f1: { f1f1: 0, f1v1: 0, f1v2: 1, v1v1: 0, v1v2: 0, v2v2: 0 },
            v1: { f1f1: 0, f1v1: 0, f1v2: 0, v1v1: 0, v1v2: 1, v2v2: 0 },
            v2: { f1f1: 0, f1v1: 0, f1v2: 0, v1v1: 0, v1v2: 0, v2v2: 1 } } };

  prob(factor_40|offspring_genotype) =
    { f1f1: { y: 1, n: 0 },
      f1v1: { y: 1, n: 0 },
      f1v2: { y: 1, n: 0 },
      v1v1: { y: 0, n: 1 },
      v1v2: { y: 0, n: 1 },
      v2v2: { y: 0, n: 1 } };

  prob(factor_41|offspring_genotype) =
    { f1f1: { y: 1, n: 0 },
      f1v1: { y: 1, n: 0 },
      f1v2: { y: 1, n: 0 },
      v1v1: { y: 0, n: 1 },
      v1v2: { y: 1, n: 0 },
      v2v2: { y: 1, n: 0 } };

  prob(factor_42|offspring_genotype) =
    { f1f1: { y: 0, n: 1 },
      f1v1: { y: 1, n: 0 },
      f1v2: { y: 0, n: 1 },
      v1v1: { y: 1, n: 0 },
      v1v2: { y: 1, n: 0 },
      v2v2: { y: 0, n: 1 } };

  prob(factor_43|offspring_genotype) =
    { f1f1: { y: 0, n: 1 },
      f1v1: { y: 1, n: 0 },
      f1v2: { y: 1, n: 0 },
      v1v1: { y: 1, n: 0 },
      v1v2: { y: 1, n: 0 },
      v2v2: { y: 1, n: 0 } };

  prob(lysis_40|factor_40) =
    { y: { 0:  0.1, 1:  1, 2:  3, 3: 4, 4: 8, 5: 14, 6: 18, 7: 20   },
      n: { 0: 20,   1: 18, 2: 14, 3: 8, 4: 4, 5:  3, 6:  1, 7:  0.1 } };

  prob(lysis_41|factor_41) =
    { y: { 0:  0.1, 1:  1, 2:  3, 3: 4, 4: 8, 5: 14, 6: 18, 7: 20   },
      n: { 0: 20,   1: 18, 2: 14, 3: 8, 4: 4, 5:  3, 6:  1, 7:  0.1 } };

  prob(lysis_42|factor_42) =
    { y: { 0:  0.1, 1:  1, 2:  3, 3: 4, 4: 8, 5: 14, 6: 18, 7: 20   },
      n: { 0: 20,   1: 18, 2: 14, 3: 8, 4: 4, 5:  3, 6:  1, 7:  0.1 } };

  prob(lysis_43|factor_43) =
    { y: { 0:  0.1, 1:  1, 2:  3, 3: 4, 4: 8, 5: 14, 6: 18, 7: 20   },
      n: { 0: 20,   1: 18, 2: 14, 3: 8, 4: 4, 5:  3, 6:  1, 7:  0.1 } };
};