@dra2020/district-analytics 10.2.6 → 10.3.0

package/dist/district-analytics.js

@@ -122,7 +122,6 @@ const analyze_1 = __webpack_require__(/*! ./analyze */ "./src/analyze.ts");
 const score_1 = __webpack_require__(/*! ./score */ "./src/score.ts");
 const results_1 = __webpack_require__(/*! ./results */ "./src/results.ts");
 const results_2 = __webpack_require__(/*! ./results */ "./src/results.ts");
-// 11-03-2020 - Added for racially polarized voting analysis
 const minority_1 = __webpack_require__(/*! ./minority */ "./src/minority.ts");
 const Poly = __importStar(__webpack_require__(/*! @dra2020/poly */ "@dra2020/poly"));
 const D = __importStar(__webpack_require__(/*! ./_data */ "./src/_data.ts"));
@@ -678,6 +677,7 @@ class Districts {
         let pacificPts = [];
         let asianPts = [];
         let nativePts = [];
+        let dataDump = []; // For testing
         // Gather [demographic %, D %] points for the selected district ID
         let i = districtID;
         // HACK - Because "this" gets ghosted inside the forEach loop below
@@ -713,37 +713,20 @@ class Districts {
                                 const pctPacific = outerThis._session.features.fieldForFeature(f, "VAP" /* VAP */, "PacC" /* PacificPop */) / totalVAP;
                                 const pctAsian = outerThis._session.features.fieldForFeature(f, "VAP" /* VAP */, "AsnC" /* AsianPop */) / totalVAP;
                                 const pctNative = outerThis._session.features.fieldForFeature(f, "VAP" /* VAP */, "NatC" /* NativePop */) / totalVAP;
-                                let bGoodPoints = true;
-                                // bGoodPoints = (Number.isNaN(pctDem)) ? false : true;
-                                // bGoodPoints = (Number.isNaN(pctWhite)) ? false : true;
-                                // bGoodPoints = (Number.isNaN(pctMinority)) ? false : true;
-                                // bGoodPoints = (Number.isNaN(pctBlack)) ? false : true;
-                                // bGoodPoints = (Number.isNaN(pctHispanic)) ? false : true;
-                                // bGoodPoints = (Number.isNaN(pctPacific)) ? false : true;
-                                // bGoodPoints = (Number.isNaN(pctAsian)) ? false : true;
-                                // bGoodPoints = (Number.isNaN(pctNative)) ? false : true;
-                                if (bGoodPoints) {
-                                    ids.push(geoID);
-                                    whitePts.push({ x: pctWhite, y: pctDem });
-                                    minorityPts.push({ x: pctMinority, y: pctDem });
-                                    blackPts.push({ x: pctBlack, y: pctDem });
-                                    hispanicPts.push({ x: pctHispanic, y: pctDem });
-                                    pacificPts.push({ x: pctPacific, y: pctDem });
-                                    asianPts.push({ x: pctAsian, y: pctDem });
-                                    nativePts.push({ x: pctNative, y: pctDem });
-                                }
-                                else {
-                                    if (bLog) {
-                                        console.log("RPV: NaN value for feature ", geoID);
-                                        console.log("* pctDem: ", pctDem);
-                                        console.log("* pctWhite: ", pctWhite);
-                                        console.log("* pctMinority: ", pctMinority);
-                                        console.log("* pctBlack: ", pctBlack);
-                                        console.log("* pctHispanic: ", pctHispanic);
-                                        console.log("* pctPacific: ", pctPacific);
-                                        console.log("* pctAsian: ", pctAsian);
-                                        console.log("* pctNative: ", pctNative);
-                                    }
+                                ids.push(geoID);
+                                whitePts.push({ x: pctWhite, y: pctDem });
+                                minorityPts.push({ x: pctMinority, y: pctDem });
+                                blackPts.push({ x: pctBlack, y: pctDem });
+                                hispanicPts.push({ x: pctHispanic, y: pctDem });
+                                pacificPts.push({ x: pctPacific, y: pctDem });
+                                asianPts.push({ x: pctAsian, y: pctDem });
+                                nativePts.push({ x: pctNative, y: pctDem });
+                                // Dump test data based on config switch
+                                const s = outerThis._session;
+                                const bDump = (('dump' in s.config) && s.config.dump) ? true : false;
+                                if (bDump) {
+                                    const row = [pctWhite, pctMinority, pctBlack, pctHispanic, pctPacific, pctAsian, pctNative, pctDem];
+                                    dataDump.push(row);
                                 }
                             }
                             else {
@@ -759,6 +742,7 @@ class Districts {
                 }
             });
             // ... but only keep points for the demographics that are going to be analyzed
+            dataDump = [];
             return {
                 ids: ids,
                 white: whitePts,
@@ -1445,7 +1429,7 @@ var __importDefault = (this && this.__importDefault) || function (mod) {
     return (mod && mod.__esModule) ? mod : { "default": mod };
 };
 Object.defineProperty(exports, "__esModule", { value: true });
-const U = __importStar(__webpack_require__(/*! ./utils */ "./src/utils.ts"));
+const RPV = __importStar(__webpack_require__(/*! @dra2020/racial-voting */ "@dra2020/racial-voting"));
 const majority_minority_json_1 = __importDefault(__webpack_require__(/*! ../static/majority-minority.json */ "./static/majority-minority.json"));
 const vra5_preclearance_json_1 = __importDefault(__webpack_require__(/*! ../static/vra5-preclearance.json */ "./static/vra5-preclearance.json"));
 // TODO - 2020: Update/revise this, when the update comes out in September:
@@ -1467,92 +1451,152 @@ function getVRASection5(s) {
     return stateVRAPre;
 }
 exports.getVRASection5 = getVRASection5;
-// 11-03-2020 - Added for racially polarized voting analysis
-// Analyze the degree of racially polarized voting for a district
+// RPV - pulled into a separate component 11-17-2020
 function doAnalyzeRacialPolarization(s, districtID, groups, bLog = false) {
-    // Either take take the minorities to analyze in as selections -or- infer them
-    // from the % VAP/CVAP for the district (not the state).
-    if (groups == undefined) {
-        const RPV_THRESHOLD = 0.35; // Using ~ the Bethune–Hill threshold
-        groups = {
-            white: true,
-            minority: false,
-            black: (s.districts.table.blackPct[districtID] > RPV_THRESHOLD) ? true : false,
-            hispanic: (s.districts.table.hispanicPct[districtID] > RPV_THRESHOLD) ? true : false,
-            pacific: (s.districts.table.pacificPct[districtID] > RPV_THRESHOLD) ? true : false,
-            asian: (s.districts.table.asianPct[districtID] > RPV_THRESHOLD) ? true : false,
-            native: (s.districts.table.nativePct[districtID] > RPV_THRESHOLD) ? true : false
-        };
-    }
+    // At least one minority group must be specified
     if (!(groups.black || groups.hispanic || groups.pacific || groups.asian || groups.native))
         return undefined;
     const points = s.districts.extractVotesByDemographic(districtID, groups, bLog);
-    if (points === undefined)
-        return undefined;
-    // Make sure there are enough points to do a regression
-    if (points.white.length <= 10)
-        return undefined;
-    // At this point, at least one single-race/ethnicity minority demographic needs
-    // to be analyzed along with white.
-    const result = {
-        ids: points.ids,
-        white: characterizeDemographicVoting(points.white),
-        minority: groups.minority ? characterizeDemographicVoting(points.minority) : undefined,
-        black: groups.black ? characterizeDemographicVoting(points.black) : undefined,
-        hispanic: groups.hispanic ? characterizeDemographicVoting(points.hispanic) : undefined,
-        pacific: groups.pacific ? characterizeDemographicVoting(points.pacific) : undefined,
-        asian: groups.asian ? characterizeDemographicVoting(points.asian) : undefined,
-        native: groups.native ? characterizeDemographicVoting(points.native) : undefined
-    };
-    return result;
+    return RPV.analyzeRacialVoting(points, districtID, groups);
 }
 exports.doAnalyzeRacialPolarization = doAnalyzeRacialPolarization;
+/* TODO - RPV: DELETE
+// 11-03-2020 - Added for racially polarized voting analysis
+// Analyze the degree of racially polarized voting for a district
+export function doAnalyzeRacialPolarization(s: AnalyticsSession, districtID: number, groups?: T.MinorityFilter, bLog: boolean = false): T.RPVAnalysis | undefined
+{
+  // Either take take the minorities to analyze in as selections -or- infer them
+  // from the % VAP/CVAP for the district (not the state).
+  if (groups == undefined)
+  {
+    const RPV_THRESHOLD = 0.35;  // Using ~ the Bethune–Hill threshold
+    groups = {
+      white: true,
+      minority: false,
+      black: (s.districts.table.blackPct[districtID] > RPV_THRESHOLD) ? true : false,
+      hispanic: (s.districts.table.hispanicPct[districtID] > RPV_THRESHOLD) ? true : false,
+      pacific: (s.districts.table.pacificPct[districtID] > RPV_THRESHOLD) ? true : false,
+      asian: (s.districts.table.asianPct[districtID] > RPV_THRESHOLD) ? true : false,
+      native: (s.districts.table.nativePct[districtID] > RPV_THRESHOLD) ? true : false
+    };
+  }
+
+  if (!(groups.black || groups.hispanic || groups.pacific || groups.asian || groups.native)) return undefined;
+
+  const points = s.districts.extractVotesByDemographic(districtID, groups, bLog);
+
+  if (points === undefined) return undefined;
+
+  // Make sure there are enough points to do a regression
+  if (points.white.length <= 10) return undefined;
+
+  // At this point, at least one single-race/ethnicity minority demographic needs
+  // to be analyzed along with white.
+
+  const result: T.RPVAnalysis = {
+    ids: points.ids,
+    white: characterizeDemographicVoting(points.white),
+    minority: groups.minority ? characterizeDemographicVoting(points.minority) : undefined,
+    black: groups.black ? characterizeDemographicVoting(points.black) : undefined,
+    hispanic: groups.hispanic ? characterizeDemographicVoting(points.hispanic) : undefined,
+    pacific: groups.pacific ? characterizeDemographicVoting(points.pacific) : undefined,
+    asian: groups.asian ? characterizeDemographicVoting(points.asian) : undefined,
+    native: groups.native ? characterizeDemographicVoting(points.native) : undefined
+  }
+
+  return result;
+}
+
 // Cloned from dra-score
-function calcProportionalDistricts(proportion, nDistricts) {
-    const roundUp = 0.0;
-    const fractional = proportion * nDistricts;
-    const integral = Math.round(fractional + roundUp);
-    return integral;
+function calcProportionalDistricts(proportion: number, nDistricts: number): number
+{
+  const roundUp = 0.0;
+  const fractional = proportion * nDistricts;
+  const integral = Math.round(fractional + roundUp);
+
+  return integral;
 }
+
 // https://trentrichardson.com/compute-linear-regressions-in-javascript.html
-function linearRegression(points) {
-    let result = {};
-    const n = points.length;
-    let sum_x = 0;
-    let sum_y = 0;
-    let sum_xy = 0;
-    let sum_xx = 0;
-    let sum_yy = 0;
-    for (let i = 0; i < n; i++) {
-        const pt = points[i];
-        sum_x += pt.x;
-        sum_y += pt.y;
-        sum_xy += (pt.x * pt.y);
-        sum_xx += (pt.x * pt.x);
-        sum_yy += (pt.y * pt.y);
-    }
-    // y = mx + b
-    result['m'] = (n * sum_xy - sum_x * sum_y) / (n * sum_xx - sum_x * sum_x);
-    result['b'] = (sum_y - result.m * sum_x) / n;
-    result['r2'] = Math.pow((n * sum_xy - sum_x * sum_y) / Math.sqrt((n * sum_xx - sum_x * sum_x) * (n * sum_yy - sum_y * sum_y)), 2);
-    return result;
+function linearRegression(points: T.Point[]): any
+{
+  let result: any = {};
+
+  const n: number = points.length;
+  let sum_x = 0;
+  let sum_y = 0;
+  let sum_xy = 0;
+  let sum_xx = 0;
+  let sum_yy = 0;
+
+  for (let i = 0; i < n; i++)
+  {
+    const pt = points[i];
+
+    sum_x += pt.x;
+    sum_y += pt.y;
+    sum_xy += (pt.x * pt.y);
+    sum_xx += (pt.x * pt.x);
+    sum_yy += (pt.y * pt.y);
+  }
+
+  // y = mx + b
+  result['m'] = (n * sum_xy - sum_x * sum_y) / (n * sum_xx - sum_x * sum_x);
+  result['b'] = (sum_y - result.m * sum_x) / n;
+  result['r2'] = Math.pow((n * sum_xy - sum_x * sum_y) / Math.sqrt((n * sum_xx - sum_x * sum_x) * (n * sum_yy - sum_y * sum_y)), 2);
+
+  return result;
 }
-function characterizeDemographicVoting(points) {
-    const { m, b, r2 } = linearRegression(points);
-    const maxVAPPct = 1.0; // 100%
-    const maxDemPct = (m * maxVAPPct) + b;
-    // NOTE - You can't trim the implied max D %, because it's a regression line!
-    // const maxDemPct: number = Math.max(0.0, Math.min(1.0, (m * maxVAPPct) + b));
-    const result = {
-        m: U.trim(m),
-        b: U.trim(b),
-        r2: U.trim(r2),
-        f1: U.trim(maxDemPct),
-        points: points
-    };
-    return result;
+
+function characterizeDemographicVoting(points: T.Point[]): T.RPVLine
+{
+  const { m, b, r2 }: any = linearRegression(points);
+  const maxVAPPct = 1.0;  // 100%
+  const maxDemPct: number = (m * maxVAPPct) + b;
+  // NOTE - You can't trim the implied max D %, because it's a regression line!
+  // const maxDemPct: number = Math.max(0.0, Math.min(1.0, (m * maxVAPPct) + b));
+
+  const result = {
+    m: U.trim(m),
+    b: U.trim(b),
+    r2: U.trim(r2),
+    sterr: U.trim(standardError(points, m, b)),
+    f1: U.trim(maxDemPct),
+    points: points
+  }
+
+  return result;
+}
+
+// TODO - RPV: Review with Stephen Ansolabehere
+// https://www.wikihow.com/Calculate-the-Standard-Error-of-Estimate
+// https://www.stat.cmu.edu/~hseltman/309/Book/chapter9.pdf
+// y = mx + b
+function standardError(points: T.Point[], m: number, b: number): number
+{
+  const n: number = points.length;
+  let sse: number = 0;
+
+  for (let i = 0; i < n; i++)
+  {
+    const pt = points[i];
+
+    const yPrime = yOnLine(pt.x, m, b);
+    const deltaY = pt.y - yPrime;
+
+    sse += Math.pow(deltaY, 2);
+  }
+
+  const sigma = Math.sqrt(sse / n);
+
+  return sigma;
 }
 
+// For interpolating points on a line
+const yOnLine = (x: number, m: number, b: number): number => { return (m * x) + b; }
+const xOnLine = (y: number, m: number, b: number): number => { return (y - b) / m; }
+*/
+
 
 /***/ }),
 
@@ -2667,6 +2711,17 @@ module.exports = require("@dra2020/poly");
 
 /***/ }),
 
+/***/ "@dra2020/racial-voting":
+/*!*****************************************!*\
+  !*** external "@dra2020/racial-voting" ***!
+  \*****************************************/
+/*! no static exports found */
+/***/ (function(module, exports) {
+
+module.exports = require("@dra2020/racial-voting");
+
+/***/ }),
+
 /***/ "@dra2020/util":
 /*!********************************!*\
   !*** external "@dra2020/util" ***!