d<-seq(0, 5, len=100)
pre<-predict(f, data.frame(x = d))
p<-exp(pre)/(1+exp(pre))
norell$y<-norell$uccess/norell$n
plot(norell$x, norell$y); lines(d, p)
d 是给出曲线横 坐标的点，pre计算预测值，p是相应预测概率，用plot函数和lines给出散点图和对应的预测曲线。
PS：f为logistic函数    norell改为数据集的名       norell$y<-norell$uccess/norell$n主要是算因变量取1的个数与总观测值的比值   
我有点说不明白吧   是从书上看到的  这是一个例子

楼主问题解决了吗，我也遇到这样的问题，能否告知怎末画的啊，不胜感激

x<-seq(-10,20,0.05)
y1<- 4.5+25/(1+exp(0.8*(3-x)))
plot(x,y1,type="l",main="四个参数four-parameter logistic")