EM algorithm

[yupinglin0301]

Hi I was curious about the EM algorithm that used in the software, would it be possible to tell us the design of the algorithm based on which equation or reference? Thanks.

Best, Yu-Ping

try: ### EM algorithm for idx_loop in range(0, 10000): ### E(num_AB|prob_AB) = 2 * num_AABB + num_AABb + num_AaBB + ### (prob_AB * (1 + prob_AB - prob_A - prob_B) * num_AbBb) / ### ((prob_A - prob_AB) * (prob_B - prob_AB) + prob_AB * (1 + prob_AB - prob_A - prob_B)) float_num_AB_estimateByEM = 2 * float(np_contigency[0, 0]) + float(np_contigency[0, 1]) + float(np_contigency[1, 0]) + (float_probability_AB * (1 + float_probability_AB - float_probability_A - float_probability_B) * float(np_contigency[1, 1])) / ((float_probability_A - float_probability_AB) * (float_probability_B - float_probability_AB) + float_probability_AB * (1 + float_probability_AB - float_probability_A - float_probability_B)) float_probability_AB_estimateByEM = float_num_AB_estimateByEM / (int_num_subject * 2) if abs(float_probability_AB_estimateByEM - float_probability_AB) < 0.0000001: break else: float_probability_AB =