I want to write a function which will perform panel first difference model with trend included.
My work so far
library(plm) fit_panel_
