如何绘制二项分布的对数似然图答案

【问题标题】：How to plot the log-likelihood of binomial distribution如何绘制二项分布的对数似然图
【发布时间】：2014-11-04 12:05:59
【问题描述】：

为了解决 (c)，我认为我需要绘制二项式分布的对数似然图。谁能帮我在R中做到这一点？数据及问题如下；

我想我需要这样的情节：

【问题讨论】：

标签： r optimization plot

【解决方案1】：

这样的事情应该可以工作：

F <- c(18,31,34,33,27,33,28,23,33,12,19,25,14,4,22,7)
M <- c(11,22,27,29,24,29,25,26,38,14,23,31,20,6,34,12)
Y <- F
N <- F + M
#a)
Y / N
#b)
sum(Y) / sum(N)
#c)
logL <- function(p) sum(log(dbinom(Y, N, p)))
#plot logL:
p.seq <- seq(0.01, 0.99, 0.01)
plot(p.seq, sapply(p.seq, logL), type="l")
#optimum:
optimize(logL, lower=0, upper=1, maximum=TRUE)

正如 Ben 所指出的（参见 cmets），数值精度可以通过使用：logL <- function(p) sum(dbinom(Y,N,p,log=TRUE)) 来提高，尤其是在 dbinom() 返回 0 但似然分数实际上只是接近 0。

【讨论】：

logL<-function(p) sum(dbinom(Y,N,p,log=TRUE)) 会稍微准确一些。
感谢 Ben，同意，但认为这种方式更能说明 OP。
JR 非常感谢。我还需要一个人情。 LogL 值出现，但 mle 值不出现。你能解决这个问题吗？
@Günal 你有来自optimize(...) 的mle，如果你想把它添加到情节中，你总是可以做一些基本的事情，比如：abline(v=optimize(logL,lower=0,upper=1,maximum=TRUE)$max)
你帮了大忙。谢谢JR