I used the data for random effect model with panel data manager then I wrote this code:
reffects <- plm(efs~sps+fe+st+maj+exp+hw,data = data.p , model = "