Create exploringPM_data.r script