RUV4 function

Remove Unwanted Variation, 4-step