stats code snippets in r