Survival of the sweetest

On receiving an advent calendar from one of our course directors, I suggested we could track each other’s chocolate consumption in a survival analysis and establish who was eating significantly more chocs. Strangely, everybody refused to take part, so I am n=1. Looking forward to being able to highlight some “shocks” on the graph.

So far, so smug

So far, so smug

If you really want to know:

library(jpeg)
today<-12
chocs<-c(5,5,7,7,7,10,10,10,12)
lc<-length(chocs)
choc.days<-rep(NA,25)
for (i in 1:25) {
choc.days[i]<-sum(chocs==i)
}
choc.surv<-25-cumsum(choc.days)
holly<-readJPEG("holly.jpg")
chocplot<-function(){
plot(1:today,choc.surv[1:today],type="s",lty=1,col="chocolate4",lwd=10,
ylim=c(0,25),xlim=c(0,25),
xlab="Day in December",ylab="Surviving chocolates",
main="Survival of chocolates in Robert's 2012 office advent calendar",
bty="n")
text(x=10,y=23,labels="Survival curve",col="chocolate4")
text(x=6.7,y=17,labels="One a day...")
text(x=18,y=20,labels="Restraint",col="red",font=3)
text(x=22,y=23,labels="Scrooge",col="green",font=3)
text(x=7,y=6,labels="Abandon",col="red",font=3)
text(x=3,y=2,labels="Nausea",col="green",font=3)
lines(x=c(1,25),y=c(25,1),lty=3)
rasterImage(holly,22,0,25,3)
}
windows()
chocplot()
jpeg("Advent_calendar.jpg")
chocplot()
dev.off(dev.cur())

Advertisements

Leave a comment

Filed under R

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s