This paper is concerned with the behavior of the homogenized coefficients associated with some random stationary ergodic medium under a Bernoulli perturbation. Introducing a new family of energy estimates that combine probability and physical spaces, we prove the analyticity of the perturbed homogenized coefficients with respect to the Bernoulli parameter. Our approach holds under the minimal assumptions of stationarity and ergodicity, both in the scalar and vector cases, and gives analytical formulas for each derivative that essentially coincide with the so-called cluster expansion used by physicists. In particular, the first term yields the celebrated (electric and elastic) Clausius-Mossotti formulas for isotropic spherical random inclusions in an isotropic reference medium. This work constitutes the first general proof of these formulas in the case of random inclusions.