From a Design Question to a Performance Target
A simulation comparison is useful only if it answers a clearly stated question. “Which sampling design performs better?” is not yet a statistical target: better could mean less bias, more precise estimates, intervals that cover the true parameter more often, or a greater chance of obtaining a usable fit. These criteria can favor different designs. The first task is therefore to define the quantity whose value would count as evidence about performance.
The earlier material on simulation workflows provides the broad sequence for comparing designs. Here the emphasis is narrower: translate the scientific or operational question into an estimand, a precisely defined target quantity. An estimand states what is to be summarized over repeated simulated studies, and under which design and scenario. It is not the result from a finite set of simulated replicates.
To make an estimand interpretable, state at least four things: the target parameter, the design being evaluated, the scenario that defines the conditions of repeated studies, and the performance quantity of interest. A fifth detail is often essential: the event or denominator used to define the quantity. This matters especially when analyses can fail or produce unusable results.
For a parameter whose true value in a specified scenario is \(\theta\), let \(\hat{\theta}\) denote the estimate produced by applying the stated analysis to one simulated study. The same design and scenario define a distribution of possible simulated datasets and analysis outcomes. Performance estimands describe features of that repeated-study distribution; they are not properties of one dataset alone.
Choose the Estimand That Matches the Question
A broad question such as “Are estimates accurate?” must be made more specific. Accuracy might refer to average signed error, typical spread, or the frequency with which an estimate lies close enough to the truth. These are related but not interchangeable.
Bias answers a question about systematic direction, not the size of typical errors. An estimator can have zero bias but substantial variability. A precision question can be represented by the standard deviation of \(\hat{\theta}\), with smaller values indicating less spread under the same conditions. Alternatively, the question may concern a practical tolerance, such as the probability that an estimate is within a specified distance of the truth. That probability requires the tolerance and scale to be stated.
For interval estimation, “Does the interval work?” usually refers to coverage: the probability that the interval procedure contains the true parameter over repeated studies. The interval procedure, including its construction method, is part of the definition. If failed analyses produce no interval, one must decide whether the target is coverage among intervals that were successfully produced or the probability of both producing an interval and covering the truth. Those answer different operational questions.
For an operational question such as “How often can this pipeline deliver an estimate we can use?”, define success as an event before counting outcomes. For example, success might require convergence, a finite estimate, and an estimate within a prespecified tolerance. The target is then the probability of that complete event across all attempted simulated studies. A looser criterion can be appropriate, but its interpretation must be stated.
State the Conditioning and Scope
A performance estimand is conditional on the conditions used to generate and analyze the simulated studies. It is therefore helpful to write the design and scenario into the question itself: for example, “What is the expected bias of the treatment-effect estimator for design A under scenario S?” This makes clear that the result does not automatically describe other designs or scenarios.
When several scenarios matter, two common targets are scenario-specific performance and an average across scenarios. A scenario-specific estimand reports performance separately under each defined condition. An average requires explicit weights. If scenario \(s\) has performance \(Q_s\) and weight \(w_s\), with weights summing to one, the weighted target is:
Equal weights answer the question “What is the average performance across these listed scenarios, treating each equally?” Weights based on a target population or decision context answer a different question. Neither choice is universally correct. The weights, scenarios, and interpretation should be specified rather than selected after seeing which design looks best.
Likewise, a comparison must name the contrast of interest. The difference in bias between two designs, the ratio of their standard deviations, and the difference in success probabilities are distinct comparisons. A design can improve one criterion while worsening another, so a single overall ranking needs an explicit decision rule or priority among criteria.
Worked Examples: Making Questions Operational
Worked Example: Bias for a Concentration Parameter
Suppose a hypothetical simulation evaluates an estimator of a population concentration parameter whose fixed true value is \(\theta=3.00\) mg/L under one specified design and scenario. Five simulated studies produce estimates of 3.1, 3.2, 3.0, 3.3, and 2.9 mg/L. Define and estimate the bias from these results.
The target estimand is the expected estimation error, \(E(\hat{\theta})-\theta\), in mg/L for this design and scenario. The five-replicate estimate of that target is the mean estimated error:
In this finite set of replicates, the average estimate is 0.100 mg/L above the true value. This is an estimate of the simulation bias, not proof that the estimator’s exact repeated-study bias is precisely 0.100 mg/L. Five replicates may give a noisy summary; the estimand remains defined by the repeated-study expectation.
Worked Example: Coverage for an Interval Procedure
In a hypothetical evaluation, 20 simulated studies each produce a 95% interval for a fixed true parameter. The true value lies inside 15 of the 20 intervals. The question is whether the interval procedure achieves its nominal coverage under this design and scenario.
The coverage estimand is the probability that the procedure’s interval contains the true value in a repeated simulated study. For the 20 replicates, estimate it by the fraction of intervals that contain the truth:
The estimated coverage is 75.0% in these replicates. It is below the nominal 95%, but 20 replicates provide limited evidence about the exact coverage probability. The denominator is appropriate here because all 20 studies produced intervals. If some replicates produced no interval, the target would need to say whether it concerns coverage conditional on an interval being produced or successful interval production and coverage jointly.
Worked Example: Probability of a Successful Fit
A hypothetical analysis pipeline is applied to 50 simulated studies. Eight fits do not converge. Of the remaining converged fits, two produce nonfinite estimates. Define success as both convergence and a finite estimate, and assume these two failure categories are disjoint as stated. Estimate the probability of success.
The estimand is the probability that one attempted study yields a converged fit with a finite estimate under this design and scenario. Because the failure categories are disjoint, 8 + 2 = 10 distinct attempts fail the criterion. Thus 50 − 10 = 40 attempts succeed, and the estimate is:
The estimated success probability is 80.0%, using all 50 attempted fits as the denominator. If the two nonfinite estimates were not known to be disjoint from the eight nonconverged fits, subtracting both counts would not be justified. The number of distinct failures would then require the overlap, or would need to be counted directly.
A Short Justification: Why Event Probabilities Are Fractions
The connection between a probability estimand and a replicate fraction follows from an indicator variable. For replicate \(r\), define \(I_r=1\) when the event of interest occurs and \(I_r=0\) otherwise. If the simulated replicates are generated under the same specified design and scenario, each indicator has expectation equal to the event probability \(q\): \(E(I_r)=1\cdot q+0\cdot(1-q)=q\). The average indicator is the observed event fraction.
Thus the fraction is a natural estimate of the probability when all \(R\) attempts contribute a defined indicator. For a success probability, each attempt has a success or failure indicator. For coverage conditional on interval production, however, attempts with no interval are excluded by definition; the resulting fraction estimates a conditional probability, not the unconditional probability of successful coverage. The denominator is part of the estimand, not a bookkeeping detail.
Diagnostics, Alternatives, and Limits
Before interpreting a performance summary, check that its operational definition matches the design question. Ask whether the target parameter and units are clear, whether the scenario is stated, whether a probability event has an unambiguous rule, and whether all relevant outcomes have a defined place in the denominator. Also check that the summary is not being interpreted beyond its scope: an estimand for one scenario does not establish performance in unexamined scenarios.
A practical comparison of estimands is useful when deciding what “better” means:
| Design question | Possible estimand | Interpretation |
|---|---|---|
| Are estimates systematically high or low? | \(E(\hat{\theta})-\theta\) | Average signed error in the parameter’s units |
| How variable are estimates? | Standard deviation of \(\hat{\theta}\) | Repeated-study spread; smaller means more precise |
| How often are estimates practically close? | \(P(|\hat{\theta}-\theta|\leq \delta)\) | Probability of falling within a stated tolerance \(\delta\) |
| How often do intervals contain the truth? | \(P(L\leq\theta\leq U)\) | Coverage probability for a specified interval procedure |
| How often is an analysis usable? | \(P(\text{defined success event})\) | Probability of meeting prespecified operational criteria |
These alternatives address different failure modes. Bias alone does not measure spread; spread alone does not reveal systematic displacement. Coverage assesses the interval procedure rather than just the point estimator. A success probability can capture computational or operational usability but says nothing about the quality of a successful estimate unless quality is included in the success event. If a design decision depends on more than one concern, report multiple estimands or specify a transparent rule for combining them.
A common mistake is to define success after reviewing which replicates failed, or to remove failed replicates without stating that the target is conditional on success. The reasoning might be: “Only converged fits can be evaluated, so exclude nonconverged fits and report the success rate among the rest.” That fraction is not an overall success rate: it conditions on an intermediate event and may even omit the failures that the success criterion is meant to count. The correction is to define the success event first and use all attempted studies in the denominator for its probability. If the intended target is conditional performance among converged fits, label it explicitly as conditional and report the convergence probability separately.
Finally, an estimand is only as informative as its definitions and the conditions under which it is evaluated. A tolerance must be scientifically meaningful and expressed in the parameter’s units. A coverage target needs a specified interval procedure. An average across scenarios needs justified weights. And a simulation result estimates the chosen target under the stated assumptions; it does not guarantee performance in scenarios that were not represented.
Check Your Understanding
For each question, identify the performance estimand or the definition needed to make it precise. Solutions follow.
- A team asks whether a clearance estimator tends to overestimate clearance. What estimand answers this question, and what are its units?
- A design produces estimates with smaller standard deviation but more bias than a competing design. Can one say it is unambiguously better? Why or why not?
- Among 100 simulated studies, 10 intervals are not produced, and 72 of the remaining 90 contain the true value. What coverage quantity does \(72/90\) estimate?
- A pipeline has 50 attempts, 8 nonconverged fits, and 2 converged fits with nonfinite estimates. If the categories are explicitly disjoint, estimate the probability of a converged fit with a finite estimate.
- Two scenarios have performance values \(Q_1=0.80\) and \(Q_2=0.60\). What is their equally weighted average, and what additional choice is needed for a population-weighted average?
Solutions. (1) Use bias, \(E(\hat{\theta})-\theta\); its units are the units of clearance, and positive values indicate overestimation on average. (2) No. The first design is more precise by standard deviation, but the second may be less biased; “better” requires a stated priority or decision rule. (3) It estimates coverage conditional on an interval being produced: \(72/90=0.800\), or 80.0%. It does not estimate the unconditional probability of producing an interval that covers. (4) Disjointness makes the number of distinct failures \(8+2=10\); success count is \(50-10=40\), so the estimated probability is \(40/50=0.800\), or 80.0%, with all attempts in the denominator. (5) The equally weighted average is \((0.80+0.60)/2=0.700\). A population-weighted average additionally requires weights reflecting the intended population or decision context.