Adjust formula for last change stability score.

The formulas for LCSS and LCSS1 are now back to the original:

LCSS1 = LCMS * LCGSS * LCLSS * LCDS1
LCSS = LCSS1 * LCDS2

A simple product turned out to be the most robust way to avoid assigning
significant scores to certain outlier cases after all.

This patch also redefines the global separation score slightly to make
it more robust (the values before the last change are now compared against
all values at or after the last change, rather than just the last change
value itself).

Finally, the default maximum durability score is adjusted from 30 to 10
(i.e. 10 seems to be sufficient to establish stability).

4 files changed, 7 insertions, 7 deletions

diff --git a/web/analysis/global.js b/web/analysis/global.js
index 1a4148d..f8f1cc1 100644
--- a/web/analysis/global.js
+++ b/web/analysis/global.js
@@ -69,7 +69,7 @@ function changeMagnitudeScore(change) {
 // ### 2 B DOCUMENTED!
 // NOTE: This function is currently duplicated elsewhere in Python!
 function qualityScore(lsd, ni, nz, nc, mdrse) {
-    var maxBadSnapshots = 30; // experimental; maybe use max durability score?
+    var maxBadSnapshots = 10; // experimental; maybe use max durability score?
     var maxSampleSize = 5;
     var maxLSD = maxBadSnapshots;
     var maxNI = maxBadSnapshots * maxSampleSize;
diff --git a/web/analysis/main.js b/web/analysis/main.js
index 30b5401..735c5d8 100644
--- a/web/analysis/main.js
+++ b/web/analysis/main.js
@@ -565,7 +565,7 @@ $(document).ready(function() {
     $("#max_rank_size option[value='10']").attr("selected", true);
     $("#diff_tol option[value='1.1']").attr("selected", true);
     $("#dur_tol_min option[value='3']").attr("selected", true);
-    $("#dur_tol_max option[value='50']").attr("selected", true);
+    $("#dur_tol_max option[value='10']").attr("selected", true);
 
     fetchContexts();
 });
diff --git a/web/analysis/tsbmbody.js b/web/analysis/tsbmbody.js
index 2c14659..8706208 100644
--- a/web/analysis/tsbmbody.js
+++ b/web/analysis/tsbmbody.js
@@ -110,7 +110,7 @@ function tooltipText_lcms() {
 function tooltipText_lcss() {
     return "Stability score for the last significant change:<br /><br />" +
         "&nbsp;&nbsp;&nbsp;&nbsp;" +
-        "(LCMS + LCGSS + LCLSS + (LCDS1 * LCDS2)) / 4<br /><br />" +
+        "LCMS * LCGSS * LCLSS * LCDS1 * LCDS2<br /><br />" +
         "The higher this score, the higher the likelihood that the last " +
         "significant change is or will become permanent.";
 }
@@ -119,10 +119,10 @@ function tooltipText_lcss1() {
     return "Stability score for the last significant change that does not " +
         "consider the history after the latter:<br /><br />" +
         "&nbsp;&nbsp;&nbsp;&nbsp;" +
-        "(LCMS + LCGSS + LCLSS + LCDS1) / 4<br /><br />" +
+        "LCMS * LCGSS * LCLSS * LCDS1<br /><br />" +
         "The higher this score, the higher the likelihood that the last " +
         "signicifant change is or will become permanent, but since LCDS2 is " +
-        "omitted from the sum, a high LCSS1 is more likely to be caused " +
+        "omitted from the product, a high LCSS1 is more likely to be caused " +
         "by an outlier than a high LCSS!";
 }
 
diff --git a/web/analysis/tsstats.js b/web/analysis/tsstats.js
index bba859f..93e58ab 100644
--- a/web/analysis/tsstats.js
+++ b/web/analysis/tsstats.js
@@ -124,10 +124,10 @@ function populatePbmTable(tableSel, data) {
                 lcds1 = stats.lc_dur1_score;
                 lcds2 = stats.lc_dur2_score;
 
-                lcss = (lcms + lcgss + lclss + (lcds1 * lcds2)) / 4.0;
+                lcss = lcms * lcgss * lclss * lcds1 * lcds2;
                 rbdy += "<td>" + lcss + "</td>";
 
-                lcss1 = (lcms + lcgss + lclss + lcds1) / 4.0;
+                lcss1 = lcms * lcgss * lclss * lcds1;
                 rbdy += "<td>" + lcss1 + "</td>";
 
                 rbdy += "<td>" + lcgss.toFixed(4) + "</td>";