curve fitting code snippets in r