我正在尝试使用正态分布的最大似然估计方法来估计参数。以下是使用mle()
的代码:
x.norm<-rnorm(100,2,1)
library(stats4)
norm<-function(mu,sigma) {
n<-100
x<-x.norm
log(sigma)+(1/2)*log(2*pi)+((x-mu)**2)/(2*sigma**2)}
est<-mle(minuslog=norm, start=list(mu=1,sigma=1))
Error in optim(start, f, method = method, hessian = TRUE, ...) : objective function in optim evaluates to length 100 not 1
library(bbmle)
est<-mle2(minuslog=norm, start=list(mu=1,sigma=1))
Error in optim(par = c(1, 1), fn = function (p) : objective function in optim evaluates to length 100 not 1
我错过了什么?
感谢您的帮助!