fixed effect code snippets in r