r - 如何在 R 中为函数添加红线

Question

对于如下图所示的函数l(p;x)=log p+3log(1-p)，

如何为 l(p;x) 添加标记值并将其 p 值标记为红色线。比如我要标记l(p;x)=-3.2及其对应的p区间。

我可以在 R 中绘制函数:

l=function(p){log(p)+3*log(1-p)}
plot(l, 0, 1)

Answer 1

这个 Base R 解决方案怎么样:

l=function(p){log(p)+3*log(1-p)}

plot(l, 0, 1)
#draw line at -3.2
abline(h=-3.2, col="red")

#find the upper and lower solutions
lower <-uniroot(function(x) {l(x)+3.2}, lower=0.01, upper=0.4)
upper <-uniroot(function(x) {l(x)+3.2}, lower=0.4, upper=1)

#plot
segments(lower$root, -20, y1=-3.2, col="blue")
segments(upper$root, -20, y1=-3.2, col="blue")

#or 
abline(v=lower$root, col="red")
abline(v=upper$root, col="red")