A general method is set up, which casts problems of electromagnetic waves in random media into a systematic formalism akin to that of supersymmetric quantum field theory, by analogy with electrons in disordered metals. The characteristic functional of the field is related to that of the medium by means of a diagrammatic expansion; for a linear medium, the vertices are the cumulants of the permittivity, conductivity and permeability. Among the auxiliary fields introduced to account for the field equations, fermionic ones can be eliminated by discarding closed loops in the diagrams. A matrix Green function relates the expectation value of the field to electric and magnetic monopole sources; its general structure and properties are reviewed. T...