Performance Metrics

The analysis module of SLM Lab calculates a set of metrics to measure the performance of an RL agent, which are logged during a run. These metrics can be computed from a time series of returns computed from periodic evaluation.
As motivation, we observe that the returns graph is commonly plotted to present qualitative visual comparison among RL algorithms, but it often contains only one quantitative result that is the final (averaged) return. We acknowledge that a single number is not sufficient to capture all the rich information in a return graph. Qualitatively, the graph visually compares if an algorithm is stronger, trains faster, is more stable, and so on. We turn these into quantitative results by designing metrics that can be compared numerically.
Prerequisites
Our metrics has two simple prerequisites. The first is the time series from the return graph that is readily available. More precisely, it should record the average return, total frames, and total optimization steps during periodic evaluations, or equivalently
R‾i, framei, optstepi, for i=1,…,N\overline{R}_i, \ frame_i, \ optstep_i, \ \text{for } i = 1, \ldots, NRi​, framei​, optstepi​, for i=1,…,N
The average return is computed by averaging over S episodes during evaluation at checkpoint i. To balance between the accuracy of the average and computation cost for evaluation, we typically use S=4, and it can be run as parallel sessions in practice.
The second prerequisite is the average return from a random policy, defined as
R‾rand=E[Rπrand]\overline{R}_{rand} = \mathbb{E}[R_{\pi^{rand}}]Rrand​=E[Rπrand​]
, also known as the random baseline. It is the expectation of the return of a policy with maximum-entropy. Since entropy is maximized for a uniform distribution, the random policy is simply a uniform distribution over the action space. Conveniently, the random policy is already implemented in any environment env using the OpenAI gym interface, and can accessed using the method env.actionspace.sample(). Using this, the random baseline can be obtained easily by running the random policy and averaging the returns over a number of episodes.
SLM Lab Metrics
Now, we propose four metrics derived from the quantitative notions from a returns graph to measure the performance of an algorithm in an environment. These are strength, efficiency (speed), stability, and consistency.
Strength
The first metric is strength, defined as
strength: str=1N∑i=0NR‾i−R‾rand\textit{strength: } str = \frac{1}{N} \sum_{i=0}^N \overline{R}_i - \overline{R}_{rand}strength: str=N1​i=0∑N​Ri​−Rrand​
This measures the average performance of an algorithm against the random baseline, which is a natural zero-point in its scale. Strength immediately tells if an algorithm is working if str > 0, or failing otherwise. Since it is simply the average return with a constant offset, it is familiar and easy to adopt to any evaluation workflow.
Note that the term inside the summation,
R‾i−R‾rand\overline{R}_i - \overline{R}_{rand}Ri​−Rrand​
is simply the local strength at a specific checkpoint i, which measures the strength at a specific frame or optimization step as
stri=strframei=stroptstepistr_i = str_{frame_i} = str_{optstep_i}stri​=strframei​​=stroptstepi​​
for example, strength at 10-million frames where i = 10M. Given a series of local strengths, the maximum and the minimum,
strmax,strminstr_{max}, str_{min}strmax​,strmin​
also capture the best and worst behaviors of an algorithm. Note that strength is simply the mean of local strengths.
The maximum strength captures the peak performance of an algorithm even in the case of regression or policy collapse. The strength, which is proportional to the area under the curve, tells if an algorithm attains its strength quickly and maintains it longer, especially when compared to the maximum strength. Additionally, to obtain a smoother plot, a moving average with a chosen window w (typically 100) can be used on top of strength; this is similar to the usual last-100 episode return average used in several papers.
With a proper zero-point, strength is relative and transitive. If algorithms A, B have strengths 20, 40 respectively, then their strength ratio B/A = 40/20 = 2, hence B is 2 times stronger than A. If an algorithm C is 3 times stronger than B, we can infer that C/A = C/B * B/A = 3 * 2 = 6, and so C is 6 times stronger than A. Strength transitivity allows a new algorithm to be compared against all the other algorithms by comparing it against only one standard algorithm of choice. This can make research easier by requiring a researcher to compare against one standard algorithm instead of many.
Efficiency
The second metric is efficiency, defined with a time unit t as
efficiency: e=∑i=0N1tistri∑i=0N1ti\textit{efficiency: } e = \frac{\sum_{i=0}^N \frac{1}{t_i} str_i}{\sum_{i=0}^N \frac{1}{t_i}}efficiency: e=∑i=0N​ti​1​∑i=0N​ti​1​stri​​
Note that it is a weighted average with 1/t as the weight. Intuitively, it measures how fast an algorithm attains its strength by assigning higher weight to strength earned early on than later. There are two ways of measuring time t that are of interest to us --- frames and optimization steps. Respectively, they yield:
sample efficiency: se=∑i=0N1frameistri∑i=0N1framei\textit{sample efficiency: } se = \frac{\sum_{i=0}^N \frac{1}{frame_i} str_i}{\sum_{i=0}^N \frac{1}{frame_i}}sample efficiency: se=∑i=0N​framei​1​∑i=0N​framei​1​stri​​
training efficiency: te=∑i=0N1optstepistri∑i=0N1optstepi    \textit{training efficiency: } te = \frac{\sum_{i=0}^N \frac{1}{optstep_i} str_i}{\sum_{i=0}^N \frac{1}{optstep_i}}training efficiency: te=∑i=0N​optstepi​1​∑i=0N​optstepi​1​stri​​
Since frames measures the amount of samples used by an agent to attain its strength, the first equation above measures sample efficiency. Similarly, since optimization steps measures the amount of gradient updates used to attain strength, the second equation captures training efficiency. Like strength, these metrics are also relative and transitive. Efficiency can also be defined point-wise by computing it at every checkpoint k simply by summing from i = 0, ..., k to yield
sek,tekse_k, te_ksek​,tek​
. Moreover, the value range of a weighted average is independent from the choice of checkpoint resolution frame_{i+1} - frame_i, although coarser resolution will cause some loss in accuracy. Typically, we use a resolution of 50,000 frames.
Stability
The third metric is stability, defined as
stability: sta=1−∣∑i=0N−1min⁡(stri+1−stri,0)∑i=0N−1stri∣\textit{stability: } sta = 1 - \left| \frac{\sum_{i=0}^{N-1} \min(str_{i+1} - str_i, 0)}{\sum_{i=0}^{N-1} str_i} \right|stability: sta=1−∣∣​∑i=0N−1​stri​∑i=0N−1​min(stri+1​−stri​,0)​∣∣​
This equation measures the average ratio of strength lost between checkpoints as an instability measure, then takes the complement of the ratio so that sta = 1 for a perfectly stable algorithm and decreases for more unstable algorithms. It also has the relative and transitive properties. Note that sta < 0 can happen when the strength drops too much, for instance during a performance collapse. Note that the definition also works for both positive and negative strengths since it is sign-independent. Stability can also be defined locally for each checkpoint i = 0, ..., N-1 as
stai=1−∣min⁡(stri+1−stri,0)stri∣sta_i = 1 - \left| \frac{\min(str_{i+1} - str_i, 0)}{str_i} \right|stai​=1−∣∣​stri​min(stri+1​−stri​,0)​∣∣​
Consistency
The fourth and the final metric is consistency, defined for multiple series of strengths str_{i,j} obtained from multiple repeated sessions j
consistency: con=1−∑i=0N2stdevj(stri,j)∑i=0Navgj(stri,j)\textit{consistency: } con = 1 - \frac{\sum_{i=0}^N 2 stdev_j(str_{i,j})}{\sum_{i=0}^N avg_j(str_{i,j})}consistency: con=1−∑i=0N​avgj​(stri,j​)∑i=0N​2stdevj​(stri,j​)​
This equation measures the average width of the error band to the average strength, then take its complement so that con = 1 for perfectly consistent sessions with 0 error band and decreases for trials with larger error bands. It also has the relative and transitive properties. The error band is defined with simple standard deviation, with the upper and lower bounds as
avgj(stri,j)+stdevj(stri,j)avgj(stri,j)−stdevj(stri,j)avg_j(str_{i,j}) + stdev_j(str_{i,j}) \\
avg_j(str_{i,j}) - stdev_j(str_{i,j})avgj​(stri,j​)+stdevj​(stri,j​)avgj​(stri,j​)−stdevj​(stri,j​)
respectively, and the error band as
(avgj(stri,j)+stdevj(stri,j))−(avgj(stri,j)−stdevj(stri,j))=2stdevj(stri,j)(avg_j(str_{i,j}) + stdev_j(str_{i,j})) - (avg_j(str_{i,j}) - stdev_j(str_{i,j})) = 2 stdev_j(str_{i,j})(avgj​(stri,j​)+stdevj​(stri,j​))−(avgj​(stri,j​)−stdevj​(stri,j​))=2stdevj​(stri,j​)
Consistency can also be defined locally for each checkpoint i = 0, ..., N as
coni=1−2stdevj(stri,j)avgj(stri,j)con_i = 1 - \frac{2 stdev_j(str_{i,j})}{avg_j(str_{i,j})}coni​=1−avgj​(stri,j​)2stdevj​(stri,j​)​
Last modified 1yr ago