compAnom function

Calculates anomalies