missing epsilon correction in R implementation to fix MassBank/MassBa…

…nk-data#248
ssmehta · Nov 14, 2023 · 030b4a5 · 030b4a5
1 parent 174fbf2
commit 030b4a5
Showing 1 changed file with 1 addition and 2 deletions.
diff --git a/splashR/R/getSplash.R b/splashR/R/getSplash.R
@@ -96,11 +96,10 @@ getBlockHist <- function(peaks, histBase, histLength, binSize) {
     summedintensities <- tapply(peaks[,2], binindex, sum)
     wrappedbinindex <- unique(binindex) %% histLength
     wrappedintensities <- tapply(summedintensities, wrappedbinindex, sum)
-    normalisedintensities <- as.integer(wrappedintensities/max(wrappedintensities)*(histBase-1))
+    normalisedintensities <- as.integer(EPS_CORRECTION + (histBase-1) * wrappedintensities / max(wrappedintensities))
 
     wrappedhist[sort(unique(wrappedbinindex))+1] <- normalisedintensities
     paste(integer2base36(wrappedhist), collapse="")
-
 }
 
 getSplash <- function(peaks) {