Value-added Plots for models estimated by plm-package in r?

If I estimated my model using the lm-function, I would be able to create Value-Added Plots typing the following code:

ols_model <- lm( log(y) ~ log(x1) + log(x2) + log(x3), data = data)avPlots(ols_model)

However, I estimated my model using the plm-function because I'm using panel data methods.

fe_model <- plm( log(y) ~ log(x1) + log(x2) + log(x3),                      model = "within",                       data = data,                       index = c("country","year"),                       effect = "twoways")

How can I create Value-Added Plots in this case?