import pandas as pd
from scipy import stats
# df = .....
df = df[(np.abs(stats.zscore(df[feature_cols])) < 3).all(axis=1)]