试图找到几何分布的MLE时出错

Question

因此，我在R中被赋予了一个任务，就是用向量中给定的值来寻找具有几何分布的随机值X的MLE：

> x<-c(4, 2, 1, 1, 4, 8, 3, 1, 2, 3, 8, 6, 2, 2, 3)

以下是MLE的功能：

> log.L<-function(p=0.3){
+ n<-length(x)
+ (sum(x)-n)*log(1-p)+n*log(p)
+ }

下面是我们的号召：

> library(stats4)
> fir<-mle(log.L)

但是，当我进行调用时，会得到以下错误：

Error in optim(start, f, method = method, hessian = TRUE, ...) : 
non-finite finite-difference value [1]

有什么想法吗？

Answer 1

问题是第一个论点是

计算负对数似然的极小函数.

因此，我们需要

log.L <- function(p = 0.3) {
  n <- length(x)
  -((sum(x) - n) * log(1 - p) + n * log(p))
}
mle(log.L)
#
# Call:
# mle(minuslogl = log.L)
#
# Coefficients:
#   p 
# 0.3