I am calculating LOD scores in genetic epidemiology where the number of non recombinants is 11 and recombinants is 1. I created a function to calculate the z scores for the following thetas. theta <- c(0,0.01,0.05,0.10,0.20,0.45,0.5) The z scores for the respective thetas in this cases are: -Inf,1.56,2.07,2.11,1.85,0.410,0.0 I attempt to plot these values using the following commands. plot(theta, zscore, xlab = "theta", ylab = "LOD Score", type = "l", ylim = c(-5,5), panel.first = grid(5,16), panel.last = abline(h=0),axes=F) axis(side = 1,at = c(0,.1,.2,.3,.4,.5)) axis(side = 2,at = c(-5,-4,-3,-2,-1,0,1,2,3,4,5)) The plot works but does not plot from the points (0,-Inf) to (0.01,1.56). As this is an asymptotic value, I'm unsure that the plot function can plot to -Inf. I am using R 1.7.1. on WinXP sp1. Is there a way to plot the asymptote? Regards Tom