add r script and shapefile