Comments (4)
install.packages('http://fgcz-ms.uzh.ch/~cpanse/rawR_0.1.3.tar.gz', repo=NULL)
cp@lilith:~/Downloads/data > cat rawR_example.R | R --no-save
R version 4.0.1 (2020-06-06) -- "See Things Now"
Copyright (C) 2020 The R Foundation for Statistical Computing
Platform: x86_64-apple-darwin17.0 (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
Natural language support but running in an English locale
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> #R
>
> # MD5 (2020_11_06_DBDI_Benzene_01.raw) = a44cd435c8775835f22fc1f76b553d70
> rawfile <- "2020_11_06_DBDI_Benzene_01.raw"
> spec <- rawR::readSpectrum(rawfile, scan = 1:100)
>
> spec[[1]]
Total Ion Current: 4870946
Scan Low Mass: 50
Scan High Mass: 250
Scan Start Time (Min): 0
Scan Number: 1
Base Peak Intensity: 1437918
Base Peak Mass: 78.04678
Scan Mode: FTMS + p NSI Full ms [50.00-250.00]
AGC: On
Micro Scan Count: 1
Scan Segment: 0
Scan Event: 0
Master Index: 0
Charge State: 1
Monoisotopic M/Z: 78.0468
Ion Injection Time (ms): 100.000
FT Resolution: 30000
MS2 Isolation Width: 0.0
Conversion Parameter B: 47557789.235
Conversion Parameter C: -2547049.695
> names(spec[[1]])
[1] "scan" "basePeak"
[3] "TIC" "massRange"
[5] "scanType" "rtinseconds"
[7] "pepmass" "centroidStream"
[9] "HasCentroidStream" "centroid.mZ"
[11] "centroid.intensity" "title"
[13] "charge" "monoisotopicMz"
[15] "mZ" "intensity"
[17] "AGC:" "Micro Scan Count:"
[19] "Ion Injection Time (ms):" "Scan Segment:"
[21] "Scan Event:" "Master Index:"
[23] "Elapsed Scan Time (sec):" "API Source CID Energy:"
[25] "Average Scan by Inst:" "Charge State:"
[27] "Monoisotopic M/Z:" "MS2 Isolation Width:"
[29] "MS3 Isolation Width:" "MS4 Isolation Width:"
[31] "MS5 Isolation Width:" "MS6 Isolation Width:"
[33] "MS7 Isolation Width:" "MS8 Isolation Width:"
[35] "MS9 Isolation Width:" "MS10 Isolation Width:"
[37] "FT Analyzer Settings:" "FT Analyzer Message:"
[39] "FT Resolution:" "Conversion Parameter I:"
[41] "Conversion Parameter A:" "Conversion Parameter B:"
[43] "Conversion Parameter C:" "Conversion Parameter D:"
[45] "Conversion Parameter E:"
>
> jpeg("Rplot.jpg")
> plot(spec[[1]]) #intensities are displayed negative
> dev.off()
null device
1
>
> spec[[1]]$"Ion Injection Time (ms):" #can't access list element by name
[1] "100.000"
>
> spec[[1]]$"Monoisotopic M/Z:"
[1] "78.0468"
>
> sessionInfo()
R version 4.0.1 (2020-06-06)
Platform: x86_64-apple-darwin17.0 (64-bit)
Running under: macOS 10.16
Matrix products: default
BLAS: /Library/Frameworks/R.framework/Versions/4.0/Resources/lib/libRblas.dylib
LAPACK: /Library/Frameworks/R.framework/Versions/4.0/Resources/lib/libRlapack.dylib
locale:
[1] C/UTF-8/C/C/C/C
attached base packages:
[1] stats graphics grDevices utils datasets methods base
loaded via a namespace (and not attached):
[1] compiler_4.0.1 rawR_0.1.3
>
cp@lilith:~/Downloads/data >
from rawrr.
So long story short: It was a bug in our code, not in the API or how it reports missing values. Shame on us! 😉 But it's fixed now! THX again @rokaempf for making us aware of the problem.
from rawrr.
Hi @rokaempf,
I would not call it a "wrong reading", because it is most likely a crude way of reporting a missing value or a null value. And our code simply didn't expect this to happen. So maybe we should "translate" this -1
to an NA
which indicates missing values in R. Or we keep the -1
(the conservative way, since that is what the API returned) and have a check in the plotting function that throws an error. Not sure which way to go yet. I would only consider the first, if Thermo confirms that our assumption is correct. I could well be that they (Thermo) are using this mechanism for all/many numeric parameters.
The other "missing parameters" seem to be truly missing. The problem with the raw format is that is has evolved organically, so maybe these key:value pairs are really not part of what the instrument logs to disc. We should check that.
Anyway, thanks for reporting this problem!
from rawrr.
Hi @rokaempf,
we found the issue that creates the -1
value and we hope to fix it soon. I will keep you posted.
Greetings,
Tobi
from rawrr.
Related Issues (20)
- Get information on gradient HOT 1
- Enhancement - Complete readIndex() function HOT 9
- Peak charges for MS1 spectras HOT 4
- Spectrum scan centroid mZ, intensity and noises values do not match HOT 2
- Error in Example: Length of "x" and "y" are not matching HOT 3
- Read noise value for profile mode mass spectra HOT 4
- Read_Spectrum - Sum Spectra
- unit should be minute / auc computation in seconds HOT 24
- validate_rawrrSpectrum 'StartTime' HOT 4
- "Error: line 1 did not have 9 elements" for readIndex() and readChromatogram() + "Error : No scan vector is provided"for readSpectrum HOT 9
- Problem executing readChromatogram inside Singularity container HOT 5
- Add a check if `input` file exists and is not empty
- Error in if (rvs != "No RAW file specified!") { : the condition has length > 1 HOT 16
- Switch to RawFileReader 5.0.93 HOT 2
- different total number of Spectra in msconvert, compomics/ThermoRawFileParser and thermofisherlsms/RawFileReader HOT 2
- Request for auc.rawrrChromatogram HOT 3
- profile mode in readSpectrum HOT 11
- speed up readIndex/readSpectrum by base::textConnection HOT 14
- rawrr::buildRawrrExe() fails HOT 6
- auc.rawrrChromatogram question
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from rawrr.