powered by
Hook to compute the combined residuals with treatment effects added back.
score_2_hook(x)
the ZLMFit