The analytical calculation of molecular dipole-moment derivatives for MCSCF wave functions is described. The formalism is based on exponential unitary transformation of the wave function and symmetric orthonormalization of the molecular orbitals. The response equations are solved using an iterative, direct technique to allow for large configuration expansions. Translational and rotational symmetries of the dipole moment are used to minimize computational costs. Sample calculations involving several thousand configurations are presented for H2O and ONF.