8  Calling MALDIquant from Clojure

A guide to using ClojisR to call R’s MALDIquant package from Clojure. This is useful for accessing MALDIquant features not yet in Ripple (peak alignment, spectrum averaging, feature matrices) and for cross-checking results against the reference implementation.

This is a companion to the NeuroKit2 comparison notebook, which demonstrates the same pattern using libpython-clj for Python interop.

Requires: R with MALDIquant installed (install.packages("MALDIquant") or ./setup_r.sh).

(ns ripple-book.maldiquant-interop
  (:require
   ;; Table processing (https://scicloj.github.io/tablecloth/):
   [tablecloth.api :as tc]
   ;; Interactive plotting via Plotly (https://scicloj.github.io/tableplot/):
   [scicloj.tableplot.v1.plotly :as plotly]
   ;; R interop (https://github.com/scicloj/clojisr):
   [clojisr.v1.r :as r :refer [r r->clj clj->r]]
   ;; Annotating kinds of visualizations (https://scicloj.github.io/kindly-noted/):
   [scicloj.kindly.v4.kind :as kind]))

Setup

ClojisR starts an R session automatically on first use. We just need to load the MALDIquant package.

(r "library(MALDIquant)")
[1] "MALDIquant" "Rserve"     "stats"      "graphics"   "grDevices" 
[6] "utils"      "datasets"   "methods"    "base"      

ClojisR Basics

ClojisR provides three core functions for R interop:

Function Purpose
(r ...) Evaluate an R expression (string or quoted form)
(r->clj x) Convert an R object to Clojure data
(clj->r x) Convert Clojure data to an R object

Quoted forms use backtick syntax. Clojure values are spliced in with tilde (~).

Evaluating R expressions

As a string:

(r "1 + 1")
[1] 2

As a quoted form (preferred — enables splicing Clojure values):

(r `(+ 1 1))
[1] 2

Converting R results to Clojure

(r->clj (r `(seq 1 5)))
[1 2 3 4 5]

Splicing Clojure values into R

Use tilde (~) to inject Clojure values into quoted R expressions:

(def my-data [10.0 20.0 30.0 40.0 50.0])
(r->clj (r `(mean ~my-data)))
[30.0]

Vectors, doubles, and strings all splice naturally. For keyword arguments in R functions, use Clojure keyword syntax:

(r->clj (r `(paste "hello" "world" :sep "-")))
["hello-world"]

Working with MALDIquant Objects

MALDIquant’s core data type is MassSpectrum — an S4 object with mass and intensity slots. We create one from Clojure vectors.

(def masses [100.0 200.0 300.0 400.0 500.0 600.0 700.0 800.0 900.0])
(def intensities [0.1 0.3 0.8 1.5 0.9 0.2 0.4 1.2 0.6])

Create a MassSpectrum object in R:

(r `(<- spec (createMassSpectrum :mass ~masses :intensity ~intensities)))
S4 class type            : MassSpectrum   
Number of m/z values     : 9              
Range of m/z values      : 100 - 900      
Range of intensity values: 1e-01 - 1.5e+00
Memory usage             : 1.227 KiB      

Extract data back to Clojure:

(def r-masses (r->clj (r `(mass spec))))
(def r-intensities (r->clj (r `(intensity spec))))
(= (mapv double masses) r-masses)
true

Single-Spectrum Preprocessing

MALDIquant ships with fiedler2009subset — 16 real MALDI-TOF mass spectra from a pancreatic cancer study (Fiedler et al. 2009). We use the first spectrum to demonstrate each preprocessing step.

(r "data('fiedler2009subset', package='MALDIquant')")
[1] "fiedler2009subset"
(def n-spectra (first (r->clj (r "length(fiedler2009subset)"))))
n-spectra
16
(def n-points (first (r->clj (r "length(fiedler2009subset[[1]])"))))
n-points
42388

Raw spectrum

(let [masses (r->clj (r "mass(fiedler2009subset[[1]])"))
      intensities (r->clj (r "intensity(fiedler2009subset[[1]])"))
      step 10
      idx (range 0 (count masses) step)]
  (-> (tc/dataset {:mass (mapv #(nth masses %) idx)
                   :intensity (mapv #(nth intensities %) idx)})
      (plotly/base {:=width 800 :=height 300 :=x :mass :=y :intensity
                    :=x-title "m/z (Da)" :=y-title "Intensity"
                    :=title "Raw spectrum (fiedler2009subset #1)"})
      (plotly/layer-line {:=mark-opacity 0.7})
      plotly/plot))

The raw spectrum shows a large baseline hump (chemical matrix signal) in the low-mass region, with analyte peaks sitting on top.

Trim to analysis range

(r "spec1 <- trim(fiedler2009subset[[1]], range=c(2000, 10000))")
S4 class type            : MassSpectrum                
Number of m/z values     : 34264                       
Range of m/z values      : 2000.137 - 9999.734         
Range of intensity values: 5 - 27518                   
Memory usage             : 411.156 KiB                 
Name                     : Pankreas_HB_L_061019_G10.M19
File                     : /data/set A - discovery leipzig/control/Pankreas_HB_L_061019_G10/0_m19/1/1SLin/fid

Variance stabilization — transformIntensity

(r "spec1 <- transformIntensity(spec1, method='sqrt')")
S4 class type            : MassSpectrum                
Number of m/z values     : 34264                       
Range of m/z values      : 2000.137 - 9999.734         
Range of intensity values: 2.236e+00 - 1.659e+02       
Memory usage             : 545 KiB                     
Name                     : Pankreas_HB_L_061019_G10.M19
File                     : /data/set A - discovery leipzig/control/Pankreas_HB_L_061019_G10/0_m19/1/1SLin/fid

Smoothing — smoothIntensity

(r "spec1 <- smoothIntensity(spec1, method='SavitzkyGolay', halfWindowSize=10)")
S4 class type            : MassSpectrum                
Number of m/z values     : 34264                       
Range of m/z values      : 2000.137 - 9999.734         
Range of intensity values: 2.91e+00 - 1.656e+02        
Memory usage             : 545 KiB                     
Name                     : Pankreas_HB_L_061019_G10.M19
File                     : /data/set A - discovery leipzig/control/Pankreas_HB_L_061019_G10/0_m19/1/1SLin/fid

Baseline removal — removeBaseline

(r "spec1 <- removeBaseline(spec1, method='SNIP', iterations=25)")
S4 class type            : MassSpectrum                
Number of m/z values     : 34264                       
Range of m/z values      : 2000.137 - 9999.734         
Range of intensity values: 0e+00 - 1.044e+02           
Memory usage             : 545 KiB                     
Name                     : Pankreas_HB_L_061019_G10.M19
File                     : /data/set A - discovery leipzig/control/Pankreas_HB_L_061019_G10/0_m19/1/1SLin/fid

Normalization — calibrateIntensity

(r "spec1 <- calibrateIntensity(spec1, method='TIC')")
S4 class type            : MassSpectrum                
Number of m/z values     : 34264                       
Range of m/z values      : 2000.137 - 9999.734         
Range of intensity values: 0e+00 - 7.579e-03           
Memory usage             : 545 KiB                     
Name                     : Pankreas_HB_L_061019_G10.M19
File                     : /data/set A - discovery leipzig/control/Pankreas_HB_L_061019_G10/0_m19/1/1SLin/fid

Verify TIC area is 1.0 (trapezoid rule):

(def real-area
  (first (r->clj (r "m <- mass(spec1)
                      i <- intensity(spec1)
                      sum(diff(m) * (i[-1] + i[-length(i)]) / 2)"))))
real-area
1.0

Peak detection — detectPeaks

(r "peaks1 <- detectPeaks(spec1, halfWindowSize=20, SNR=2, method='MAD')")
S4 class type            : MassPeaks                   
Number of m/z values     : 252                         
Range of m/z values      : 2012.127 - 9952.986         
Range of intensity values: 9.657e-05 - 7.579e-03       
Range of snr values      : 2.003 - 157.225             
Memory usage             : 15.523 KiB                  
Name                     : Pankreas_HB_L_061019_G10.M19
File                     : /data/set A - discovery leipzig/control/Pankreas_HB_L_061019_G10/0_m19/1/1SLin/fid
(def n-real-peaks (first (r->clj (r "length(peaks1)"))))
n-real-peaks
252

Visualizing the result

(let [spec-m (r->clj (r "mass(spec1)"))
      spec-i (r->clj (r "intensity(spec1)"))
      peak-m (r->clj (r "mass(peaks1)"))
      peak-i (r->clj (r "intensity(peaks1)"))
      step 5
      idx (range 0 (count spec-m) step)
      sub-m (mapv #(nth spec-m %) idx)
      sub-i (mapv #(nth spec-i %) idx)
      ns (count sub-m)
      np (count peak-m)]
  (-> (tc/dataset {:mass (vec (concat sub-m peak-m))
                   :intensity (vec (concat sub-i peak-i))
                   :series (vec (concat (repeat ns "spectrum")
                                        (repeat np "peaks")))})
      (plotly/base {:=width 800 :=height 400 :=x :mass :=y :intensity :=color :series
                    :=x-title "m/z (Da)" :=y-title "Normalized intensity"
                    :=title "Preprocessed spectrum with detected peaks"})
      (plotly/layer-line {:=mark-opacity 0.4})
      (plotly/layer-point {:=mark-size 4 :=mark-opacity 0.8})
      plotly/plot))

Multi-Sample Workflow

So far we have processed a single spectrum. In a real experiment you have many spectra from different samples. MALDIquant provides batch functions that operate on lists of spectra — this is the workflow described in Gibb & Strimmer (2016), Sections 11.3.3–11.3.7.

Step 1: Batch preprocessing

All single-spectrum functions accept a list of spectra:

(r "spectra <- fiedler2009subset
spectra <- trim(spectra, range=c(2000, 10000))
spectra <- transformIntensity(spectra, method='sqrt')
spectra <- smoothIntensity(spectra, method='SavitzkyGolay', halfWindowSize=10)
spectra <- removeBaseline(spectra, method='SNIP', iterations=25)
spectra <- calibrateIntensity(spectra, method='TIC')")
[[1]]
S4 class type            : MassSpectrum                
Number of m/z values     : 34264                       
Range of m/z values      : 2000.137 - 9999.734         
Range of intensity values: 0e+00 - 7.579e-03           
Memory usage             : 545 KiB                     
Name                     : Pankreas_HB_L_061019_G10.M19
File                     : /data/set A - discovery leipzig/control/Pankreas_HB_L_061019_G10/0_m19/1/1SLin/fid

[[2]]
S4 class type            : MassSpectrum                
Number of m/z values     : 34264                       
Range of m/z values      : 2000.137 - 9999.734         
Range of intensity values: 0e+00 - 7.353e-03           
Memory usage             : 545 KiB                     
Name                     : Pankreas_HB_L_061019_G10.M20
File                     : /data/set A - discovery leipzig/control/Pankreas_HB_L_061019_G10/0_m20/1/1SLin/fid

[[3]]
S4 class type            : MassSpectrum               
Number of m/z values     : 34264                      
Range of m/z values      : 2000.137 - 9999.734        
Range of intensity values: 0e+00 - 6.992e-03          
Memory usage             : 545 KiB                    
Name                     : Pankreas_HB_L_061019_H7.O14
File                     : /data/set A - discovery leipzig/control/Pankreas_HB_L_061019_H7/0_o14/1/1SLin/fid

[[4]]
S4 class type            : MassSpectrum               
Number of m/z values     : 34264                      
Range of m/z values      : 2000.137 - 9999.734        
Range of intensity values: 0e+00 - 6.641e-03          
Memory usage             : 545 KiB                    
Name                     : Pankreas_HB_L_061019_H7.P13
File                     : /data/set A - discovery leipzig/control/Pankreas_HB_L_061019_H7/0_p13/1/1SLin/fid

[[5]]
S4 class type            : MassSpectrum                
Number of m/z values     : 34264                       
Range of m/z values      : 2000.137 - 9999.734         
Range of intensity values: 0e+00 - 6.213e-03           
Memory usage             : 545 KiB                     
Name                     : Pankreas_HB_L_061019_F10.L19
File                     : /data/set A - discovery leipzig/tumor/Pankreas_HB_L_061019_F10/0_l19/1/1SLin/fid

[[6]]
S4 class type            : MassSpectrum                
Number of m/z values     : 34264                       
Range of m/z values      : 2000.137 - 9999.734         
Range of intensity values: 0e+00 - 6.144e-03           
Memory usage             : 545 KiB                     
Name                     : Pankreas_HB_L_061019_F10.L20
File                     : /data/set A - discovery leipzig/tumor/Pankreas_HB_L_061019_F10/0_l20/1/1SLin/fid

[[7]]
S4 class type            : MassSpectrum               
Number of m/z values     : 34264                      
Range of m/z values      : 2000.137 - 9999.734        
Range of intensity values: 0e+00 - 5.975e-03          
Memory usage             : 545 KiB                    
Name                     : Pankreas_HB_L_061019_F9.L17
File                     : /data/set A - discovery leipzig/tumor/Pankreas_HB_L_061019_F9/0_l17/1/1SLin/fid

[[8]]
S4 class type            : MassSpectrum               
Number of m/z values     : 34264                      
Range of m/z values      : 2000.137 - 9999.734        
Range of intensity values: 0e+00 - 6.002e-03          
Memory usage             : 545 KiB                    
Name                     : Pankreas_HB_L_061019_F9.L18
File                     : /data/set A - discovery leipzig/tumor/Pankreas_HB_L_061019_F9/0_l18/1/1SLin/fid

[[9]]
S4 class type            : MassSpectrum               
Number of m/z values     : 34264                      
Range of m/z values      : 2000.137 - 9999.734        
Range of intensity values: 0e+00 - 6.703e-03          
Memory usage             : 545.008 KiB                
Name                     : Pankreas_HB_L_061019_A6.A11
File                     : /data/set B - discovery heidelberg/control/Pankreas_HB_L_061019_A6/0_a11/1/1SLin/fid

[[10]]
S4 class type            : MassSpectrum               
Number of m/z values     : 34264                      
Range of m/z values      : 2000.137 - 9999.734        
Range of intensity values: 0e+00 - 6.895e-03          
Memory usage             : 545.008 KiB                
Name                     : Pankreas_HB_L_061019_A6.A12
File                     : /data/set B - discovery heidelberg/control/Pankreas_HB_L_061019_A6/0_a12/1/1SLin/fid

[[11]]
S4 class type            : MassSpectrum               
Number of m/z values     : 34264                      
Range of m/z values      : 2000.137 - 9999.734        
Range of intensity values: 0e+00 - 5.557e-03          
Memory usage             : 545.008 KiB                
Name                     : Pankreas_HB_L_061019_A8.A15
File                     : /data/set B - discovery heidelberg/control/Pankreas_HB_L_061019_A8/0_a15/1/1SLin/fid

[[12]]
S4 class type            : MassSpectrum               
Number of m/z values     : 34264                      
Range of m/z values      : 2000.137 - 9999.734        
Range of intensity values: 0e+00 - 5.85e-03           
Memory usage             : 545.008 KiB                
Name                     : Pankreas_HB_L_061019_A8.A16
File                     : /data/set B - discovery heidelberg/control/Pankreas_HB_L_061019_A8/0_a16/1/1SLin/fid

[[13]]
S4 class type            : MassSpectrum              
Number of m/z values     : 34264                     
Range of m/z values      : 2000.137 - 9999.734       
Range of intensity values: 0e+00 - 4.145e-03         
Memory usage             : 545.008 KiB               
Name                     : Pankreas_HB_L_061019_C4.F7
File                     : /data/set B - discovery heidelberg/tumor/Pankreas_HB_L_061019_C4/0_f7/1/1SLin/fid

[[14]]
S4 class type            : MassSpectrum              
Number of m/z values     : 34264                     
Range of m/z values      : 2000.137 - 9999.734       
Range of intensity values: 0e+00 - 4.126e-03         
Memory usage             : 545.008 KiB               
Name                     : Pankreas_HB_L_061019_C4.F8
File                     : /data/set B - discovery heidelberg/tumor/Pankreas_HB_L_061019_C4/0_f8/1/1SLin/fid

[[15]]
S4 class type            : MassSpectrum               
Number of m/z values     : 34264                      
Range of m/z values      : 2000.137 - 9999.734        
Range of intensity values: 0e+00 - 4.285e-03          
Memory usage             : 545.008 KiB                
Name                     : Pankreas_HB_L_061019_D9.G17
File                     : /data/set B - discovery heidelberg/tumor/Pankreas_HB_L_061019_D9/0_g17/1/1SLin/fid

[[16]]
S4 class type            : MassSpectrum               
Number of m/z values     : 34264                      
Range of m/z values      : 2000.137 - 9999.734        
Range of intensity values: 0e+00 - 4.019e-03          
Memory usage             : 545.008 KiB                
Name                     : Pankreas_HB_L_061019_D9.G18
File                     : /data/set B - discovery heidelberg/tumor/Pankreas_HB_L_061019_D9/0_g18/1/1SLin/fid

(def n-preprocessed (first (r->clj (r "length(spectra)"))))
n-preprocessed
16

Step 2: Peak detection on all spectra

(r "all_peaks <- detectPeaks(spectra, halfWindowSize=20, SNR=2, method='MAD')")
[[1]]
S4 class type            : MassPeaks                   
Number of m/z values     : 252                         
Range of m/z values      : 2012.127 - 9952.986         
Range of intensity values: 9.657e-05 - 7.579e-03       
Range of snr values      : 2.003 - 157.225             
Memory usage             : 15.523 KiB                  
Name                     : Pankreas_HB_L_061019_G10.M19
File                     : /data/set A - discovery leipzig/control/Pankreas_HB_L_061019_G10/0_m19/1/1SLin/fid

[[2]]
S4 class type            : MassPeaks                   
Number of m/z values     : 243                         
Range of m/z values      : 2002.301 - 9860.139         
Range of intensity values: 9.209e-05 - 7.353e-03       
Range of snr values      : 2.012 - 160.65              
Memory usage             : 15.312 KiB                  
Name                     : Pankreas_HB_L_061019_G10.M20
File                     : /data/set A - discovery leipzig/control/Pankreas_HB_L_061019_G10/0_m20/1/1SLin/fid

[[3]]
S4 class type            : MassPeaks                  
Number of m/z values     : 268                        
Range of m/z values      : 2013.284 - 9954.596        
Range of intensity values: 9.873e-05 - 6.992e-03      
Range of snr values      : 2.002 - 141.751            
Memory usage             : 15.898 KiB                 
Name                     : Pankreas_HB_L_061019_H7.O14
File                     : /data/set A - discovery leipzig/control/Pankreas_HB_L_061019_H7/0_o14/1/1SLin/fid

[[4]]
S4 class type            : MassPeaks                  
Number of m/z values     : 265                        
Range of m/z values      : 2018.063 - 9965.227        
Range of intensity values: 9.431e-05 - 6.641e-03      
Range of snr values      : 2.007 - 141.312            
Memory usage             : 15.828 KiB                 
Name                     : Pankreas_HB_L_061019_H7.P13
File                     : /data/set A - discovery leipzig/control/Pankreas_HB_L_061019_H7/0_p13/1/1SLin/fid

[[5]]
S4 class type            : MassPeaks                   
Number of m/z values     : 218                         
Range of m/z values      : 2011.114 - 9717.072         
Range of intensity values: 9.39e-05 - 6.213e-03        
Range of snr values      : 2.007 - 132.815             
Memory usage             : 14.727 KiB                  
Name                     : Pankreas_HB_L_061019_F10.L19
File                     : /data/set A - discovery leipzig/tumor/Pankreas_HB_L_061019_F10/0_l19/1/1SLin/fid

[[6]]
S4 class type            : MassPeaks                   
Number of m/z values     : 249                         
Range of m/z values      : 2002.301 - 9977.152         
Range of intensity values: 9.989e-05 - 6.144e-03       
Range of snr values      : 2 - 123.044                 
Memory usage             : 15.453 KiB                  
Name                     : Pankreas_HB_L_061019_F10.L20
File                     : /data/set A - discovery leipzig/tumor/Pankreas_HB_L_061019_F10/0_l20/1/1SLin/fid

[[7]]
S4 class type            : MassPeaks                  
Number of m/z values     : 231                        
Range of m/z values      : 2024.008 - 9954.596        
Range of intensity values: 9.836e-05 - 5.975e-03      
Range of snr values      : 2.003 - 121.669            
Memory usage             : 15.031 KiB                 
Name                     : Pankreas_HB_L_061019_F9.L17
File                     : /data/set A - discovery leipzig/tumor/Pankreas_HB_L_061019_F9/0_l17/1/1SLin/fid

[[8]]
S4 class type            : MassPeaks                  
Number of m/z values     : 221                        
Range of m/z values      : 2020.672 - 9819.476        
Range of intensity values: 9.469e-05 - 6.002e-03      
Range of snr values      : 2.001 - 126.825            
Memory usage             : 14.797 KiB                 
Name                     : Pankreas_HB_L_061019_F9.L18
File                     : /data/set A - discovery leipzig/tumor/Pankreas_HB_L_061019_F9/0_l18/1/1SLin/fid

[[9]]
S4 class type            : MassPeaks                  
Number of m/z values     : 226                        
Range of m/z values      : 2009.524 - 9928.206        
Range of intensity values: 9.046e-05 - 6.703e-03      
Range of snr values      : 2.003 - 148.4              
Memory usage             : 14.922 KiB                 
Name                     : Pankreas_HB_L_061019_A6.A11
File                     : /data/set B - discovery heidelberg/control/Pankreas_HB_L_061019_A6/0_a11/1/1SLin/fid

[[10]]
S4 class type            : MassPeaks                  
Number of m/z values     : 232                        
Range of m/z values      : 2020.672 - 9798.697        
Range of intensity values: 8.948e-05 - 6.895e-03      
Range of snr values      : 2.001 - 154.204            
Memory usage             : 15.062 KiB                 
Name                     : Pankreas_HB_L_061019_A6.A12
File                     : /data/set B - discovery heidelberg/control/Pankreas_HB_L_061019_A6/0_a12/1/1SLin/fid

[[11]]
S4 class type            : MassPeaks                  
Number of m/z values     : 203                        
Range of m/z values      : 2024.733 - 9987.795        
Range of intensity values: 9.276e-05 - 5.557e-03      
Range of snr values      : 2.025 - 121.319            
Memory usage             : 14.383 KiB                 
Name                     : Pankreas_HB_L_061019_A8.A15
File                     : /data/set B - discovery heidelberg/control/Pankreas_HB_L_061019_A8/0_a15/1/1SLin/fid

[[12]]
S4 class type            : MassPeaks                  
Number of m/z values     : 205                        
Range of m/z values      : 2008.945 - 9969.094        
Range of intensity values: 8.707e-05 - 5.85e-03       
Range of snr values      : 2.011 - 135.137            
Memory usage             : 14.43 KiB                  
Name                     : Pankreas_HB_L_061019_A8.A16
File                     : /data/set B - discovery heidelberg/control/Pankreas_HB_L_061019_A8/0_a16/1/1SLin/fid

[[13]]
S4 class type            : MassPeaks                 
Number of m/z values     : 243                       
Range of m/z values      : 2022.992 - 9977.152       
Range of intensity values: 1.061e-04 - 4.145e-03     
Range of snr values      : 2.001 - 78.191            
Memory usage             : 15.32 KiB                 
Name                     : Pankreas_HB_L_061019_C4.F7
File                     : /data/set B - discovery heidelberg/tumor/Pankreas_HB_L_061019_C4/0_f7/1/1SLin/fid

[[14]]
S4 class type            : MassPeaks                 
Number of m/z values     : 220                       
Range of m/z values      : 2012.416 - 9978.442       
Range of intensity values: 1.029e-04 - 4.126e-03     
Range of snr values      : 2.009 - 80.584            
Memory usage             : 14.781 KiB                
Name                     : Pankreas_HB_L_061019_C4.F8
File                     : /data/set B - discovery heidelberg/tumor/Pankreas_HB_L_061019_C4/0_f8/1/1SLin/fid

[[15]]
S4 class type            : MassPeaks                  
Number of m/z values     : 233                        
Range of m/z values      : 2009.524 - 9958.139        
Range of intensity values: 1.05e-04 - 4.285e-03       
Range of snr values      : 2 - 81.63                  
Memory usage             : 15.086 KiB                 
Name                     : Pankreas_HB_L_061019_D9.G17
File                     : /data/set B - discovery heidelberg/tumor/Pankreas_HB_L_061019_D9/0_g17/1/1SLin/fid

[[16]]
S4 class type            : MassPeaks                  
Number of m/z values     : 249                        
Range of m/z values      : 2007.355 - 9956.529        
Range of intensity values: 1.026e-04 - 4.019e-03      
Range of snr values      : 2.004 - 78.473             
Memory usage             : 15.461 KiB                 
Name                     : Pankreas_HB_L_061019_D9.G18
File                     : /data/set B - discovery heidelberg/tumor/Pankreas_HB_L_061019_D9/0_g18/1/1SLin/fid

(def peaks-per-spectrum
  (r->clj (r "sapply(all_peaks, length)")))
peaks-per-spectrum
[252 243 268 265 218 249 231 221 226 232 203 205 243 220 233 249]

Step 3: Peak alignment

Peaks from different spectra may be slightly shifted in m/z due to instrument drift. alignSpectra uses a reference set of peaks to compute a warping function and applies it to the spectra.

(r "spectra_aligned <- alignSpectra(spectra, halfWindowSize=20, SNR=2)")
[[1]]
S4 class type            : MassSpectrum                
Number of m/z values     : 34264                       
Range of m/z values      : 2000.224 - 10000.114        
Range of intensity values: 0e+00 - 7.579e-03           
Memory usage             : 545 KiB                     
Name                     : Pankreas_HB_L_061019_G10.M19
File                     : /data/set A - discovery leipzig/control/Pankreas_HB_L_061019_G10/0_m19/1/1SLin/fid

[[2]]
S4 class type            : MassSpectrum                
Number of m/z values     : 34264                       
Range of m/z values      : 2000.14 - 10001.921         
Range of intensity values: 0e+00 - 7.353e-03           
Memory usage             : 545 KiB                     
Name                     : Pankreas_HB_L_061019_G10.M20
File                     : /data/set A - discovery leipzig/control/Pankreas_HB_L_061019_G10/0_m20/1/1SLin/fid

[[3]]
S4 class type            : MassSpectrum               
Number of m/z values     : 34264                      
Range of m/z values      : 2000.453 - 10001.949       
Range of intensity values: 0e+00 - 6.992e-03          
Memory usage             : 545 KiB                    
Name                     : Pankreas_HB_L_061019_H7.O14
File                     : /data/set A - discovery leipzig/control/Pankreas_HB_L_061019_H7/0_o14/1/1SLin/fid

[[4]]
S4 class type            : MassSpectrum               
Number of m/z values     : 34264                      
Range of m/z values      : 2001.03 - 10005.194        
Range of intensity values: 0e+00 - 6.641e-03          
Memory usage             : 545 KiB                    
Name                     : Pankreas_HB_L_061019_H7.P13
File                     : /data/set A - discovery leipzig/control/Pankreas_HB_L_061019_H7/0_p13/1/1SLin/fid

[[5]]
S4 class type            : MassSpectrum                
Number of m/z values     : 34264                       
Range of m/z values      : 1999.912 - 9998.049         
Range of intensity values: 0e+00 - 6.213e-03           
Memory usage             : 545 KiB                     
Name                     : Pankreas_HB_L_061019_F10.L19
File                     : /data/set A - discovery leipzig/tumor/Pankreas_HB_L_061019_F10/0_l19/1/1SLin/fid

[[6]]
S4 class type            : MassSpectrum                
Number of m/z values     : 34264                       
Range of m/z values      : 2000.215 - 9999.715         
Range of intensity values: 0e+00 - 6.144e-03           
Memory usage             : 545 KiB                     
Name                     : Pankreas_HB_L_061019_F10.L20
File                     : /data/set A - discovery leipzig/tumor/Pankreas_HB_L_061019_F10/0_l20/1/1SLin/fid

[[7]]
S4 class type            : MassSpectrum               
Number of m/z values     : 34264                      
Range of m/z values      : 1999.79 - 9997.317         
Range of intensity values: 0e+00 - 5.975e-03          
Memory usage             : 545 KiB                    
Name                     : Pankreas_HB_L_061019_F9.L17
File                     : /data/set A - discovery leipzig/tumor/Pankreas_HB_L_061019_F9/0_l17/1/1SLin/fid

[[8]]
S4 class type            : MassSpectrum               
Number of m/z values     : 34264                      
Range of m/z values      : 1999.792 - 9998.879        
Range of intensity values: 0e+00 - 6.002e-03          
Memory usage             : 545 KiB                    
Name                     : Pankreas_HB_L_061019_F9.L18
File                     : /data/set A - discovery leipzig/tumor/Pankreas_HB_L_061019_F9/0_l18/1/1SLin/fid

[[9]]
S4 class type            : MassSpectrum               
Number of m/z values     : 34264                      
Range of m/z values      : 2000.616 - 10002.042       
Range of intensity values: 0e+00 - 6.703e-03          
Memory usage             : 545.008 KiB                
Name                     : Pankreas_HB_L_061019_A6.A11
File                     : /data/set B - discovery heidelberg/control/Pankreas_HB_L_061019_A6/0_a11/1/1SLin/fid

[[10]]
S4 class type            : MassSpectrum               
Number of m/z values     : 34264                      
Range of m/z values      : 2000.727 - 10001.581       
Range of intensity values: 0e+00 - 6.895e-03          
Memory usage             : 545.008 KiB                
Name                     : Pankreas_HB_L_061019_A6.A12
File                     : /data/set B - discovery heidelberg/control/Pankreas_HB_L_061019_A6/0_a12/1/1SLin/fid

[[11]]
S4 class type            : MassSpectrum               
Number of m/z values     : 34264                      
Range of m/z values      : 2000.633 - 10001.254       
Range of intensity values: 0e+00 - 5.557e-03          
Memory usage             : 545.008 KiB                
Name                     : Pankreas_HB_L_061019_A8.A15
File                     : /data/set B - discovery heidelberg/control/Pankreas_HB_L_061019_A8/0_a15/1/1SLin/fid

[[12]]
S4 class type            : MassSpectrum               
Number of m/z values     : 34264                      
Range of m/z values      : 2000.6 - 10000.147         
Range of intensity values: 0e+00 - 5.85e-03           
Memory usage             : 545.008 KiB                
Name                     : Pankreas_HB_L_061019_A8.A16
File                     : /data/set B - discovery heidelberg/control/Pankreas_HB_L_061019_A8/0_a16/1/1SLin/fid

[[13]]
S4 class type            : MassSpectrum              
Number of m/z values     : 34264                     
Range of m/z values      : 1998.955 - 9994.75        
Range of intensity values: 0e+00 - 4.145e-03         
Memory usage             : 545.008 KiB               
Name                     : Pankreas_HB_L_061019_C4.F7
File                     : /data/set B - discovery heidelberg/tumor/Pankreas_HB_L_061019_C4/0_f7/1/1SLin/fid

[[14]]
S4 class type            : MassSpectrum              
Number of m/z values     : 34264                     
Range of m/z values      : 1999.234 - 9995.479       
Range of intensity values: 0e+00 - 4.126e-03         
Memory usage             : 545.008 KiB               
Name                     : Pankreas_HB_L_061019_C4.F8
File                     : /data/set B - discovery heidelberg/tumor/Pankreas_HB_L_061019_C4/0_f8/1/1SLin/fid

[[15]]
S4 class type            : MassSpectrum               
Number of m/z values     : 34264                      
Range of m/z values      : 1999.754 - 9996.686        
Range of intensity values: 0e+00 - 4.285e-03          
Memory usage             : 545.008 KiB                
Name                     : Pankreas_HB_L_061019_D9.G17
File                     : /data/set B - discovery heidelberg/tumor/Pankreas_HB_L_061019_D9/0_g17/1/1SLin/fid

[[16]]
S4 class type            : MassSpectrum               
Number of m/z values     : 34264                      
Range of m/z values      : 1999.721 - 9999.414        
Range of intensity values: 0e+00 - 4.019e-03          
Memory usage             : 545.008 KiB                
Name                     : Pankreas_HB_L_061019_D9.G18
File                     : /data/set B - discovery heidelberg/tumor/Pankreas_HB_L_061019_D9/0_g18/1/1SLin/fid

(r "peaks_aligned <- detectPeaks(spectra_aligned, halfWindowSize=20, SNR=2, method='MAD')")
[[1]]
S4 class type            : MassPeaks                   
Number of m/z values     : 252                         
Range of m/z values      : 2012.214 - 9953.366         
Range of intensity values: 9.657e-05 - 7.579e-03       
Range of snr values      : 2.003 - 157.225             
Memory usage             : 15.523 KiB                  
Name                     : Pankreas_HB_L_061019_G10.M19
File                     : /data/set A - discovery leipzig/control/Pankreas_HB_L_061019_G10/0_m19/1/1SLin/fid

[[2]]
S4 class type            : MassPeaks                   
Number of m/z values     : 243                         
Range of m/z values      : 2002.304 - 9862.327         
Range of intensity values: 9.209e-05 - 7.353e-03       
Range of snr values      : 2.012 - 160.65              
Memory usage             : 15.312 KiB                  
Name                     : Pankreas_HB_L_061019_G10.M20
File                     : /data/set A - discovery leipzig/control/Pankreas_HB_L_061019_G10/0_m20/1/1SLin/fid

[[3]]
S4 class type            : MassPeaks                  
Number of m/z values     : 268                        
Range of m/z values      : 2013.601 - 9956.811        
Range of intensity values: 9.873e-05 - 6.992e-03      
Range of snr values      : 2.002 - 141.751            
Memory usage             : 15.898 KiB                 
Name                     : Pankreas_HB_L_061019_H7.O14
File                     : /data/set A - discovery leipzig/control/Pankreas_HB_L_061019_H7/0_o14/1/1SLin/fid

[[4]]
S4 class type            : MassPeaks                  
Number of m/z values     : 265                        
Range of m/z values      : 2018.956 - 9970.686        
Range of intensity values: 9.431e-05 - 6.641e-03      
Range of snr values      : 2.007 - 141.312            
Memory usage             : 15.828 KiB                 
Name                     : Pankreas_HB_L_061019_H7.P13
File                     : /data/set A - discovery leipzig/control/Pankreas_HB_L_061019_H7/0_p13/1/1SLin/fid

[[5]]
S4 class type            : MassPeaks                   
Number of m/z values     : 218                         
Range of m/z values      : 2010.89 - 9715.387          
Range of intensity values: 9.39e-05 - 6.213e-03        
Range of snr values      : 2.007 - 132.815             
Memory usage             : 14.727 KiB                  
Name                     : Pankreas_HB_L_061019_F10.L19
File                     : /data/set A - discovery leipzig/tumor/Pankreas_HB_L_061019_F10/0_l19/1/1SLin/fid

[[6]]
S4 class type            : MassPeaks                   
Number of m/z values     : 249                         
Range of m/z values      : 2002.379 - 9977.133         
Range of intensity values: 9.989e-05 - 6.144e-03       
Range of snr values      : 2 - 123.044                 
Memory usage             : 15.453 KiB                  
Name                     : Pankreas_HB_L_061019_F10.L20
File                     : /data/set A - discovery leipzig/tumor/Pankreas_HB_L_061019_F10/0_l20/1/1SLin/fid

[[7]]
S4 class type            : MassPeaks                  
Number of m/z values     : 231                        
Range of m/z values      : 2023.662 - 9952.179        
Range of intensity values: 9.836e-05 - 5.975e-03      
Range of snr values      : 2.003 - 121.669            
Memory usage             : 15.031 KiB                 
Name                     : Pankreas_HB_L_061019_F9.L17
File                     : /data/set A - discovery leipzig/tumor/Pankreas_HB_L_061019_F9/0_l17/1/1SLin/fid

[[8]]
S4 class type            : MassPeaks                  
Number of m/z values     : 221                        
Range of m/z values      : 2020.327 - 9818.621        
Range of intensity values: 9.469e-05 - 6.002e-03      
Range of snr values      : 2.001 - 126.825            
Memory usage             : 14.797 KiB                 
Name                     : Pankreas_HB_L_061019_F9.L18
File                     : /data/set A - discovery leipzig/tumor/Pankreas_HB_L_061019_F9/0_l18/1/1SLin/fid

[[9]]
S4 class type            : MassPeaks                  
Number of m/z values     : 226                        
Range of m/z values      : 2010.002 - 9930.513        
Range of intensity values: 9.046e-05 - 6.703e-03      
Range of snr values      : 2.003 - 148.4              
Memory usage             : 14.922 KiB                 
Name                     : Pankreas_HB_L_061019_A6.A11
File                     : /data/set B - discovery heidelberg/control/Pankreas_HB_L_061019_A6/0_a11/1/1SLin/fid

[[10]]
S4 class type            : MassPeaks                  
Number of m/z values     : 232                        
Range of m/z values      : 2021.263 - 9800.543        
Range of intensity values: 8.948e-05 - 6.895e-03      
Range of snr values      : 2.001 - 154.204            
Memory usage             : 15.062 KiB                 
Name                     : Pankreas_HB_L_061019_A6.A12
File                     : /data/set B - discovery heidelberg/control/Pankreas_HB_L_061019_A6/0_a12/1/1SLin/fid

[[11]]
S4 class type            : MassPeaks                  
Number of m/z values     : 203                        
Range of m/z values      : 2025.23 - 9989.315         
Range of intensity values: 9.276e-05 - 5.557e-03      
Range of snr values      : 2.025 - 121.319            
Memory usage             : 14.383 KiB                 
Name                     : Pankreas_HB_L_061019_A8.A15
File                     : /data/set B - discovery heidelberg/control/Pankreas_HB_L_061019_A8/0_a15/1/1SLin/fid

[[12]]
S4 class type            : MassPeaks                  
Number of m/z values     : 205                        
Range of m/z values      : 2009.408 - 9969.506        
Range of intensity values: 8.707e-05 - 5.85e-03       
Range of snr values      : 2.011 - 135.137            
Memory usage             : 14.43 KiB                  
Name                     : Pankreas_HB_L_061019_A8.A16
File                     : /data/set B - discovery heidelberg/control/Pankreas_HB_L_061019_A8/0_a16/1/1SLin/fid

[[13]]
S4 class type            : MassPeaks                 
Number of m/z values     : 243                       
Range of m/z values      : 2021.81 - 9972.168        
Range of intensity values: 1.061e-04 - 4.145e-03     
Range of snr values      : 2.001 - 78.191            
Memory usage             : 15.32 KiB                 
Name                     : Pankreas_HB_L_061019_C4.F7
File                     : /data/set B - discovery heidelberg/tumor/Pankreas_HB_L_061019_C4/0_f7/1/1SLin/fid

[[14]]
S4 class type            : MassPeaks                 
Number of m/z values     : 220                       
Range of m/z values      : 2011.513 - 9974.187       
Range of intensity values: 1.029e-04 - 4.126e-03     
Range of snr values      : 2.009 - 80.584            
Memory usage             : 14.781 KiB                
Name                     : Pankreas_HB_L_061019_C4.F8
File                     : /data/set B - discovery heidelberg/tumor/Pankreas_HB_L_061019_C4/0_f8/1/1SLin/fid

[[15]]
S4 class type            : MassPeaks                  
Number of m/z values     : 233                        
Range of m/z values      : 2009.141 - 9955.091        
Range of intensity values: 1.05e-04 - 4.285e-03       
Range of snr values      : 2 - 81.63                  
Memory usage             : 15.086 KiB                 
Name                     : Pankreas_HB_L_061019_D9.G17
File                     : /data/set B - discovery heidelberg/tumor/Pankreas_HB_L_061019_D9/0_g17/1/1SLin/fid

[[16]]
S4 class type            : MassPeaks                  
Number of m/z values     : 249                        
Range of m/z values      : 2006.94 - 9956.209         
Range of intensity values: 1.026e-04 - 4.019e-03      
Range of snr values      : 2.004 - 78.473             
Memory usage             : 15.461 KiB                 
Name                     : Pankreas_HB_L_061019_D9.G18
File                     : /data/set B - discovery heidelberg/tumor/Pankreas_HB_L_061019_D9/0_g18/1/1SLin/fid

Step 4: Peak binning

Even after alignment, exact m/z values differ slightly between spectra. binPeaks clusters nearby m/z values into shared bins.

(r "peaks_binned <- binPeaks(peaks_aligned)")
[[1]]
S4 class type            : MassPeaks                   
Number of m/z values     : 252                         
Range of m/z values      : 2010.464 - 9954.731         
Range of intensity values: 9.657e-05 - 7.579e-03       
Range of snr values      : 2.003 - 157.225             
Memory usage             : 15.523 KiB                  
Name                     : Pankreas_HB_L_061019_G10.M19
File                     : /data/set A - discovery leipzig/control/Pankreas_HB_L_061019_G10/0_m19/1/1SLin/fid

[[2]]
S4 class type            : MassPeaks                   
Number of m/z values     : 243                         
Range of m/z values      : 2002.342 - 9868.06          
Range of intensity values: 9.209e-05 - 7.353e-03       
Range of snr values      : 2.012 - 160.65              
Memory usage             : 15.312 KiB                  
Name                     : Pankreas_HB_L_061019_G10.M20
File                     : /data/set A - discovery leipzig/control/Pankreas_HB_L_061019_G10/0_m20/1/1SLin/fid

[[3]]
S4 class type            : MassPeaks                  
Number of m/z values     : 268                        
Range of m/z values      : 2010.464 - 9954.731        
Range of intensity values: 9.873e-05 - 6.992e-03      
Range of snr values      : 2.002 - 141.751            
Memory usage             : 15.898 KiB                 
Name                     : Pankreas_HB_L_061019_H7.O14
File                     : /data/set A - discovery leipzig/control/Pankreas_HB_L_061019_H7/0_o14/1/1SLin/fid

[[4]]
S4 class type            : MassPeaks                  
Number of m/z values     : 265                        
Range of m/z values      : 2019.388 - 9973.725        
Range of intensity values: 9.431e-05 - 6.641e-03      
Range of snr values      : 2.007 - 141.312            
Memory usage             : 15.828 KiB                 
Name                     : Pankreas_HB_L_061019_H7.P13
File                     : /data/set A - discovery leipzig/control/Pankreas_HB_L_061019_H7/0_p13/1/1SLin/fid

[[5]]
S4 class type            : MassPeaks                   
Number of m/z values     : 218                         
Range of m/z values      : 2010.464 - 9715.852         
Range of intensity values: 9.39e-05 - 6.213e-03        
Range of snr values      : 2.007 - 132.815             
Memory usage             : 14.727 KiB                  
Name                     : Pankreas_HB_L_061019_F10.L19
File                     : /data/set A - discovery leipzig/tumor/Pankreas_HB_L_061019_F10/0_l19/1/1SLin/fid

[[6]]
S4 class type            : MassPeaks                   
Number of m/z values     : 249                         
Range of m/z values      : 2002.342 - 9973.725         
Range of intensity values: 9.989e-05 - 6.144e-03       
Range of snr values      : 2 - 123.044                 
Memory usage             : 15.453 KiB                  
Name                     : Pankreas_HB_L_061019_F10.L20
File                     : /data/set A - discovery leipzig/tumor/Pankreas_HB_L_061019_F10/0_l20/1/1SLin/fid

[[7]]
S4 class type            : MassPeaks                  
Number of m/z values     : 231                        
Range of m/z values      : 2023.298 - 9954.731        
Range of intensity values: 9.836e-05 - 5.975e-03      
Range of snr values      : 2.003 - 121.669            
Memory usage             : 15.031 KiB                 
Name                     : Pankreas_HB_L_061019_F9.L17
File                     : /data/set A - discovery leipzig/tumor/Pankreas_HB_L_061019_F9/0_l17/1/1SLin/fid

[[8]]
S4 class type            : MassPeaks                  
Number of m/z values     : 221                        
Range of m/z values      : 2019.388 - 9816.199        
Range of intensity values: 9.469e-05 - 6.002e-03      
Range of snr values      : 2.001 - 126.825            
Memory usage             : 14.797 KiB                 
Name                     : Pankreas_HB_L_061019_F9.L18
File                     : /data/set A - discovery leipzig/tumor/Pankreas_HB_L_061019_F9/0_l18/1/1SLin/fid

[[9]]
S4 class type            : MassPeaks                  
Number of m/z values     : 226                        
Range of m/z values      : 2010.464 - 9925.968        
Range of intensity values: 9.046e-05 - 6.703e-03      
Range of snr values      : 2.003 - 148.4              
Memory usage             : 14.922 KiB                 
Name                     : Pankreas_HB_L_061019_A6.A11
File                     : /data/set B - discovery heidelberg/control/Pankreas_HB_L_061019_A6/0_a11/1/1SLin/fid

[[10]]
S4 class type            : MassPeaks                  
Number of m/z values     : 232                        
Range of m/z values      : 2023.298 - 9796.211        
Range of intensity values: 8.948e-05 - 6.895e-03      
Range of snr values      : 2.001 - 154.204            
Memory usage             : 15.062 KiB                 
Name                     : Pankreas_HB_L_061019_A6.A12
File                     : /data/set B - discovery heidelberg/control/Pankreas_HB_L_061019_A6/0_a12/1/1SLin/fid

[[11]]
S4 class type            : MassPeaks                  
Number of m/z values     : 203                        
Range of m/z values      : 2023.298 - 9989.315        
Range of intensity values: 9.276e-05 - 5.557e-03      
Range of snr values      : 2.025 - 121.319            
Memory usage             : 14.383 KiB                 
Name                     : Pankreas_HB_L_061019_A8.A15
File                     : /data/set B - discovery heidelberg/control/Pankreas_HB_L_061019_A8/0_a15/1/1SLin/fid

[[12]]
S4 class type            : MassPeaks                  
Number of m/z values     : 205                        
Range of m/z values      : 2010.464 - 9973.725        
Range of intensity values: 8.707e-05 - 5.85e-03       
Range of snr values      : 2.011 - 135.137            
Memory usage             : 14.43 KiB                  
Name                     : Pankreas_HB_L_061019_A8.A16
File                     : /data/set B - discovery heidelberg/control/Pankreas_HB_L_061019_A8/0_a16/1/1SLin/fid

[[13]]
S4 class type            : MassPeaks                 
Number of m/z values     : 243                       
Range of m/z values      : 2023.298 - 9973.725       
Range of intensity values: 1.061e-04 - 4.145e-03     
Range of snr values      : 2.001 - 78.191            
Memory usage             : 15.32 KiB                 
Name                     : Pankreas_HB_L_061019_C4.F7
File                     : /data/set B - discovery heidelberg/tumor/Pankreas_HB_L_061019_C4/0_f7/1/1SLin/fid

[[14]]
S4 class type            : MassPeaks                 
Number of m/z values     : 220                       
Range of m/z values      : 2010.464 - 9973.725       
Range of intensity values: 1.029e-04 - 4.126e-03     
Range of snr values      : 2.009 - 80.584            
Memory usage             : 14.781 KiB                
Name                     : Pankreas_HB_L_061019_C4.F8
File                     : /data/set B - discovery heidelberg/tumor/Pankreas_HB_L_061019_C4/0_f8/1/1SLin/fid

[[15]]
S4 class type            : MassPeaks                  
Number of m/z values     : 233                        
Range of m/z values      : 2010.464 - 9954.731        
Range of intensity values: 1.05e-04 - 4.285e-03       
Range of snr values      : 2 - 81.63                  
Memory usage             : 15.086 KiB                 
Name                     : Pankreas_HB_L_061019_D9.G17
File                     : /data/set B - discovery heidelberg/tumor/Pankreas_HB_L_061019_D9/0_g17/1/1SLin/fid

[[16]]
S4 class type            : MassPeaks                  
Number of m/z values     : 249                        
Range of m/z values      : 2010.464 - 9954.731        
Range of intensity values: 1.026e-04 - 4.019e-03      
Range of snr values      : 2.004 - 78.473             
Memory usage             : 15.461 KiB                 
Name                     : Pankreas_HB_L_061019_D9.G18
File                     : /data/set B - discovery heidelberg/tumor/Pankreas_HB_L_061019_D9/0_g18/1/1SLin/fid

Step 5: Peak filtering

filterPeaks keeps only peaks that appear in at least a minimum fraction of spectra. With minFrequency=0.25, a peak must appear in at least 4 of 16 spectra.

(r "peaks_filtered <- filterPeaks(peaks_binned, minFrequency=0.25)")
[[1]]
S4 class type            : MassPeaks                   
Number of m/z values     : 225                         
Range of m/z values      : 2010.464 - 9954.731         
Range of intensity values: 9.711e-05 - 7.579e-03       
Range of snr values      : 2.014 - 157.225             
Memory usage             : 14.891 KiB                  
Name                     : Pankreas_HB_L_061019_G10.M19
File                     : /data/set A - discovery leipzig/control/Pankreas_HB_L_061019_G10/0_m19/1/1SLin/fid

[[2]]
S4 class type            : MassPeaks                   
Number of m/z values     : 218                         
Range of m/z values      : 2023.298 - 9868.06          
Range of intensity values: 9.209e-05 - 7.353e-03       
Range of snr values      : 2.012 - 160.65              
Memory usage             : 14.727 KiB                  
Name                     : Pankreas_HB_L_061019_G10.M20
File                     : /data/set A - discovery leipzig/control/Pankreas_HB_L_061019_G10/0_m20/1/1SLin/fid

[[3]]
S4 class type            : MassPeaks                  
Number of m/z values     : 233                        
Range of m/z values      : 2010.464 - 9954.731        
Range of intensity values: 9.873e-05 - 6.992e-03      
Range of snr values      : 2.002 - 141.751            
Memory usage             : 15.078 KiB                 
Name                     : Pankreas_HB_L_061019_H7.O14
File                     : /data/set A - discovery leipzig/control/Pankreas_HB_L_061019_H7/0_o14/1/1SLin/fid

[[4]]
S4 class type            : MassPeaks                  
Number of m/z values     : 230                        
Range of m/z values      : 2019.388 - 9973.725        
Range of intensity values: 9.431e-05 - 6.641e-03      
Range of snr values      : 2.007 - 141.312            
Memory usage             : 15.008 KiB                 
Name                     : Pankreas_HB_L_061019_H7.P13
File                     : /data/set A - discovery leipzig/control/Pankreas_HB_L_061019_H7/0_p13/1/1SLin/fid

[[5]]
S4 class type            : MassPeaks                   
Number of m/z values     : 205                         
Range of m/z values      : 2010.464 - 9715.852         
Range of intensity values: 9.39e-05 - 6.213e-03        
Range of snr values      : 2.007 - 132.815             
Memory usage             : 14.422 KiB                  
Name                     : Pankreas_HB_L_061019_F10.L19
File                     : /data/set A - discovery leipzig/tumor/Pankreas_HB_L_061019_F10/0_l19/1/1SLin/fid

[[6]]
S4 class type            : MassPeaks                   
Number of m/z values     : 210                         
Range of m/z values      : 2023.298 - 9973.725         
Range of intensity values: 9.989e-05 - 6.144e-03       
Range of snr values      : 2 - 123.044                 
Memory usage             : 14.539 KiB                  
Name                     : Pankreas_HB_L_061019_F10.L20
File                     : /data/set A - discovery leipzig/tumor/Pankreas_HB_L_061019_F10/0_l20/1/1SLin/fid

[[7]]
S4 class type            : MassPeaks                  
Number of m/z values     : 213                        
Range of m/z values      : 2023.298 - 9954.731        
Range of intensity values: 9.836e-05 - 5.975e-03      
Range of snr values      : 2.003 - 121.669            
Memory usage             : 14.609 KiB                 
Name                     : Pankreas_HB_L_061019_F9.L17
File                     : /data/set A - discovery leipzig/tumor/Pankreas_HB_L_061019_F9/0_l17/1/1SLin/fid

[[8]]
S4 class type            : MassPeaks                  
Number of m/z values     : 205                        
Range of m/z values      : 2019.388 - 9816.199        
Range of intensity values: 9.469e-05 - 6.002e-03      
Range of snr values      : 2.001 - 126.825            
Memory usage             : 14.422 KiB                 
Name                     : Pankreas_HB_L_061019_F9.L18
File                     : /data/set A - discovery leipzig/tumor/Pankreas_HB_L_061019_F9/0_l18/1/1SLin/fid

[[9]]
S4 class type            : MassPeaks                  
Number of m/z values     : 208                        
Range of m/z values      : 2010.464 - 9925.968        
Range of intensity values: 9.046e-05 - 6.703e-03      
Range of snr values      : 2.003 - 148.4              
Memory usage             : 14.5 KiB                   
Name                     : Pankreas_HB_L_061019_A6.A11
File                     : /data/set B - discovery heidelberg/control/Pankreas_HB_L_061019_A6/0_a11/1/1SLin/fid

[[10]]
S4 class type            : MassPeaks                  
Number of m/z values     : 214                        
Range of m/z values      : 2023.298 - 9758.985        
Range of intensity values: 8.948e-05 - 6.895e-03      
Range of snr values      : 2.001 - 154.204            
Memory usage             : 14.641 KiB                 
Name                     : Pankreas_HB_L_061019_A6.A12
File                     : /data/set B - discovery heidelberg/control/Pankreas_HB_L_061019_A6/0_a12/1/1SLin/fid

[[11]]
S4 class type            : MassPeaks                  
Number of m/z values     : 194                        
Range of m/z values      : 2023.298 - 9973.725        
Range of intensity values: 9.276e-05 - 5.557e-03      
Range of snr values      : 2.025 - 121.319            
Memory usage             : 14.172 KiB                 
Name                     : Pankreas_HB_L_061019_A8.A15
File                     : /data/set B - discovery heidelberg/control/Pankreas_HB_L_061019_A8/0_a15/1/1SLin/fid

[[12]]
S4 class type            : MassPeaks                  
Number of m/z values     : 194                        
Range of m/z values      : 2010.464 - 9973.725        
Range of intensity values: 8.707e-05 - 5.85e-03       
Range of snr values      : 2.011 - 135.137            
Memory usage             : 14.172 KiB                 
Name                     : Pankreas_HB_L_061019_A8.A16
File                     : /data/set B - discovery heidelberg/control/Pankreas_HB_L_061019_A8/0_a16/1/1SLin/fid

[[13]]
S4 class type            : MassPeaks                 
Number of m/z values     : 214                       
Range of m/z values      : 2023.298 - 9973.725       
Range of intensity values: 1.061e-04 - 4.145e-03     
Range of snr values      : 2.001 - 78.191            
Memory usage             : 14.641 KiB                
Name                     : Pankreas_HB_L_061019_C4.F7
File                     : /data/set B - discovery heidelberg/tumor/Pankreas_HB_L_061019_C4/0_f7/1/1SLin/fid

[[14]]
S4 class type            : MassPeaks                 
Number of m/z values     : 195                       
Range of m/z values      : 2010.464 - 9973.725       
Range of intensity values: 1.03e-04 - 4.126e-03      
Range of snr values      : 2.011 - 80.584            
Memory usage             : 14.195 KiB                
Name                     : Pankreas_HB_L_061019_C4.F8
File                     : /data/set B - discovery heidelberg/tumor/Pankreas_HB_L_061019_C4/0_f8/1/1SLin/fid

[[15]]
S4 class type            : MassPeaks                  
Number of m/z values     : 207                        
Range of m/z values      : 2010.464 - 9954.731        
Range of intensity values: 1.052e-04 - 4.285e-03      
Range of snr values      : 2.004 - 81.63              
Memory usage             : 14.477 KiB                 
Name                     : Pankreas_HB_L_061019_D9.G17
File                     : /data/set B - discovery heidelberg/tumor/Pankreas_HB_L_061019_D9/0_g17/1/1SLin/fid

[[16]]
S4 class type            : MassPeaks                  
Number of m/z values     : 213                        
Range of m/z values      : 2010.464 - 9954.731        
Range of intensity values: 1.026e-04 - 4.019e-03      
Range of snr values      : 2.004 - 78.473             
Memory usage             : 14.617 KiB                 
Name                     : Pankreas_HB_L_061019_D9.G18
File                     : /data/set B - discovery heidelberg/tumor/Pankreas_HB_L_061019_D9/0_g18/1/1SLin/fid

(def n-filtered-peaks
  (r->clj (r "sapply(peaks_filtered, length)")))
n-filtered-peaks
[225 218 233 230 205 210 213 205 208 214 194 194 214 195 207 213]

Step 6: Feature matrix

intensityMatrix converts the aligned, binned, filtered peaks into a samples-by-features matrix suitable for machine learning.

(r "fm <- intensityMatrix(peaks_filtered, spectra_aligned)")
      2010.46372181388 2019.38802650132 2023.2977876444 2038.72956667587
 [1,]     1.210506e-04     1.993417e-04    1.070342e-04     7.300137e-05
 [2,]     4.400536e-05     3.790488e-05    1.025702e-04     1.047463e-04
 [3,]     1.155278e-04     1.448812e-04    2.626509e-05     9.281477e-05
 [4,]     3.132117e-05     1.544475e-04    6.939479e-05     7.105812e-05
 [5,]     1.044336e-04     1.630882e-05    2.115012e-04     2.371837e-05
 [6,]     2.504093e-05     6.697287e-05    1.941454e-04     5.363917e-05
 [7,]     7.886063e-05     6.570303e-05    1.825488e-04     9.492398e-05
 [8,]     3.341655e-05     1.394128e-04    1.247749e-04     4.045832e-05
 [9,]     1.074566e-04     6.510684e-05    1.265661e-04     5.672352e-05
[10,]     6.208333e-05     9.906698e-05    1.228416e-04     7.211646e-05
[11,]     4.089187e-05     9.318478e-05    2.124132e-04     1.422758e-04
[12,]     9.587934e-05     8.995848e-05    1.950347e-04     1.056014e-04
[13,]     9.775870e-06     3.692000e-04    8.278236e-04     2.269902e-04
[14,]     1.202038e-04     3.645247e-04    8.898072e-04     1.527356e-04
[15,]     1.067410e-04     1.903003e-04    7.910257e-04     1.263589e-04
[16,]     2.563730e-04     2.718498e-04    7.406616e-04     1.256374e-04
      2045.65282117107 2067.61159972805 2072.04781721356 2081.62304745686
 [1,]     0.0001342010     0.0001425297     8.445436e-05     0.0005918102
 [2,]     0.0001402854     0.0002319625     1.721385e-04     0.0006515445
 [3,]     0.0001099131     0.0001758044     1.085525e-04     0.0008590112
 [4,]     0.0001621597     0.0001824497     1.082504e-04     0.0009218033
 [5,]     0.0002368342     0.0002051817     1.821955e-04     0.0006923273
 [6,]     0.0001907535     0.0001579290     2.161255e-04     0.0007315509
 [7,]     0.0001876065     0.0001862506     8.549708e-05     0.0004805097
 [8,]     0.0001601275     0.0002114925     9.081995e-05     0.0004360625
 [9,]     0.0002233786     0.0001618136     1.179099e-04     0.0003716667
[10,]     0.0002480737     0.0001736907     1.562085e-04     0.0003607110
[11,]     0.0002183129     0.0001794301     1.368684e-04     0.0008051578
[12,]     0.0002086279     0.0001439165     1.509857e-04     0.0007614646
[13,]     0.0001823328     0.0002658753     9.440444e-05     0.0014384219
[14,]     0.0002556911     0.0002067906     3.187510e-05     0.0016660713
[15,]     0.0001630247     0.0002604932     4.044431e-05     0.0023278542
[16,]     0.0002095761     0.0002488383     7.037197e-05     0.0019315328
      2092.4017847307 2105.13340864733 2113.53888094399 2123.28263737811
 [1,]    1.185319e-04     0.0004186724     1.628411e-05     9.859979e-05
 [2,]    2.042644e-04     0.0004713007     1.086596e-04     6.004790e-05
 [3,]    8.327090e-05     0.0005463139     6.833207e-05     6.850584e-05
 [4,]    7.617156e-05     0.0005623556     1.275649e-04     2.613654e-05
 [5,]    3.738036e-04     0.0007856701     1.734118e-04     3.550277e-05
 [6,]    3.497580e-04     0.0007212871     2.574773e-04     6.003622e-05
 [7,]    2.668241e-04     0.0003932749     1.183426e-04     1.953994e-05
 [8,]    1.604252e-04     0.0005785555     1.776502e-04     4.870442e-05
 [9,]    1.794671e-04     0.0007017110     1.006928e-04     7.668267e-05
[10,]    1.024838e-04     0.0007434170     9.958325e-05     6.163105e-05
[11,]    1.862946e-04     0.0008241476     1.617020e-04     4.131826e-05
[12,]    3.194351e-04     0.0007634470     1.438166e-04     8.794021e-05
[13,]    4.503198e-04     0.0005155775     2.353579e-05     1.549241e-04
[14,]    6.458137e-04     0.0006397791     3.365175e-05     7.578327e-05
[15,]    6.579865e-04     0.0004945284     1.904985e-04     1.854927e-04
[16,]    3.212661e-04     0.0004231703     1.533510e-04     7.832647e-05
      2134.16642829054 2141.35288165096 2153.28194753832 2159.37044610478
 [1,]     4.726874e-05     1.839750e-04     1.124081e-04     1.379477e-04
 [2,]     2.361132e-05     1.031389e-04     1.400458e-04     1.138354e-04
 [3,]     5.414960e-05     7.138999e-05     1.215775e-04     1.209412e-04
 [4,]     7.707720e-05     5.468827e-05     1.373768e-04     1.389748e-04
 [5,]     1.642870e-04     4.204271e-05     9.377908e-05     4.054494e-04
 [6,]     1.274355e-04     7.631858e-05     3.387959e-05     4.567983e-04
 [7,]     1.319369e-04     9.896772e-05     1.453400e-04     1.525286e-04
 [8,]     1.213553e-04     1.051518e-04     1.644955e-04     1.222403e-04
 [9,]     2.078657e-04     1.432580e-04     1.987542e-04     1.610447e-04
[10,]     1.302583e-04     4.216021e-05     1.283282e-04     9.673127e-05
[11,]     1.192442e-04     7.995527e-05     6.898373e-04     1.667249e-04
[12,]     1.209645e-04     4.580054e-05     6.300736e-04     2.445718e-04
[13,]     1.796098e-04     1.483155e-04     1.355487e-04     7.561801e-05
[14,]     1.531785e-04     1.917656e-04     2.318144e-04     1.730311e-04
[15,]     1.725470e-04     3.863328e-05     5.249275e-04     6.606523e-05
[16,]     1.799566e-04     1.118773e-04     3.680441e-04     4.797371e-05
      2163.25287420978 2178.82247406275 2190.67003646009 2194.26201501391
 [1,]     1.502802e-04     1.863214e-05     2.594915e-04     1.911162e-04
 [2,]     1.689144e-04     1.223341e-04     2.195439e-04     1.644015e-04
 [3,]     1.301144e-04     1.720309e-05     7.482642e-05     5.419167e-05
 [4,]     1.574197e-04     9.869539e-05     1.636902e-04     1.001798e-04
 [5,]     3.431557e-04     1.400267e-04     2.046047e-04     1.594414e-04
 [6,]     3.449122e-04     1.264888e-04     3.206671e-04     2.012030e-04
 [7,]     1.821983e-04     1.180233e-04     1.169750e-04     1.497595e-04
 [8,]     1.692560e-04     8.661295e-05     1.373337e-04     1.266002e-04
 [9,]     1.291682e-04     1.422056e-04     7.360140e-05     1.968361e-05
[10,]     1.014572e-04     1.290877e-04     9.730383e-05     8.839511e-06
[11,]     1.592804e-04     1.408048e-04     1.784765e-04     6.904852e-05
[12,]     2.506036e-04     1.767715e-04     2.185821e-04     1.436901e-04
[13,]     1.038445e-04     3.549528e-04     1.773862e-04     2.008118e-04
[14,]     4.514342e-05     3.685949e-04     1.277698e-04     1.952820e-04
[15,]     1.621268e-04     1.189537e-04     6.529312e-05     1.010653e-04
[16,]     6.842879e-05     1.080535e-04     1.013197e-04     3.495206e-05
      2209.90566739532 2230.59383637624 2245.13341811979 2251.27533270642
 [1,]     0.0005026185     1.133986e-04     1.172736e-04     1.371885e-04
 [2,]     0.0005730879     7.728257e-06     1.139034e-04     1.420245e-04
 [3,]     0.0005817462     3.257954e-05     1.078990e-04     1.275623e-04
 [4,]     0.0007156735     9.932051e-05     1.138395e-04     2.433278e-05
 [5,]     0.0003700478     4.589567e-06     9.613422e-05     1.483900e-04
 [6,]     0.0003360071     4.537446e-05     1.772963e-04     1.227221e-04
 [7,]     0.0001671910     1.502636e-04     1.292791e-05     8.084299e-05
 [8,]     0.0001847546     1.213602e-04     1.545467e-05     7.416765e-05
 [9,]     0.0001886906     2.951223e-05     9.450751e-05     1.610857e-04
[10,]     0.0001576984     1.022401e-04     9.150112e-05     1.472656e-04
[11,]     0.0003855563     6.512682e-05     1.220892e-04     1.246931e-04
[12,]     0.0004315023     1.076251e-04     1.556182e-04     1.553644e-04
[13,]     0.0008125446     2.838138e-04     2.277801e-04     2.280742e-04
[14,]     0.0010096752     4.136386e-04     2.810983e-04     2.153144e-04
[15,]     0.0013348338     2.183419e-04     4.598612e-05     2.766289e-04
[16,]     0.0010992321     2.665814e-04     5.787841e-06     2.082191e-04
      2264.90793102928 2279.20283071085 2294.35901377343 2306.622848006
 [1,]     1.829270e-04     2.457157e-04     0.0004374749   1.243566e-04
 [2,]     2.092193e-04     2.283668e-04     0.0004924751   6.211185e-05
 [3,]     1.454510e-04     1.942717e-04     0.0002252753   1.214029e-04
 [4,]     1.417238e-04     2.542556e-04     0.0003357565   1.135370e-04
 [5,]     1.505699e-04     3.238407e-04     0.0002605278   1.097341e-04
 [6,]     1.779764e-04     2.760804e-04     0.0003215382   6.530919e-05
 [7,]     4.823964e-05     2.601626e-04     0.0004651447   6.546152e-05
 [8,]     9.540997e-05     2.327038e-04     0.0003345341   4.937287e-05
 [9,]     5.419122e-05     2.315878e-04     0.0003199438   6.215578e-05
[10,]     5.824273e-05     3.594388e-04     0.0002817886   1.003158e-04
[11,]     5.602852e-05     3.368095e-04     0.0002992864   1.031659e-04
[12,]     8.904675e-05     3.009937e-04     0.0003458662   4.789239e-05
[13,]     1.956681e-04     1.934717e-04     0.0002218202   1.516982e-04
[14,]     1.177938e-04     2.502679e-04     0.0002060572   1.516835e-04
[15,]     2.225413e-04     8.721711e-05     0.0003775174   2.190758e-05
[16,]     1.375694e-04     5.988392e-05     0.0003663471   2.254469e-05
      2318.54348347902 2325.93176963302 2337.12821512602 2347.7175129187
 [1,]     1.238909e-04     2.138919e-04     1.166121e-04    0.0001706054
 [2,]     6.808896e-05     2.742421e-04     1.019963e-04    0.0001599298
 [3,]     1.228191e-04     1.557572e-04     7.142836e-05    0.0001861974
 [4,]     4.608879e-05     2.193549e-04     1.745118e-04    0.0001796775
 [5,]     1.080953e-04     2.619737e-04     1.293430e-04    0.0005923731
 [6,]     1.939910e-05     2.565609e-04     9.988878e-05    0.0006462369
 [7,]     3.517646e-05     2.028612e-04     1.293507e-04    0.0003095121
 [8,]     9.925096e-05     1.267505e-04     1.366139e-04    0.0001898630
 [9,]     7.105093e-05     1.356901e-04     1.080640e-04    0.0002098793
[10,]     4.811016e-05     1.761726e-04     1.463797e-04    0.0001935491
[11,]     1.187857e-04     2.133390e-04     8.071830e-05    0.0003500119
[12,]     5.119687e-05     1.483440e-04     9.062143e-05    0.0004387959
[13,]     9.525553e-05     1.727064e-04     1.622828e-04    0.0002118830
[14,]     6.049727e-05     1.803115e-04     1.441695e-04    0.0001804813
[15,]     3.650622e-05     7.288604e-05     1.372503e-04    0.0002041231
[16,]     8.021492e-05     1.493321e-04     1.105451e-04    0.0002046691
      2354.24766595096 2368.9430574178 2379.27591342701 2393.05103731591
 [1,]     1.192660e-04    1.091742e-04     0.0001888534     1.130613e-04
 [2,]     7.583024e-05    1.205953e-04     0.0001475808     5.027920e-05
 [3,]     1.512519e-04    1.188988e-04     0.0002159035     1.071039e-04
 [4,]     1.062182e-04    3.567815e-05     0.0002092932     4.308708e-05
 [5,]     3.529362e-04    1.849871e-04     0.0001687255     2.065127e-05
 [6,]     2.378844e-04    2.097522e-04     0.0001975256     6.666997e-05
 [7,]     1.856675e-04    1.846877e-04     0.0001474264     3.221204e-05
 [8,]     1.597400e-04    1.020625e-04     0.0001808097     6.353401e-05
 [9,]     8.770636e-05    4.546645e-05     0.0003056160     4.454428e-05
[10,]     6.079276e-05    3.021889e-05     0.0003084305     3.154424e-05
[11,]     1.688449e-04    1.817090e-05     0.0003045908     5.766195e-05
[12,]     1.004627e-04    1.942487e-05     0.0002807963     1.441520e-05
[13,]     1.115136e-04    1.259230e-04     0.0002768048     6.396481e-05
[14,]     6.962496e-05    1.452142e-05     0.0002534977     1.055760e-04
[15,]     2.427989e-04    1.152666e-04     0.0003018159     1.365914e-04
[16,]     2.096202e-04    3.320767e-05     0.0002961587     1.490199e-04
      2403.46795853612 2409.22571493166 2423.58654051572 2428.26413301708
 [1,]     1.577363e-05     5.954957e-05     0.0003329138     0.0001162420
 [2,]     5.268437e-05     1.110753e-04     0.0003983258     0.0002016261
 [3,]     5.258884e-05     1.094004e-04     0.0002584966     0.0002135025
 [4,]     8.332537e-05     1.405406e-04     0.0002801275     0.0001235750
 [5,]     3.236710e-05     6.058568e-05     0.0008880456     0.0002970524
 [6,]     1.219210e-04     1.181719e-04     0.0007940647     0.0002802373
 [7,]     1.034128e-04     7.019435e-05     0.0005659396     0.0002735074
 [8,]     1.035426e-04     1.042858e-04     0.0005382558     0.0002139753
 [9,]     1.347162e-04     1.032616e-04     0.0005838508     0.0002144275
[10,]     9.014908e-05     9.271911e-05     0.0004796576     0.0002228970
[11,]     1.044312e-04     3.727586e-05     0.0010702982     0.0003735978
[12,]     8.807097e-05     1.052361e-04     0.0011106388     0.0003948748
[13,]     6.639507e-05     1.905727e-04     0.0002338446     0.0005767909
[14,]     6.870946e-05     1.185889e-04     0.0003232034     0.0005955444
[15,]     1.391586e-05     9.403466e-05     0.0003855334     0.0002736294
[16,]     7.051330e-05     1.727922e-04     0.0002954989     0.0003036939
      2449.52823053379 2467.0020493968 2481.94695020369 2495.00406132144
 [1,]     1.741170e-04    0.0002245977     1.231931e-04     0.0003557956
 [2,]     1.478056e-04    0.0002737704     6.956824e-05     0.0004749649
 [3,]     1.376724e-04    0.0003369703     9.498596e-06     0.0003419785
 [4,]     1.516775e-04    0.0003534002     7.201056e-05     0.0003743264
 [5,]     1.596141e-04    0.0001532583     1.451102e-04     0.0014182648
 [6,]     1.251367e-04    0.0001684937     5.022206e-05     0.0014322568
 [7,]     8.287858e-05    0.0003500193     1.583734e-04     0.0005858222
 [8,]     1.937461e-04    0.0003210217     1.065091e-04     0.0005766024
 [9,]     1.234161e-04    0.0004939895     3.610891e-05     0.0006040369
[10,]     1.705453e-04    0.0004344638     3.539286e-05     0.0005652330
[11,]     9.311752e-05    0.0003451684     4.440599e-05     0.0011138608
[12,]     1.677007e-04    0.0003131997     7.342917e-05     0.0012423658
[13,]     8.239925e-05    0.0002619775     3.452373e-05     0.0003628592
[14,]     3.494458e-05    0.0002890333     1.475881e-06     0.0003364380
[15,]     1.102017e-04    0.0007833182     1.108998e-04     0.0001663123
[16,]     6.243737e-05    0.0006928464     1.199072e-04     0.0001793832
      2509.80866216721 2522.31009107036 2535.61788550317 2545.10431985288
 [1,]     0.0002341934     1.262669e-04     1.074606e-04     1.467599e-04
 [2,]     0.0002272636     9.983892e-05     1.038475e-04     1.317387e-04
 [3,]     0.0001753728     1.282657e-04     1.209328e-04     8.850545e-05
 [4,]     0.0002231120     1.456878e-04     2.685878e-06     1.432198e-04
 [5,]     0.0004603191     3.626883e-05     5.409371e-05     1.466076e-04
 [6,]     0.0003343744     1.081736e-04     6.740275e-05     1.157324e-04
 [7,]     0.0002777326     1.214569e-04     1.077940e-04     1.584884e-04
 [8,]     0.0002816667     1.120383e-04     1.830927e-05     1.259729e-04
 [9,]     0.0002314807     1.005364e-04     9.663434e-05     2.016582e-04
[10,]     0.0001741831     5.575098e-05     1.005965e-04     2.213752e-04
[11,]     0.0002396932     9.402163e-05     7.512046e-05     2.416285e-04
[12,]     0.0002355167     3.874879e-05     3.260740e-05     3.226106e-04
[13,]     0.0001374061     1.287237e-04     1.300188e-04     6.491970e-05
[14,]     0.0001215154     1.145191e-04     2.595531e-05     7.403464e-05
[15,]     0.0002538274     8.293095e-05     6.538947e-05     1.297435e-04
[16,]     0.0001608631     1.217651e-04     3.667331e-05     1.118590e-04
      2553.94520319165 2562.81146954903 2581.17152095104 2588.17914452778
 [1,]     0.0016921621     0.0001805253     1.314602e-04     1.509396e-04
 [2,]     0.0017005445     0.0002199767     3.633624e-05     7.585383e-05
 [3,]     0.0015372080     0.0002135825     5.841212e-05     1.368010e-04
 [4,]     0.0016134724     0.0002075952     1.031470e-04     2.293849e-04
 [5,]     0.0007664925     0.0006464653     7.721300e-05     1.379954e-04
 [6,]     0.0007095102     0.0005983284     3.392696e-05     5.350343e-05
 [7,]     0.0011134868     0.0006775430     1.054959e-04     7.715240e-05
 [8,]     0.0011610707     0.0004946099     1.244043e-04     1.070357e-04
 [9,]     0.0019012602     0.0002488017     5.861198e-05     6.165463e-05
[10,]     0.0017702499     0.0002564462     3.933615e-05     1.161954e-04
[11,]     0.0018411659     0.0004106344     6.799838e-05     8.354504e-05
[12,]     0.0019052157     0.0004545514     6.327832e-05     2.158814e-05
[13,]     0.0006958463     0.0004419241     1.691146e-04     1.018141e-04
[14,]     0.0007826321     0.0003345508     1.769179e-04     2.134458e-04
[15,]     0.0022623817     0.0001861984     2.351255e-04     1.319842e-04
[16,]     0.0020081991     0.0001106463     1.936966e-04     3.777592e-05
      2604.23549575005 2610.98333444184 2616.83787567359 2623.80629073153
 [1,]     0.0003053761     6.919066e-05     6.223112e-05     1.640415e-04
 [2,]     0.0002909495     9.160758e-05     9.854179e-05     1.521220e-04
 [3,]     0.0002437793     7.189441e-05     5.907487e-05     1.769425e-04
 [4,]     0.0002533557     1.463797e-04     4.259987e-05     2.034289e-04
 [5,]     0.0002686150     3.480534e-04     1.853330e-04     5.273761e-05
 [6,]     0.0002636389     3.430407e-04     1.528201e-04     4.243997e-05
 [7,]     0.0001109944     2.198244e-04     1.998431e-04     1.317537e-04
 [8,]     0.0000988639     2.570901e-04     1.264263e-04     8.261296e-05
 [9,]     0.0002924740     1.258040e-04     5.560219e-05     1.941957e-04
[10,]     0.0002456050     9.031550e-05     1.854846e-05     1.363088e-04
[11,]     0.0002267997     1.321494e-04     1.438200e-04     6.845722e-05
[12,]     0.0002716446     4.866834e-05     1.032435e-04     4.316043e-05
[13,]     0.0005517672     6.310992e-05     3.517649e-05     8.150204e-05
[14,]     0.0006729497     1.879187e-04     3.883435e-05     1.029886e-04
[15,]     0.0003541311     1.292561e-04     7.773488e-05     5.793062e-05
[16,]     0.0003525418     1.293469e-04     6.250734e-05     6.514368e-05
      2632.54298462564 2643.16607949533 2660.07717852527 2672.76734152112
 [1,]     4.426755e-08     0.0002661511      0.003537573     0.0006656668
 [2,]     1.171183e-04     0.0002763753      0.003808063     0.0008060158
 [3,]     1.210669e-04     0.0001670622      0.002101817     0.0006908562
 [4,]     5.917626e-05     0.0001114647      0.002225474     0.0007127513
 [5,]     7.570244e-05     0.0002847497      0.004521747     0.0006062187
 [6,]     1.053248e-04     0.0003115975      0.004720231     0.0005654903
 [7,]     1.181953e-04     0.0002782829      0.004565412     0.0006458589
 [8,]     4.803618e-05     0.0001705004      0.003585267     0.0005417207
 [9,]     9.126657e-05     0.0001588126      0.002254183     0.0009522634
[10,]     2.423906e-05     0.0001451571      0.001908104     0.0008062360
[11,]     1.258460e-05     0.0002679257      0.004490451     0.0008466136
[12,]     2.934194e-05     0.0003934217      0.005149722     0.0009027641
[13,]     8.973417e-08     0.0002490969      0.001350831     0.0010783927
[14,]     6.556419e-07     0.0003322327      0.001467521     0.0011202152
[15,]     3.741848e-06     0.0002211354      0.003521867     0.0009909815
[16,]     2.578881e-05     0.0002514576      0.002720156     0.0008715605
      2682.19711838672 2686.81292226065 2698.68476215796 2709.03904326631
 [1,]     2.616348e-05     5.076770e-05     1.089463e-04     1.831892e-05
 [2,]     1.020577e-04     4.169381e-05     8.042229e-05     3.193942e-05
 [3,]     1.454209e-04     1.882853e-04     1.114359e-04     5.368874e-05
 [4,]     1.332495e-04     1.171913e-04     1.317865e-04     5.685112e-05
 [5,]     1.072447e-04     5.755269e-05     6.108763e-05     6.372412e-06
 [6,]     1.069654e-04     1.097021e-04     1.247434e-04     9.001408e-05
 [7,]     9.916510e-05     1.590203e-04     1.392243e-04     1.152726e-04
 [8,]     1.246572e-04     1.264697e-04     5.868602e-05     8.821805e-05
 [9,]     1.339363e-04     1.352573e-05     7.742228e-05     1.258196e-04
[10,]     7.196892e-05     2.986024e-05     2.565616e-05     1.263072e-04
[11,]     3.897824e-04     2.556863e-04     6.398274e-05     4.431737e-05
[12,]     4.301002e-04     2.703387e-04     9.411577e-05     4.838735e-05
[13,]     7.818987e-05     3.690253e-05     1.240300e-04     1.959866e-04
[14,]     5.005926e-05     1.567375e-08     8.868192e-06     2.279342e-04
[15,]     1.338315e-04     2.089598e-04     5.207338e-05     1.549203e-04
[16,]     2.099077e-04     2.276442e-04     2.479785e-05     1.135057e-04
      2714.06202511835 2723.83906564336 2730.78972177116 2739.43658146684
 [1,]     0.0002209574     1.360015e-04     0.0001969174     2.507704e-04
 [2,]     0.0002525267     1.151346e-04     0.0001582766     2.503589e-04
 [3,]     0.0002126467     9.872852e-05     0.0001537007     1.394096e-04
 [4,]     0.0002314461     1.156227e-04     0.0001735538     1.941606e-04
 [5,]     0.0003137415     3.201888e-05     0.0001986661     1.944254e-04
 [6,]     0.0002263105     3.041794e-05     0.0001870889     2.206798e-04
 [7,]     0.0003192300     1.927063e-04     0.0001978740     1.883629e-04
 [8,]     0.0002929959     8.804252e-05     0.0001458328     2.322899e-04
 [9,]     0.0001738626     1.426799e-04     0.0001491316     2.506701e-04
[10,]     0.0001256693     1.012868e-04     0.0001720564     2.681279e-04
[11,]     0.0002377925     1.413901e-04     0.0004186885     2.287749e-04
[12,]     0.0002919226     2.048972e-04     0.0004492546     1.915118e-04
[13,]     0.0001491596     7.334903e-05     0.0002973193     1.874436e-04
[14,]     0.0001291227     1.670369e-04     0.0003074093     7.109311e-05
[15,]     0.0002438243     1.479155e-04     0.0003669193     4.084168e-04
[16,]     0.0002182037     1.255267e-04     0.0003111488     3.477880e-04
      2754.06948240275 2769.29399575416 2789.91447076973 2808.20313365643
 [1,]     0.0002938698     0.0032741805     0.0005628564     9.784434e-05
 [2,]     0.0002688179     0.0033374330     0.0005129844     7.524510e-05
 [3,]     0.0002299486     0.0023804920     0.0003150285     1.254255e-04
 [4,]     0.0002144061     0.0022632759     0.0002878156     1.160694e-04
 [5,]     0.0008152410     0.0010387203     0.0007271850     2.351779e-05
 [6,]     0.0007629857     0.0011083043     0.0006513702     7.142778e-05
 [7,]     0.0001821824     0.0035670205     0.0007428231     6.465765e-05
 [8,]     0.0002450439     0.0034683981     0.0008220129     6.758818e-05
 [9,]     0.0001248770     0.0038225341     0.0002284599     5.500768e-05
[10,]     0.0001034747     0.0036779890     0.0001604716     2.947301e-06
[11,]     0.0001145298     0.0030947390     0.0010437568     2.029041e-05
[12,]     0.0001110695     0.0032766441     0.0010278558     2.558524e-05
[13,]     0.0012199264     0.0008810005     0.0003208050     1.123841e-04
[14,]     0.0012980289     0.0009271560     0.0002885045     9.760301e-05
[15,]     0.0006592308     0.0034752589     0.0004987017     5.892828e-05
[16,]     0.0006248008     0.0032207349     0.0005122422     4.583278e-05
      2819.46276457694 2830.95253222116 2845.17603302829 2862.36503198417
 [1,]     1.518008e-04     1.054776e-04     0.0001911377     0.0020345473
 [2,]     1.320455e-04     3.854839e-05     0.0002041531     0.0021108712
 [3,]     1.014307e-04     6.119191e-05     0.0001686516     0.0012463829
 [4,]     1.461146e-04     9.738782e-06     0.0001566293     0.0013074987
 [5,]     1.482120e-04     1.921800e-04     0.0002360857     0.0016446294
 [6,]     1.133939e-04     1.547522e-04     0.0002111171     0.0017051932
 [7,]     1.166905e-04     2.454200e-04     0.0002405309     0.0017828221
 [8,]     5.362285e-05     2.029887e-04     0.0001986093     0.0017520747
 [9,]     5.897938e-05     1.274370e-04     0.0002322812     0.0013932054
[10,]     1.494727e-04     5.600089e-05     0.0002541643     0.0013325602
[11,]     7.930960e-05     1.932031e-04     0.0002095002     0.0015948762
[12,]     3.941698e-05     9.838696e-05     0.0001985670     0.0017725358
[13,]     1.773081e-04     6.274176e-05     0.0001906422     0.0004674001
[14,]     6.189232e-05     4.468748e-05     0.0002292164     0.0005526141
[15,]     1.062438e-04     2.164913e-04     0.0004583430     0.0016546644
[16,]     1.120741e-04     2.326766e-04     0.0004060958     0.0013142013
      2876.95950979815 2887.74682900504 2902.46246124011 2915.28322026525
 [1,]     1.085623e-04     1.870847e-04     1.633247e-04     1.675985e-04
 [2,]     3.934910e-05     2.153735e-04     2.169954e-04     2.078474e-04
 [3,]     1.143256e-04     1.831010e-04     1.152528e-04     1.108577e-04
 [4,]     1.649391e-04     1.940446e-04     1.671966e-04     1.469708e-04
 [5,]     4.025379e-05     1.687393e-04     2.216160e-04     1.648064e-04
 [6,]     7.918250e-07     1.547637e-04     1.944685e-04     1.056152e-04
 [7,]     5.423956e-05     5.291189e-05     1.288727e-04     6.057524e-05
 [8,]     9.129195e-05     1.546019e-04     2.225966e-04     1.344756e-04
 [9,]     1.010083e-04     1.883435e-04     1.146404e-04     1.372236e-04
[10,]     1.531355e-05     1.876426e-04     1.778163e-04     1.314390e-04
[11,]     8.675392e-05     1.823885e-04     2.279829e-04     1.941893e-04
[12,]     4.349605e-05     1.627686e-04     1.432405e-04     1.724459e-04
[13,]     2.831985e-05     1.609034e-04     8.987887e-05     9.389974e-05
[14,]     2.050682e-05     1.262065e-04     2.916265e-05     8.419264e-05
[15,]     9.605152e-06     1.415825e-04     1.007252e-04     1.090214e-04
[16,]     2.900238e-05     1.920104e-04     4.917274e-05     2.054325e-04
      2932.33770504227 2952.54779084813 2972.57761432583 2978.5174561921
 [1,]      0.003938789     0.0027973920     8.298786e-05    2.161455e-04
 [2,]      0.003851041     0.0027495412     9.734718e-05    2.846605e-04
 [3,]      0.003655281     0.0026152556     1.670927e-04    1.236085e-04
 [4,]      0.003487755     0.0026051759     1.035926e-04    2.267035e-04
 [5,]      0.002835269     0.0023102773     9.788339e-05    2.088281e-04
 [6,]      0.002709161     0.0022215284     1.444758e-05    2.283804e-04
 [7,]      0.002883580     0.0018874737     4.942958e-05    7.065404e-05
 [8,]      0.002961211     0.0022646446     1.346102e-04    1.312481e-04
 [9,]      0.003268424     0.0024338989     1.005113e-04    7.919477e-05
[10,]      0.003238496     0.0024733365     9.989191e-05    1.212249e-04
[11,]      0.002777936     0.0021314957     1.210403e-04    8.882970e-05
[12,]      0.002884061     0.0019337926     1.302882e-04    1.412509e-04
[13,]      0.002325408     0.0009732434     8.077474e-05    3.149833e-04
[14,]      0.002371288     0.0010912111     7.461827e-05    3.373581e-04
[15,]      0.004284750     0.0018152710     1.016260e-04    3.183310e-04
[16,]      0.004018639     0.0016576365     5.436542e-05    3.611261e-04
      2990.57378388783 3002.09428002708 3016.7275646013 3035.63652710832
 [1,]     0.0004424794     3.680037e-05    1.035137e-04     1.541661e-05
 [2,]     0.0003980588     1.444157e-04    7.454831e-05     4.670326e-05
 [3,]     0.0003605540     1.182972e-04    1.173917e-04     4.888562e-05
 [4,]     0.0004507546     9.701265e-05    9.887601e-05     1.198595e-04
 [5,]     0.0004221049     2.434265e-04    4.470081e-05     7.054840e-05
 [6,]     0.0005068402     2.272165e-04    1.481157e-04     2.269475e-05
 [7,]     0.0003468949     1.444292e-04    9.918232e-05     1.164186e-04
 [8,]     0.0003944192     1.740683e-04    9.543320e-05     1.535679e-05
 [9,]     0.0006523188     1.086926e-04    1.176932e-04     1.893876e-05
[10,]     0.0005851036     1.133242e-04    1.030414e-04     1.313868e-05
[11,]     0.0008728312     1.715446e-04    2.364791e-04     3.111390e-05
[12,]     0.0009643621     1.612642e-04    1.647322e-04     4.203681e-05
[13,]     0.0001437295     7.437895e-05    1.812347e-04     1.221769e-04
[14,]     0.0002025958     1.044672e-04    1.622769e-04     3.267958e-05
[15,]     0.0005430047     2.384499e-06    2.154336e-04     7.302449e-05
[16,]     0.0004234705     4.542415e-05    2.687996e-04     1.100427e-04
      3045.08697912476 3059.82037561276 3075.7859028279 3087.66973033219
 [1,]     0.0002128342     1.214230e-04    5.164359e-05     9.713146e-05
 [2,]     0.0002413207     1.326898e-04    1.518626e-04     9.513298e-05
 [3,]     0.0002335237     1.669515e-04    4.132943e-06     1.179071e-04
 [4,]     0.0002434758     1.328135e-04    9.565711e-05     8.275215e-05
 [5,]     0.0001658738     1.179488e-04    9.714517e-05     1.792807e-04
 [6,]     0.0001992451     7.301297e-05    7.981919e-05     1.289830e-04
 [7,]     0.0001397244     8.296167e-05    1.489837e-05     1.351988e-04
 [8,]     0.0001919628     6.087040e-05    0.000000e+00     8.773319e-05
 [9,]     0.0002191422     9.088242e-05    2.272506e-05     1.224713e-04
[10,]     0.0001837115     1.338026e-04    3.781776e-05     9.422613e-05
[11,]     0.0001397604     1.829637e-04    5.651764e-05     1.163002e-04
[12,]     0.0001446239     1.093318e-04    7.168068e-06     9.608191e-05
[13,]     0.0001232200     8.305125e-05    1.087574e-04     1.467263e-04
[14,]     0.0001053160     4.844503e-05    1.191845e-04     1.386740e-04
[15,]     0.0002125417     1.352597e-04    4.856194e-05     3.104567e-05
[16,]     0.0002182373     6.014725e-05    9.292666e-05     4.329613e-05
      3096.10898387045 3104.49234242311 3117.84333653732 3143.32635463077
 [1,]     0.0002920887     1.161782e-04     0.0002166238     0.0002702262
 [2,]     0.0003423394     1.572251e-04     0.0002739520     0.0002755580
 [3,]     0.0004768277     1.012784e-04     0.0002468459     0.0002297675
 [4,]     0.0005428548     1.218804e-04     0.0002136075     0.0002236279
 [5,]     0.0003978804     1.485268e-04     0.0002156888     0.0002528043
 [6,]     0.0002746296     1.253898e-04     0.0001433648     0.0002898431
 [7,]     0.0002297750     1.069505e-04     0.0001204387     0.0002667662
 [8,]     0.0003920035     8.742062e-05     0.0001876731     0.0002286207
 [9,]     0.0004202568     1.485762e-04     0.0002305567     0.0003436511
[10,]     0.0004458420     1.154436e-04     0.0002362219     0.0002668743
[11,]     0.0002461071     1.319038e-04     0.0002312322     0.0003451925
[12,]     0.0001451143     1.536237e-04     0.0002109391     0.0004045032
[13,]     0.0003295599     1.655150e-04     0.0002055697     0.0002034376
[14,]     0.0004183712     1.265574e-04     0.0002319091     0.0001916089
[15,]     0.0002260938     2.100762e-04     0.0005421599     0.0002216822
[16,]     0.0001808643     1.735409e-04     0.0004726154     0.0002124413
      3158.62751162253 3175.40845486099 3191.72121064019 3207.90056647982
 [1,]     0.0002051099     1.686318e-04      0.005304437     1.459014e-04
 [2,]     0.0001875835     1.862909e-04      0.005188222     1.382578e-04
 [3,]     0.0002364519     1.597982e-04      0.004730564     1.087619e-04
 [4,]     0.0003104078     1.592646e-04      0.004525649     1.797575e-04
 [5,]     0.0002119085     1.779859e-04      0.004155132     3.194098e-05
 [6,]     0.0002669684     1.119208e-04      0.004050205     7.007577e-05
 [7,]     0.0001636922     1.698545e-04      0.004202110     3.106053e-05
 [8,]     0.0001911078     1.709600e-04      0.004288897     8.815848e-05
 [9,]     0.0001097417     1.734571e-04      0.005069016     1.504947e-04
[10,]     0.0001925065     1.734951e-04      0.005180290     1.190843e-04
[11,]     0.0002327538     1.811956e-04      0.004075057     1.332962e-04
[12,]     0.0002168523     9.614572e-05      0.004162858     1.423376e-04
[13,]     0.0005456134     6.651359e-05      0.002382513     1.876028e-04
[14,]     0.0006249138     1.059363e-04      0.002435514     1.215845e-04
[15,]     0.0003252258     2.828152e-05      0.003113090     1.453585e-04
[16,]     0.0003137503     1.073036e-04      0.003020691     8.057664e-05
      3217.0855720029 3231.96956065605 3240.80749728491 3262.80152524051
 [1,]    2.028216e-04     1.334515e-04      0.002234515      0.007579381
 [2,]    1.956495e-04     1.369544e-04      0.002219740      0.007352572
 [3,]    1.719491e-04     1.271249e-04      0.001725685      0.006991991
 [4,]    2.778992e-04     1.500083e-04      0.001772405      0.006641271
 [5,]    2.606786e-04     1.500280e-04      0.002225312      0.006213153
 [6,]    1.839817e-04     9.587363e-05      0.002127383      0.006143984
 [7,]    1.042224e-04     7.021253e-05      0.002190770      0.005975066
 [8,]    2.417218e-04     4.942335e-05      0.002182222      0.006002429
 [9,]    2.229729e-04     6.023166e-05      0.002122166      0.006702627
[10,]    2.174623e-04     1.092359e-04      0.002015544      0.006895259
[11,]    7.790827e-05     6.117395e-05      0.003398792      0.005557398
[12,]    9.500797e-05     8.029861e-05      0.003672830      0.005849704
[13,]    3.104555e-04     2.157795e-05      0.001031333      0.003120228
[14,]    3.558467e-04     8.315373e-05      0.001118675      0.003309616
[15,]    1.281633e-04     4.845635e-05      0.002698777      0.003102633
[16,]    1.260668e-04     0.000000e+00      0.002487427      0.002969711
      3279.24494158417 3294.33551981694 3303.23597838704 3316.02792610326
 [1,]     0.0003359546     1.178626e-04     1.508510e-04     0.0003108147
 [2,]     0.0002698837     1.604642e-04     2.002365e-04     0.0004567797
 [3,]     0.0003081351     6.342827e-05     2.668230e-04     0.0003539394
 [4,]     0.0003210338     1.254242e-04     1.887002e-04     0.0005981134
 [5,]     0.0002261317     1.130603e-04     1.053269e-04     0.0005940681
 [6,]     0.0002236354     5.069033e-05     1.395685e-04     0.0004723192
 [7,]     0.0002043225     1.170829e-04     1.411418e-04     0.0003233411
 [8,]     0.0001851841     7.657641e-05     1.268335e-04     0.0004140523
 [9,]     0.0002049039     6.473771e-05     1.293844e-04     0.0005287597
[10,]     0.0002618303     9.161647e-05     9.425853e-05     0.0004755743
[11,]     0.0002451950     1.461385e-04     1.864819e-04     0.0003331456
[12,]     0.0001562907     1.563313e-04     1.572905e-04     0.0002614848
[13,]     0.0001144251     1.093855e-04     1.073594e-04     0.0006857765
[14,]     0.0002074661     5.573348e-05     2.224333e-05     0.0007656296
[15,]     0.0002262508     9.153663e-05     1.206989e-04     0.0003835202
[16,]     0.0001998820     1.449401e-04     1.741464e-04     0.0004561342
      3325.27556213481 3334.41156997893 3365.63706580383 3377.6603556791
 [1,]     7.378432e-05     1.479906e-04     0.0002516389    0.0004933789
 [2,]     1.734369e-04     1.236188e-04     0.0002800737    0.0005177855
 [3,]     2.854168e-04     5.390627e-05     0.0002285605    0.0004619338
 [4,]     4.386177e-05     1.283184e-04     0.0002253612    0.0004603760
 [5,]     1.322831e-04     1.351438e-04     0.0002858200    0.0003920655
 [6,]     1.040929e-04     2.009602e-04     0.0002617051    0.0003043497
 [7,]     1.748600e-04     8.832457e-05     0.0002774224    0.0003006107
 [8,]     1.247374e-04     1.050520e-04     0.0003173476    0.0003764848
 [9,]     2.322615e-04     6.219847e-05     0.0002585099    0.0005214265
[10,]     9.633249e-05     9.436199e-05     0.0002961662    0.0005450427
[11,]     1.050002e-04     7.442936e-05     0.0002635710    0.0003596776
[12,]     1.196900e-04     2.849773e-05     0.0002699502    0.0003552182
[13,]     2.909603e-06     2.247062e-04     0.0002832186    0.0002140715
[14,]     5.756695e-05     1.962065e-04     0.0002787968    0.0003173934
[15,]     7.709932e-05     2.142030e-04     0.0002178271    0.0002661428
[16,]     1.119729e-04     2.336847e-04     0.0002370202    0.0002385382
      3408.43059823469 3413.84696290897 3419.96586728355 3434.40827852602
 [1,]     1.868905e-04     7.428708e-05     1.530092e-04     2.807565e-04
 [2,]     1.709129e-04     1.381483e-04     1.582002e-04     1.913403e-04
 [3,]     1.521975e-04     1.082721e-04     1.643666e-04     1.790011e-04
 [4,]     1.563481e-04     3.619751e-05     1.688414e-04     1.740910e-04
 [5,]     2.011409e-04     1.406972e-04     1.350798e-04     2.007371e-04
 [6,]     1.766116e-04     8.402382e-05     1.351538e-04     2.072723e-04
 [7,]     1.293129e-04     1.766041e-04     1.416681e-04     9.835969e-05
 [8,]     1.558671e-04     9.414690e-05     5.868462e-05     2.106584e-04
 [9,]     1.618691e-04     6.027151e-05     9.479298e-05     9.731384e-05
[10,]     1.229264e-04     9.635213e-05     6.678300e-05     1.547620e-04
[11,]     1.057917e-04     1.134005e-04     5.714679e-05     1.287933e-04
[12,]     1.070273e-04     9.114206e-05     6.682017e-05     1.218741e-04
[13,]     1.272100e-04     1.116087e-04     1.571343e-04     1.302262e-04
[14,]     1.721920e-04     1.178688e-04     4.939824e-05     1.359855e-04
[15,]     5.908076e-05     1.777447e-04     1.126835e-04     6.342584e-05
[16,]     1.337623e-04     8.072569e-05     1.041765e-04     7.912484e-05
      3448.20196997724 3466.1202754243 3489.13005197377 3506.14680512151
 [1,]     0.0009496084    1.197715e-04     6.130098e-05     5.701097e-05
 [2,]     0.0009608224    1.146746e-04     5.852842e-05     4.758720e-05
 [3,]     0.0008009206    1.280688e-04     1.244533e-04     6.315811e-05
 [4,]     0.0007960359    1.136885e-04     1.051758e-04     1.600311e-04
 [5,]     0.0007664039    1.047418e-04     2.539723e-05     1.244876e-04
 [6,]     0.0008226790    9.469580e-05     1.576648e-05     1.270146e-04
 [7,]     0.0006747301    1.258394e-04     6.472584e-05     1.389105e-04
 [8,]     0.0007722087    6.581344e-05     4.690050e-05     1.137101e-04
 [9,]     0.0007870383    6.010385e-05     5.166088e-05     1.405005e-04
[10,]     0.0008159099    9.913773e-05     5.157139e-05     1.945955e-04
[11,]     0.0006605950    2.674704e-05     1.191449e-04     2.454994e-04
[12,]     0.0006474472    9.687109e-06     2.447636e-06     2.076993e-04
[13,]     0.0003925269    1.012863e-06     2.630986e-05     5.789062e-05
[14,]     0.0004047804    5.868777e-05     6.580916e-05     3.425513e-05
[15,]     0.0003788736    3.847914e-05     5.940983e-05     1.638882e-04
[16,]     0.0003681003    1.303537e-04     1.104330e-04     2.024041e-04
      3524.49712185398 3540.24266495582 3568.70957577374 3585.14781658734
 [1,]     0.0001904684     0.0004064625     2.022796e-05     5.082702e-05
 [2,]     0.0002021671     0.0004623092     1.084645e-04     1.099476e-04
 [3,]     0.0002283624     0.0004058314     7.038467e-05     4.157850e-05
 [4,]     0.0002314173     0.0003460820     3.815141e-05     5.454676e-05
 [5,]     0.0002494383     0.0004187000     2.803500e-05     1.294046e-04
 [6,]     0.0002135500     0.0004589310     1.014466e-04     1.697214e-04
 [7,]     0.0001251579     0.0002978926     2.241576e-05     6.544272e-05
 [8,]     0.0001589518     0.0002229461     6.451793e-05     1.076786e-04
 [9,]     0.0002700401     0.0001552558     8.701697e-05     1.314897e-04
[10,]     0.0002218577     0.0001659572     5.130904e-05     1.413023e-04
[11,]     0.0002709470     0.0003452944     5.256874e-05     3.613177e-04
[12,]     0.0003007792     0.0003533584     1.133948e-04     3.672112e-04
[13,]     0.0005451610     0.0008420193     7.392150e-05     2.770261e-05
[14,]     0.0005688845     0.0008441930     6.414824e-05     7.601811e-05
[15,]     0.0005003630     0.0005416031     1.434548e-04     3.526336e-04
[16,]     0.0004832466     0.0004999740     1.221708e-04     2.294353e-04
      3605.54707105848 3625.18037536826 3634.95540515964 3676.09975876223
 [1,]     5.062396e-05     6.610228e-05     2.104911e-04     1.391591e-04
 [2,]     4.705888e-05     5.947652e-05     2.623146e-04     9.866208e-05
 [3,]     1.172156e-04     1.465909e-04     3.040396e-04     1.262924e-04
 [4,]     1.540309e-04     1.129245e-04     2.899335e-04     1.502591e-04
 [5,]     1.044126e-04     1.174632e-04     2.392446e-04     1.280696e-04
 [6,]     1.042188e-04     1.316073e-04     2.887691e-04     1.568984e-04
 [7,]     1.258334e-04     8.300264e-05     2.579447e-04     1.412899e-04
 [8,]     3.014669e-05     6.834376e-05     2.111586e-04     1.853218e-04
 [9,]     9.751056e-05     6.902880e-05     3.118498e-04     2.880321e-05
[10,]     2.383677e-05     5.113986e-05     2.471992e-04     5.302859e-05
[11,]     1.163733e-04     3.563910e-05     2.638889e-04     2.083050e-05
[12,]     4.179218e-05     5.515167e-05     2.219615e-04     5.935285e-05
[13,]     1.603659e-04     8.554956e-05     5.202247e-05     1.122045e-04
[14,]     6.390851e-06     4.150248e-05     2.351909e-04     1.014515e-04
[15,]     1.129564e-04     1.232148e-04     1.950012e-04     2.892333e-05
[16,]     1.623202e-04     7.272245e-05     2.367470e-04     3.039698e-05
      3680.95667927439 3704.58015419771 3719.56244085583 3734.74097421538
 [1,]     6.530373e-05     3.302681e-06     6.072070e-05     1.148080e-04
 [2,]     1.021884e-04     2.510782e-05     6.333805e-05     1.371578e-04
 [3,]     8.711543e-05     1.199399e-04     1.065247e-04     9.997326e-05
 [4,]     1.309923e-04     9.687881e-05     3.935396e-05     6.615163e-05
 [5,]     1.508803e-04     5.483369e-05     1.053696e-04     1.842816e-04
 [6,]     1.362887e-04     1.297767e-04     9.813203e-08     1.363037e-04
 [7,]     1.785284e-04     5.042255e-05     1.223617e-04     1.018161e-04
 [8,]     1.391204e-04     1.281269e-04     6.583517e-05     1.350507e-04
 [9,]     9.264524e-05     2.339953e-05     6.742895e-05     7.033556e-05
[10,]     9.376860e-05     5.028467e-05     1.074614e-05     1.215826e-04
[11,]     1.420212e-04     1.053244e-07     3.298597e-05     1.477730e-04
[12,]     1.014226e-04     1.336812e-05     1.789032e-05     6.811064e-05
[13,]     1.976343e-04     6.227420e-05     1.245260e-04     1.175157e-04
[14,]     2.003606e-04     3.806239e-05     5.644300e-05     8.628692e-05
[15,]     1.243771e-04     1.275335e-04     1.442423e-05     1.107336e-04
[16,]     1.271938e-04     4.606270e-05     4.081496e-05     1.230661e-04
      3744.04677948767 3763.75996298824 3768.64436784783 3782.78260788104
 [1,]     6.175389e-05     1.390584e-04     7.871198e-05     1.192336e-04
 [2,]     4.799663e-05     1.055474e-04     4.981512e-05     1.158782e-04
 [3,]     1.431203e-04     1.377450e-04     7.649917e-05     1.166793e-04
 [4,]     6.273323e-05     1.013141e-04     3.612966e-05     1.401536e-04
 [5,]     7.250670e-05     1.846388e-04     8.224700e-05     7.328214e-05
 [6,]     1.385823e-04     2.005389e-04     7.143812e-05     6.122090e-05
 [7,]     1.334105e-04     2.044975e-04     6.830346e-05     3.010282e-05
 [8,]     1.559986e-04     2.773041e-04     4.157186e-05     6.755296e-05
 [9,]     4.804812e-05     1.242018e-04     1.079457e-04     1.099434e-04
[10,]     5.916591e-05     4.521693e-06     9.325401e-05     2.587247e-05
[11,]     4.685937e-05     1.690949e-04     5.463263e-04     2.966046e-05
[12,]     2.259568e-05     1.612026e-04     5.971921e-04     4.488422e-07
[13,]     1.406128e-05     1.731710e-04     1.411085e-04     7.383870e-05
[14,]     1.166280e-04     1.435916e-04     9.770658e-05     1.456862e-04
[15,]     6.131930e-05     1.363875e-04     2.375100e-04     1.089468e-04
[16,]     3.706763e-05     1.886615e-04     2.431045e-04     5.275963e-05
      3795.35706752575 3816.79231944162 3860.96809262121 3874.16235898778
 [1,]     1.186627e-04     0.0001849383     3.650602e-05     0.0004060549
 [2,]     6.595890e-07     0.0001704351     7.880031e-05     0.0003951679
 [3,]     4.333439e-05     0.0002440809     1.938117e-04     0.0003547629
 [4,]     5.407767e-05     0.0001938839     1.422483e-04     0.0003277201
 [5,]     1.169078e-04     0.0001563735     1.080394e-04     0.0002949292
 [6,]     5.067673e-05     0.0002116955     1.598707e-04     0.0002326361
 [7,]     6.355461e-05     0.0001358603     1.236985e-04     0.0004136348
 [8,]     5.014514e-05     0.0001999945     1.376331e-04     0.0004360146
 [9,]     1.616670e-04     0.0001336715     9.503125e-05     0.0001967385
[10,]     1.042670e-04     0.0001656033     1.354050e-04     0.0002434770
[11,]     1.162500e-04     0.0001185548     5.728244e-05     0.0002183300
[12,]     1.956570e-04     0.0001004099     8.903313e-05     0.0001822049
[13,]     1.398099e-04     0.0003100427     1.627900e-04     0.0001498235
[14,]     1.689729e-04     0.0003408862     1.419398e-04     0.0002393552
[15,]     1.154293e-04     0.0002503599     1.171160e-04     0.0001616363
[16,]     1.076449e-04     0.0002772144     1.289764e-04     0.0001281068
      3882.93414759827 3903.2878538192 3909.62570592724 3918.98433918382
 [1,]      0.002186050    9.833565e-05     1.214554e-04     1.806672e-04
 [2,]      0.002216933    1.308252e-04     2.136385e-04     1.646736e-04
 [3,]      0.002357195    1.916364e-04     1.603698e-04     1.737438e-04
 [4,]      0.002650211    1.852947e-04     2.866344e-04     2.199919e-04
 [5,]      0.001911640    1.340526e-04     1.478486e-04     1.278609e-04
 [6,]      0.001870409    5.884042e-05     1.498566e-04     1.208507e-04
 [7,]      0.002073204    1.182780e-04     2.030171e-04     1.852587e-04
 [8,]      0.002448330    1.316905e-04     1.928584e-04     1.079577e-04
 [9,]      0.002022032    7.801032e-05     1.756376e-04     1.562381e-04
[10,]      0.002029527    1.121256e-04     1.510718e-04     1.621489e-04
[11,]      0.001585821    1.048136e-04     1.144387e-04     1.293065e-04
[12,]      0.001338594    5.138494e-05     6.451248e-05     1.415866e-04
[13,]      0.002083500    1.826684e-04     1.862952e-04     1.372231e-04
[14,]      0.002219554    2.038109e-04     1.604955e-04     5.506598e-05
[15,]      0.001631029    1.309391e-04     1.302800e-04     1.417423e-04
[16,]      0.001534938    1.067808e-04     6.604414e-05     7.607880e-05
      3934.21297348999 3952.56089161257 3961.83127480521 3993.46801555899
 [1,]     0.0001301017     1.275965e-04     0.0002017873     1.060107e-04
 [2,]     0.0001230168     1.440596e-04     0.0002126868     5.189753e-05
 [3,]     0.0001634939     1.467818e-04     0.0002137947     4.077726e-05
 [4,]     0.0001272958     1.345605e-04     0.0002108300     1.054277e-04
 [5,]     0.0002613717     1.615447e-04     0.0002387653     1.139835e-04
 [6,]     0.0002251300     1.309166e-04     0.0001884305     6.743462e-05
 [7,]     0.0001401680     1.252061e-04     0.0001998188     1.078735e-04
 [8,]     0.0001621973     1.601234e-04     0.0002080579     2.442139e-05
 [9,]     0.0002956121     2.446977e-05     0.0001718788     1.115523e-04
[10,]     0.0002848260     1.289514e-04     0.0001638940     4.421211e-06
[11,]     0.0003609043     1.263062e-04     0.0001379955     5.792570e-05
[12,]     0.0004039233     9.858818e-05     0.0001128373     4.476945e-05
[13,]     0.0004581724     3.707576e-05     0.0002194368     8.680932e-05
[14,]     0.0004945893     5.901910e-05     0.0002186345     7.109668e-05
[15,]     0.0002979041     1.070107e-04     0.0001746811     1.269116e-04
[16,]     0.0003655538     1.545080e-04     0.0001534663     1.027210e-04
      4019.24209316322 4036.89868748243 4053.87685443975 4063.52812679879
 [1,]     1.729008e-04     1.274492e-04     0.0004049548     4.246586e-05
 [2,]     1.520991e-04     1.100461e-04     0.0003873266     7.270726e-05
 [3,]     1.390519e-04     1.218555e-04     0.0003878726     1.503034e-04
 [4,]     1.494777e-04     4.191038e-05     0.0003819637     8.747854e-05
 [5,]     1.113079e-04     1.035390e-04     0.0009290330     1.199887e-04
 [6,]     1.096462e-04     8.976160e-05     0.0008973565     1.111993e-04
 [7,]     1.722787e-04     1.059603e-04     0.0005860883     7.317429e-05
 [8,]     1.165627e-04     1.165947e-04     0.0007515631     3.394467e-05
 [9,]     1.225738e-04     1.005960e-04     0.0011776317     1.393768e-04
[10,]     9.289061e-05     1.179187e-04     0.0012134568     1.252144e-04
[11,]     1.400956e-04     1.482056e-04     0.0014081468     1.794307e-04
[12,]     1.561446e-04     1.178325e-04     0.0015858733     1.661441e-04
[13,]     1.068830e-04     1.556608e-04     0.0015987611     7.280709e-04
[14,]     1.161956e-04     1.133329e-04     0.0015195681     7.958920e-04
[15,]     5.370693e-05     1.588932e-04     0.0008824790     4.335262e-04
[16,]     1.591526e-04     4.809712e-05     0.0008663292     3.891836e-04
      4071.56137672196 4091.28327473639 4110.14070456692 4123.03360546748
 [1,]     0.0004315673      0.001248159     3.050394e-07     0.0001619631
 [2,]     0.0003922276      0.001311819     6.768872e-06     0.0001847762
 [3,]     0.0004103585      0.001363917     6.976036e-05     0.0002167023
 [4,]     0.0004994564      0.001331316     3.768849e-05     0.0001798320
 [5,]     0.0003518900      0.001659261     1.318641e-05     0.0003850546
 [6,]     0.0003313661      0.001621391     1.167492e-04     0.0004434966
 [7,]     0.0003204449      0.001266535     3.240964e-05     0.0002104265
 [8,]     0.0004180778      0.001560510     3.729001e-05     0.0002796693
 [9,]     0.0002936314      0.001468808     1.110110e-04     0.0001566239
[10,]     0.0002990908      0.001600597     1.017793e-04     0.0001911188
[11,]     0.0001493632      0.001761646     1.634931e-04     0.0002036648
[12,]     0.0001564508      0.001789221     8.964431e-05     0.0002134345
[13,]     0.0003428179      0.001420248     1.660093e-04     0.0002212043
[14,]     0.0003538670      0.001508371     1.636111e-04     0.0002122252
[15,]     0.0002178054      0.001419965     1.277095e-04     0.0001109314
[16,]     0.0002172062      0.001386046     1.309375e-04     0.0001059146
      4145.67241202184 4152.80101824032 4166.88592458118 4173.34919307908
 [1,]     1.051130e-04     3.799127e-05     0.0001297208     8.296044e-05
 [2,]     9.898931e-05     3.906466e-05     0.0001634838     1.304177e-04
 [3,]     8.930942e-05     1.113179e-04     0.0001491960     5.709919e-05
 [4,]     1.079913e-04     9.682765e-05     0.0001199826     1.262513e-04
 [5,]     1.454244e-04     9.860001e-05     0.0001615580     1.350094e-04
 [6,]     1.682114e-04     9.504287e-07     0.0002125679     5.715154e-05
 [7,]     7.697788e-05     6.754024e-06     0.0001479840     1.168246e-04
 [8,]     7.688057e-05     3.143200e-05     0.0001709451     9.335072e-05
 [9,]     9.365469e-05     5.043100e-05     0.0001879037     9.220178e-05
[10,]     5.794988e-05     3.858389e-05     0.0001566531     9.586917e-05
[11,]     1.127768e-04     1.165856e-04     0.0001586421     6.657914e-05
[12,]     1.060285e-04     6.413264e-05     0.0002356264     1.522895e-04
[13,]     4.985150e-06     9.385538e-05     0.0001924734     4.747097e-05
[14,]     4.823436e-05     6.081305e-05     0.0001164482     1.080565e-04
[15,]     3.432946e-05     7.457832e-05     0.0001448963     1.512696e-04
[16,]     2.061961e-05     1.257671e-04     0.0001063621     7.873616e-05
      4192.54093829617 4209.96272693358 4235.38885389076 4249.60616500285
 [1,]     0.0002062974      0.003340931     8.214948e-05     0.0001855307
 [2,]     0.0002609073      0.003267795     5.936543e-05     0.0001848799
 [3,]     0.0002715186      0.003573256     1.293700e-04     0.0001656909
 [4,]     0.0003045061      0.003475206     1.021090e-04     0.0001733511
 [5,]     0.0003542119      0.004576164     1.364678e-04     0.0002361554
 [6,]     0.0004028670      0.004366179     7.451039e-05     0.0001367199
 [7,]     0.0003209115      0.003787964     1.023481e-04     0.0001789170
 [8,]     0.0003298610      0.004153320     1.152590e-04     0.0001848074
 [9,]     0.0003606409      0.004273747     1.861509e-04     0.0002139574
[10,]     0.0003627109      0.004413050     1.235926e-04     0.0002860834
[11,]     0.0004011702      0.005017185     1.689452e-04     0.0002234005
[12,]     0.0004552599      0.005169507     1.459362e-04     0.0002967366
[13,]     0.0003103667      0.004145229     9.291744e-05     0.0001913713
[14,]     0.0002793309      0.004126192     6.901553e-05     0.0002504001
[15,]     0.0002313901      0.002932236     1.344968e-04     0.0001557128
[16,]     0.0002236560      0.003169481     3.802088e-05     0.0001848384
      4266.14514653073 4301.88531578493 4337.88827949057 4354.64856011389
 [1,]     0.0006212864     0.0002638961     1.442479e-04     2.713724e-05
 [2,]     0.0006176165     0.0002455455     1.158295e-04     1.450040e-04
 [3,]     0.0006470418     0.0003291351     1.584050e-05     1.317227e-04
 [4,]     0.0006111320     0.0002607807     3.569428e-05     8.265893e-05
 [5,]     0.0008588819     0.0002140757     2.834893e-05     1.421530e-04
 [6,]     0.0008047625     0.0002201217     5.597821e-05     1.798516e-04
 [7,]     0.0007163447     0.0001466654     5.689017e-05     1.219495e-04
 [8,]     0.0008011721     0.0001808546     1.257300e-04     1.447289e-04
 [9,]     0.0009859699     0.0003045707     1.041537e-04     1.105084e-04
[10,]     0.0009887460     0.0003422769     1.445543e-04     2.350784e-05
[11,]     0.0010030804     0.0005578741     2.220888e-05     5.687060e-05
[12,]     0.0010709045     0.0004294555     6.631237e-05     1.041568e-04
[13,]     0.0007360542     0.0006006867     1.851698e-04     1.360283e-04
[14,]     0.0007816228     0.0006476959     1.514094e-04     6.717997e-05
[15,]     0.0007115186     0.0006194366     1.268443e-04     1.344461e-04
[16,]     0.0007773354     0.0006041530     7.494501e-05     1.348470e-04
      4363.0248725225 4382.47913715518 4397.03853101479 4416.70956743285
 [1,]    0.0001631619     5.639219e-05     1.575863e-04     1.426639e-04
 [2,]    0.0001535885     2.997672e-05     1.117858e-04     5.628198e-05
 [3,]    0.0001680826     1.149663e-04     1.987687e-04     1.105343e-04
 [4,]    0.0001822221     1.089702e-04     2.283673e-04     1.453620e-04
 [5,]    0.0002161427     1.121762e-04     1.198794e-04     1.183578e-04
 [6,]    0.0002177409     1.323331e-04     1.301748e-04     1.107147e-04
 [7,]    0.0001592930     9.676335e-05     5.699920e-05     6.417149e-05
 [8,]    0.0001878740     4.066980e-05     4.948589e-05     7.202582e-05
 [9,]    0.0002522432     1.114405e-04     2.097229e-04     9.776404e-05
[10,]    0.0001868669     9.259516e-05     9.980906e-05     6.574791e-05
[11,]    0.0001837841     9.414772e-05     1.728099e-04     1.333051e-04
[12,]    0.0002617401     6.595108e-05     1.148037e-04     1.098130e-04
[13,]    0.0001210266     1.011450e-04     1.668421e-04     3.985624e-05
[14,]    0.0000980900     8.521306e-05     1.410860e-04     4.274617e-05
[15,]    0.0001663975     7.176971e-05     2.161804e-04     1.052339e-04
[16,]    0.0001459555     5.219433e-05     2.832833e-04     3.786212e-05
      4430.78313009185 4452.02787263112 4466.98888410858 4493.8217253974
 [1,]     1.573797e-04     3.945038e-07     0.0004253456    6.552370e-05
 [2,]     2.009879e-04     1.110822e-05     0.0004059392    3.773356e-05
 [3,]     1.796283e-04     1.064449e-04     0.0003960915    1.365259e-04
 [4,]     2.951253e-05     9.802872e-05     0.0005167025    1.159648e-04
 [5,]     1.314152e-04     2.291525e-06     0.0005671780    3.860534e-05
 [6,]     9.365619e-05     1.465245e-04     0.0004292988    5.632119e-05
 [7,]     1.054540e-04     4.891080e-06     0.0004821108    9.924032e-05
 [8,]     1.197527e-04     8.131052e-06     0.0006242013    8.613390e-05
 [9,]     1.446653e-04     2.702852e-05     0.0004311018    9.688777e-06
[10,]     1.366352e-04     3.165854e-05     0.0004260966    6.414382e-05
[11,]     9.830214e-05     1.797474e-05     0.0006439478    8.493693e-05
[12,]     5.351820e-05     1.339331e-05     0.0005268177    6.032343e-05
[13,]     1.938554e-04     1.324993e-04     0.0020367307    2.518461e-04
[14,]     2.739004e-04     7.666419e-05     0.0019219363    2.739447e-04
[15,]     8.551294e-05     2.081462e-05     0.0013901127    2.717453e-04
[16,]     9.198962e-05     3.212655e-05     0.0015549979    2.475648e-04
      4530.40799218225 4567.34533543643 4586.41456787549 4599.51839509017
 [1,]     0.0002560050     1.154731e-04     1.014072e-04     6.456217e-05
 [2,]     0.0002762637     1.121143e-04     1.170088e-04     1.155984e-04
 [3,]     0.0002117529     4.741229e-05     1.034922e-04     6.314246e-05
 [4,]     0.0003379185     9.484982e-05     1.426106e-04     2.966808e-05
 [5,]     0.0003783846     1.178855e-04     8.148882e-05     3.049817e-05
 [6,]     0.0002275866     6.022317e-05     1.476112e-04     1.029085e-04
 [7,]     0.0003858762     7.788413e-05     7.225992e-05     3.557589e-05
 [8,]     0.0004309206     6.054134e-05     1.140382e-04     6.951104e-05
 [9,]     0.0003356630     1.001824e-04     1.277319e-04     4.062625e-05
[10,]     0.0002865310     5.406944e-05     1.845048e-05     6.860488e-05
[11,]     0.0001660770     1.058581e-04     4.300209e-05     2.751060e-05
[12,]     0.0001778229     7.673094e-05     4.275101e-05     1.204313e-05
[13,]     0.0003347139     1.668659e-04     9.300480e-05     1.060601e-04
[14,]     0.0003833485     1.525124e-04     4.034429e-05     5.611396e-05
[15,]     0.0001681661     8.466315e-05     8.395453e-06     1.096369e-04
[16,]     0.0002172601     1.433792e-04     1.157772e-04     1.466037e-04
      4613.85025288703 4644.42111389705 4663.77171872873 4671.40987352484
 [1,]     3.442870e-05      0.002416619     1.337205e-04     0.0002423466
 [2,]     5.345844e-05      0.002233803     6.413270e-05     0.0002386332
 [3,]     5.466911e-05      0.003001030     1.215940e-04     0.0003203479
 [4,]     1.067677e-04      0.003135096     1.840003e-04     0.0004512870
 [5,]     6.241661e-05      0.002795479     1.261539e-04     0.0003346127
 [6,]     5.036556e-05      0.002674422     1.202560e-04     0.0002559675
 [7,]     9.265473e-06      0.002715929     9.214965e-05     0.0002493553
 [8,]     1.184045e-05      0.003053412     1.090292e-04     0.0003945425
 [9,]     5.507491e-05      0.002954519     4.827728e-05     0.0002915746
[10,]     5.423223e-05      0.002861514     6.145328e-05     0.0002112297
[11,]     1.964640e-04      0.002001516     2.360224e-05     0.0001391996
[12,]     2.161946e-04      0.001684884     6.547548e-05     0.0001353210
[13,]     2.579148e-05      0.002764311     6.522620e-05     0.0002427110
[14,]     4.891264e-08      0.002867245     2.150910e-04     0.0003278548
[15,]     1.622985e-04      0.001797896     2.467230e-05     0.0001955162
[16,]     1.746566e-04      0.001789026     3.727405e-05     0.0002037430
      4696.96495414691 4711.60113524674 4717.14836388862 4738.66026830518
 [1,]     1.134666e-04     7.722719e-05     5.305972e-05     1.292614e-06
 [2,]     9.661748e-05     9.867651e-05     1.019636e-04     1.040465e-04
 [3,]     6.485004e-05     6.790479e-05     6.664116e-05     1.857448e-05
 [4,]     1.248065e-04     1.042563e-04     1.932370e-04     7.300556e-05
 [5,]     5.007947e-05     8.843450e-05     1.951268e-04     4.635917e-05
 [6,]     1.193600e-04     2.249024e-04     1.837367e-04     4.575025e-05
 [7,]     8.055708e-05     1.069018e-04     1.649969e-04     1.272748e-04
 [8,]     1.335887e-04     8.866566e-05     1.470143e-04     2.893482e-05
 [9,]     2.309045e-05     7.629486e-05     1.650403e-04     9.513921e-05
[10,]     1.124551e-04     1.367141e-04     1.154319e-04     1.028136e-04
[11,]     1.137066e-05     1.858406e-04     1.343718e-04     4.757099e-06
[12,]     4.087844e-06     1.212634e-04     1.476578e-04     6.168010e-06
[13,]     1.413714e-04     2.403836e-04     1.876590e-04     4.148402e-05
[14,]     4.715460e-06     1.638243e-04     2.116437e-04     6.361136e-05
[15,]     6.131375e-05     1.225064e-04     1.830333e-04     2.666932e-05
[16,]     5.585543e-05     2.037284e-04     1.373318e-04     3.811995e-05
      4756.70119020648 4787.35841210002 4818.73389558634 4858.58348616161
 [1,]     2.065905e-04     0.0004265728     0.0001473230     1.336936e-04
 [2,]     1.384903e-04     0.0004324996     0.0001985432     1.022186e-04
 [3,]     1.715065e-04     0.0003755936     0.0001465495     1.079806e-04
 [4,]     1.598141e-04     0.0003215047     0.0004325609     3.641924e-05
 [5,]     1.771767e-04     0.0003823927     0.0004494675     2.252833e-05
 [6,]     1.910047e-04     0.0003008218     0.0002140275     5.279958e-05
 [7,]     1.390322e-04     0.0002557415     0.0002444113     5.412298e-05
 [8,]     1.178294e-04     0.0002908686     0.0002375233     6.210567e-05
 [9,]     9.990369e-05     0.0005997571     0.0001777974     9.046192e-05
[10,]     1.462710e-05     0.0005558499     0.0001759821     9.131179e-05
[11,]     1.530482e-04     0.0007527810     0.0001271095     1.270332e-04
[12,]     2.120331e-04     0.0007795038     0.0001442019     6.111219e-05
[13,]     2.725709e-04     0.0009243312     0.0002653166     7.236452e-06
[14,]     3.563547e-04     0.0009003499     0.0003645268     4.870185e-05
[15,]     2.844437e-04     0.0012901235     0.0003345602     1.125097e-04
[16,]     2.855070e-04     0.0013071277     0.0003342047     1.133711e-04
      4889.08531319992 4915.91856017224 4929.08511037096 4963.95577399317
 [1,]     4.730173e-05     1.228498e-04     9.711190e-05     0.0007464950
 [2,]     9.299445e-05     5.964462e-05     7.687649e-05     0.0008112504
 [3,]     1.233842e-04     1.188421e-04     1.198838e-04     0.0005821682
 [4,]     2.888570e-05     1.670499e-04     1.495746e-05     0.0006188056
 [5,]     1.140674e-04     5.919726e-05     2.237563e-05     0.0014726618
 [6,]     6.953635e-05     1.036329e-04     1.157589e-04     0.0013530412
 [7,]     5.110512e-05     7.230649e-05     2.155537e-05     0.0012710971
 [8,]     1.139270e-04     9.468978e-05     2.199854e-05     0.0014057738
 [9,]     9.171834e-05     8.422090e-05     1.434972e-05     0.0008160756
[10,]     1.223505e-04     1.030202e-04     5.517171e-05     0.0007204076
[11,]     2.007431e-04     1.746969e-04     6.686083e-05     0.0008027316
[12,]     1.795232e-04     1.549454e-04     3.709127e-05     0.0008187075
[13,]     5.549608e-05     1.177636e-04     1.151284e-04     0.0007284584
[14,]     1.094618e-04     4.578925e-05     7.408778e-05     0.0007216042
[15,]     1.564122e-04     1.485310e-04     1.604559e-05     0.0014141203
[16,]     1.391571e-04     1.628555e-04     4.523585e-05     0.0013805370
      5003.96096248945 5025.24385627375 5045.64515476984 5113.13827310569
 [1,]     0.0003160602     7.153395e-05     0.0002084518     3.888811e-05
 [2,]     0.0003382652     2.016356e-05     0.0001988173     1.445277e-05
 [3,]     0.0004139160     1.063553e-05     0.0002598668     4.666709e-05
 [4,]     0.0003531564     1.464887e-04     0.0002108283     2.838529e-05
 [5,]     0.0004342921     1.021718e-04     0.0002397619     3.893629e-05
 [6,]     0.0003350467     1.251788e-04     0.0001649994     2.607667e-05
 [7,]     0.0003651336     9.413900e-05     0.0001583693     1.057563e-04
 [8,]     0.0003396573     2.905094e-05     0.0001553697     5.485447e-05
 [9,]     0.0002117377     9.586800e-05     0.0001226091     1.739057e-05
[10,]     0.0001673359     5.406218e-05     0.0001202710     1.207982e-05
[11,]     0.0003495602     1.058840e-04     0.0001662464     2.896603e-05
[12,]     0.0003053438     5.777545e-05     0.0002509664     1.162910e-04
[13,]     0.0004708190     4.291380e-05     0.0003164240     3.657083e-05
[14,]     0.0005055644     4.107280e-05     0.0003756503     4.915318e-05
[15,]     0.0007480489     1.282625e-04     0.0003157454     1.212699e-04
[16,]     0.0006189057     1.088443e-04     0.0003472505     1.081799e-04
      5131.49700360209 5159.75511489671 5204.14945235999 5224.97647985704
 [1,]     0.0008902415     1.303576e-04     2.449517e-05     1.004579e-04
 [2,]     0.0009078894     7.215000e-05     5.154321e-05     2.529736e-05
 [3,]     0.0008199384     1.036763e-04     5.791249e-05     2.749484e-05
 [4,]     0.0008360045     2.623608e-05     6.324520e-05     5.219899e-05
 [5,]     0.0004220085     5.946637e-05     5.050723e-05     3.088832e-05
 [6,]     0.0003478683     1.285054e-04     1.074796e-04     1.002731e-04
 [7,]     0.0003712600     6.981937e-05     1.154554e-04     3.307403e-05
 [8,]     0.0004843937     2.474820e-05     2.338220e-05     1.069121e-04
 [9,]     0.0003552378     1.034292e-06     9.712629e-05     1.085629e-04
[10,]     0.0003872204     4.944523e-05     1.421452e-05     3.245429e-05
[11,]     0.0005580732     4.383828e-05     3.849311e-05     1.402465e-05
[12,]     0.0004829251     4.386437e-05     5.691485e-05     3.030819e-05
[13,]     0.0003877630     5.312174e-05     1.076383e-04     1.165500e-04
[14,]     0.0003840961     1.165976e-04     1.537879e-04     1.113487e-04
[15,]     0.0001416581     2.602767e-05     1.519257e-05     4.005618e-05
[16,]     0.0001477282     1.068927e-04     1.432862e-04     1.227795e-04
      5248.49109786014 5265.25380412582 5292.67680160368 5320.12492890514
 [1,]     0.0002864990     1.839745e-04     0.0001710560     0.0002516085
 [2,]     0.0002666611     1.270923e-04     0.0001877078     0.0002136896
 [3,]     0.0002094497     1.594493e-04     0.0002079273     0.0002529062
 [4,]     0.0002440496     1.840488e-04     0.0001966472     0.0002097772
 [5,]     0.0002228024     5.184934e-05     0.0001760094     0.0001674137
 [6,]     0.0002010659     1.091541e-04     0.0001699895     0.0001167281
 [7,]     0.0002467292     1.796473e-04     0.0001522196     0.0001514988
 [8,]     0.0002144708     1.486963e-04     0.0001381053     0.0002076439
 [9,]     0.0003208054     1.030196e-04     0.0002249978     0.0002764442
[10,]     0.0002781961     1.518114e-04     0.0002043579     0.0002512945
[11,]     0.0002702800     7.068135e-05     0.0001660049     0.0002095762
[12,]     0.0002845496     3.472338e-05     0.0001719048     0.0002794310
[13,]     0.0002010190     3.252291e-05     0.0002775708     0.0002343803
[14,]     0.0002296774     2.173065e-05     0.0002505213     0.0002234294
[15,]     0.0002905230     7.622332e-05     0.0002339665     0.0001855118
[16,]     0.0001988635     5.131116e-05     0.0002457504     0.0001978019
      5336.94475143711 5376.66957046922 5391.09177499478 5438.42327020182
 [1,]      0.003088365     7.412615e-05     0.0003059230     1.169348e-04
 [2,]      0.003065289     1.203826e-04     0.0003923409     6.207572e-05
 [3,]      0.003112527     1.742231e-04     0.0003262418     9.905747e-05
 [4,]      0.003091151     9.431327e-05     0.0003237742     9.632139e-05
 [5,]      0.002728466     1.178914e-04     0.0002543969     7.781234e-05
 [6,]      0.002620679     1.304766e-04     0.0002823979     5.215966e-05
 [7,]      0.002474524     6.041251e-05     0.0002078120     1.061193e-04
 [8,]      0.002591986     6.775103e-05     0.0002018662     9.321149e-05
 [9,]      0.003802145     9.664386e-05     0.0003652520     7.503079e-05
[10,]      0.003743696     1.025231e-04     0.0003149559     1.061013e-04
[11,]      0.003017052     1.014798e-04     0.0002736128     1.072001e-04
[12,]      0.002993227     6.568304e-05     0.0002531195     3.851955e-05
[13,]      0.003402453     1.326990e-04     0.0002532959     6.014682e-05
[14,]      0.003401297     1.358545e-04     0.0002835555     1.247122e-04
[15,]      0.002683888     5.289219e-05     0.0002160900     3.506794e-05
[16,]      0.002693494     1.180440e-04     0.0002883953     1.136037e-04
      5463.82459392421 5481.72111516643 5507.57173049439 5523.0695388462
 [1,]     1.141899e-04     0.0001424261     1.175940e-04    0.0002180331
 [2,]     7.458436e-05     0.0001129288     7.909682e-05    0.0001980973
 [3,]     5.209135e-05     0.0001800506     6.667297e-05    0.0002191504
 [4,]     1.049603e-04     0.0001458221     5.380398e-05    0.0002219989
 [5,]     4.874269e-05     0.0001437502     1.011012e-05    0.0002011310
 [6,]     4.695315e-05     0.0001299746     1.322313e-04    0.0001513589
 [7,]     6.579896e-05     0.0001243286     3.332696e-05    0.0001970752
 [8,]     3.819511e-05     0.0001235606     9.901175e-05    0.0001953691
 [9,]     6.619565e-05     0.0001637736     9.353465e-05    0.0003196639
[10,]     1.075040e-04     0.0001564527     1.056511e-04    0.0003145029
[11,]     4.738493e-05     0.0001202939     8.437161e-05    0.0001989539
[12,]     4.048843e-05     0.0001035968     1.129633e-05    0.0002721987
[13,]     6.875094e-05     0.0001608439     1.167298e-04    0.0002252092
[14,]     2.816343e-06     0.0001836105     5.959595e-05    0.0002600959
[15,]     7.329065e-05     0.0001455964     6.682414e-05    0.0002008362
[16,]     1.425727e-04     0.0001372106     4.510102e-05    0.0002545408
      5549.99183805551 5592.82362141132 5631.53828626015 5664.26693631609
 [1,]     6.072025e-05     1.034941e-05     3.653581e-05     1.374662e-04
 [2,]     3.221543e-05     1.565244e-05     4.508084e-05     9.823708e-05
 [3,]     5.349529e-05     1.569785e-04     1.963747e-04     1.263315e-04
 [4,]     7.694213e-05     2.750544e-05     1.902159e-04     3.962835e-05
 [5,]     1.078201e-04     3.756612e-06     9.917728e-05     1.100780e-04
 [6,]     2.913729e-05     4.866559e-06     4.226820e-06     6.551013e-05
 [7,]     1.398286e-06     1.417377e-04     1.354807e-04     1.471130e-06
 [8,]     5.469243e-05     1.023815e-04     1.633270e-04     3.439579e-05
 [9,]     7.653008e-05     4.427816e-05     1.624820e-05     6.420076e-05
[10,]     9.414495e-05     1.110837e-04     9.791351e-05     8.213728e-05
[11,]     4.805521e-05     4.783788e-05     1.235032e-06     2.792383e-05
[12,]     5.506959e-05     5.181036e-05     3.019065e-08     7.390743e-06
[13,]     1.148129e-04     6.523943e-05     8.636458e-05     1.074139e-04
[14,]     1.044478e-04     6.943598e-05     1.034062e-04     2.349505e-06
[15,]     4.437924e-06     1.377060e-04     1.109228e-04     8.978584e-05
[16,]     7.212177e-06     1.454254e-04     1.045077e-04     3.810004e-05
      5709.0677056475 5754.39113904892 5805.84744428392 5861.03605679657
 [1,]    1.055637e-04     2.460789e-04     0.0004492554     0.0004080457
 [2,]    4.968648e-05     1.924180e-04     0.0004371165     0.0004462756
 [3,]    6.024291e-05     4.798983e-04     0.0005093875     0.0003894447
 [4,]    9.612556e-05     4.313482e-04     0.0004706251     0.0004141037
 [5,]    4.110160e-05     1.732603e-04     0.0005879437     0.0003555838
 [6,]    3.473669e-05     1.941304e-04     0.0005725246     0.0003634321
 [7,]    4.696550e-05     3.097730e-04     0.0005303274     0.0002958095
 [8,]    3.069433e-05     4.428773e-04     0.0006390917     0.0003426473
 [9,]    6.816180e-05     5.706048e-05     0.0003848927     0.0003368813
[10,]    3.739020e-05     6.607506e-05     0.0004153139     0.0003222925
[11,]    9.971565e-05     8.549246e-05     0.0005682714     0.0002766608
[12,]    1.501981e-05     1.050588e-04     0.0005777440     0.0003264111
[13,]    3.050030e-05     5.430088e-05     0.0002156724     0.0001894542
[14,]    1.108707e-04     1.281992e-04     0.0002118438     0.0002117597
[15,]    6.058445e-05     1.413767e-04     0.0002724563     0.0002064774
[16,]    5.477728e-07     1.432513e-04     0.0002790576     0.0002967624
      5888.01996610497 5904.8416229427 5944.59390105404 5958.67584536279
 [1,]     0.0004547650     0.005961930     1.981813e-04     0.0009281123
 [2,]     0.0004524012     0.005698564     2.434852e-04     0.0010641328
 [3,]     0.0004668910     0.005346926     2.534688e-04     0.0008177657
 [4,]     0.0004948327     0.005236728     2.434380e-04     0.0008861159
 [5,]     0.0004282888     0.005135802     1.930042e-04     0.0008221709
 [6,]     0.0004175543     0.005071343     2.673975e-04     0.0008315228
 [7,]     0.0004088752     0.004914612     1.425687e-04     0.0007042711
 [8,]     0.0003949093     0.004944022     2.527877e-04     0.0007693195
 [9,]     0.0004564545     0.005273369     1.329654e-04     0.0006947347
[10,]     0.0004649319     0.005268427     1.255526e-04     0.0006416639
[11,]     0.0004335448     0.005017247     1.575745e-04     0.0006810367
[12,]     0.0004617116     0.005170046     1.156504e-04     0.0006429693
[13,]     0.0002767731     0.002788945     8.131470e-05     0.0003082577
[14,]     0.0001943959     0.002782726     1.330554e-04     0.0003146466
[15,]     0.0003612513     0.004041620     1.109738e-04     0.0004398736
[16,]     0.0003441454     0.003983635     9.058133e-05     0.0006161386
      6004.03860366915 6049.85815913964 6075.11159149524 6090.33276680834
 [1,]     1.153510e-04     0.0003785172     2.238513e-04     0.0006593598
 [2,]     9.665426e-05     0.0003173530     1.572495e-04     0.0006864656
 [3,]     8.036462e-05     0.0002869883     1.583097e-04     0.0006448303
 [4,]     7.854874e-05     0.0002779845     1.472055e-04     0.0005964953
 [5,]     1.045598e-04     0.0002393113     1.476879e-04     0.0005510108
 [6,]     8.320169e-05     0.0003155621     8.005096e-05     0.0005181508
 [7,]     1.762998e-05     0.0002141392     1.101920e-04     0.0004835484
 [8,]     7.397199e-05     0.0003244369     1.175081e-04     0.0005032730
 [9,]     9.086050e-05     0.0002587501     1.087210e-04     0.0005379042
[10,]     9.512121e-05     0.0002826322     9.656866e-05     0.0005110448
[11,]     4.103469e-05     0.0002311208     4.834890e-05     0.0004682753
[12,]     8.491378e-05     0.0002325906     9.389735e-05     0.0004581811
[13,]     1.089984e-04     0.0001515710     7.657924e-05     0.0003785363
[14,]     3.366854e-05     0.0001341181     8.575283e-05     0.0003542324
[15,]     4.390011e-05     0.0001913691     8.043958e-05     0.0004778903
[16,]     1.232767e-04     0.0002180969     1.228129e-04     0.0005093181
      6142.79597906585 6175.07811529631 6191.47316724826 6226.72066206662
 [1,]     7.020125e-05     0.0002069752     9.801551e-05     8.135877e-05
 [2,]     7.415967e-05     0.0002436456     1.007914e-04     7.935729e-05
 [3,]     4.840170e-05     0.0001976049     7.245768e-05     2.273270e-05
 [4,]     1.459019e-04     0.0002010387     1.251557e-04     5.860394e-05
 [5,]     5.320548e-05     0.0001924022     1.349882e-04     9.527173e-05
 [6,]     1.833119e-05     0.0002609434     1.648488e-04     2.260335e-04
 [7,]     2.386622e-05     0.0002411817     2.061023e-04     4.398201e-05
 [8,]     4.714987e-05     0.0003214705     2.026059e-04     5.655032e-05
 [9,]     5.346320e-05     0.0001131746     7.393088e-05     2.630813e-04
[10,]     3.039991e-05     0.0001627799     1.159018e-04     2.826383e-04
[11,]     1.048098e-04     0.0001426536     6.423709e-05     2.242881e-04
[12,]     4.946802e-05     0.0002130328     4.082552e-05     1.651957e-04
[13,]     5.769845e-05     0.0002016778     1.618216e-04     3.756717e-04
[14,]     4.287050e-05     0.0002160538     7.744993e-05     3.606440e-04
[15,]     1.083584e-04     0.0002856136     1.289735e-04     1.097467e-04
[16,]     1.260865e-04     0.0002302208     1.086466e-04     1.113976e-04
      6271.80687063463 6305.28147764681 6332.67940029196 6357.33152142188
 [1,]     1.128043e-04     0.0002251132     1.528091e-04     5.875258e-05
 [2,]     1.291672e-04     0.0002006182     9.504794e-05     9.786330e-06
 [3,]     1.050896e-04     0.0002048426     1.568513e-04     1.206276e-04
 [4,]     1.107552e-04     0.0001866841     1.443333e-04     1.060619e-04
 [5,]     2.488616e-05     0.0002656893     1.830946e-04     7.151027e-09
 [6,]     9.223732e-05     0.0002761589     1.263732e-04     2.480843e-05
 [7,]     4.592758e-05     0.0001041438     1.035577e-04     5.442823e-05
 [8,]     3.309860e-05     0.0002692963     1.495533e-04     9.840658e-05
 [9,]     5.431232e-05     0.0003011839     1.189237e-04     1.511010e-05
[10,]     3.160365e-05     0.0003983494     1.489200e-04     1.540542e-05
[11,]     1.045750e-04     0.0002134747     1.269868e-04     1.679852e-05
[12,]     1.046321e-04     0.0001707981     1.100373e-04     1.387951e-05
[13,]     2.382335e-05     0.0003120392     2.356604e-04     3.071771e-05
[14,]     7.339941e-06     0.0003120533     2.159969e-04     4.565393e-05
[15,]     5.429242e-05     0.0001503627     1.871124e-04     1.051966e-04
[16,]     7.123923e-05     0.0001089884     1.897225e-04     5.615154e-05
      6378.45261795257 6389.36063457517 6433.13858751605 6486.53350672773
 [1,]     1.301784e-04     5.027073e-05     0.0002875094     1.096733e-04
 [2,]     1.625996e-04     9.013871e-05     0.0005058485     6.034037e-05
 [3,]     1.972537e-04     7.683909e-05     0.0003389775     5.738553e-05
 [4,]     2.732307e-04     7.027988e-05     0.0006052795     1.039860e-04
 [5,]     2.417591e-04     9.771302e-05     0.0008792028     1.092225e-04
 [6,]     2.044766e-04     5.858657e-05     0.0004987679     1.272353e-05
 [7,]     1.729251e-04     3.213528e-05     0.0005274600     8.799362e-06
 [8,]     1.394125e-04     5.077735e-05     0.0005551410     1.080415e-04
 [9,]     9.240146e-05     1.077915e-04     0.0006410953     9.465412e-05
[10,]     1.281367e-04     9.017524e-05     0.0005646639     6.025551e-05
[11,]     1.545984e-04     6.838163e-05     0.0003842195     1.044600e-04
[12,]     1.773386e-04     7.899794e-05     0.0004480768     5.718496e-05
[13,]     3.263548e-04     1.264670e-04     0.0009481186     9.334212e-05
[14,]     2.637229e-04     5.545414e-05     0.0009069153     1.060290e-04
[15,]     2.955306e-04     9.217470e-05     0.0004624737     3.488932e-05
[16,]     2.491993e-04     2.474105e-04     0.0007182943     6.171207e-05
      6492.81796708918 6518.51214477903 6561.92037763536 6590.74358448274
 [1,]     4.134182e-05     1.079134e-04     2.092105e-04     2.816122e-05
 [2,]     2.357097e-05     1.162127e-04     1.677960e-04     1.360242e-04
 [3,]     6.050035e-05     8.835301e-05     1.131995e-04     1.082724e-04
 [4,]     1.856077e-05     8.794843e-06     1.627599e-04     3.759542e-05
 [5,]     2.552999e-05     4.483058e-05     8.151089e-05     1.853908e-06
 [6,]     1.026655e-04     1.208267e-04     1.010256e-04     8.789257e-06
 [7,]     1.079610e-04     4.756280e-05     5.470451e-05     2.382437e-05
 [8,]     5.606862e-05     7.093019e-05     5.999767e-05     3.544726e-05
 [9,]     1.038217e-04     2.900487e-05     4.318043e-05     5.094130e-05
[10,]     9.523246e-05     1.032006e-04     1.187446e-04     9.965902e-05
[11,]     5.998805e-05     4.894216e-05     1.358290e-04     3.688652e-06
[12,]     2.609810e-05     9.133048e-05     8.515754e-05     3.016925e-05
[13,]     1.633679e-06     1.202000e-04     1.434181e-06     2.811470e-05
[14,]     4.194263e-05     1.455352e-05     1.085482e-04     3.088651e-05
[15,]     4.138816e-05     1.374183e-04     5.553831e-05     1.413084e-05
[16,]     5.796295e-05     1.026419e-04     1.422984e-04     1.336751e-04
      6614.65757361826 6631.48716897254 6669.4795205433 6685.71154030979
 [1,]     2.940697e-05     0.0004498604    1.610271e-04     2.566649e-05
 [2,]     1.123380e-04     0.0007652165    1.414527e-04     7.581960e-05
 [3,]     1.334930e-04     0.0007436718    2.717369e-04     1.344982e-05
 [4,]     5.817424e-05     0.0010769472    2.768227e-04     1.312845e-04
 [5,]     8.976433e-05     0.0013115991    2.428615e-04     1.141849e-04
 [6,]     1.057280e-04     0.0010604741    1.811256e-04     7.000242e-05
 [7,]     9.145385e-05     0.0008809762    1.181128e-04     2.355760e-05
 [8,]     9.602798e-05     0.0008711227    1.303649e-04     7.451948e-05
 [9,]     1.150870e-04     0.0012156924    8.010051e-05     8.630502e-05
[10,]     9.629682e-05     0.0012655953    9.726752e-05     3.212519e-05
[11,]     6.549798e-05     0.0008584471    1.539036e-04     5.772284e-05
[12,]     8.434251e-05     0.0009081684    1.658972e-04     2.685761e-05
[13,]     8.324379e-05     0.0016140486    3.456865e-04     1.078310e-04
[14,]     1.047018e-04     0.0015761002    3.282078e-04     1.368786e-04
[15,]     6.176425e-05     0.0007870085    3.274833e-04     4.805059e-05
[16,]     1.321066e-04     0.0010626166    3.310250e-04     1.062870e-04
      6749.08869581467 6805.5111901082 6831.8438152041 6852.73944871908
 [1,]     1.079157e-04    1.016284e-04    6.601326e-05     1.312511e-04
 [2,]     1.185029e-04    6.693779e-05    1.036583e-04     5.127593e-05
 [3,]     7.021816e-05    1.100846e-04    1.161709e-05     2.592794e-05
 [4,]     1.009103e-04    5.610479e-05    4.721990e-05     8.091887e-05
 [5,]     3.283028e-05    6.525752e-05    2.494973e-05     3.251717e-05
 [6,]     7.599736e-05    1.218206e-04    4.008118e-05     8.954416e-07
 [7,]     3.254328e-05    1.501050e-04    9.847702e-05     7.468508e-05
 [8,]     1.539596e-05    1.080285e-04    2.529631e-05     1.266667e-05
 [9,]     6.764899e-06    7.122361e-05    2.987660e-05     1.705231e-05
[10,]     1.438152e-05    6.600934e-05    9.066433e-05     8.947785e-05
[11,]     2.132605e-05    8.823520e-05    4.072842e-05     1.218349e-05
[12,]     2.511285e-05    3.701180e-05    2.271225e-05     6.203578e-05
[13,]     2.086894e-05    5.506680e-05    1.130062e-04     1.326931e-04
[14,]     1.260630e-04    1.108075e-04    7.283107e-05     1.242931e-04
[15,]     7.149981e-05    1.942024e-04    3.990972e-06     8.465653e-05
[16,]     1.091570e-04    1.318107e-04    3.805372e-05     1.685827e-04
      6909.41489438275 6961.4734152769 7008.65679889363 7068.14341465598
 [1,]     6.629931e-05    6.241293e-05     1.637883e-04     1.208702e-05
 [2,]     9.939505e-05    4.048453e-05     1.135172e-04     6.841276e-05
 [3,]     1.182052e-04    1.038518e-04     1.707857e-04     7.992120e-05
 [4,]     3.672116e-05    4.395376e-05     1.701485e-04     9.587097e-05
 [5,]     7.744652e-05    1.033517e-04     1.405041e-04     1.093285e-04
 [6,]     5.897693e-05    1.126095e-04     1.771875e-04     6.306313e-05
 [7,]     2.214074e-05    1.275872e-04     2.100562e-04     2.933138e-05
 [8,]     6.554725e-05    7.454758e-05     1.942245e-04     4.616960e-05
 [9,]     1.233229e-04    6.140232e-05     2.186439e-04     5.067098e-06
[10,]     7.734256e-05    3.983239e-05     2.319773e-04     3.363262e-05
[11,]     2.126035e-05    3.741937e-05     1.116694e-04     2.010432e-05
[12,]     4.188521e-06    6.453737e-05     1.148219e-04     8.990686e-05
[13,]     4.166699e-05    7.545316e-05     1.134935e-04     1.434165e-04
[14,]     7.669534e-05    8.693573e-05     1.707377e-04     3.426320e-05
[15,]     1.377714e-04    1.356772e-04     1.257391e-04     4.065765e-05
[16,]     6.294362e-05    5.563039e-05     2.470849e-05     5.082218e-05
      7131.32716337575 7163.24208666278 7191.42574725989 7222.08538000553
 [1,]     1.084978e-05     4.392270e-05     1.503617e-04     2.902761e-05
 [2,]     1.199057e-05     1.045038e-04     1.108743e-04     7.088013e-05
 [3,]     5.969969e-05     1.269347e-04     4.922389e-05     1.181955e-04
 [4,]     5.326225e-05     1.531344e-05     1.148864e-04     1.183524e-04
 [5,]     9.931855e-05     5.546058e-07     3.016322e-05     9.659441e-05
 [6,]     1.524581e-04     1.251345e-04     1.326838e-04     4.545740e-05
 [7,]     1.650450e-05     4.816040e-05     1.406574e-04     1.389983e-04
 [8,]     2.602647e-06     2.876670e-05     1.218063e-04     1.125665e-04
 [9,]     1.317781e-05     9.004612e-06     8.415447e-05     1.010809e-04
[10,]     1.244733e-04     2.753979e-06     1.455903e-04     1.479262e-04
[11,]     2.233008e-05     4.541203e-06     4.589024e-05     1.993664e-05
[12,]     3.392529e-05     1.436143e-05     6.387459e-05     4.269144e-05
[13,]     1.085533e-04     1.333312e-04     7.728702e-05     1.441576e-04
[14,]     2.474203e-05     1.427176e-05     1.105084e-04     3.296566e-05
[15,]     5.068791e-05     4.538850e-05     1.198806e-04     1.678114e-05
[16,]     3.973852e-05     2.010423e-05     7.478769e-05     1.275841e-05
      7291.2700162216 7307.41203973464 7343.96223714406 7399.45219041599
 [1,]    4.609579e-05     7.651926e-06     9.185712e-06     7.888216e-05
 [2,]    4.040541e-05     9.644380e-05     9.962778e-05     4.607467e-05
 [3,]    6.112981e-05     5.695981e-05     2.456599e-05     1.008821e-04
 [4,]    8.847325e-06     1.728233e-05     3.286364e-05     1.214840e-04
 [5,]    2.161887e-05     1.926670e-05     7.936128e-06     2.276290e-05
 [6,]    6.323335e-05     7.874446e-05     1.076436e-04     1.115122e-04
 [7,]    1.075889e-04     1.237592e-04     1.034669e-04     1.024789e-04
 [8,]    9.641824e-05     5.466362e-05     1.441528e-05     1.062981e-04
 [9,]    2.926623e-05     4.335950e-05     3.461585e-05     2.904699e-05
[10,]    4.548727e-05     1.034732e-04     7.236914e-07     4.495927e-05
[11,]    2.653983e-05     1.291660e-04     1.835333e-05     2.967808e-05
[12,]    6.354471e-05     9.344829e-05     2.580892e-05     2.011028e-05
[13,]    1.311628e-04     2.759542e-05     5.513123e-05     4.659317e-05
[14,]    1.175867e-04     1.112440e-05     4.501580e-05     3.843095e-05
[15,]    1.301433e-04     5.476191e-05     1.078994e-04     2.567753e-05
[16,]    5.388976e-05     6.261572e-05     1.901182e-05     9.403609e-05
      7432.80326677328 7471.77342742474 7550.61399083004 7565.97150499574
 [1,]     9.990255e-05     0.0001772832     1.751484e-05     0.0003009970
 [2,]     3.004521e-05     0.0001304203     4.608869e-05     0.0002647176
 [3,]     1.104437e-04     0.0001211745     1.049788e-04     0.0002790347
 [4,]     8.611638e-05     0.0001762308     1.082145e-04     0.0002560815
 [5,]     5.414968e-05     0.0001069742     6.713514e-05     0.0002135309
 [6,]     8.797558e-06     0.0001615774     6.421848e-05     0.0001778226
 [7,]     1.766797e-04     0.0001656311     1.163916e-04     0.0003250583
 [8,]     1.113850e-04     0.0001343157     5.789451e-05     0.0002647018
 [9,]     1.078039e-04     0.0001397171     5.575685e-05     0.0002883723
[10,]     4.293981e-05     0.0001217358     5.804654e-05     0.0002801147
[11,]     2.653015e-05     0.0001921106     6.593081e-05     0.0001746955
[12,]     3.222120e-05     0.0001587055     6.068329e-05     0.0002011087
[13,]     1.087923e-04     0.0002023428     1.153664e-04     0.0001937267
[14,]     5.251365e-05     0.0001471685     8.084032e-05     0.0001821217
[15,]     1.595134e-05     0.0001550476     7.930254e-05     0.0001882411
[16,]     5.769536e-05     0.0001780699     2.920069e-05     0.0001397298
      7636.79659628074 7721.88431895238 7766.71168655629 7806.77244003716
 [1,]     1.643237e-04     1.287151e-04      0.002968064     0.0001106211
 [2,]     1.181747e-04     1.029749e-04      0.002882008     0.0001526789
 [3,]     1.757314e-04     1.688218e-04      0.003073595     0.0002516052
 [4,]     1.307811e-04     2.259811e-04      0.003086729     0.0002104736
 [5,]     2.093552e-05     7.663029e-05      0.002621178     0.0001414057
 [6,]     3.442294e-05     1.456289e-04      0.002659649     0.0001749863
 [7,]     1.265345e-04     8.939988e-05      0.003091128     0.0001344565
 [8,]     1.420287e-04     1.217931e-04      0.003017838     0.0002007781
 [9,]     6.213151e-05     1.246100e-04      0.002813987     0.0001642806
[10,]     9.091061e-05     1.316255e-04      0.002866035     0.0001205672
[11,]     4.387280e-05     1.336578e-04      0.002415155     0.0001081745
[12,]     7.919068e-05     9.583185e-05      0.002454852     0.0001521287
[13,]     6.714663e-05     1.240525e-04      0.002934838     0.0002673567
[14,]     6.685217e-05     2.067548e-04      0.002848654     0.0003461664
[15,]     1.363078e-04     6.148549e-05      0.002324636     0.0001000702
[16,]     6.366999e-05     1.264455e-04      0.002345571     0.0001146530
      7820.65441315211 7836.69091874102 7875.9707895942 7909.14811077727
 [1,]     0.0002597252     1.295522e-04    3.267987e-05     1.690771e-04
 [2,]     0.0003844805     1.342340e-05    8.050937e-05     1.722287e-04
 [3,]     0.0002523752     1.730217e-04    1.475380e-04     8.561630e-05
 [4,]     0.0004074782     2.791810e-04    9.669218e-05     1.544113e-04
 [5,]     0.0002627476     6.086698e-05    9.596121e-05     5.155144e-05
 [6,]     0.0002652001     1.214140e-05    7.551828e-05     5.951604e-05
 [7,]     0.0003552379     1.761099e-04    2.568717e-05     1.074149e-04
 [8,]     0.0003471571     1.531449e-04    3.895225e-05     1.620029e-04
 [9,]     0.0002980332     1.099802e-04    6.360401e-05     1.096458e-04
[10,]     0.0002015069     1.217157e-04    1.022684e-05     1.286970e-04
[11,]     0.0002267941     1.117095e-04    1.251255e-07     1.132053e-04
[12,]     0.0002379248     6.658337e-05    1.463022e-06     1.205715e-04
[13,]     0.0002674417     1.460158e-05    1.069393e-04     9.238864e-05
[14,]     0.0002433285     4.271296e-05    4.582096e-05     1.351628e-04
[15,]     0.0002085901     6.562039e-05    8.676326e-05     1.505656e-04
[16,]     0.0002867898     1.262529e-04    6.909160e-05     4.306220e-05
      7924.8597573429 8038.51347193325 8092.12356717488 8129.58689059504
 [1,]    0.0004185041     4.069078e-05     1.157076e-04     0.0002499391
 [2,]    0.0004409984     1.703559e-05     2.249677e-05     0.0002320492
 [3,]    0.0003758891     1.185696e-04     1.397059e-04     0.0002109783
 [4,]    0.0003429366     1.246885e-04     6.536798e-06     0.0001979859
 [5,]    0.0003798360     5.184517e-05     3.927754e-06     0.0002107596
 [6,]    0.0003182167     2.663849e-05     3.634790e-05     0.0002371271
 [7,]    0.0003805914     5.876265e-05     1.392002e-04     0.0002041612
 [8,]    0.0003939195     4.640422e-05     1.283389e-05     0.0002485390
 [9,]    0.0003715028     1.606648e-05     2.312496e-05     0.0002848179
[10,]    0.0003476720     9.058281e-05     6.939779e-05     0.0003158082
[11,]    0.0002647681     3.382740e-05     8.530922e-06     0.0003310458
[12,]    0.0002932638     1.664930e-05     9.879380e-05     0.0003336474
[13,]    0.0003950839     1.151150e-04     1.306107e-05     0.0010441444
[14,]    0.0003900028     6.338580e-05     1.555888e-06     0.0010256885
[15,]    0.0003262674     8.968171e-05     1.171714e-04     0.0006674996
[16,]    0.0003006409     1.331753e-04     1.700616e-05     0.0007034744
      8143.62691281769 8194.45551490086 8217.086521642 8262.10426631392
 [1,]     0.0008607159     1.323838e-04   1.349245e-05     2.414642e-05
 [2,]     0.0008641804     1.250418e-04   4.866079e-05     2.907355e-05
 [3,]     0.0009191043     1.297581e-04   1.063763e-04     7.424332e-05
 [4,]     0.0009470598     1.067116e-04   4.287885e-05     9.987789e-05
 [5,]     0.0007924824     1.836492e-05   6.408673e-05     1.612654e-05
 [6,]     0.0006626809     6.422021e-05   1.545045e-04     3.125874e-05
 [7,]     0.0009058698     1.986640e-04   5.931171e-05     6.181019e-05
 [8,]     0.0008968294     1.148749e-04   7.925167e-05     1.401205e-05
 [9,]     0.0007875340     9.355185e-05   3.870683e-05     1.009886e-04
[10,]     0.0007227464     4.403453e-05   9.045822e-05     4.939221e-05
[11,]     0.0006056406     3.181333e-05   1.239586e-04     2.729462e-05
[12,]     0.0006101023     1.000193e-04   1.283305e-04     6.351354e-05
[13,]     0.0005620446     1.220120e-04   1.154177e-04     3.544594e-05
[14,]     0.0005575915     1.146108e-04   1.177903e-04     1.349444e-04
[15,]     0.0005164784     2.150600e-05   1.162818e-05     1.079866e-04
[16,]     0.0005020924     1.454068e-04   2.069485e-05     1.087812e-04
      8301.56510152212 8364.22970743122 8429.66401181617 8452.85272629876
 [1,]     2.483909e-05     1.230419e-04     1.379903e-04     1.208035e-04
 [2,]     9.337451e-05     9.863060e-05     9.209457e-05     1.676654e-05
 [3,]     1.178081e-04     5.203644e-05     3.316052e-05     1.120195e-04
 [4,]     4.861959e-05     5.139208e-05     3.333383e-05     1.364822e-05
 [5,]     2.433694e-06     7.435424e-05     4.511768e-05     6.109158e-05
 [6,]     1.083122e-04     1.086488e-04     1.126001e-04     3.910381e-05
 [7,]     3.530870e-05     2.969570e-05     6.188665e-05     5.577509e-05
 [8,]     1.055143e-05     4.660576e-05     1.804524e-05     2.486845e-05
 [9,]     2.630510e-05     9.439161e-05     9.611358e-05     9.274145e-05
[10,]     1.244388e-05     1.525178e-04     4.290026e-05     3.984753e-05
[11,]     0.000000e+00     4.826470e-05     9.381053e-05     2.466070e-05
[12,]     2.672827e-05     3.129349e-06     1.696664e-05     3.945723e-05
[13,]     1.503475e-04     5.745236e-05     1.074065e-04     4.127970e-05
[14,]     1.062609e-04     2.223510e-05     3.733655e-05     1.726768e-05
[15,]     1.102107e-04     1.573953e-05     7.431757e-05     6.571752e-05
[16,]     8.587509e-05     1.231620e-04     5.887002e-05     1.058435e-04
      8500.58327853541 8534.32599080284 8566.37795359577 8604.21921683076
 [1,]     3.786741e-05     1.612851e-04     6.313794e-05     0.0002136551
 [2,]     5.572230e-06     1.772115e-04     1.335274e-05     0.0002592179
 [3,]     1.016294e-04     1.225093e-04     6.051870e-05     0.0001761719
 [4,]     9.434943e-05     1.124738e-04     7.852631e-05     0.0001971061
 [5,]     2.774850e-05     1.483205e-05     1.246020e-04     0.0002379248
 [6,]     1.071139e-05     4.523342e-05     7.760886e-05     0.0001643412
 [7,]     2.107922e-05     6.546473e-05     1.164215e-04     0.0001136347
 [8,]     2.842873e-05     5.046010e-05     1.570622e-04     0.0001627434
 [9,]     1.547122e-05     1.006203e-04     1.059783e-04     0.0002723230
[10,]     2.487094e-06     1.378100e-05     2.758014e-05     0.0002745756
[11,]     1.028475e-04     1.117535e-05     4.752411e-05     0.0007612668
[12,]     2.115813e-06     5.993451e-05     2.202301e-05     0.0006333327
[13,]     6.258645e-05     1.619435e-05     7.354832e-05     0.0006705469
[14,]     5.660246e-05     1.132055e-05     8.271655e-05     0.0006160834
[15,]     1.181576e-04     1.131051e-05     6.929155e-05     0.0007411999
[16,]     1.270004e-04     1.085217e-04     3.680096e-05     0.0007363214
      8657.56901717083 8663.53353909783 8677.05207560643 8737.82217095891
 [1,]     1.176087e-04     6.794434e-05     2.045784e-04     5.733209e-05
 [2,]     7.517349e-05     6.695462e-05     1.797768e-04     7.736686e-05
 [3,]     2.728258e-05     1.116808e-04     1.087675e-04     1.001629e-04
 [4,]     7.303911e-06     5.827326e-06     6.559095e-05     6.701545e-06
 [5,]     5.220064e-05     8.652095e-05     9.457578e-05     3.996625e-05
 [6,]     7.668091e-05     3.112761e-05     1.660761e-05     9.225594e-06
 [7,]     9.301811e-05     1.170494e-04     2.035856e-04     4.015792e-05
 [8,]     7.030515e-05     4.565005e-05     2.077512e-04     1.249788e-04
 [9,]     1.128732e-04     3.606667e-06     1.671677e-04     2.119269e-05
[10,]     3.060497e-05     4.550543e-05     1.938467e-04     4.814652e-05
[11,]     1.511416e-04     4.785447e-05     3.193479e-05     2.533779e-05
[12,]     1.029394e-04     1.072566e-04     5.029100e-05     2.090522e-05
[13,]     1.161288e-04     5.134152e-05     1.649780e-04     1.297642e-04
[14,]     2.412191e-05     5.711057e-05     2.212845e-04     1.619469e-05
[15,]     6.959899e-05     1.099390e-04     7.880490e-05     1.250659e-04
[16,]     1.323218e-04     9.982136e-05     1.196239e-04     2.954459e-05
      8757.21688745984 8767.25578963908 8810.00218944477 8864.01135962225
 [1,]     5.060848e-05     5.302185e-05     1.558217e-05     2.269900e-04
 [2,]     6.205004e-06     5.179442e-05     9.220775e-05     2.284466e-04
 [3,]     1.486301e-05     4.665740e-05     4.154166e-05     1.112106e-04
 [4,]     5.579412e-05     1.167011e-04     4.251063e-06     1.475148e-04
 [5,]     2.693291e-05     1.386366e-04     9.564101e-05     1.891748e-04
 [6,]     1.284572e-04     4.808189e-05     1.049798e-04     1.681095e-04
 [7,]     3.217572e-05     2.733263e-05     7.868452e-05     2.659381e-04
 [8,]     1.200830e-04     1.370487e-04     6.225700e-05     2.644205e-04
 [9,]     9.244492e-05     8.297438e-05     6.264980e-05     2.395601e-04
[10,]     6.685708e-05     1.190246e-04     5.290747e-05     2.546597e-04
[11,]     3.005688e-05     1.157718e-04     2.701677e-05     4.865177e-05
[12,]     2.500448e-05     4.487476e-05     1.365002e-05     9.390601e-05
[13,]     1.104853e-04     7.746345e-05     7.666811e-05     2.382909e-04
[14,]     3.346089e-05     3.173301e-05     6.253918e-05     2.311732e-04
[15,]     1.200015e-04     3.177759e-05     1.071868e-04     4.886903e-05
[16,]     7.353126e-05     3.000126e-05     1.317125e-04     6.309932e-05
      8935.25253689238 8989.09691156978 9062.59283961026 9131.45920094808
 [1,]     0.0003761294     1.190625e-04     0.0002677583     1.136724e-04
 [2,]     0.0003381886     9.449349e-05     0.0002425881     1.114333e-04
 [3,]     0.0003639052     1.023860e-04     0.0001416579     6.421829e-05
 [4,]     0.0004036331     1.347503e-04     0.0001977291     1.061837e-04
 [5,]     0.0003501925     3.057944e-05     0.0001639960     1.979559e-04
 [6,]     0.0003546520     6.257605e-05     0.0001236520     1.176646e-04
 [7,]     0.0004589074     7.613468e-05     0.0004068735     8.012923e-05
 [8,]     0.0004713617     1.122904e-04     0.0003378264     9.902114e-05
 [9,]     0.0004199782     6.025571e-05     0.0003334265     7.702748e-05
[10,]     0.0003957223     3.597648e-05     0.0003574385     9.906295e-05
[11,]     0.0006934210     1.017105e-04     0.0001618980     1.926058e-04
[12,]     0.0005972795     1.300323e-04     0.0001581274     2.309636e-04
[13,]     0.0015857888     2.942831e-04     0.0002347753     2.338028e-04
[14,]     0.0014289032     2.658309e-04     0.0002301146     2.118090e-04
[15,]     0.0011782602     2.462047e-04     0.0001336956     1.710283e-04
[16,]     0.0013211702     3.128011e-04     0.0001037638     2.094104e-04
      9136.94358480636 9176.13651228705 9290.92755913792 9345.09979718134
 [1,]     7.820678e-05     5.382989e-05      0.001912067     0.0003198169
 [2,]     1.146135e-04     1.041424e-04      0.001902960     0.0004079370
 [3,]     1.215120e-04     1.311041e-04      0.002409656     0.0004121462
 [4,]     9.867801e-05     1.601087e-04      0.002424606     0.0004784966
 [5,]     1.192410e-04     9.390285e-05      0.002232488     0.0004030589
 [6,]     1.658873e-04     1.077221e-04      0.002260217     0.0004363428
 [7,]     1.595055e-05     1.259783e-04      0.002470011     0.0004347837
 [8,]     5.810650e-05     7.735012e-05      0.002327020     0.0004144292
 [9,]     1.035681e-04     1.499039e-04      0.002407357     0.0004056656
[10,]     1.372139e-04     1.274449e-04      0.002420663     0.0003150365
[11,]     1.474431e-04     6.222727e-05      0.001865506     0.0002810194
[12,]     1.725843e-04     4.418732e-05      0.001876134     0.0002775032
[13,]     2.305488e-04     1.119799e-04      0.002313272     0.0003409817
[14,]     2.274177e-04     1.639701e-04      0.002249897     0.0002885116
[15,]     2.312486e-04     6.740981e-05      0.001612229     0.0002496030
[16,]     2.440957e-04     3.826013e-05      0.001534821     0.0003151860
      9425.91517823556 9438.43733716691 9478.95784114679 9525.01614140679
 [1,]     0.0001075953     1.242442e-04     1.070771e-04     1.097545e-04
 [2,]     0.0001868120     1.326399e-04     2.457778e-05     5.263377e-05
 [3,]     0.0001996740     1.530838e-04     1.336458e-04     1.002961e-04
 [4,]     0.0001011570     1.773814e-04     7.942590e-05     2.952309e-05
 [5,]     0.0001982479     9.864598e-05     6.464284e-05     5.838870e-05
 [6,]     0.0002143551     6.528834e-05     5.178883e-06     1.087077e-05
 [7,]     0.0001186221     1.180920e-04     5.974827e-05     4.176670e-05
 [8,]     0.0001291386     1.623924e-04     4.351741e-05     9.599229e-05
 [9,]     0.0001670179     1.417093e-04     1.397082e-04     7.614912e-06
[10,]     0.0001885729     1.009507e-04     8.596630e-06     8.061077e-05
[11,]     0.0002659734     9.996632e-05     1.060916e-04     6.674046e-05
[12,]     0.0002309277     1.159410e-04     1.061692e-04     4.728000e-05
[13,]     0.0003286658     1.084962e-04     1.218729e-04     3.365245e-05
[14,]     0.0001931007     1.573012e-04     9.230907e-07     4.973866e-05
[15,]     0.0002287929     5.940442e-05     5.147401e-05     1.070110e-04
[16,]     0.0002910252     1.661708e-04     1.100604e-04     1.073678e-04
      9590.43015419575 9629.07937520637 9715.85186727084 9724.35761206149
 [1,]     1.023977e-04     1.669892e-05     1.062896e-04     5.050348e-05
 [2,]     2.870870e-05     2.415084e-05     1.012562e-04     7.786326e-05
 [3,]     6.945030e-05     9.438526e-07     1.521285e-04     2.896603e-05
 [4,]     1.009321e-04     3.930046e-05     5.373022e-05     6.573548e-05
 [5,]     1.598246e-05     2.561729e-05     9.536226e-05     3.871777e-05
 [6,]     2.083652e-05     7.967534e-05     1.084318e-04     1.059579e-04
 [7,]     1.166538e-04     1.094850e-04     9.987854e-05     1.372494e-04
 [8,]     1.099063e-04     1.038141e-04     4.826107e-05     5.146403e-05
 [9,]     1.011720e-05     1.006503e-04     7.867312e-05     5.290968e-05
[10,]     3.809215e-05     1.091992e-04     7.941984e-05     9.066691e-05
[11,]     3.046586e-05     4.878608e-05     1.078194e-04     9.018278e-08
[12,]     1.153790e-05     9.398728e-05     1.014597e-04     7.793442e-05
[13,]     1.207546e-04     5.400128e-05     1.650936e-04     3.307070e-05
[14,]     3.363755e-05     3.846311e-05     9.003498e-05     3.105949e-05
[15,]     4.633976e-05     3.094512e-07     8.797136e-05     1.194676e-04
[16,]     6.499309e-07     3.048933e-05     5.785986e-05     6.867118e-05
      9758.98484948898 9778.58875327006 9816.19928572776 9868.0596360162
 [1,]     9.930708e-05     4.146051e-05     7.267820e-06    2.259877e-05
 [2,]     2.988199e-05     1.073717e-04     8.501099e-06    1.049191e-04
 [3,]     1.169614e-05     3.755046e-05     1.381443e-04    9.387456e-05
 [4,]     2.357303e-05     1.194799e-04     5.398691e-05    1.036863e-04
 [5,]     4.120060e-06     2.066768e-05     7.151736e-05    3.017069e-05
 [6,]     9.993914e-05     1.440047e-04     2.984839e-05    1.005554e-04
 [7,]     1.078357e-04     4.158405e-05     8.862955e-05    9.998097e-05
 [8,]     3.404820e-05     1.771820e-05     9.486378e-05    1.498082e-05
 [9,]     4.784168e-05     9.505246e-05     9.453302e-06    6.172306e-07
[10,]     9.193471e-05     8.216938e-05     1.574864e-05    3.292314e-06
[11,]     4.019197e-05     3.537123e-05     9.276336e-05    2.902876e-05
[12,]     1.438380e-05     5.411481e-05     9.709943e-06    6.647868e-07
[13,]     1.546480e-05     9.782903e-06     3.677364e-06    1.086279e-04
[14,]     8.758532e-06     5.572453e-05     4.820779e-05    5.757439e-05
[15,]     2.270904e-05     5.937699e-06     6.938062e-05    1.145165e-04
[16,]     6.676996e-06     6.737977e-05     1.114246e-04    7.289682e-07
      9885.32063493821 9925.96821730658 9954.731082253 9973.72548001882
 [1,]     9.872940e-05     9.778104e-05   1.268935e-04     3.655475e-05
 [2,]     3.550631e-05     3.760143e-05   6.908382e-06     2.337041e-05
 [3,]     4.982446e-05     1.111388e-04   1.313879e-04     3.940600e-05
 [4,]     2.475594e-05     9.702897e-05   3.584827e-06     1.010448e-04
 [5,]     3.200797e-05     1.155256e-05   6.371528e-05     6.154948e-05
 [6,]     1.003540e-04     5.605564e-05   8.255832e-05     1.340360e-04
 [7,]     1.025027e-04     6.626735e-05   1.451279e-04     9.513360e-05
 [8,]     1.866566e-05     4.526266e-05   6.601159e-05     5.462779e-05
 [9,]     7.994335e-06     9.199696e-05   1.253754e-05     4.112972e-05
[10,]     5.602521e-05     1.825694e-05   3.643667e-05     2.247582e-05
[11,]     4.681403e-05     1.681284e-05   3.646940e-05     9.619929e-05
[12,]     3.403309e-05     5.396253e-05   3.627615e-05     8.706792e-05
[13,]     1.087972e-04     1.227843e-04   5.098598e-05     1.204490e-04
[14,]     7.184304e-05     1.978984e-05   5.369319e-05     1.037623e-04
[15,]     2.761308e-05     1.065435e-04   1.216641e-04     2.497984e-05
[16,]     3.127558e-06     3.814242e-05   1.048949e-04     7.129167e-07
attr(,"mass")
  [1] 2010.464 2019.388 2023.298 2038.730 2045.653 2067.612 2072.048 2081.623
  [9] 2092.402 2105.133 2113.539 2123.283 2134.166 2141.353 2153.282 2159.370
 [17] 2163.253 2178.822 2190.670 2194.262 2209.906 2230.594 2245.133 2251.275
 [25] 2264.908 2279.203 2294.359 2306.623 2318.543 2325.932 2337.128 2347.718
 [33] 2354.248 2368.943 2379.276 2393.051 2403.468 2409.226 2423.587 2428.264
 [41] 2449.528 2467.002 2481.947 2495.004 2509.809 2522.310 2535.618 2545.104
 [49] 2553.945 2562.811 2581.172 2588.179 2604.235 2610.983 2616.838 2623.806
 [57] 2632.543 2643.166 2660.077 2672.767 2682.197 2686.813 2698.685 2709.039
 [65] 2714.062 2723.839 2730.790 2739.437 2754.069 2769.294 2789.914 2808.203
 [73] 2819.463 2830.953 2845.176 2862.365 2876.960 2887.747 2902.462 2915.283
 [81] 2932.338 2952.548 2972.578 2978.517 2990.574 3002.094 3016.728 3035.637
 [89] 3045.087 3059.820 3075.786 3087.670 3096.109 3104.492 3117.843 3143.326
 [97] 3158.628 3175.408 3191.721 3207.901 3217.086 3231.970 3240.807 3262.802
[105] 3279.245 3294.336 3303.236 3316.028 3325.276 3334.412 3365.637 3377.660
[113] 3408.431 3413.847 3419.966 3434.408 3448.202 3466.120 3489.130 3506.147
[121] 3524.497 3540.243 3568.710 3585.148 3605.547 3625.180 3634.955 3676.100
[129] 3680.957 3704.580 3719.562 3734.741 3744.047 3763.760 3768.644 3782.783
[137] 3795.357 3816.792 3860.968 3874.162 3882.934 3903.288 3909.626 3918.984
[145] 3934.213 3952.561 3961.831 3993.468 4019.242 4036.899 4053.877 4063.528
[153] 4071.561 4091.283 4110.141 4123.034 4145.672 4152.801 4166.886 4173.349
[161] 4192.541 4209.963 4235.389 4249.606 4266.145 4301.885 4337.888 4354.649
[169] 4363.025 4382.479 4397.039 4416.710 4430.783 4452.028 4466.989 4493.822
[177] 4530.408 4567.345 4586.415 4599.518 4613.850 4644.421 4663.772 4671.410
[185] 4696.965 4711.601 4717.148 4738.660 4756.701 4787.358 4818.734 4858.583
[193] 4889.085 4915.919 4929.085 4963.956 5003.961 5025.244 5045.645 5113.138
[201] 5131.497 5159.755 5204.149 5224.976 5248.491 5265.254 5292.677 5320.125
[209] 5336.945 5376.670 5391.092 5438.423 5463.825 5481.721 5507.572 5523.070
[217] 5549.992 5592.824 5631.538 5664.267 5709.068 5754.391 5805.847 5861.036
[225] 5888.020 5904.842 5944.594 5958.676 6004.039 6049.858 6075.112 6090.333
[233] 6142.796 6175.078 6191.473 6226.721 6271.807 6305.281 6332.679 6357.332
[241] 6378.453 6389.361 6433.139 6486.534 6492.818 6518.512 6561.920 6590.744
[249] 6614.658 6631.487 6669.480 6685.712 6749.089 6805.511 6831.844 6852.739
[257] 6909.415 6961.473 7008.657 7068.143 7131.327 7163.242 7191.426 7222.085
[265] 7291.270 7307.412 7343.962 7399.452 7432.803 7471.773 7550.614 7565.972
[273] 7636.797 7721.884 7766.712 7806.772 7820.654 7836.691 7875.971 7909.148
[281] 7924.860 8038.513 8092.124 8129.587 8143.627 8194.456 8217.087 8262.104
[289] 8301.565 8364.230 8429.664 8452.853 8500.583 8534.326 8566.378 8604.219
[297] 8657.569 8663.534 8677.052 8737.822 8757.217 8767.256 8810.002 8864.011
[305] 8935.253 8989.097 9062.593 9131.459 9136.944 9176.137 9290.928 9345.100
[313] 9425.915 9438.437 9478.958 9525.016 9590.430 9629.079 9715.852 9724.358
[321] 9758.985 9778.589 9816.199 9868.060 9885.321 9925.968 9954.731 9973.725
(def fm-dims (r->clj (r "dim(fm)")))
fm-dims
[16 328]

The feature matrix has 16 rows (one per spectrum) and hundreds of columns (one per shared peak position).

(let [col-masses (r->clj (r "as.numeric(colnames(fm))"))
      n-cols (count col-masses)]
  {:n-samples (first fm-dims)
   :n-features (second fm-dims)
   :mass-range [(first col-masses) (last col-masses)]})
{:n-samples 16,
 :n-features 328,
 :mass-range [2010.46372181388 9973.72548001882]}

Summary

ClojisR makes it straightforward to call MALDIquant from Clojure. The key patterns are:

  • (r ...) with backtick-quoted forms and tilde splicing
  • (r->clj ...) to bring R results back to Clojure vectors
  • R string expressions for complex nested calls

For multi-sample workflows (alignment, binning, feature matrices), MALDIquant can be called directly via ClojisR as shown above.

MALDIquant function Purpose Key parameters
createMassSpectrum Create spectrum object mass, intensity
trim Restrict mass range range
transformIntensity Variance stabilization method (“sqrt”)
smoothIntensity Noise reduction method, halfWindowSize
removeBaseline Background subtraction method (“SNIP”), iterations
calibrateIntensity Area normalization method (“TIC”)
detectPeaks Find peaks halfWindowSize, SNR, method
alignSpectra Correct m/z drift halfWindowSize, SNR
binPeaks Cluster peaks across spectra (tolerance)
filterPeaks Remove infrequent peaks minFrequency
intensityMatrix Build sample x feature matrix peaks, spectra
source: notebooks/ripple_book/maldiquant_interop.clj