robust statistics code snippets in r