我正在使用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
函数吗?
谢谢。