@sjcrh/proteinpaint-server 2.145.0 → 2.145.1-0

package/dataset/termdb.test.js

@@ -314,7 +314,8 @@ function termdb_test_default() {
       },
       rnaseqGeneCount: {
         storage_type: "HDF5",
-        file: "files/hg38/TermdbTest/rnaseq/TermdbTest.geneCounts.h5"
+        newformat: true,
+        file: "files/hg38/TermdbTest/rnaseq/TermdbTest.geneCounts.new.h5"
       },
       singleCell: {
         samples: {

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@sjcrh/proteinpaint-server",
-  "version": "2.145.0",
+  "version": "2.145.1-0",
   "type": "module",
   "description": "a genomics visualization tool for exploring a cohort's genotype and phenotype data",
   "main": "src/app.js",
@@ -62,10 +62,10 @@
   "dependencies": {
     "@sjcrh/augen": "2.143.0",
     "@sjcrh/proteinpaint-python": "2.144.0",
-    "@sjcrh/proteinpaint-r": "2.137.2-0",
-    "@sjcrh/proteinpaint-rust": "2.145.0",
-    "@sjcrh/proteinpaint-shared": "2.143.0",
-    "@sjcrh/proteinpaint-types": "2.145.0",
+    "@sjcrh/proteinpaint-r": "2.145.1-0",
+    "@sjcrh/proteinpaint-rust": "2.145.1-0",
+    "@sjcrh/proteinpaint-shared": "2.145.1-0",
+    "@sjcrh/proteinpaint-types": "2.145.1-0",
     "@types/express": "^5.0.0",
     "@types/express-session": "^1.18.1",
     "better-sqlite3": "^9.4.1",

package/routes/termdb.DE.js

@@ -204,7 +204,12 @@ async function run_DE(param, ds, term_results, term_results2) {
   const sample_size_limit = 8;
   if (group1names.length <= sample_size_limit && group2names.length <= sample_size_limit || param.method == "edgeR" || param.method == "limma") {
     const time12 = (/* @__PURE__ */ new Date()).valueOf();
-    const result2 = JSON.parse(await run_R("edge.R", JSON.stringify(expression_input)));
+    let result2;
+    if (q.newformat) {
+      result2 = JSON.parse(await run_R("edge_newh5.R", JSON.stringify(expression_input)));
+    } else {
+      result2 = JSON.parse(await run_R("edge.R", JSON.stringify(expression_input)));
+    }
     mayLog("Time taken to run edgeR:", formatElapsedTime(Date.now() - time12));
     param.method = "edgeR";
     const ql_imagePath = path.join(serverconfig.cachedir, result2.edgeR_ql_image_name[0]);
@@ -254,15 +259,30 @@ async function validate_query_rnaseqGeneCount(ds) {
     if (ds.queries.rnaseqGeneCount.storage_type == "text") {
       samples = (await get_header_txt(q.file, null)).split("	").slice(4);
     } else if (ds.queries.rnaseqGeneCount.storage_type == "HDF5") {
-      const get_samples_from_hdf5 = {
-        input_file: q.file,
-        data_type: "get_samples"
-      };
-      const time1 = (/* @__PURE__ */ new Date()).valueOf();
-      const result = await run_rust("DEanalysis", JSON.stringify(get_samples_from_hdf5));
-      const time2 = (/* @__PURE__ */ new Date()).valueOf();
-      mayLog("Time taken to query gene expression:", time2 - time1, "ms");
-      samples = result.split(",");
+      if (q.newformat) {
+        const get_samples_from_hdf5 = {
+          hdf5_file: q.file,
+          validate: true
+        };
+        const time1 = (/* @__PURE__ */ new Date()).valueOf();
+        const result = await run_rust("readH5", JSON.stringify(get_samples_from_hdf5));
+        const time2 = (/* @__PURE__ */ new Date()).valueOf();
+        mayLog("Time taken to query gene expression:", time2 - time1, "ms");
+        const vr = JSON.parse(result);
+        if (vr.status !== "success") throw vr.message;
+        if (!Array.isArray(vr.samples)) throw "HDF5 file has no samples, please check file.";
+        samples = vr.samples;
+      } else {
+        const get_samples_from_hdf5 = {
+          input_file: q.file,
+          data_type: "get_samples"
+        };
+        const time1 = (/* @__PURE__ */ new Date()).valueOf();
+        const result = await run_rust("DEanalysis", JSON.stringify(get_samples_from_hdf5));
+        const time2 = (/* @__PURE__ */ new Date()).valueOf();
+        mayLog("Time taken to query gene expression:", time2 - time1, "ms");
+        samples = result.split(",");
+      }
     } else throw "unknown storage type:" + ds.queries.rnaseqGeneCount.storage_type;
     q.allSampleSet = new Set(samples);
     const unknownSamples = [];

package/routes/termdb.boxplot.js

@@ -140,7 +140,7 @@ function setUncomputableValues(values) {
 function parseValues(q, data, sampleType, isLog, overlayTerm, divideTerm) {
   const chart2plot2values = /* @__PURE__ */ new Map();
   const uncomputableValues = {};
-  let absMin = null, absMax = null;
+  let absMin = Infinity, absMax = -Infinity;
   for (const val of Object.values(data.samples)) {
     const value = val[q.tw.$id];
     if (!Number.isFinite(value?.value)) continue;
@@ -175,8 +175,8 @@ function parseValues(q, data, sampleType, isLog, overlayTerm, divideTerm) {
     if (!plot2values.has(plot)) plot2values.set(plot, []);
     const values = plot2values.get(plot);
     values.push(value.value);
-    if (absMin === null || value.value < absMin) absMin = value.value;
-    if (absMax === null || value.value > absMax) absMax = value.value;
+    if (value.value < absMin) absMin = value.value;
+    if (value.value > absMax) absMax = value.value;
   }
   return { absMax, absMin, chart2plot2values, uncomputableValues };
 }