statistical modeling code snippets in r