Add some R analysis code and plot