Coffee and Productivity

[This article was first published on Freakonometrics » R-english, and kindly contributed to R-bloggers]. (You can report issue about the content on this page here)
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.

On Twitter, I was asked if there were serious research papers published on coffee consumption and labour productivity. There are some papers on coffee breaks and productivity, e.g. Productivity Through Coffee Breaks, but I could not find anything on coffee consumptions. Since I could not find any dataset with personal consumption (maybe I should start keeping tracks of my own consumption to run a study) I tried to find data for national consumption instead (even if we know that both are – clearly – not equivalent)

  • last year, Sabine published on a dataset with consumption of coffee, per country (and per unhabitants),
  • on we can find a dataset with GDP per hour worked for some countries (which can be seen as a common  measure of the productivity of a country)

If we merge those two datasets, we get

> base=read.table(
+ "",
+ header=TRUE,sep=";",dec=",")
> b=base[!$GDP.PPP),]
> plot(b[,3],b[,4],xlab="Coffee Consumption",
+ ylab="GDP per hour worked")
> text(b[,3],b[,4]+1.6,b[,1],cex=.6)
> library(splines)
> X=b[,3]
> Y=b[,4]
> B=data.frame(X,Y)
> reg=glm(Y~bs(X),data=B)
> y=predict(reg,newdata=data.frame(
+ X=seq(0,10,by=.1)))
> lines(seq(0,10,by=.1),y,col="red")

To leave a comment for the author, please follow the link and comment on their blog: Freakonometrics » R-english. offers daily e-mail updates about R news and tutorials about learning R and many other topics. Click here if you're looking to post or find an R/data-science job.
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.

Never miss an update!
Subscribe to R-bloggers to receive
e-mails with the latest R posts.
(You will not see this message again.)

Click here to close (This popup will not appear again)