analysing random effect model and plotting panel data random effect

Question

I used the following code for a random effect panel data model:

reffects <- plm(efs~sps+fe+st+maj+exp+hw,data = data.p , model = "random", index = c("ma