AI, Python, Heredity, Joint Probability, Bayesian Networks.
one_gene in powerset(names): for two_genes in powerset(names - one_gene): # Update probabilities with new joint probability p = joint_probability(people, one_gene, two_genes, have_trait) update(probabilities, one_gene, two_genes, have_trait, p) # Ensure probabilities sum to 1