plot(0,xlim=c(0,4),ylim=c(0,4),xlab="Age",ylab="Calendar Year",type="n",
main="A Lexis diagram",sub="Numerical values from Halley's table")
segments(0,0,4,4)
segments(0,1,3,4)
segments(0,2,2,4)
segments(0,3,1,4)
segments(1,0,1,4)
segments(2,0,2,4)
segments(3,0,3,4)
segments(0,1,1,1)
segments(0,2,2,2)
segments(1,3,3,3)
segments(1,1,2,2,lwd=3)
segments(1,1,1,2,lwd=3)
segments(2,2,2,3,lwd=3)
segments(1,2,2,3,lwd=3)
segments(0,0.8,1.5,2.3,lty=5)
points(1.5,2.3,pch=16)
text(0,1.5,"1238")
text(0.25,1.75,"238")
text(0.5,2.0,"1000")
text(0.75,2.25,"110")
text(1,2.5,"890")
text(1.25,2.75,"35")
text(1.5,3,"855")
text(1.75,3.25,"35")
text(2,3.5,"820")
text(2.25,3.75,"22")
text(2.5,4,"798")
text(0.1,0.5,expression(l[0]))
text(0.25,0.75,expression(d[0]^b))
text(0.5,0.95,expression(L[0]))
text(0.75,1.25,expression(d[0]^a))
text(1.1,1.5,expression(l[1]))
text(1.25,1.75,expression(d[1]^a))
text(1.5,1.95,expression(L[1]))
text(1.75,2.25,expression(d[2]^b))
text(2.1,2.5,expression(l[2]))
text(2.25,2.75,expression(d[2]^a))
text(2.5,2.95,expression(L[2]))
text(2.75,3.25,expression(d[3]^b))
text(3.1,3.5,expression(l[3]))
text(3.25,3.75,expression(l[0]))
text(3.5,3.95,expression(L[3]))