From Simulated Exposures to PTA
A simulated exposure distribution is useful for dose justification only when it is summarized against a clearly defined target. Probability of target attainment (PTA) is the proportion of subjects in the target population whose exposure meets a prespecified exposure criterion under a specified regimen. For example, the criterion might be an AUC at or above a stated threshold. The criterion must be stated before the simulated results are classified; changing it after reviewing the exposures changes the question being answered.
The calculation is straightforward, but its meaning depends on the denominator. Every subject counted in the denominator must have a usable exposure result for the specified metric and regimen. A value that fails the exposure criterion is still a valid result and belongs in the denominator. A missing or failed output is different: it does not establish whether that subject attained the criterion.
This estimate describes the simulated cohort for the model, population representation, regimen, and exposure metric used. When the virtual subjects are independent draws from the intended population distribution, the fraction is also a Monte Carlo estimate of the population probability of meeting the criterion. It is not a guarantee for every patient, nor does it by itself establish clinical benefit. The result is conditional on the model and simulation setup.
Specify the Criterion Before Counting
A criterion must identify the exposure metric, any relevant time or dosing condition, the threshold, and the comparison rule. The rule matters at equality. “At least 50” means \(E_i\geq 50\), so a value exactly equal to 50 qualifies. “Greater than 50” means \(E_i>50\), so equality does not qualify. These definitions are not interchangeable.
Classify the unrounded exposure values whenever possible. Rounding a value before comparison can change whether it lies on a threshold. For instance, an underlying value of 49.96 may display as 50.0 when rounded to one decimal place. If the criterion is \(E_i\geq 50\), classification using the displayed value would incorrectly count it as attaining the target.
A practical calculation starts by documenting the rule and then making one classification per subject. Each valid result receives either an attainment indicator of 1 or a nonattainment indicator of 0. Keeping these classifications explicit makes it possible to check the count and to trace an unexpected result back to individual records.
Name the exposure metric, threshold, and comparison rule, such as AUC at least 50 mg·h/L.
Confirm which subject-level results are valid and which are missing or failed. Do not silently convert an unavailable result into a zero.
Count valid results meeting the criterion, then divide by the number of valid results evaluated.
Give the numerator, denominator, percentage, criterion, and number of missing or failed outputs.
Worked Example: An Inclusive AUC Threshold
In a hypothetical simulation, 12 virtual subjects have valid AUC results, in mg·h/L: 40, 50, 52, 61, 49, 55, 48, 70, 50, 44, 58, and 39. The prespecified criterion is AUC \(\geq 50\) mg·h/L. Calculate PTA.
The values meeting the inclusive threshold are 50, 52, 61, 55, 70, 50, and 58 mg·h/L. The two results exactly equal to 50 qualify because the rule is \(\geq\), not \(>\). Thus, 7 of the 12 valid results meet the criterion.
Calculation:
For this hypothetical simulated cohort, estimated PTA is 58.3% for AUC at least 50 mg·h/L. The units belong to the threshold and exposure values; PTA itself is a proportion or percentage and has no exposure units.
Worked Example: A Strict Threshold and Boundary Values
A hypothetical simulation produces 10 valid \(C_{\max}\) results, in mg/L: 9.8, 10.0, 10.1, 11.4, 9.9, 10.0, 12.0, 8.7, 10.6, and 10.0. The criterion is \(C_{\max}>10.0\) mg/L. Find PTA.
The qualifying results are 10.1, 11.4, 12.0, and 10.6 mg/L. The three results equal to 10.0 mg/L do not qualify: the criterion is strictly greater than 10.0, not greater than or equal to 10.0. Therefore, 4 of 10 valid results meet the criterion.
Using an inclusive comparison would instead count the three boundary values and yield \(7/10=70.00\%\). That is a different estimand, not a harmless formatting choice. The comparison rule must match the stated clinical or pharmacological target.
Why the Fraction Estimates a Population Probability
The arithmetic fraction has a probability interpretation when the simulation is designed to represent the target population. Suppose each simulated subject is an independent draw from that population under the specified regimen, and the exposure criterion is fixed. Let \(Z_i\) equal 1 when the draw meets the criterion and 0 otherwise. Each \(Z_i\) is then a Bernoulli random variable with success probability \(p\), where \(p\) is the population probability of meeting that criterion.
The sample fraction is the average of these indicators. Its expectation follows directly from linearity of expectation:
Thus, under the stated sampling conditions and with complete valid outcomes, the simulated fraction is an unbiased Monte Carlo estimate of \(p\). This justification does not mean that a particular finite fraction exactly equals the population probability. It explains why the fraction is an appropriate estimator. Quantifying its Monte Carlo error is a separate question.
Worked Example: Missing Exposure Outputs
A hypothetical run requests results for 20 virtual subjects. Eighteen subjects have valid AUC outputs. Of those 18, 11 meet the prespecified criterion of AUC \(\geq 50\) mg·h/L; two subjects have no usable AUC because their simulation outputs failed. Calculate the fraction among evaluable results and explain its scope.
The valid-result denominator is 18, not 20. The two failed outputs cannot be classified as either attaining or not attaining the AUC criterion. The observed fraction among valid results is
The defensible description is: 11 of 18 subjects with valid outputs met the criterion, a PTA estimate of 61.1% among evaluable simulated results; 2 of 20 requested outputs were unavailable. Dividing by 20 gives \(11/20=55.00\%\), which silently treats both failures as nonattainment and is not justified by the information given.
The fraction among the 18 valid results is not automatically an unbiased estimate for the intended full cohort. If failures are related to exposure or to characteristics that affect exposure, the valid results may not represent the subjects with missing outputs. Investigate the failures and, where possible, correct the simulation or obtain valid results. If they cannot be resolved, report the missingness and qualify the interpretation rather than assigning an invented exposure.
Checks, Alternatives, and Failure Modes
Before reporting PTA, verify that the simulated records match the question being asked. There should be one exposure result per virtual subject for the chosen metric and regimen. Repeated records for a subject must not accidentally count that subject more than once. Confirm units and threshold units agree, and verify that the rule applied in code matches the written criterion. A mismatch between, for example, mg/L and µg/L can change the classifications by orders of magnitude.
Keep missing, invalid, and valid nonattaining results distinct. A valid exposure below a lower target is a nonattainment and contributes 0 to the indicator while remaining in the denominator. A missing value, a numerical failure, or a non-finite result cannot be classified without an explicit, justified rule. A true zero can be a valid exposure in some simulation settings; it should not be removed just because it is zero. Conversely, a failed calculation should not be recorded as zero merely to make the data table complete.
The empirical fraction is the direct approach when the simulation already provides subject-level exposures. A model-based probability calculation could be an alternative if the exposure distribution and criterion permit a valid analytic calculation. A deterministic prediction for a typical subject is not a substitute for PTA: it does not count the fraction of a population distribution meeting the target. The alternatives answer different questions and may rely on different assumptions.
PTA also has limitations. It summarizes attainment relative to one specified criterion and does not describe how far nonattaining subjects fall below it, nor does it establish a clinical outcome. Its interpretation depends on the target definition, the population represented in simulation, the model, and the regimen. A precise fraction cannot compensate for an inappropriate criterion or unrepresentative virtual cohort.
Check Your Understanding
Answer each question, then compare your reasoning with the solutions below.
- For the criterion AUC \(\geq 30\) mg·h/L, hypothetical results are 28, 30, 31, 29, and 35 mg·h/L. What is the PTA?
- For the criterion \(C_{\max}>8\) mg/L, hypothetical valid results are 7.5, 8.0, 8.2, and 9.1 mg/L. What is the PTA, and how is the boundary value treated?
- A simulation requests 50 subjects. Forty-eight have valid results, and 36 of those meet the criterion. What fraction should be reported among evaluable results, and what additional information should accompany it?
- Why can a subject with a valid exposure below a target not be handled in the same way as a subject whose exposure output is missing?
Solutions. (1) Three of five results are at least 30, so PTA is \(3/5=60.00\%\). The value exactly equal to 30 qualifies. (2) Two of four results are strictly greater than 8, so PTA is \(2/4=50.00\%\); the value equal to 8 does not qualify under \(>\). (3) Report \(36/48=75.00\%\) among evaluable results, and state that 2 of the 50 requested outputs were missing or invalid. Investigate the failures because the evaluable fraction may not represent all 50 subjects. (4) The valid below-target exposure is known not to meet the criterion, so it is a valid nonattainment and belongs in the denominator. A missing output does not reveal whether its subject attained the criterion and must not be assigned a status without justification.