How can I use lm() with poly() and ggeffect() to plot model effects when variable in quadratic term is not being plotted

Question

I am trying to plot model effects using the ggeffects package. Usually this is no problem. However, I have a model in which there is an interaction between a quadratic term