We applied finite difference time domain (FDTD) algorithm to the study of field and intensity correlations in random media. Close to the onset of Anderson localization, we observe deviations of the correlation functions, in both shape and magnitude, from those predicted by the diffusion theory. Physical implications of the observed phenomena are discussed.