I am trying to figure out how to write a routine that would compute a two way fixed effects regression model of the form: y=Xβ+u
where u=ci+λ
