Add r script
Update README