Learning sources of variability from high-dimensional observational studies