Determination of the optimal number of strata for propensity score subclassification