#current_params = np.array ( [(mu_x, mu_z, sigma_x, sigma_z, rho), # params 1ère loi normale # (mu_x, mu_z, sigma_x, sigma_z, rho)] ) # params 2ème loi normale current_params = np.array([[ 3.28778309, 69.89705882, 1.13927121, 13.56996002, 0. ], [ 3.68778309, 71.89705882, 1.13927121, 13.56996002, 0. ]]) # current_weights = np.array ( [ pi_0, pi_1 ] ) current_weights = np.array ( [ 0.5, 0.5 ] ) T = Q_i ( data, current_params, current_weights )