delta.s.estimate function

Calculates model-based or robust residual treatment effect