Dear Katarzyna Julia Marosz and Firoz Ahmad, thank you for your sharing. I know how to plot interaction effect in LM model. Those packages seem not be able to plot interaction effect in SPATIAL REGRESSION MODELS which use lagsarlm() and errorsarlm() functions rather than lm().
Have you found the answer to your question? I've encountered the same problem with the calculation and visualization of the interation effect in spatial models.