#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))