我如何在与ggplot中Y中给定值相对应的x值处添加垂直线

Question

library(polynom)
set.seed(12345)
x<-as.numeric(seq(1,5,0.01))
lp<-rnorm(401,-0.7,1)-4*x+0.9*x^2
link_lp <- exp(lp)/(1 + exp(lp))
y<-ifelse((runif(401) < link_lp),0,1)
df<-data.frame(y=y,x=x)
f<-glm(y~poly(x,degree=2,raw=TRUE),family="binomial")
ymark<-c(0.2,0.3,0.35,0.4,0.45,0.5,0.6)
ggplot(df, aes(x=x, y=y))+
  geom_jitter(size=1, alpha=0.2,height=0.05)+
  stat_smooth(method="glm",
     formula=y~poly(x,degree=2,raw=TRUE),
     method.args =list(family = "binomial"),
     colour="blue", size=1.5)+
  xlab("x")+
  ylab("Probability of event")+theme_bw()+
  geom_hline(yintercept=ymark,col="red")+
  scale_y_continuous(breaks=c(0,ymark,1))

figure 1

接下来我要添加垂直线对应水平红线又添加了一条geom_vline()

ggplot(df, aes(x=x, y=y))+
   geom_jitter(size=1, alpha=0.2,height=0.05)+
   stat_smooth(method="glm",
     formula=y~poly(x,degree=2,raw=TRUE),
     method.args =list(family = "binomial"),
     colour="blue", size=1.5)+
xlab("x")+
ylab("Probability of event")+theme_bw()+
geom_hline(yintercept=ymark,col="red")+
scale_y_continuous(breaks=c(0,ymark,1))+
geom_vline(xintercept=solve(polynomial(coef = coef(f)),
   b=log(ymark/(1-ymark))),col="green")

但是报错，请问如何更改我的代码？

Answer 1

solve returns 复数。如果你只是想要真实的部分，请在其周围放置 Re。

ggplot(df, aes(x=x, y=y))+
  geom_jitter(size=1, alpha=0.2,height=0.05)+
  stat_smooth(method="glm",
              formula=y~poly(x,degree=2,raw=TRUE),
              method.args =list(family = "binomial"),
              colour="blue", size=1.5)+
  xlab("x")+
  ylab("Probability of event")+theme_bw()+
  geom_hline(yintercept=ymark,col="red")+
  scale_y_continuous(breaks=c(0,ymark,1))+
  geom_vline(xintercept=Re(solve(polynomial(coef = coef(f)),
                                         b=log(ymark/(1-ymark)))),col="green")

我如何在与ggplot中Y中给定值相对应的x值处添加垂直线

how may I add vertical lines at x values corresponding to given values in Y in ggplot

r

ggplot2

polynomials