Recreate some SAS procedures in R
