This paper describes a sample partitioning approach to retain or reject samples from an initial distribution of stability maps using milling test results. The stability maps are calculated using distributions of uncertain modal parameters that represent the tool tip frequency response functions and cutting force model coefficients. Test points for sample partitioning are selected using either (1) the combination of spindle speed and mean axial depth from the available samples that provides the high material removal rate, or (2) a spindle speed based on the chatter frequency and mean axial depth at that spindle speed. The latter is selected when an unstable (chatter) result is obtained from a test. Because the stability model input parameters are also partitioned using the test results, their uncertainty is reduced using a limited number of tests and the milling stability model accuracy is increased. A case study is provided to evaluate the algorithm.