lm() code snippets in r

fit a fixed effect model in r