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

Question

为了解决(c),我想我需要一个二项分布的对数似然图.有人可以帮我在R做吗？数据和问题如下;

我想我需要这种情节:

Answer 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)

正如本(注释)所指出的那样,通过使用以下方法可以提高数值的准确性:logL <- function(p) sum(dbinom(Y,N,p,log=TRUE))相反,尤其是在dbinom()返回0但是似然得分实际上接近0的情况下,它可以"拯救"你.