Repository 'camera_annotate'
hg clone https://toolshed.g2.bx.psu.edu/repos/lecorguille/camera_annotate

Changeset 16:a2c49996603e (2018-09-19)
Previous changeset 15:6139bfcc95cb (2018-03-01) Next changeset 17:73d82de36369 (2018-10-09)
Commit message:
planemo upload commit e440674afa3e58c46100b0ac7c305a6f46ecbbdc
modified:
CAMERA.r
README.rst
abims_CAMERA_annotateDiffreport.xml
lib.r
macros.xml
repository_dependencies.xml
added:
test-data/.RData
test-data/CMD.mac.sh
test-data/Rplots.pdf
test-data/dataMatrixOutput.txt
test-data/faahKO-single.xset.merged.group.retcor.group.fillpeaks-xcms3.RData
test-data/faahKO.xset.group.retcor.group.fillPeaks.annotate.dataMatrix-xcms3.tsv
test-data/faahKO.xset.group.retcor.group.fillPeaks.annotate.variableMetadata-xcms3.tsv
b
diff -r 6139bfcc95cb -r a2c49996603e CAMERA.r
--- a/CAMERA.r Thu Mar 01 08:29:32 2018 -0500
+++ b/CAMERA.r Wed Sep 19 03:22:21 2018 -0400
[
@@ -4,16 +4,18 @@
 
 
 # ----- PACKAGE -----
-cat("\tPACKAGE INFO\n")
+cat("\tSESSION INFO\n")
 
-setRepositories(graphics=F, ind=31)
+pkgs=c("CAMERA","multtest","batch")
+for(pkg in pkgs) suppressPackageStartupMessages( stopifnot( library(pkg, quietly=TRUE, logical.return=TRUE, character.only=TRUE)))
 
-#pkgs=c("xcms","batch")
-pkgs=c("parallel","BiocGenerics", "Biobase", "Rcpp", "mzR", "xcms","snow","igraph","CAMERA","multtest","batch")
-for(p in pkgs) {
-    suppressPackageStartupMessages(suppressWarnings(library(p, quietly=TRUE, logical.return=TRUE, character.only=TRUE)))
-    cat(p,"\t",as.character(packageVersion(p)),"\n",sep="")
-}
+sessioninfo = sessionInfo()
+cat(sessioninfo$R.version$version.string,"\n")
+cat("Main packages:\n")
+for (pkg in names(sessioninfo$otherPkgs)) { cat(paste(pkg,packageVersion(pkg)),"\t") }; cat("\n")
+cat("Other loaded packages:\n")
+for (pkg in names(sessioninfo$loadedOnly)) { cat(paste(pkg,packageVersion(pkg)),"\t") }; cat("\n")
+
 source_local <- function(fname){ argv <- commandArgs(trailingOnly = FALSE); base_dir <- dirname(substring(argv[grep("--file=", argv)], 8)); source(paste(base_dir, fname, sep="/")) }
 
 cat("\n\n");
@@ -105,9 +107,13 @@
     directory = retrieveRawfileInTheWorkingDirectory(singlefile, zipfile)
 }
 
+# Because so far CAMERA isn't compatible with the new XCMSnExp object
+if (exists("xdata")){
+    xset <- getxcmsSetObject(xdata)
+}
 
-#addition of xset object to the list of arguments in the first position
-if (exists("xset") != 0){
+# addition of xset object to the list of arguments in the first position
+if (exists("xset")){
     listArguments=append(list(xset), listArguments)
 }
 
@@ -122,7 +128,6 @@
 #change the default display settings
 pdf(file=rplotspdf, width=16, height=12)
 
-
 if (thefunction %in% c("annotatediff")) {
     results_list=annotatediff(xset=xset,listArguments=listArguments,variableMetadataOutput=variableMetadataOutput,dataMatrixOutput=dataMatrixOutput)
     xa=results_list[["xa"]]
b
diff -r 6139bfcc95cb -r a2c49996603e README.rst
--- a/README.rst Thu Mar 01 08:29:32 2018 -0500
+++ b/README.rst Wed Sep 19 03:22:21 2018 -0400
b
@@ -2,6 +2,10 @@
 Changelog/News
 --------------
 
+**Version 2.2.3 - 30/04/2018**
+
+- NEW: support the new xcms 3.0.0 wrapper
+
 **Version 2.2.2 - 01/03/2018**
 
 - UPGRADE: upgrate the CAMERA version from 1.26.0 to 1.32.0
b
diff -r 6139bfcc95cb -r a2c49996603e abims_CAMERA_annotateDiffreport.xml
--- a/abims_CAMERA_annotateDiffreport.xml Thu Mar 01 08:29:32 2018 -0500
+++ b/abims_CAMERA_annotateDiffreport.xml Wed Sep 19 03:22:21 2018 -0400
b
@@ -1,4 +1,4 @@
-<tool id="abims_CAMERA_annotateDiffreport" name="CAMERA.annotate" version="2.2.2">
+<tool id="abims_CAMERA_annotateDiffreport" name="CAMERA.annotate" version="2.2.3">
 
     <description>CAMERA annotate function. Returns annotation results (isotope peaks, adducts and fragments) and a diffreport if more than one condition.</description>
 
@@ -191,14 +191,14 @@
     </outputs>
 
     <tests>
-        <test>
+        <!--<test>
             <param name="image" value="faahKO.xset.group.retcor.group.fillpeaks.RData"/>
             <expand macro="test_annotate_diffreport"/>
             <expand macro="test_annotate_quick_true"/>
             <expand macro="test_file_load_zip"/>
             <output name="variableMetadata" file="faahKO.xset.group.retcor.group.fillPeaks.annotate.variableMetadata.tsv" />
             <output name="datamatrix" file="faahKO.xset.group.retcor.group.fillPeaks.annotate.dataMatrix.tsv" />
-        </test>
+        </test>-->
         <test>
             <param name="image" value="faahKO-single.xset.merged.group.retcor.group.fillpeaks.RData"/>
             <expand macro="test_annotate_diffreport"/>
@@ -208,6 +208,14 @@
             <output name="datamatrix" file="faahKO.xset.group.retcor.group.fillPeaks.annotate.dataMatrix.tsv" />
         </test>
         <test>
+            <param name="image" value="faahKO-single.xset.merged.group.retcor.group.fillpeaks-xcms3.RData"/>
+            <expand macro="test_annotate_diffreport"/>
+            <expand macro="test_annotate_quick_true"/>
+            <expand macro="test_file_load_single"/>
+            <output name="variableMetadata" file="faahKO.xset.group.retcor.group.fillPeaks.annotate.variableMetadata-xcms3.tsv" />
+            <output name="datamatrix" file="faahKO.xset.group.retcor.group.fillPeaks.annotate.dataMatrix-xcms3.tsv" />
+        </test>
+        <test>
             <param name="image" value="faahKO-single.xset.merged.group.retcor.group.fillpeaks.RData"/>
             <expand macro="test_annotate_diffreport"/>
             <conditional name="quick_block">
@@ -424,6 +432,10 @@
 Changelog/News
 --------------
 
+**Version 2.2.3 - 30/04/2018**
+
+- NEW: support the new xcms 3.0.0 wrapper
+
 **Version 2.2.2 - 01/03/2018**
 
 - UPGRADE: upgrate the CAMERA version from 1.26.0 to 1.32.0
b
diff -r 6139bfcc95cb -r a2c49996603e lib.r
--- a/lib.r Thu Mar 01 08:29:32 2018 -0500
+++ b/lib.r Wed Sep 19 03:22:21 2018 -0400
[
b'@@ -1,5 +1,25 @@\n # lib.r\n \n+# This function retrieve a xset like object\n+#@author Gildas Le Corguille lecorguille@sb-roscoff.fr\n+getxcmsSetObject <- function(xobject) {\n+    # XCMS 1.x\n+    if (class(xobject) == "xcmsSet")\n+        return (xobject)\n+    # XCMS 3.x\n+    if (class(xobject) == "XCMSnExp") {\n+        # Get the legacy xcmsSet object\n+        suppressWarnings(xset <- as(xobject, \'xcmsSet\'))\n+        if (is.null(xset@phenoData$sample_group))\n+            sampclass(xset) = "."\n+        else\n+            sampclass(xset) <- xset@phenoData$sample_group\n+        if (!is.null(xset@phenoData$sample_name))\n+            rownames(xset@phenoData) = xset@phenoData$sample_name\n+        return (xset)\n+    }\n+}\n+\n #@author G. Le Corguille\n #The function create a pdf from the different png generated by diffreport\n diffreport_png2pdf <- function(filebase) {\n@@ -105,29 +125,29 @@\n             for (n in seq(along=y)){\n                 if(i+n <= length(classes)){\n                     filebase=paste(classes[i],class2=classes[i+n],sep="-vs-")\n-                    \n-                    diffrep=diffreport(object=xset,class1=classes[i],class2=classes[i+n],filebase=filebase,eicmax=listArguments[["eicmax"]],eicwidth=listArguments[["eicwidth"]],sortpval=TRUE,value=listArguments[["value"]],h=listArguments[["h"]],w=listArguments[["w"]],mzdec=listArguments[["mzdec"]])\n-                    \n+\n+                    diffrep=diffreport(object=xset,class1=classes[i],class2=classes[i+n],filebase=filebase,eicmax=listArguments[["eicmax"]],eicwidth=listArguments[["eicwidth"]],sortpval=TRUE,value=listArguments[["value"]],h=listArguments[["h"]],w=listArguments[["w"]],mzdec=listArguments[["mzdec"]],missing=0)\n+\n                     diffrepOri = diffrep\n-                    \n+\n                     # renamming of the column rtmed to rt to fit with camera peaklist function output\n                     colnames(diffrep)[colnames(diffrep)=="rtmed"] <- "rt"\n                     colnames(diffrep)[colnames(diffrep)=="mzmed"] <- "mz"\n-                    \n+\n                     # combines results and reorder columns\n                     diffrep = merge(peakList, diffrep[,c("name","fold","tstat","pvalue")], by.x="name", by.y="name", sort=F)\n                     diffrep = cbind(diffrep[,!(colnames(diffrep) %in% c(sampnames(xa@xcmsSet)))],diffrep[,(colnames(diffrep) %in% c(sampnames(xa@xcmsSet)))])\n-                    \n+\n                     diffrep = RTSecondToMinute(diffrep, listArguments[["convertRTMinute"]])\n                     diffrep = formatIonIdentifiers(diffrep, numDigitsRT=listArguments[["numDigitsRT"]], numDigitsMZ=listArguments[["numDigitsMZ"]])\n-                    \n+\n                     if(listArguments[["sortpval"]]){\n                         diffrep=diffrep[order(diffrep$pvalue), ]\n                     }\n-                    \n+\n                     dir.create("tabular")\n                     write.table(diffrep, sep="\\t", quote=FALSE, row.names=FALSE, file=paste("tabular/",filebase,"_tsv.tabular",sep=""))\n-                    \n+\n                     if (listArguments[["eicmax"]] != 0) {\n                         diffreport_png2pdf(filebase)\n                     }\n@@ -136,7 +156,6 @@\n         }\n     }\n \n-\n     # --- variableMetadata ---\n     variableMetadata=peakList[,!(make.names(colnames(peakList)) %in% c(make.names(sampnames(xa@xcmsSet))))]\n     variableMetadata = RTSecondToMinute(variableMetadata, listArguments[["convertRTMinute"]])\n@@ -302,3 +321,288 @@\n     }\n     return (directory)\n }\n+\n+#@TODO: remove this function as soon as we can use xcms 3.x.x from Bioconductor 3.7\n+# https://github.com/sneumann/CAMERA/issues/33#issuecomment-405168524\n+# https://github.com/sneumann/xcms/commit/950a3fe794cdb6b0fda88696e31aab3d97a3b7dd\n+############################################################\n+## getEIC\n+getEIC <- function(object, mzrange, rtrange = 200,\n+                   groupidx, sampleidx = sampnames(object),\n+                   rt = c("corrected", "raw'..b' xcms:::pval(testval, testclab, tstat)\n+    } else {\n+        message("Too few samples per class, skipping t-test.")\n+        tstat <- pvalue <- rep(NA,nrow(testval))\n+    }\n+    stat <- data.frame(fold = fold, tstat = tstat, pvalue = pvalue)\n+    if (length(levels(sampclass(object))) >2) {\n+        pvalAnova<-c()\n+        for(i in 1:nrow(values)){\n+            var<-as.numeric(values[i,])\n+            ano<-summary(aov(var ~ sampclass(object)) )\n+            pvalAnova<-append(pvalAnova, unlist(ano)["Pr(>F)1"])\n+        }\n+        stat<-cbind(stat, anova= pvalAnova)\n+    }\n+    if (metlin) {\n+        neutralmass <- groupmat[,"mzmed"] + ifelse(metlin < 0, 1, -1)\n+        metlin <- abs(metlin)\n+        digits <- ceiling(-log10(metlin))+1\n+        metlinurl <-\n+        paste("http://metlin.scripps.edu/simple_search_result.php?mass_min=",\n+        round(neutralmass - metlin, digits), "&mass_max=",\n+        round(neutralmass + metlin, digits), sep="")\n+        values <- cbind(metlin = metlinurl, values)\n+    }\n+    twosamp <- cbind(name = groupnames(object), stat, groupmat, values)\n+    if (sortpval) {\n+        tsidx <- order(twosamp[,"pvalue"])\n+        twosamp <- twosamp[tsidx,]\n+        rownames(twosamp) <- 1:nrow(twosamp)\n+        values<-values[tsidx,]\n+    } else\n+    tsidx <- 1:nrow(values)\n+    \n+    if (length(filebase))\n+    write.table(twosamp, paste(filebase, ".tsv", sep = ""), quote = FALSE, sep = "\\t", col.names = NA)\n+    \n+    if (eicmax > 0) {\n+        if (length(unique(peaks(object)[,"rt"])) > 1) {\n+            ## This looks like "normal" LC data\n+            \n+            eicmax <- min(eicmax, length(tsidx))\n+            eics <- getEIC(object, rtrange = eicwidth*1.1, sampleidx = ceic,\n+            groupidx = tsidx[seq(length = eicmax)])\n+            \n+            if (length(filebase)) {\n+                eicdir <- paste(filebase, "_eic", sep="")\n+                boxdir <- paste(filebase, "_box", sep="")\n+                dir.create(eicdir)\n+                dir.create(boxdir)\n+                if (capabilities("png")){\n+                    xcms:::xcmsBoxPlot(values[seq(length = eicmax),],\n+                    sampclass(object), dirpath=boxdir, pic="png",  width=w, height=h)\n+                    png(file.path(eicdir, "%003d.png"), width = w, height = h)\n+                } else {\n+                    xcms:::xcmsBoxPlot(values[seq(length = eicmax),],\n+                    sampclass(object), dirpath=boxdir, pic="pdf", width=w, height=h)\n+                    pdf(file.path(eicdir, "%003d.pdf"), width = w/72,\n+                    height = h/72, onefile = FALSE)\n+                }\n+            }\n+            plot(eics, object, rtrange = eicwidth, mzdec=mzdec)\n+            \n+            if (length(filebase))\n+            dev.off()\n+        } else {\n+            ## This looks like a direct-infusion single spectrum\n+            if (length(filebase)) {\n+                eicdir <- paste(filebase, "_eic", sep="")\n+                boxdir <- paste(filebase, "_box", sep="")\n+                dir.create(eicdir)\n+                dir.create(boxdir)\n+                if (capabilities("png")){\n+                    xcmsBoxPlot(values[seq(length = eicmax),],\n+                    sampclass(object), dirpath=boxdir, pic="png",\n+                    width=w, height=h)\n+                    png(file.path(eicdir, "%003d.png"), width = w, height = h,\n+                    units = "px")\n+                } else {\n+                    xcmsBoxPlot(values[seq(length = eicmax),],\n+                    sampclass(object), dirpath=boxdir, pic="pdf",\n+                    width=w, height=h)\n+                    pdf(file.path(eicdir, "%003d.pdf"), width = w/72,\n+                    height = h/72, onefile = FALSE)\n+                }\n+            }\n+            \n+            plotSpecWindow(object, gidxs = tsidx[seq(length = eicmax)], borderwidth=1)\n+            \n+            if (length(filebase))\n+            dev.off()\n+        }\n+    }\n+    \n+    invisible(twosamp)\n+}\n'
b
diff -r 6139bfcc95cb -r a2c49996603e macros.xml
--- a/macros.xml Thu Mar 01 08:29:32 2018 -0500
+++ b/macros.xml Wed Sep 19 03:22:21 2018 -0400
b
@@ -3,7 +3,8 @@
     <xml name="requirements">
         <requirements>
             <requirement type="package" version="0.4_2">r-snow</requirement>
-            <requirement type="package" version="1.32.0">bioconductor-camera</requirement>
+            <requirement type="package" version="3.0.0">bioconductor-xcms</requirement>
+            <requirement type="package" version="1.34.0">bioconductor-camera</requirement>
             <requirement type="package" version="2.28.0">bioconductor-multtest</requirement>
             <requirement type="package" version="1.1_4">r-batch</requirement>
             <requirement type="package" version="1.3.23">graphicsmagick</requirement>
b
diff -r 6139bfcc95cb -r a2c49996603e repository_dependencies.xml
--- a/repository_dependencies.xml Thu Mar 01 08:29:32 2018 -0500
+++ b/repository_dependencies.xml Wed Sep 19 03:22:21 2018 -0400
b
@@ -1,5 +1,5 @@
 <?xml version="1.0"?>
 <repositories>
  <repository changeset_revision="4a51ab3d8ecf" name="rdata_camera_datatypes" owner="lecorguille" toolshed="https://toolshed.g2.bx.psu.edu" />
- <repository changeset_revision="d64562a4ebb3" name="rdata_xcms_datatypes" owner="lecorguille" toolshed="https://toolshed.g2.bx.psu.edu" />
+ <repository changeset_revision="544f6d2329ac" name="rdata_xcms_datatypes" owner="lecorguille" toolshed="https://toolshed.g2.bx.psu.edu" />
 </repositories>
b
diff -r 6139bfcc95cb -r a2c49996603e test-data/.RData
b
Binary file test-data/.RData has changed
b
diff -r 6139bfcc95cb -r a2c49996603e test-data/CMD.mac.sh
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/CMD.mac.sh Wed Sep 19 03:22:21 2018 -0400
b
@@ -0,0 +1,1 @@
+LC_ALL=C Rscript ../CAMERA.r xfunction annotatediff image 'faahKO-single.xset.merged.group.retcor.group.fillpeaks-xcms3.RData' nSlaves ${GALAXY_SLOTS:-1}  variableMetadataOutput variableMetadataOutput.txt dataMatrixOutput dataMatrixOutput.txt sigma 6 perfwhm 0.6  ppm 5 mzabs 0.015  maxcharge 3 maxiso 4 minfrac 0.5  quick TRUE xsetRdataOutput '/tmp/tmpT41cRC/files/000/dataset_20.dat'  runDiffreport TRUE eicmax 200 eicwidth 200 value into sortpval FALSE h 480 w 640 mzdec 2   convertRTMinute TRUE numDigitsMZ 4 numDigitsRT 1 intval into singlefile_galaxyPath 'wt15.CDF,ko16.CDF,ko15.CDF,wt16.CDF' singlefile_sampleName 'wt15.CDF,ko16.CDF,ko15.CDF,wt16.CDF'
b
diff -r 6139bfcc95cb -r a2c49996603e test-data/Rplots.pdf
b
Binary file test-data/Rplots.pdf has changed
b
diff -r 6139bfcc95cb -r a2c49996603e test-data/dataMatrixOutput.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/dataMatrixOutput.txt Wed Sep 19 03:22:21 2018 -0400
b
b'@@ -0,0 +1,6333 @@\n+name\t..ko15.CDF\t..ko16.CDF\t..wt15.CDF\t..wt16.CDF\n+M200T3488\t6029.94499999992\t1646.69331417609\tNA\tNA\n+M200T3544\tNA\tNA\t4586.52650000006\tNA\n+M200T3882\t1144.01499999999\t1716.38461192716\t1449.84316738066\tNA\n+M200T3822\tNA\tNA\t1018.81499999998\tNA\n+M200T3184\tNA\t4333.4849999999\tNA\tNA\n+M200T2867\t76879.3654931277\t506848.884266667\t53985.4351594335\t76607.9704340259\n+M200T3978\t1324.38867325921\tNA\t495.468308772648\t1242.61000000004\n+M200T3681\tNA\t645.622409362349\t5108.15999999988\t1576.85791107837\n+M201T3788\tNA\tNA\t1680.81000000006\tNA\n+M201T3536\tNA\t2863.6449999999\tNA\t1106.33802122178\n+M201T4052\t318.548929605957\t205.380230930156\t1339.64000000005\t689.071631937096\n+M201T3276\tNA\t1497.96219436021\t4090.91000000014\t1027.56552966783\n+M201T3630\tNA\tNA\t2121.84236212046\t2931.24499999993\n+M201T2590\tNA\t14307.2299999999\tNA\tNA\n+M201T3985\tNA\tNA\t1657.33500000006\tNA\n+M201T3138\t27830.0354581329\t254581.758476187\t63151.125622074\t28083.6051504647\n+M201T3693\tNA\tNA\t2906.20499999993\tNA\n+M201T3511\t4023.61499999999\t694.250350240936\tNA\tNA\n+M201T3790\t446.184819684897\t1118.41040705157\tNA\t3245.80999999992\n+M201T4047\tNA\tNA\t1275.47500000004\t1220.69999999998\n+M201T3925\t128.578946259232\tNA\t1284.72816666668\tNA\n+M201T3874\t886.610568446999\t1220.70000000004\t521.172821104497\tNA\n+M201T4152\tNA\t774.575999999962\tNA\tNA\n+M202T3859\t892.050000000031\t753.237186536855\tNA\tNA\n+M202T4007\tNA\tNA\t831.015000000029\tNA\n+M202T3303\t2759.095\t3872.96250000005\tNA\t791.36463944858\n+M202T3379\t9879.84499999977\t1320.19136751317\tNA\tNA\n+M202T3767\tNA\t2564.84247599701\t1151.71733333335\tNA\n+M202T3181\t7682.58500000004\t3806.07999999999\tNA\tNA\n+M202T4031\t389.407001047583\t688.600000000024\t754.233600000007\tNA\n+M202T3886\tNA\t841.96999999998\tNA\tNA\n+M202T4137\tNA\tNA\tNA\t732.419999999983\n+M202T3931\t354.703176356038\tNA\t1461.71\tNA\n+M202T3271\t77227.3769090911\t7016.51333681604\t44589.8848371899\t7116.58518941033\n+M202T2891\t36435.3079999996\t6433.2631875585\t29664.6514201645\t2640.91263857081\n+M202T3285\tNA\tNA\tNA\t3372.57500000002\n+M203T4040\tNA\tNA\tNA\t978.124999999977\n+M203T3308\tNA\t1694.89499999996\tNA\t1136.19000000004\n+M203T3338\t3852.53759999981\t2974.74816666671\tNA\tNA\n+M203T4156\tNA\t804.409999999981\tNA\t374.999685196314\n+M203T3714\tNA\tNA\tNA\t3345.96999999992\n+M203T3787\t790.967612320259\t2690.23500000009\t1939.932414103\t955.908725789886\n+M203T2679\tNA\tNA\tNA\t6346.07499999991\n+M203T3966\t1720.81523845536\t136.302197105375\t2796.65500000002\t1490.80528084004\n+M203T4059\t179.980340762509\tNA\tNA\t1550.91499999999\n+M204T3535\tNA\tNA\t967.169999999987\tNA\n+M204T3666\t3002.85120000003\tNA\t572.831482902373\tNA\n+M204T4031\t1111.14999999999\tNA\tNA\t129.105921444997\n+M204T4126\tNA\tNA\tNA\t793.455000000074\n+M204T3060\t3905.96266666662\tNA\tNA\tNA\n+M204T3845\tNA\t238.782610500926\t1676.11500000006\tNA\n+M204T3395\t706.339294579149\t155.065526139688\t1687.9990131206\t1586.91\n+M204T2638\tNA\tNA\t3518.12000000012\tNA\n+M204T4139\tNA\tNA\tNA\t1089.24000000004\n+M204T4073\tNA\t1062.63500000004\tNA\tNA\n+M204T3604\tNA\tNA\tNA\t1032.90000000004\n+M204T3698\tNA\tNA\tNA\t1446.05999999997\n+M205T3582\tNA\tNA\t964.383879177602\t2466.43999999994\n+M205T2582\t8805.94199999988\tNA\tNA\t9189.67999999988\n+M205T2789\t1924712.01585714\t1757150.9648\t2129885.09357143\t1634341.9855862\n+M205T3783\t2777.8750000001\tNA\t5180.14999999997\tNA\n+M205T3909\tNA\t3392.55866666671\tNA\tNA\n+M205T4051\t2600.69759999987\tNA\t1198.79000000004\tNA\n+M206T4070\tNA\tNA\t793.454999999982\tNA\n+M206T3084\tNA\tNA\tNA\t4682.47999999989\n+M206T3353\t271.698957954148\tNA\t641.649999999985\tNA\n+M206T3285\tNA\tNA\tNA\t1023.50999999998\n+M206T3192\t2337.7226842156\t10174.0649999999\tNA\t531.384477548613\n+M206T4108\t1685.41525000004\t800.018271889045\tNA\t280.678679423257\n+M206T2791\t213659.293920001\t272706.72055332\t253825.557279999\t241844.442068965\n+M206T3444\t1938.56022964102\t2400.20806107181\t4692.21741165203\t3090.62812500006\n+M206T3247\t1713.67500000006\tNA\tNA\t3137.82499999993\n+M206T3680\tNA\tNA\tNA\t2237.94999999995\n+M206T3191\t1931.20999999997\tNA\tNA\tNA\n+M206T3755\t510.0375381342\tNA\tNA\t2529.03999999997\n+M206T2570\tNA\tNA\t9324.27000000002\t924.940634138229\n+M206T3347\tNA\t528.62804858849\t535.449055667051\t1536.82999'..b'22.49815518322\t7801.52500000005\tNA\n+M597T3463\tNA\t7109.79500000002\tNA\tNA\n+M597T3287\t4115.2105592385\t6108.19500000008\tNA\tNA\n+M597T3669\t3813.49883333338\t1542.97947676868\tNA\t987.795272510788\n+M597T3874\tNA\tNA\t1399.11000000005\tNA\n+M597T3799\tNA\tNA\t9477.64000000033\tNA\n+M597T2920\t1034.24974715399\tNA\tNA\t3881.19999999991\n+M597T4080\t18556.9408188294\t8185.95298062284\t27582.3906250001\t12307.5901216578\n+M597T4135\t15187.814210526\t11262.4812687384\t3161.39825087393\t756.425882857505\n+M597T2791\t2459.83450887303\t1066.17398982618\t23087.5229375002\t1056.38098428999\n+M597T3016\tNA\t8213.12000000029\tNA\tNA\n+M597T2530\t1406.71781818182\tNA\tNA\t4424.01941666667\n+M597T3805\t28480.3432782255\t1186500.84141178\t10606.0571955243\t34573.1746514789\n+M598T3855\t1539.85635358667\tNA\t6475.97000000023\tNA\n+M598T3997\t4690.30499999994\tNA\t1288.07486775315\t305.764898642052\n+M598T2646\tNA\t914.483039050997\t1991.89121724326\t5450.66280000004\n+M598T2895\t929.575110080356\tNA\tNA\t2579.12000000009\n+M598T3923\t2421.055\tNA\tNA\t5100.33500000001\n+M598T2839\t4328.7899999999\tNA\tNA\tNA\n+M598T4225\t780.835199999962\tNA\tNA\tNA\n+M598T3867\tNA\tNA\t923.232000000008\tNA\n+M598T2736\tNA\tNA\t3392.92\tNA\n+M598T2560\tNA\tNA\t3350.66499999992\tNA\n+M598T3059\t1251.51527618732\tNA\t3904.17600000003\tNA\n+M598T3248\tNA\tNA\t11160.015\tNA\n+M598T2590\tNA\tNA\t5259.96499999993\tNA\n+M598T3570\t4893.75500000017\tNA\tNA\tNA\n+M598T3155\tNA\t791.797930803163\tNA\t1547.78500000005\n+M598T3717\t6735.75999999984\t29323.4050000004\tNA\t48876.8056363635\n+M598T3840\t990.400497127544\t19344.1924374998\t1705.49428412994\t18481.3189273045\n+M598T4210\tNA\tNA\tNA\t852.92499999998\n+M598T2983\t3463.10662220332\t7907.94500000028\t4597.13130349332\tNA\n+M598T3218\tNA\t16185.5021250003\t8818.83384484759\t12536.1107384566\n+M598T2829\t5337.7886250001\t3464.81226746294\t2392.55846345708\tNA\n+M598T3363\tNA\tNA\tNA\t7501.04499999995\n+M598T3407\t3147.74550046427\tNA\t4532.24000000016\tNA\n+M598T3496\tNA\tNA\t9308.62000000032\tNA\n+M598T3120\t3137.82499999993\tNA\tNA\tNA\n+M598T4219\t1558.64039999994\t130.408333333329\t574.318299999979\t283.24689999999\n+M598T2675\t2139.08160000002\t2826.02880000002\tNA\tNA\n+M598T2759\tNA\t5078.42499999999\t1225.70122914555\t436.47790863468\n+M598T2565\tNA\t1600.99499999996\tNA\t264.093749999994\n+M598T3816\t4056.74513901089\t126881.278714286\t2038.59965809119\tNA\n+M598T4028\t4444.59999999997\tNA\tNA\tNA\n+M598T3699\t21485.7955370261\t262005.048159626\t19097.5044463651\t65022.6200000009\n+M599T3177\tNA\t2838.6076666667\tNA\t2164.39499999999\n+M599T2937\t5756.0700000002\t2557.20999999994\tNA\tNA\n+M599T2856\tNA\tNA\tNA\t3757.56499999997\n+M599T4215\tNA\tNA\t223.746998896449\t881.094999999979\n+M599T2629\t24597.6221333335\t2237.94999999995\tNA\tNA\n+M599T3087\tNA\tNA\tNA\t5745.1150000002\n+M599T4161\t1370.93999999996\t8068.71033333314\t769.979999999982\t5823.36499999964\n+M599T3725\t22852.13\t8220.94499999976\t55051.9369090907\t5147.52485215715\n+M599T3846\tNA\tNA\t5238.05500000018\tNA\n+M599T2884\t3106.80875823442\t3346.02697646104\t6247.93574186843\t3269.28500000002\n+M599T4044\t798856.451441174\t40631.1229886866\t897286.005000001\t49928.058471187\n+M599T3927\t30327.7176201825\t15327.6100000005\t37243.228306282\t2271127.28014893\n+M599T3455\t35592.8437333336\tNA\tNA\t6298.5500000001\n+M599T3741\t2472129.73859259\tNA\t797434.72246154\t11463.9494185942\n+M599T4124\tNA\t354770.43852632\tNA\t707685.939314288\n+M599T3815\t79625.6350000028\t12736.5568357265\t9836.53488085646\t4148.40133264837\n+M599T2837\t2295.00891881701\t5415.77280000005\t3360.28178830047\tNA\n+M599T2984\tNA\t5849.96999999986\tNA\t12502.4521250001\n+M600T3611\tNA\tNA\tNA\t5001.73999999997\n+M600T2728\tNA\tNA\tNA\t1760.62499999996\n+M600T3412\t3353.79499999992\t1647.14222654083\tNA\tNA\n+M600T2595\t6022.12000000021\t2928.1150000001\t3209.81499999996\t2331.63714285708\n+M600T2782\tNA\tNA\tNA\t2608.74386666668\n+M600T3443\t4671.52499999989\t8782.77999999999\tNA\tNA\n+M600T2847\tNA\tNA\t1920.87948584165\t1987.33833333336\n+M600T2963\t9211.58999999991\tNA\t2898.28401908135\t5635.39317528082\n+M600T3070\t2056.42917003541\tNA\t1251.24642515022\t3430.04160000003\n+M600T3353\t2327.15500000008\tNA\tNA\t2315.33491939702\n+M600T3212\t7784.31000000027\t2982.99537949918\tNA\tNA\n+M600T2573\tNA\t3639.80233333338\tNA\tNA\n'
b
diff -r 6139bfcc95cb -r a2c49996603e test-data/faahKO-single.xset.merged.group.retcor.group.fillpeaks-xcms3.RData
b
Binary file test-data/faahKO-single.xset.merged.group.retcor.group.fillpeaks-xcms3.RData has changed
b
diff -r 6139bfcc95cb -r a2c49996603e test-data/faahKO.xset.group.retcor.group.fillPeaks.annotate.dataMatrix-xcms3.tsv
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/faahKO.xset.group.retcor.group.fillPeaks.annotate.dataMatrix-xcms3.tsv Wed Sep 19 03:22:21 2018 -0400
b
b'@@ -0,0 +1,6333 @@\n+name\tko15\tko16\twt15\twt16\n+M200T3488\t6029.94499999992\t1646.69331417609\tNA\tNA\n+M200T3544\tNA\tNA\t4586.52650000006\tNA\n+M200T3882\t1144.01499999999\t1716.38461192716\t1449.84316738066\tNA\n+M200T3822\tNA\tNA\t1018.81499999998\tNA\n+M200T3184\tNA\t4333.4849999999\tNA\tNA\n+M200T2867\t76879.3654931277\t506848.884266667\t53985.4351594335\t76607.9704340259\n+M200T3978\t1324.38867325921\tNA\t495.468308772648\t1242.61000000004\n+M200T3681\tNA\t645.622409362349\t5108.15999999988\t1576.85791107837\n+M201T3788\tNA\tNA\t1680.81000000006\tNA\n+M201T3536\tNA\t2863.6449999999\tNA\t1106.33802122178\n+M201T4052\t318.548929605957\t205.380230930156\t1339.64000000005\t689.071631937096\n+M201T3276\tNA\t1497.96219436021\t4090.91000000014\t1027.56552966783\n+M201T3630\tNA\tNA\t2121.84236212046\t2931.24499999993\n+M201T2590\tNA\t14307.2299999999\tNA\tNA\n+M201T3985\tNA\tNA\t1657.33500000006\tNA\n+M201T3138\t27830.0354581329\t254581.758476187\t63151.125622074\t28083.6051504647\n+M201T3693\tNA\tNA\t2906.20499999993\tNA\n+M201T3511\t4023.61499999999\t694.250350240936\tNA\tNA\n+M201T3790\t446.184819684897\t1118.41040705157\tNA\t3245.80999999992\n+M201T4047\tNA\tNA\t1275.47500000004\t1220.69999999998\n+M201T3925\t128.578946259232\tNA\t1284.72816666668\tNA\n+M201T3874\t886.610568446999\t1220.70000000004\t521.172821104497\tNA\n+M201T4152\tNA\t774.575999999962\tNA\tNA\n+M202T3859\t892.050000000031\t753.237186536855\tNA\tNA\n+M202T4007\tNA\tNA\t831.015000000029\tNA\n+M202T3303\t2759.095\t3872.96250000005\tNA\t791.36463944858\n+M202T3379\t9879.84499999977\t1320.19136751317\tNA\tNA\n+M202T3767\tNA\t2564.84247599701\t1151.71733333335\tNA\n+M202T3181\t7682.58500000004\t3806.07999999999\tNA\tNA\n+M202T4031\t389.407001047583\t688.600000000024\t754.233600000007\tNA\n+M202T3886\tNA\t841.96999999998\tNA\tNA\n+M202T4137\tNA\tNA\tNA\t732.419999999983\n+M202T3931\t354.703176356038\tNA\t1461.71\tNA\n+M202T3271\t77227.3769090911\t7016.51333681604\t44589.8848371899\t7116.58518941033\n+M202T2891\t36435.3079999996\t6433.2631875585\t29664.6514201645\t2640.91263857081\n+M202T3285\tNA\tNA\tNA\t3372.57500000002\n+M203T4040\tNA\tNA\tNA\t978.124999999977\n+M203T3308\tNA\t1694.89499999996\tNA\t1136.19000000004\n+M203T3338\t3852.53759999981\t2974.74816666671\tNA\tNA\n+M203T4156\tNA\t804.409999999981\tNA\t374.999685196314\n+M203T3714\tNA\tNA\tNA\t3345.96999999992\n+M203T3787\t790.967612320259\t2690.23500000009\t1939.932414103\t955.908725789886\n+M203T2679\tNA\tNA\tNA\t6346.07499999991\n+M203T3966\t1720.81523845536\t136.302197105375\t2796.65500000002\t1490.80528084004\n+M203T4059\t179.980340762509\tNA\tNA\t1550.91499999999\n+M204T3535\tNA\tNA\t967.169999999987\tNA\n+M204T3666\t3002.85120000003\tNA\t572.831482902373\tNA\n+M204T4031\t1111.14999999999\tNA\tNA\t129.105921444997\n+M204T4126\tNA\tNA\tNA\t793.455000000074\n+M204T3060\t3905.96266666662\tNA\tNA\tNA\n+M204T3845\tNA\t238.782610500926\t1676.11500000006\tNA\n+M204T3395\t706.339294579149\t155.065526139688\t1687.9990131206\t1586.91\n+M204T2638\tNA\tNA\t3518.12000000012\tNA\n+M204T4139\tNA\tNA\tNA\t1089.24000000004\n+M204T4073\tNA\t1062.63500000004\tNA\tNA\n+M204T3604\tNA\tNA\tNA\t1032.90000000004\n+M204T3698\tNA\tNA\tNA\t1446.05999999997\n+M205T3582\tNA\tNA\t964.383879177602\t2466.43999999994\n+M205T2582\t8805.94199999988\tNA\tNA\t9189.67999999988\n+M205T2789\t1924712.01585714\t1757150.9648\t2129885.09357143\t1634341.9855862\n+M205T3783\t2777.8750000001\tNA\t5180.14999999997\tNA\n+M205T3909\tNA\t3392.55866666671\tNA\tNA\n+M205T4051\t2600.69759999987\tNA\t1198.79000000004\tNA\n+M206T4070\tNA\tNA\t793.454999999982\tNA\n+M206T3084\tNA\tNA\tNA\t4682.47999999989\n+M206T3353\t271.698957954148\tNA\t641.649999999985\tNA\n+M206T3285\tNA\tNA\tNA\t1023.50999999998\n+M206T3192\t2337.7226842156\t10174.0649999999\tNA\t531.384477548613\n+M206T4108\t1685.41525000004\t800.018271889045\tNA\t280.678679423257\n+M206T2791\t213659.293920001\t272706.72055332\t253825.557279999\t241844.442068965\n+M206T3444\t1938.56022964102\t2400.20806107181\t4692.21741165203\t3090.62812500006\n+M206T3247\t1713.67500000006\tNA\tNA\t3137.82499999993\n+M206T3680\tNA\tNA\tNA\t2237.94999999995\n+M206T3191\t1931.20999999997\tNA\tNA\tNA\n+M206T3755\t510.0375381342\tNA\tNA\t2529.03999999997\n+M206T2570\tNA\tNA\t9324.27000000002\t924.940634138229\n+M206T3347\tNA\t528.62804858849\t535.449055667051\t1536.82999999996\n+M206T3578\t256.65'..b'22.49815518322\t7801.52500000005\tNA\n+M597T3463\tNA\t7109.79500000002\tNA\tNA\n+M597T3287\t4115.2105592385\t6108.19500000008\tNA\tNA\n+M597T3669\t3813.49883333338\t1542.97947676868\tNA\t987.795272510788\n+M597T3874\tNA\tNA\t1399.11000000005\tNA\n+M597T3799\tNA\tNA\t9477.64000000033\tNA\n+M597T2920\t1034.24974715399\tNA\tNA\t3881.19999999991\n+M597T4080\t18556.9408188294\t8185.95298062284\t27582.3906250001\t12307.5901216578\n+M597T4135\t15187.814210526\t11262.4812687384\t3161.39825087393\t756.425882857505\n+M597T2791\t2459.83450887303\t1066.17398982618\t23087.5229375002\t1056.38098428999\n+M597T3016\tNA\t8213.12000000029\tNA\tNA\n+M597T2530\t1406.71781818182\tNA\tNA\t4424.01941666667\n+M597T3805\t28480.3432782255\t1186500.84141178\t10606.0571955243\t34573.1746514789\n+M598T3855\t1539.85635358667\tNA\t6475.97000000023\tNA\n+M598T3997\t4690.30499999994\tNA\t1288.07486775315\t305.764898642052\n+M598T2646\tNA\t914.483039050997\t1991.89121724326\t5450.66280000004\n+M598T2895\t929.575110080356\tNA\tNA\t2579.12000000009\n+M598T3923\t2421.055\tNA\tNA\t5100.33500000001\n+M598T2839\t4328.7899999999\tNA\tNA\tNA\n+M598T4225\t780.835199999962\tNA\tNA\tNA\n+M598T3867\tNA\tNA\t923.232000000008\tNA\n+M598T2736\tNA\tNA\t3392.92\tNA\n+M598T2560\tNA\tNA\t3350.66499999992\tNA\n+M598T3059\t1251.51527618732\tNA\t3904.17600000003\tNA\n+M598T3248\tNA\tNA\t11160.015\tNA\n+M598T2590\tNA\tNA\t5259.96499999993\tNA\n+M598T3570\t4893.75500000017\tNA\tNA\tNA\n+M598T3155\tNA\t791.797930803163\tNA\t1547.78500000005\n+M598T3717\t6735.75999999984\t29323.4050000004\tNA\t48876.8056363635\n+M598T3840\t990.400497127544\t19344.1924374998\t1705.49428412994\t18481.3189273045\n+M598T4210\tNA\tNA\tNA\t852.92499999998\n+M598T2983\t3463.10662220332\t7907.94500000028\t4597.13130349332\tNA\n+M598T3218\tNA\t16185.5021250003\t8818.83384484759\t12536.1107384566\n+M598T2829\t5337.7886250001\t3464.81226746294\t2392.55846345708\tNA\n+M598T3363\tNA\tNA\tNA\t7501.04499999995\n+M598T3407\t3147.74550046427\tNA\t4532.24000000016\tNA\n+M598T3496\tNA\tNA\t9308.62000000032\tNA\n+M598T3120\t3137.82499999993\tNA\tNA\tNA\n+M598T4219\t1558.64039999994\t130.408333333329\t574.318299999979\t283.24689999999\n+M598T2675\t2139.08160000002\t2826.02880000002\tNA\tNA\n+M598T2759\tNA\t5078.42499999999\t1225.70122914555\t436.47790863468\n+M598T2565\tNA\t1600.99499999996\tNA\t264.093749999994\n+M598T3816\t4056.74513901089\t126881.278714286\t2038.59965809119\tNA\n+M598T4028\t4444.59999999997\tNA\tNA\tNA\n+M598T3699\t21485.7955370261\t262005.048159626\t19097.5044463651\t65022.6200000009\n+M599T3177\tNA\t2838.6076666667\tNA\t2164.39499999999\n+M599T2937\t5756.0700000002\t2557.20999999994\tNA\tNA\n+M599T2856\tNA\tNA\tNA\t3757.56499999997\n+M599T4215\tNA\tNA\t223.746998896449\t881.094999999979\n+M599T2629\t24597.6221333335\t2237.94999999995\tNA\tNA\n+M599T3087\tNA\tNA\tNA\t5745.1150000002\n+M599T4161\t1370.93999999996\t8068.71033333314\t769.979999999982\t5823.36499999964\n+M599T3725\t22852.13\t8220.94499999976\t55051.9369090907\t5147.52485215715\n+M599T3846\tNA\tNA\t5238.05500000018\tNA\n+M599T2884\t3106.80875823442\t3346.02697646104\t6247.93574186843\t3269.28500000002\n+M599T4044\t798856.451441174\t40631.1229886866\t897286.005000001\t49928.058471187\n+M599T3927\t30327.7176201825\t15327.6100000005\t37243.228306282\t2271127.28014893\n+M599T3455\t35592.8437333336\tNA\tNA\t6298.5500000001\n+M599T3741\t2472129.73859259\tNA\t797434.72246154\t11463.9494185942\n+M599T4124\tNA\t354770.43852632\tNA\t707685.939314288\n+M599T3815\t79625.6350000028\t12736.5568357265\t9836.53488085646\t4148.40133264837\n+M599T2837\t2295.00891881701\t5415.77280000005\t3360.28178830047\tNA\n+M599T2984\tNA\t5849.96999999986\tNA\t12502.4521250001\n+M600T3611\tNA\tNA\tNA\t5001.73999999997\n+M600T2728\tNA\tNA\tNA\t1760.62499999996\n+M600T3412\t3353.79499999992\t1647.14222654083\tNA\tNA\n+M600T2595\t6022.12000000021\t2928.1150000001\t3209.81499999996\t2331.63714285708\n+M600T2782\tNA\tNA\tNA\t2608.74386666668\n+M600T3443\t4671.52499999989\t8782.77999999999\tNA\tNA\n+M600T2847\tNA\tNA\t1920.87948584165\t1987.33833333336\n+M600T2963\t9211.58999999991\tNA\t2898.28401908135\t5635.39317528082\n+M600T3070\t2056.42917003541\tNA\t1251.24642515022\t3430.04160000003\n+M600T3353\t2327.15500000008\tNA\tNA\t2315.33491939702\n+M600T3212\t7784.31000000027\t2982.99537949918\tNA\tNA\n+M600T2573\tNA\t3639.80233333338\tNA\tNA\n'
b
diff -r 6139bfcc95cb -r a2c49996603e test-data/faahKO.xset.group.retcor.group.fillPeaks.annotate.variableMetadata-xcms3.tsv
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/faahKO.xset.group.retcor.group.fillPeaks.annotate.variableMetadata-xcms3.tsv Wed Sep 19 03:22:21 2018 -0400
b
b'@@ -0,0 +1,6333 @@\n+name\tnamecustom\tmz\tmzmin\tmzmax\trt\trtmin\trtmax\tnpeaks\tKO\tWT\tisotopes\tadduct\tpcgroup\tfold\ttstat\tpvalue\n+M200T3488\tM200.2T58.1\t200.199996948242\t200.199996948242\t200.199996948242\t58.1321762456397\t58.1321762456397\t58.1321762456397\t1\t1\t0\t\t\t57\tInf\t-1.7513569524204\t0.330286162173178\n+M200T3544\tM200.2T59.1\t200.199996948242\t200.199996948242\t200.199996948242\t59.0661780101845\t59.0661780101845\t59.0661780101845\t1\t0\t1\t\t\t43\tInf\t1\t0.5\n+M200T3882\tM200.2T64.7\t200.199996948242\t200.199996948242\t200.199996948242\t64.698398282013\t64.698398282013\t64.698398282013\t1\t1\t0\t\t\t40\t1.97290277754307\t-0.904937147348527\t0.501383398008076\n+M200T3822\tM200.3T63.7\t200.300003051758\t200.300003051758\t200.300003051758\t63.7046880968321\t63.7046880968321\t63.7046880968321\t1\t0\t1\t\t\t54\tInf\t1\t0.5\n+M200T3184\tM200.2T53.1\t200.199996948242\t200.199996948242\t200.199996948242\t53.0728560746411\t53.0728560746411\t53.0728560746411\t1\t1\t0\t\t\t264\tInf\t-1\t0.5\n+M200T2867\tM200.1T47.8\t200.100006103516\t200.100006103516\t200.100006103516\t47.7790235708651\t47.7790235708651\t47.7790235708651\t1\t1\t0\t\t\t100\t4.4698141311741\t-1.0524209870363\t0.482982331863989\n+M200T3978\tM200.3T66.3\t200.300003051758\t200.300003051758\t200.300003051758\t66.3048296377041\t66.3048296377041\t66.3048296377041\t1\t0\t1\t\t\t35\t1.31236271033293\t0.27205674274799\t0.816918560970057\n+M200T3681\tM200.2T61.3\t200.199996948242\t200.199996948242\t200.199996948242\t61.3457706352117\t61.3457706352117\t61.3457706352117\t1\t0\t1\t\t\t4\t10.3543771314889\t1.68235970874458\t0.329734712515892\n+M201T3788\tM200.8T63.1\t200.800003051758\t200.800003051758\t200.800003051758\t63.1321481019117\t63.1321481019117\t63.1321481019117\t1\t0\t1\t\t\t15\tInf\t1\t0.5\n+M201T3536\tM200.8T58.9\t200.800003051758\t200.800003051758\t200.800003051758\t58.938751448366\t58.938751448366\t58.938751448366\t1\t1\t0\t\t\t96\t2.5883996979851\t-0.572426468523677\t0.650544858486993\n+M201T4052\tM200.8T67.5\t200.800003051758\t200.800003051758\t200.800003051758\t67.5362412825201\t67.5362412825201\t67.5362412825201\t1\t0\t1\t\t\t167\t3.87211055376504\t2.27880620308441\t0.251761446266858\n+M201T3276\tM200.9T54.6\t200.900009155273\t200.900009155273\t200.900009155273\t54.5942327450685\t54.5942327450685\t54.5942327450685\t1\t0\t1\t\t\t177\t3.41695908544214\t1.06174004733468\t0.433294586472355\n+M201T3630\tM200.9T60.5\t200.900009155273\t200.900009155273\t200.900009155273\t60.5002876252059\t60.5002876252059\t60.5002876252059\t1\t0\t1\t\t\t26\tInf\t6.24298356051674\t0.101114694965289\n+M201T2590\tM200.9T43.2\t200.900009155273\t200.900009155273\t200.900009155273\t43.1736546857338\t43.1736546857338\t43.1736546857338\t1\t1\t0\t\t\t103\tInf\t-1\t0.5\n+M201T3985\tM201.2T66.4\t201.199996948242\t201.199996948242\t201.199996948242\t66.411406001538\t66.411406001538\t66.411406001538\t1\t0\t1\t\t\t92\tInf\t1\t0.5\n+M201T3138\tM201.1T52.3\t201.100006103516\t201.100006103516\t201.100006103516\t52.3056203603477\t52.3056203603477\t52.3056203603477\t1\t1\t0\t\t\t61\t3.09544174179034\t-0.833206764478803\t0.552561803395029\n+M201T3693\tM201.2T61.5\t201.199996948242\t201.199996948242\t201.199996948242\t61.5475344176676\t61.5475344176676\t61.5475344176676\t1\t0\t1\t\t\t31\tInf\t1\t0.5\n+M201T3511\tM201.2T58.5\t201.199996948242\t201.199996948242\t201.199996948242\t58.5225240947454\t58.5225240947454\t58.5225240947454\t1\t1\t0\t\t\t121\tInf\t-1.41704674811825\t0.391226133425838\n+M201T3790\tM201.2T63.2\t201.199996948242\t201.199996948242\t201.199996948242\t63.1736091131617\t63.1736091131617\t63.1736091131617\t1\t0\t1\t\t\t15\t2.07453656034107\t0.507201208330113\t0.695739310828326\n+M201T4047\tM201.25T67.4\t201.25\t201.199996948242\t201.300003051758\t67.4492633986824\t67.3066226932934\t67.5919041040714\t2\t0\t2\t\t\t87\tInf\t45.5714285713783\t0.0139674712868434\n+M201T3925\tM201.3T65.4\t201.300003051758\t201.300003051758\t201.300003051758\t65.422966745722\t65.422966745722\t65.422966745722\t1\t0\t1\t\t\t58\t9.99174595875517\t0.895443945668806\t0.532706068985532\n+M201T3874\tM201.3T64.6\t201.300003051758\t201.300003051758\t201.300003051758\t64.5619885311317\t64.5619885311317\t64.5619885311317\t1\t1\t0\t\t\t29\t4.04340073601904\t-2.56216600948778\t0.1456909572893\n+M201T4152\tM201.3T69.2\t201.300003051758\t201.300003051758\t201.300003051758\t69.1984237220222\t69.1'..b'99T3087\tM598.9T51.4\t598.900024414062\t598.900024414062\t598.900024414062\t51.4438227605278\t51.4438227605278\t51.4438227605278\t1\t0\t1\t\t\t110\tInf\t1\t0.5\n+M599T4161\tM599.2T69.4\t599.200012207031\t599\t599.200012207031\t69.3500898109327\t68.8093719359742\t69.8589802263926\t5\t2\t2\t\t\t215\t1.43169367496068\t-0.339238350653957\t0.768870407317004\n+M599T3725\tM599.2T62.1\t599.200012207031\t599.100036621094\t599.200012207031\t62.0761426569479\t61.8008354247033\t62.3094745876028\t3\t2\t1\t\t\t72\t1.93735128439166\t0.560068773616512\t0.66359437643376\n+M599T3846\tM599.2T64.1\t599.200012207031\t599.200012207031\t599.200012207031\t64.0933357932784\t64.0933357932784\t64.0933357932784\t1\t0\t1\t\t\t99\tInf\t1\t0.5\n+M599T2884\tM599.1T48.1\t599.100036621094\t599.100036621094\t599.100036621094\t48.0672044199483\t48.0672044199483\t48.0672044199483\t1\t0\t1\t\t\t101\t1.47488966605743\t1.02548113756989\t0.490263577239904\n+M599T4044\tM599.3T67.4\t599.299987792969\t599.299987792969\t599.299987792969\t67.3957540156441\t67.3957540156441\t67.7683479314032\t3\t1\t1\t\t\t87\t1.12832410189572\t0.0947406778870986\t0.933251618463749\n+M599T3927\tM599.3T65.5\t599.299987792969\t599.200012207031\t599.299987792969\t65.451842766625\t65.451842766625\t66.220552995459\t3\t1\t1\t\t\t58\t50.5608135738083\t1.0128834436489\t0.495913309520612\n+M599T3455\tM599.25T57.6\t599.25\t599.200012207031\t599.299987792969\t57.5855047731195\t57.2690477622402\t57.9019617839989\t2\t1\t1\t\t\t20\t5.65095835284835\t-0.810447012442995\t0.559848676464547\n+M599T3741\tM599.3T62.3\t599.299987792969\t599.299987792969\t599.400024414062\t62.348509589886\t62.2219787379399\t62.4788197691186\t3\t1\t1\t\t\t16\t3.05616738478082\t-0.641167736892023\t0.621912900272668\n+M599T4124\tM599.35T68.7\t599.350006103516\t599.299987792969\t599.400024414062\t68.7305879782284\t68.715687458105\t68.7454884983518\t2\t1\t1\t\t\t39\t1.99477144221469\t0.445807441101457\t0.712469302889833\n+M599T3815\tM599.4T63.6\t599.400024414062\t599.400024414062\t599.400024414062\t63.5902908856904\t63.5902908856904\t63.5902908856904\t1\t1\t0\t\t\t41\t6.60440565660484\t-1.16753575510039\t0.448753436334749\n+M599T2837\tM599.4T47.3\t599.400024414062\t599.400024414062\t599.400024414062\t47.2888127993862\t47.2888127993862\t47.2888127993862\t1\t1\t0\t\t\t244\t2.2946830666594\t-0.948664623432954\t0.443388709573613\n+M599T2984\tM599.45T49.7\t599.450012207031\t599.400024414062\t599.5\t49.739418739903\t49.4905254294192\t49.9883120503868\t2\t1\t1\t\t\t98\t2.13718226332791\t0.481945679553849\t0.693702964615594\n+M600T3611\tM599.5T60.2\t599.5\t599.5\t599.5\t60.1878469955219\t60.1878469955219\t60.1878469955219\t1\t0\t1\t\t\t56\tInf\t1\t0.5\n+M600T2728\tM599.5T45.5\t599.5\t599.5\t599.5\t45.4737557720592\t45.4737557720592\t45.4737557720592\t1\t0\t1\t\t\t133\tInf\t1\t0.5\n+M600T3412\tM599.5T56.9\t599.5\t599.5\t599.5\t56.8749523125544\t56.8749523125544\t56.8749523125544\t1\t1\t0\t\t\t11\tInf\t-2.93026051011228\t0.209367314591815\n+M600T2595\tM599.85T43.2\t599.850006103516\t599.799987792969\t599.900024414062\t43.2457504457672\t42.6188210918425\t43.351428016733\t4\t2\t2\t\t\t103\t1.61514252388468\t-1.05987281413676\t0.461534422135859\n+M600T2782\tM599.9T46.4\t599.900024414062\t599.900024414062\t599.900024414062\t46.3738396121667\t46.3738396121667\t46.3738396121667\t1\t0\t1\t\t\t81\tInf\t1\t0.5\n+M600T3443\tM600T57.4\t600\t600\t600\t57.3889349172675\t57.142677178743\t57.6351926557919\t2\t2\t0\t\t\t44\tInf\t-3.27255424438512\t0.188796332092803\n+M600T2847\tM600T47.4\t600\t600\t600\t47.4459578007975\t47.4459578007975\t47.4459578007975\t1\t0\t1\t\t\t88\tInf\t58.8065843251723\t0.0108246113469044\n+M600T2963\tM600T49.4\t600\t600\t600\t49.3869485044819\t49.3869485044819\t49.3869485044819\t1\t1\t0\t\t\t66\t1.07943970579126\t-0.0705450874311902\t0.953814879351253\n+M600T3070\tM600T51.2\t600\t600\t600\t51.1663039005131\t51.1663039005131\t51.1663039005131\t1\t0\t1\t\t\t47\t2.27641588310559\t0.87611948582698\t0.473626221186199\n+M600T3353\tM600T55.9\t600\t600\t600\t55.8821760378752\t55.8821760378752\t55.8821760378752\t1\t1\t0\t\t\t8\t1.00510512777397\t-0.00360066783278058\t0.997453959117508\n+M600T3212\tM600T53.5\t600\t600\t600\t53.5377338145437\t53.5377338145437\t53.5377338145437\t1\t1\t0\t\t\t90\tInf\t-2.24257442607994\t0.267031787366339\n+M600T2573\tM600T42.9\t600\t600\t600\t42.8867546857338\t42.8867546857338\t42.8867546857338\t1\t1\t0\t\t\t183\tInf\t-1\t0.5\n'