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 =