We report a new stochastic model for studying the polarization-dependent gain in Raman fiber amplifiers with randomly varying birefringence as a function of the pump state of polarization and the polarization mode dispersion parameter. Our theoretical results agree with the previously obtained experimental data.