Я делаю оценку максимального правдоподобия с использованием функции 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
?
Заранее спасибо.