marginal-MPHstar-method function

Marginal method for MPHstar class