A method was developed for studying mixing of cohesive pharmaceutical mixtures. A combination of accurate sampling and NIR spectroscopic analysis was developed as a suitable method to determine homogenization of magnesium stearate as a function of blending variables. A typical pharmaceutical blend containing a ratio 35:64:1 lactose, avicel, and magnesium stearate was used as a model system. The method accounted for variability of the concentration of magnesium stearate as well as variability of the excipients. Levels of magnesium stearate as low as 0.05% could be resolved by the method, and showed a predicting confidence interval above 98%.