zscore = (f - f.mean()) / f.std()