Calculate time used for each possession.

Note that this counts from the pull or turnover, not from the first catch/take
(since we don't have that information for the opposing team, and we would
like it to stay consistent between teams).

diff --git a/ultimate.js b/ultimate.js
index 1fbb3d15f84e0dcb2d493fe17a076ea0697485fd..f782e754bd376ad4bad3781cc35e683a64657a4f 100644 (file)
--- a/ultimate.js
+++ b/ultimate.js
@@ -119,6 +119,9 @@ function add_3cell_ci(tr, ci) {
        } else {
                text = ci.val.toFixed(2);
        }
+       if (ci.unit !== undefined) {
+               text += ' '+ ci.unit;
+       }
        let element = add_3cell(tr, text);
        let to_x = (val) => { return ci_width * (val - ci.min) / (ci.max - ci.min); };
 
@@ -242,6 +245,10 @@ function calc_stats(json, filters) {
 
                'touches_for_turnover': 0,
                'touches_for_goal': 0,
+               'time_for_turnover': [],
+               'time_for_goal': [],
+               'their_time_for_turnover': [],
+               'their_time_for_goal': [],
        };
        let globals = players['globals'];
 
@@ -270,6 +277,8 @@ function calc_stats(json, filters) {
                let we_lost_disc = false;
                let they_lost_disc = false;
                let touches_this_possession = 0;
+               let possession_started = null;
+               let extra_ms_this_possession = 0;  // Used at stoppages.
 
                // The last used formations of the given kind, if any; they may be reused
                // when the point starts, if nothing else is set.
@@ -507,17 +516,26 @@ function calc_stats(json, filters) {
                                if (keep && type !== 'goal' && !(type === 'set_defense' && last_pull_was_ours === null)) {
                                        ++globals.turnovers_lost;
                                        globals.touches_for_turnover += touches_this_possession;
+                                       if (possession_started != null) globals.time_for_turnover.push((t - possession_started) + extra_ms_this_possession);
                                } else if (keep && type === 'goal') {
                                        globals.touches_for_goal += touches_this_possession;
+                                       if (possession_started != null) globals.time_for_goal.push((t - possession_started) + extra_ms_this_possession);
                                }
                                touches_this_possession = 0;
+                               possession_started = null;
+                               extra_ms_this_possession = 0;
                        } else if (type === 'set_offense' || type === 'their_goal' || type === 'their_throwaway' || type === 'defense' || type === 'interception') {
                                offense = true;
                                they_lost_disc = true;
                                if (keep && type !== 'their_goal' && !(type === 'set_offense' && last_pull_was_ours === null)) {
                                        ++globals.turnovers_won;
+                                       if (possession_started != null) globals.their_time_for_turnover.push((t - possession_started) + extra_ms_this_possession);
+                               } else if (keep && type === 'their_goal') {
+                                       if (possession_started != null) globals.their_time_for_goal.push((t - possession_started) + extra_ms_this_possession);
                                }
                                touches_this_possession = 0;
+                               possession_started = null;
+                               extra_ms_this_possession = 0;
                        }
                        if (type === 'goal' || type === 'their_goal') {
                                between_points = true;
@@ -544,6 +562,19 @@ function calc_stats(json, filters) {
                                }
                                live_since = t;
                        }
+                       // The rules for possessions are slightly different; we want to start at pulls,
+                       // not previous goals.
+                       if ((last_offense !== offense && type !== 'goal' && type !== 'their_goal' && live_since !== null) ||
+                           type === 'pull' || type === 'their_pull') {
+                               possession_started = t;
+                       } else if (type === 'stoppage') {
+                               if (possession_started !== null) {
+                                       extra_ms_this_possession = t - possession_started;
+                                       possession_started = null;
+                               }
+                       } else if (type === 'restart' && live_since !== null && !between_points) {
+                               possession_started = t;
+                       }
 
                        if (type === 'pull') {
                                last_pull_was_ours = true;
@@ -927,6 +958,49 @@ function make_poisson_ci_large(val, num, z)
        };
 }
 
+// We don't really know what to expect for possessions; it appears roughly
+// lognormal, but my head doesn't work well enough to extract the CI for
+// the estimated sample mean, and now my head spins around whether we do
+// this correctly even for Poisson :-) But we do a simple (non-BCa) bootstrap
+// here, which should give us what we want.
+//
+// FIXME: z is ignored.
+function make_ci_duration(vals, z)
+{
+       // Bootstrap.
+       let sample_means = [];
+       for (let i = 0; i < 300; ++i) {
+               let sum = 0.0;
+               for (let j = 0; j < vals.length; ++j) {
+                       sum += vals[Math.floor(Math.random() * vals.length)];
+               }
+               sample_means.push(sum / vals.length);
+       }
+       sample_means.sort();
+
+       // Interpolating here would probably be better, but OK.
+       let low  = sample_means[Math.floor(0.025 * sample_means.length)];
+       let high = sample_means[Math.floor(0.975 * sample_means.length)];
+
+       // Find the point estimate of the mean.
+       let sum = 0.0;
+       for (const v of vals) {
+               sum += v;
+       }
+       let avg = sum / vals.length;
+
+       return {
+               'val': avg / 1000.0,
+               'lower_ci': low / 1000.0,
+               'upper_ci': high / 1000.0,
+               'min': 0.0,
+               'max': 120.0,
+               'desired': 0.0,
+               'inverted': false,
+               'unit': 'sec',
+       };
+}
+
 function make_table_general(players) {
        let rows = [];
        {
@@ -1037,7 +1111,6 @@ function make_table_team_wide(players) {
                        if (q === 'globals') continue;
                        goals += p.goals;
                }
-
                let row = document.createElement('tr');
                let name = add_3cell(row, 'Touches per possession (all)', 'name');
                let touches = globals.touches_for_turnover + globals.touches_for_goal;
@@ -1065,6 +1138,35 @@ function make_table_team_wide(players) {
                rows.push(row);
        }
 
+       // Time. Here we have (roughly) for both.
+       {
+               let row = document.createElement('tr');
+               let name = add_3cell(row, 'Time per possession (all)', 'name');
+               let time = globals.time_for_turnover.concat(globals.time_for_goal);
+               let their_time = globals.their_time_for_turnover.concat(globals.their_time_for_goal);
+               add_3cell(row, '');
+               add_3cell_ci(row, make_ci_duration(time, z));
+               add_3cell(row, '');
+               add_3cell_ci(row, make_ci_duration(their_time, z));
+               rows.push(row);
+
+               row = document.createElement('tr');
+               add_3cell(row, 'Time per possession (goals)', 'name');
+               add_3cell(row, '');
+               add_3cell_ci(row, make_ci_duration(globals.time_for_goal, z));
+               add_3cell(row, '');
+               add_3cell_ci(row, make_ci_duration(globals.their_time_for_goal, z));
+               rows.push(row);
+
+               row = document.createElement('tr');
+               add_3cell(row, 'Time per possession (turnovers)', 'name');
+               add_3cell(row, '');
+               add_3cell_ci(row, make_ci_duration(globals.time_for_turnover, z));
+               add_3cell(row, '');
+               add_3cell_ci(row, make_ci_duration(globals.their_time_for_turnover, z));
+               rows.push(row);
+       }
+
        return rows;
 }