Оценка максимального правдоподобия в R нечисловой ошибке аргумента

Question

Это моя функция правдоподобия для отрицательного биномиального распределения. Что я делаю не так, чтобы получить эту ошибку?

• Ошибка в k + z: нечисловой аргумент бинарного оператора

Код R

Answer 1

Вы установили z к вектору, но когда вы звоните optim ты говоришь z=data, data является встроенным объектом R (функция), так k+z в итоге пытается добавить функцию к числу. Я думаю, что вы хотели z=z в вашем optim вызов...