Aggregate coeffcients from linear regression with interactions [r]

Question

I am trying to write a function that would take linear regression model, Design matrix X and would return marginal effects of our coefficients.

Marginal effects are g