enh: added group statistcs calculations

Author: akeshavan

dmriprepViewer/src/components/Bucket.vue

@@ -4,6 +4,11 @@
       <div class="col-md-3">
         <div class="container">
         <h3>{{bucket}} ({{manifestEntries.length}})</h3>
+          <div class="mb-3">
+            <b-button v-if="!statsReady" @click="getAllReports">
+              Compute Statistics
+            </b-button>
+          </div>
           <b-nav vertical pills class="w-100">
             <!-- <b-nav-item active>Active</b-nav-item> -->
             <b-nav-item v-for="(subject, index) in manifestEntries"
@@ -18,6 +23,9 @@
         <h1 v-if="manifestEntries.length">
           {{manifestEntries[currentReportIdx].split('/')[0]}}
         </h1>
+        <div v-if="statsReady">
+          <GroupStats :data="allReports" :individual="currentReport.eddy_quad"/>
+        </div>
         <div v-if="ready">
           <report
             :reportProp="currentReport"
@@ -35,6 +43,7 @@
 import axios from 'axios';
 import _ from 'lodash';
 import Report from './Report';
+import GroupStats from './GroupStats';
 
 export default {
   name: 'bucket',
@@ -44,12 +53,33 @@ export default {
       currentReportIdx: 0,
       currentReport: {},
       ready: false,
+      allReports: [],
+      statsReady: false,
     };
   },
   components: {
     Report,
+    GroupStats,
   },
   methods: {
+    getReport(r) {
+      return axios.get(`https://s3-us-west-2.amazonaws.com/${this.bucket}/${r}`);
+    },
+    /**
+    *
+    */
+    async getAllReports() {
+      this.statsReady = false;
+      const reports = await _.map(this.manifestEntries, m => this.getReport(m));
+      _.map(reports, (r) => {
+        r.then((resp) => {
+          if (resp.data.eddy_quad) {
+            this.allReports.push(resp.data.eddy_quad);
+          }
+        });
+      });
+      this.statsReady = true;
+    },
     /**
     * XML parser for pubmed query returns.
     */

dmriprepViewer/src/components/GroupStats.vue

@@ -0,0 +1,94 @@
+<template>
+  <div>
+    <div class="mb-3">
+      <p class="lead">
+        Comparison to Group Eddy_Quad Statistics
+      </p>
+      <p>Absolute Motion (z-score)</p>
+      <b-progress class="w-50 mx-auto">
+        <b-progress-bar :value="scaleZ(individual_abs_mot_z)"
+         :variant="getColor(0, individual_abs_mot_z)">
+          z = {{individual_abs_mot_z}}
+        </b-progress-bar>
+      </b-progress>
+    </div>
+    <div class="mb-3">
+      <p>Relative Motion (z-score)</p>
+      <b-progress class="w-50 mx-auto">
+        <b-progress-bar :value="scaleZ(individual_rel_mot_z)"
+         :variant="getColor(0, individual_rel_mot_z)">
+          z = {{individual_rel_mot_z}}
+        </b-progress-bar>
+      </b-progress>
+    </div>
+
+  </div>
+</template>
+
+<script>
+// import _ from 'lodash';
+const d3 = require('d3');
+
+export default {
+  name: 'GroupStats',
+  props: ['data', 'individual'],
+  data() {
+    return {
+
+    };
+  },
+  methods: {
+    scaleZ(val) {
+      const scaler = d3.scaleLinear().range([0, 100]).domain([-3, 3]);
+      return scaler(val);
+    },
+    getColor(direction, zScore) {
+      if (direction) {
+        console.log(zScore);
+        // positive is good
+        if (zScore <= -1) {
+          return 'danger';
+        } else if (zScore >= 1) {
+          return 'success';
+        }
+      }
+      console.log('here', zScore);
+      if (zScore <= -1) {
+        return 'success';
+      } else if (zScore >= 1) {
+        return 'danger';
+      }
+      return 'primary';
+    },
+  },
+  computed: {
+    mean_group_abs_mot() {
+      return d3.mean(this.data, d => d.qc_mot_abs);
+    },
+    std_group_abs_mot() {
+      return d3.deviation(this.data, d => d.qc_mot_abs);
+    },
+    individual_abs_mot_z() {
+      if (this.individual) {
+        return (this.individual.qc_mot_abs - this.mean_group_abs_mot) / this.std_group_abs_mot;
+      }
+      return null;
+    },
+    mean_group_rel_mot() {
+      return d3.mean(this.data, d => d.qc_mot_rel);
+    },
+    std_group_rel_mot() {
+      return d3.deviation(this.data, d => d.qc_mot_rel);
+    },
+    individual_rel_mot_z() {
+      if (this.individual) {
+        return (this.individual.qc_mot_rel - this.mean_group_rel_mot) / this.std_group_rel_mot;
+      }
+      return null;
+    },
+  },
+};
+</script>
+
+<style>
+</style>