2014.11.23 R code explanation added.