对于R的新手来说这很简单但很困难。如何在没有mu和sigma的情况下获取 mu 和 sigma 的值。我希望以矩阵,表格或矢量形式的估计值,如下标题我希望的
LL <- function(mu, sigma) {
R = suppressWarnings(dnorm(x, mu, sigma))
-sum(log(R))
}
library(stats4)
mle(LL, start = list(mu = 1, sigma=1))@coef
mu sigma
2.550000 1.417393
# I wish
1 2
[1,] 1.1871891 4.2896111