使用R中的mle2()函数计算正态分布的参数估计

时间:2013-07-01 10:19:16

标签: r

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

我错过了什么?

感谢您的帮助!

0 个答案:

没有答案