我正在尝试使用正态分布的最大似然估计方法来估计参数。下面是我使用的代码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
然后我尝试使用 mle2():
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
我错过了什么??
谢谢您的帮助!