I have been working on a project to write a function for a two-way fixed effects model. I am familiar with the plm function which is what I am trying to replicate at a simpl