Add R-script SurfStats.R