I want to create a function which will perform within panel regression :
My work so far
library(plm) fit_panel_lr <- function(y, x)