x=seq(20,70,.5) w =-4.5173+(0.04590*x)+(0.00686*140)+(-0.00694*90)+ (0.00631*200)+(-0.07400*69)+(0.02014*200) w1=-4.5173+(0.04590*x)+(0.00686*140)+(-0.00694*90)+ (0.00631*300)+(-0.07400*69)+(0.02014*200) y=exp(w)/(1+exp(w)) y1=exp(w1)/(1+exp(w1)) plot(x,y1,type="l",xlim=c(20,70),ylim=c(0,.5), ylab="Fitted",xlab="Age",lty=2) lines(x,y,type="l",lty=1) legend("topleft",c("Chol","300","200"),lty=c(NA,2,1))