Я не уверен, что я делаю не так с этим кодом.Я пытаюсь реализовать алгоритм максимизации ожидания (скриншот псевдокода ниже), и я думаю, что я не минимизирую бета-версию правильно.Любые мысли о том, как я могу это исправить?
Это ссылка на псевдокод здесь> EM Algo
Ниже приведен код, который я написал для дублирования псевдокода *[007]
b_1 = np.argmin((np.sum(J_1[data_y] - (J_1[data_x])*b_1)**2)**.5)
b_2 = np.argmin((np.sum(J_2[data_y] - (J_2[data_x])*b_2)**2)**.5)
Я получаю эту ошибку в строках выше
only integer scalar arrays can be converted to a scalar index
В конечном счете, я хочу найти радиус вокруг истинных бета-версий того, как далеко от первоначального предположения Бета-версийможет быть из истинных бета-версий, где они будут сходиться к истинным бета-версиям за 500 итераций