How to add a prediction band to the plot of a mixed effect model?

I have created a plot for my lme model using plot_model():

m_lme = lme(Valency ~ logTFI, random = ~1|H, data = data1, method = "ML")
plot_model(A_lme) +
geom_point(size = 3)

It gives me the following plot:

Does anybody know how to add a prediction band to the plot? I would be immensely grateful for your help!