Although NMR spectroscopic techniques coupled with multivariate statistics can yield much useful information for classifying biological samples based on metabolic profiles, biomarker identification remains a time-consuming and complex procedure involving separation methods, two-dimensional NMR, and other spectroscopic tools. We present a new approach to aid complex biomixture analysis that combines diffusion ordered (DO) NMR spectroscopy with statistical total correlation spectroscopy (STOCSY) and demonstrate its application in the characterization of urinary biomarkers and enhanced information recovery from plasma NMR spectra. This method relies on calculation and display of the covariance of signal intensities from the various nuclei on the same molecule across a series of spectra collected under different pulsed field gradient conditions that differentially attenuate the signal intensities according to translational molecular diffusion rates. We term this statistical diffusion-ordered spectroscopy (S-DOSY). We also have developed a new visualization tool in which the apparent diffusion coefficients from DO spectra are projected onto a 1D NMR spectrum (diffusion-ordered projection spectroscopy, DOPY). Both methods either alone or in combination have the potential for general applications to any complex mixture analysis where the sample contains compounds with a range of diffusion coefficients.