Add the R script and codebook and update README