我现在想通过R计算半数致死量，数据如下：
conc           nlive
1000000      0
100000        1
10000          8

conc列是毒物浓度，nlive为存活的动物数。每种剂量的动物总数为10只，总共三个浓度。

conc<-c(1000000,100000,10000)
dose<-log10(conc)
numlive<-c(0,1,8)
total<-cbind(numlive,10-numlive)
fit<-glm(total~dose,family = binomial)
library(MASS)
tmp<-dose.p(fit,p=0.5)
logld50<-tmp[1]
SE<-attr(tmp,"SE")
ld50<-10**logld50
ci<-10**(logld50+c(-1.96,1.96)*SE)
logld50
ld50
ci

参考了文献后设计了以上代码计算ld50，及其95%置信区间。但我想画一个以剂量为横坐标，以死亡率为纵坐标的“S”型曲线。这该如何实现呢？