Why Compare Designs on a Pareto Front?
A simulation study can show that one sampling design gives narrower confidence intervals, while another achieves coverage closer to its target with fewer patient visits or less computation. These criteria answer different questions. Unless there is a defensible basis for assigning each a common value, combining them into one score can conceal consequential trade-offs.
A Pareto comparison keeps the criteria separate. It identifies designs for which no alternative is at least as good on every criterion and strictly better on at least one. Those designs form a shortlist of trade-offs, not a final ranking: choosing among them still requires a decision about which trade-offs are acceptable.
For example, a study might minimize confidence-interval width, absolute coverage error, patient-visit burden, and computation time. Confidence-interval width is a precision measure, in the same units as the estimated parameter. Coverage error can be defined as the absolute difference between simulated coverage and nominal coverage, reported in percentage points. Patient visits and computation time are burdens. Each metric must have a clear direction before dominance can be assessed.
The choice of coverage objective needs particular care. Maximizing coverage without qualification could reward intervals that are unnecessarily wide. If the goal is coverage near a nominal level, use distance from that level as an objective. If the study instead has a minimum acceptable coverage requirement, coverage can be treated as a constraint. These encode different decision questions and should not be silently interchanged.
How to Determine Dominance
With \(K\) objectives to minimize, write a design's values as \(f(d)=(f_1(d),\ldots,f_K(d))\). Design \(a\) dominates design \(b\) precisely when \(f_k(a)\leq f_k(b)\) for every \(k\), and \(f_j(a)<f_j(b)\) for at least one objective \(j\). If one design is better on one objective but worse on another, neither dominates the other. A tie on one or several objectives does not prevent dominance if the candidate is strictly better on at least one remaining objective and worse on none.
This is a coordinate-wise comparison, not a comparison of overall distance from an imagined ideal design. Units do not need to match: a width in milligrams per liter can be compared with a burden in patient-visits because each is compared only with the same objective for another design. However, all candidates must use the same definition, units, scenario aggregation, and direction for each objective.
For a finite candidate list, a direct pairwise procedure is sufficient: compare each design with every other design, mark a candidate if any alternative dominates it, then retain the unmarked candidates. This requires no distributional assumption about the metrics. It does require that the candidate set and objective definitions be fixed for that comparison.
Worked Example: Finding the Pareto Set Across Four Objectives
Suppose four hypothetical designs are evaluated using median confidence-interval width, absolute coverage error from a 95% target, patient visits, and computation time. Smaller is preferable for every column. The values below are illustrative, not results from a real study.
| Design | Interval width (mg/L) | Coverage error (percentage points) | Patient visits | Computation (hours) |
|---|---|---|---|---|
| A | 0.20 | 2.0 | 120 | 3.0 |
| B | 0.17 | 3.0 | 150 | 2.5 |
| C | 0.23 | 1.0 | 100 | 2.0 |
| D | 0.21 | 2.5 | 130 | 4.0 |
Compare A with D. A has narrower intervals (\(0.20<0.21\) mg/L), smaller coverage error (\(2.0<2.5\) percentage points), fewer visits (\(120<130\)), and less computation (\(3.0<4.0\) hours). Thus A is better on every objective and dominates D.
B is not dominated by A: B has narrower intervals and less computation, but its coverage error and patient-visit burden are greater. C is also not dominated by A: C has smaller coverage error, fewer visits, and less computation, while A has narrower intervals. Comparing B and C likewise gives a trade-off: B has narrower intervals, but C is better on the other three objectives. Therefore A, B, and C are nondominated, while D is dominated. The Pareto set is \(\{A,B,C\}\); it does not say which of those three should be chosen.
Constraints and Objective Choices
Sometimes a criterion represents a minimum standard rather than a benefit to trade against other benefits. For instance, a decision process might require simulated coverage to be at least 93%. In that case, first identify feasible designs that satisfy the requirement, and then compare the feasible designs on the remaining objectives. A design that fails a required standard should not re-enter the shortlist merely because it is excellent on precision or cost.
Define each metric, its units, how it is summarized across scenarios, and whether smaller or larger is preferred.
Mark a design infeasible if it violates a stated requirement, such as a minimum coverage threshold or a maximum acceptable visit burden.
Use coordinate-wise dominance on the remaining objectives and retain the nondominated candidates.
State which criteria distinguish the retained candidates; do not imply that nondominance determines a unique winner.
Worked Example: Applying a Coverage Requirement
Consider four hypothetical designs. The requirement is empirical coverage of at least 93%; among designs that meet it, the objectives are to minimize interval width, patient visits, and computation time.
| Design | Coverage | Interval width (mg/L) | Patient visits | Computation (hours) |
|---|---|---|---|---|
| F | 94% | 0.20 | 120 | 2.4 |
| G | 96% | 0.18 | 140 | 4.0 |
| H | 93.5% | 0.24 | 100 | 1.8 |
| I | 95% | 0.22 | 130 | 3.0 |
All four meet the 93% coverage requirement, so all are feasible. Compare F with I on the three objectives: \(0.20<0.22\) mg/L, \(120<130\) visits, and \(2.4<3.0\) hours. F dominates I, so I is removed.
F, G, and H do not dominate one another. G has the narrowest intervals, but more visits and computation than F. H has the fewest visits and lowest computation, but the widest intervals. F lies between these alternatives on the displayed criteria. The feasible Pareto set is therefore \(\{F,G,H\}\). If the requirement had instead been 94%, H would be infeasible at 93.5%, despite its low burden. A threshold is a decision rule, not a claim that 93.5% coverage is statistically indistinguishable from 94%.
A criterion can be represented either as an objective or as a constraint, but the choice changes the question. Treating coverage error as an objective permits trade-offs between coverage and other metrics. Treating coverage as a threshold rules out any design below the threshold, regardless of how favorable its other metrics are. The threshold should therefore have a substantive justification and be stated in advance where possible.
Trade-Offs, Transformations, and Uncertainty
The numerical scale of an objective does not itself determine its influence in a dominance comparison. If a lower-is-better objective \(f_k\) is replaced by a strictly increasing transformation \(g_k(f_k)\), its ordering is unchanged: \(f_k(a)\leq f_k(b)\) exactly when \(g_k(f_k(a))\leq g_k(f_k(b))\). Consequently, changing an objective from hours to minutes, or from a fraction to a percentage, leaves the Pareto set unchanged if the transformation is applied consistently. Reversing the preference direction does not: it changes the decision problem.
Unlike a weighted sum, Pareto dominance does not require assigning a numerical exchange rate between, for example, one milligram per liter of interval width and ten patient visits. Its cost is that it may leave several candidates, sometimes many, on the Pareto front. A weighted sum can give an overall ranking when defensible weights are available, but its result depends on those weights and may fail to select some nondominated candidates in a discrete or nonconvex set. A lexicographic rule is another alternative: it ranks criteria by priority and considers the next criterion only when candidates tie on the higher-priority one. That is explicit, but it also encodes a strong priority ordering. These approaches answer different questions; none makes trade-offs disappear.
The front is conditional on the selected metrics, scenario summaries, constraints, and candidate designs. Leaving out computation time can make a computationally burdensome design appear attractive. Summarizing performance across scenarios can conceal poor subgroup performance; a separate subgroup requirement or objective may be needed if it matters to the decision. A Pareto-efficient design can still be clinically unacceptable, operationally impractical, or based on a poorly performing analysis pipeline.
Simulation-based metrics are estimates. When two designs are close, Monte Carlo uncertainty may make their apparent ordering unstable. A point-estimate comparison can label one design as dominating another even when the estimated difference is small relative to its uncertainty. Use the uncertainty assessments from the simulation comparison to qualify borderline claims; do not treat a tiny numerical advantage as decisive simply because it is positive. Also check that convergence failures and other estimation failures have not been omitted in a way that makes a design's performance look artificially favorable.
Worked Example: A Dominated Candidate and a Robust Shortlist
Suppose three metrics are minimized: interval width, absolute coverage error, and computation time. For a fourth candidate, suppose the hypothetical values are: J, \(0.40\) mg/L, 1.5 percentage points, and 5 hours; K, \(0.44\) mg/L, 1.0 percentage point, and 3 hours; L, \(0.38\) mg/L, 2.5 percentage points, and 4 hours; M, \(0.45\) mg/L, 2.0 percentage points, and 6 hours.
J dominates M because \(0.40<0.45\) mg/L, \(1.5<2.0\) percentage points, and \(5<6\) hours. K does not dominate J: K has lower coverage error and computation time, but a wider interval (\(0.44>0.40\) mg/L). L does not dominate J either: L has narrower intervals, but higher coverage error and more computation. Similar pairwise comparisons show K and L are not dominated by the other candidates: K trades wider intervals for lower coverage error and computation; L trades a narrower interval for higher coverage error and computation. The point-estimate Pareto set is \(\{J,K,L\}\); M is excluded.
This calculation identifies which candidates merit further decision attention, not a definitive winner. If J's interval-width advantage over K is comparable to simulation uncertainty, the practical distinction may be weak. The appropriate report is that the point estimates put both on the front, with the precision comparison uncertain, rather than claiming that J is unequivocally superior.
Diagnostics and Common Mistakes
- Mixing directions. A larger coverage error is worse, while higher coverage is not automatically better when the target is nominal coverage. Convert the stated goal into a clear objective or constraint before comparing.
- Using inconsistent definitions. Comparing median interval width for one design with mean width for another, or using different scenario weights, makes the coordinates non-comparable.
- Calling every trade-off a tie. Nondominated designs are not equal on all criteria; they are simply not ordered by dominance. Report the criteria on which each candidate gains or loses.
- Choosing the front member with the most favorable single metric. This discards the other objectives without justification. If one metric is decisive, state it as a constraint or priority rather than presenting it as a consequence of Pareto analysis.
- Over-interpreting tiny differences. Rounding can create or erase apparent dominance, and Monte Carlo uncertainty can make near-ties unstable. Compare using adequate precision and qualify conclusions supported only by small estimated gaps.
- Assuming the front guarantees acceptability. Nondominance is relative to the objectives and candidates considered; it does not certify adequate coverage, acceptable burden, or reliable estimation unless those conditions are explicitly assessed.
Check Your Understanding
Answer each question, then compare your reasoning with the solutions.
- Design P has objective values \((0.30, 2, 80)\), and Q has \((0.32, 2, 90)\), with all objectives minimized. Does P dominate Q? Does Q dominate P?
- Why might absolute error from nominal coverage be a more suitable objective than maximizing coverage?
- Three feasible designs have objective values R \((0.25, 1, 100)\), S \((0.23, 2, 120)\), and T \((0.27, 0.5, 90)\), with all objectives minimized. Which designs are nondominated?
- A design has excellent interval width and low computation, but fails a prespecified minimum-coverage requirement. Should it be included in the feasible Pareto comparison? Explain.
Solutions. (1) P dominates Q: it has a smaller first value (\(0.30<0.32\)), ties on the second, and has a smaller third value (\(80<90\)). Q does not dominate P. (2) Coverage is targeted, not maximized without limit: a higher-than-nominal coverage rate may be achieved with unnecessarily wide intervals. Absolute error measures closeness to the target; a minimum-coverage rule instead expresses a feasibility requirement. (3) All three are nondominated. S has the best first value, but R and T have better values on other objectives; R and T also each trade off at least one objective against the other. (4) No. It is infeasible under the stated rule and should be excluded before Pareto comparisons among feasible designs. Its favorable values do not override the requirement.