0

我正在尝试使用正态分布的最大似然估计方法来估计参数。下面是我使用的代码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

我错过了什么??

谢谢您的帮助!

4

0 回答 0