#jan 14 graphclick

# exemestane
tx=read.table("exemestane.txt",head=T,col.names=c("t","ci"));tx$ci=tx$ci/100
y = c(seq(0,5,1))
t.n = c(2285,1902,1468,980,464,77);
n = length(tx$ci); n
S = 1 - tx$ci
ratio = S[2:n] / S[1:(n-1)] ; ratio
n.at.risk = 1 / (1 - ratio); 
plot(tx$t[2:n],n.at.risk,pch=19,cex=1, xlab="Time (years)", ylim=c(0,2500), xlim=c(0,5),ylab="No. At Risk")
points(y,t.n,col=2,cex=3,pch=1)
points(2,2350,col=2,cex=3,pch=1)
text(2,2350,"Reported numbers at time 0, 1, 2, ... ",col=2,
  adj=c(-0.1,0.5) )

points(2,2150,col=1,cex=1,pch=19)
text(2,2150,"Derived numbers at times of events",col=1, 
  adj=c(-0.1,0.5))