fixed effects code snippets in r