通过数据框计算MLE,以获取不同的观察结果

时间:2018-10-19 21:24:20

标签: r normal-distribution mle

我正在使用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))

我使用前6次观察估计了mu的参数。我想为任何mu估计datl = dat[1:i,]的{​​{1}}。我该怎么办?我可以在i=2,3,3,..30中使用apply函数吗?

谢谢。

0 个答案:

没有答案