Предположим, у меня есть функция логарифмического правдоподобия из этой картинки: введите описание изображения здесь
Я хочу сделать MLE для v (t) = v1, v2. Если бы я сделал nll2, я получил бы ошибку:
> nll2<-function(v1,v2){
+ b=dbinom(mix(x2543b1),mix(x2543b),
+ exp(v1*mix(fee1)+v2*mix(fee2))/(1+exp(v1*mix(fee1)+v2*mix(fee2))),
+ log = T)
+ -sum(b)
+
+ }
> library(stats4)
> res2 = mle(minuslog=nll2,start=list(v1=0.2,v2=0.2),
+ method = "BFGS")
Error in optim(start, f, method = method, hessian = TRUE, ...) :
initial value in 'vmmin' is not finite
Я получил
"Error in optim(start, f, method = method, hessian = TRUE, ...) :
initial value in 'vmmin' is not finite"
Как я могу преодолеть это?