Exercise 4: Quantifying probability of detection#
Goal#
Understand how MAPIT can be used to evaluate probability of detection. So far, exercises have focused on exploring qualitative changes that are induced by changes in measurement uncertainty or material loss. This section will focus on developing quantitative metrics for safeguards performance.
Caution
This exercise assumes that you are familiar with MAPIT and can perform tasks discussed in the previous exercises; exercise 1, exercise 2 and exercise 3.
Problem setup#
Start this exercise by launching MAPIT, selecting the
Normal
dataset.This tutorial starts by assuming you are at the main MAPIT interface
Select all the checkboxes for available statistical tests
A total of five (5) checkboxes should be checked
MUF
Cumulative MUF
Sigma MUF
SITMUF
Page’s test on SITMUF
Next, configure the required parameters in the statistical box using the same parameters from exercise 1 and 2:
MBP
: 416Iterations
: 50If running on a lower performance device, try running with
Iterations
: 20
Analysis Element/Index
: UTemporal Offset
: Empty
Set the measurement errors by pressing the Select Errors button. Choose 0.5% for all random and systematic errors.
Run MAPIT by pressing the Run button.
Note
Safeguards statistical tests are commonly tuned to have an average false alarm probability of 5% per year. The next steps will use MAPIT functionality to determine this threshold for Page’s trend test on SITMUF.
Caution
The data provided in these examples are 6480 hours long, which translate to about 270 days (a reasonable estimate of an operational year). Therefore, one iteration is a single year of simulated operation. When using your own data, use caution in determining statistical thresholds if your datasets have different lengths of time.
Determining statistical thresholds#
Plot Page’s trend test on SITMUF
This is labeled as
U Page's test on SITMUF
in MAPITPlotting a quantity is necessary to use the threshold functionality
Note
Plots in this excercise are using 300 iterations to have better statistics, so if less iterations are being run the results may differ slightly from those shown here. In large iteration datasets, MAPIT only plots 15 iterations to ensure the plotting window is not overloaded.
Use the
Statistical Thresholds
box to determine an appropriate threshold (i.e. 5% FAP)This is performed by entering values in the
Enter Threshold
boxMAPIT returns the number of iterations that exceed this threshold
MAPIT operates on the entire dataset, not just the iterations that are plotted. This can be important in cases where many iterations are calculated and only a few are plotted.
Your threshold should be near 0.5
Due to the randomness of the calculations, your threshold might be slightly different
Increasing the number of iterations can help obtain a more precise estimate of the threshold
The uncertainty in the threshold itself should roughly decrease with sqrt(iterations).
Evaluating probability of detection#
Load the
Abrupt
dataset.
Select the same options from steps 2, 3, and 4 then run MAPIT.
Select all the options in the Tests/Uranium box
Select a value of 0.5% for all random and systematic errors
MBP
: 416Iterations
: 50If running on a lower performance device, try running with
Iterations
: 20
Analysis Element/Index
: UTemporal Offset = empty
: Empty
Plot Page’s trend test on U SITMUF
This is labeled as
U Page's test on SITMUF
in MAPIT
Enter in the previously determined threshold to determine the probability of detection for this material loss.
The value should be approximately 50-60%
Note
In this exercise, we changed the measurement uncertainties for all locations to be the same value. However, in practice, different measurement technologies are deployed based on safeguards need. On your own, try changing just some high impact measurement locations to 0.5% while leaving others at higher levels to see if you can reach similar results.
Summary#
🎉 You have finished the fourth MAPIT exercise. 🎉
This exercise introduced the the SITMUF test in MAPIT along with the capability to set threshold and evaluate probability of detection. In this exercise you:
Set a threshold based on a false alarm probability
Evaluated probabilities of detection
Explored the SITMUF transform
Explored Page’s trend test on SITMUF
For an optional challenge, continue to exercise 5.