Вычислить MLE через фрейм данных для другого сегмента наблюдений - PullRequest
0 голосов
/ 20 октября 2018

Я делаю оценку максимального правдоподобия с использованием функции mle в r.Вот мой r код.

library(stats4)
set.seed(650)
y = c(rnorm(15, 1, 1), rnorm(15,3,1))
dat = data.frame(y)
LLl <- function(mu) {
  datl = dat[1:6,]
  R = dnorm(datl , mu, 1)
  -sum(log(R))

}
mle(LLl, start = list(mu=1))

Я оценил параметр для mu для использования первых 6 наблюдений.Я хочу оценить mu для datl = dat[1:i,] для любого i=2,3,3,..30.Как я могу сделать это?Могу ли я использовать функцию apply в r?

Заранее спасибо.

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...