Issue/317 si stat

R_work/include/statistics.R

@@ -117,6 +117,7 @@ calcRMSE    = function(d){ return( sqrt(calcMSE(d)) ); }
 calcESTDEV    = function(d){ return( calcStdDev( calcME(d) * d$total, calcMSE(d) * d$total, d$total) ); }
 calcBCMSE    = function(d){ return( calcMSE(d) - (d$fbar - d$obar)^2 ); }
 calcBCRMSE    = function(d){ return( sqrt(calcBCMSE(d)) ); }
+calcSI    = function(d){return( calcRMSE(d)/ calcOBAR(d));}
 
 # SSVAR stat calculations
 calcSSVAR_FBAR    = function(d){ return( calcFBAR(d) ); }

java/edu/ucar/metviewer/MVUtil.java

@@ -546,6 +546,7 @@ public class MVUtil {
     statsCnt.put("MSESS", new String[]{"bc", SL1L2});
     statsCnt.put("RMSFA", new String[]{"bc", SAL1L2});
     statsCnt.put("RMSOA", new String[]{"bc", SAL1L2});
+    statsCnt.put("SI", new String[]{"bc", SL1L2});
   }
 
   static {

sql/mv_mysql.sql

@@ -395,6 +395,10 @@ CREATE TABLE line_data_cnt
     anom_corr_uncntr_bcl DOUBLE DEFAULT -9999,
     anom_corr_uncntr_bcu DOUBLE DEFAULT -9999,
 
+    si                   DOUBLE DEFAULT -9999,
+    si_bcl               DOUBLE DEFAULT -9999,
+    si_bcu               DOUBLE DEFAULT -9999,
+
 
     CONSTRAINT line_data_cnt_data_file_id_pk
         FOREIGN KEY (data_file_id)

webapp/metviewer/js/metviewer_common.js

@@ -202,6 +202,7 @@ value_to_desc_map['FGBAR'] = 'Mean of absolute value of forecast gradients';
 value_to_desc_map['MGBAR'] = 'Mean of maximum of absolute values of forecast and observed gradients';
 value_to_desc_map['EGBAR'] = 'Mean of absolute value of forecast minus observed gradients';
 value_to_desc_map['S1'] = 'S1 score';
+value_to_desc_map['SI'] = 'Scatter Index';
 value_to_desc_map['S1_OG'] = 'S1 score with respect to observed gradient';
 value_to_desc_map['FGOG_RATIO'] = 'Ratio of forecast and observed gradients';
 value_to_desc_map['SSVAR_SPREAD'] = 'RHIST_SPREAD is preferred';