wolffd@0: function engine = enter_evidence(engine, pos, neg)
wolffd@0: % ENTER_EVIDENCE Add evidence to the QMR network
wolffd@0: % engine = enter_evidence(engine, pos, neg)
wolffd@0: %
wolffd@0: % pos = list of leaves that have positive observations
wolffd@0: % neg = list of leaves that have negative observations
wolffd@0: 
wolffd@0: % Extract params for the observed findings
wolffd@0: obs = myunion(pos, neg);
wolffd@0: %inhibit_obs = engine.inhibit(obs, :);
wolffd@0: inhibit_obs = engine.inhibit(:,obs)';
wolffd@0: leak_obs = engine.leak(obs);
wolffd@0: 
wolffd@0: % Find what nodes correspond to the original observed leaves
wolffd@0: pos2 = find_equiv_posns(pos, obs);
wolffd@0: neg2 = find_equiv_posns(neg, obs);
wolffd@0: engine.post = quickscore(pos2, neg2, inhibit_obs, engine.prior, leak_obs); 
wolffd@0: %engine.post = C_quickscore(pos2, neg2, inhibit_obs, engine.prior, leak_obs); 
wolffd@0: