R dataset and code