Calculating average marginal/incremental effects using recycled predictions in R

Question

I am fitting a variety of models for health services use and cost data (heavy zero mass, right skew). For one of my categorical predictors (disease), I would li