If you are looking for a conceptual explanation, here you go. Approximate a dielectric material as a large collection of tiny spheres packed together. Assume each sphere is small enough that the applied electric field across its interior is constant. Then each sphere responds independently to the distinct applied electric field its feels at its center. It is a fairly straightforward calculation to show in general that a dielectric sphere placed in a previously uniform electric field will respond as if it were a perfect electric dipole at its center pointing in the direction of the applied field, i.e. the material response field is the same as that due to a aligned dipole. Finding the total material response field of an extended object then just involves adding up (integrating) all the induced dipole fields of all the little spheres. The material's response field therefore ends up being essentially the polarization of the material, which by definition is the average dipole density.