Если вы говорите о функции mle()
из пакета stats4
, то, насколько я знаю, у вас не может быть векторных параметров.Что вы можете сделать, это обернуть вашу функцию LL()
в другую функцию, которая передается в mle()
.Например,
LL1 <- function(position1, position2, position3, position4,
position5, position6, position7, position8,
weight1, weight2,
beta1, beta2,
prob_s1, prob_s2)
LL(position = c(position1, position2, position3, position4,
position5, position6, position7, position8),
weight = c(weight1, weight2),
beta = c(beta1, beta2),
prob_s = c(prob_s1, prob_s2))
fit <- mle(LL1, start = unlist(inits))