XCMS包裝 - 保留時間

Question

是否有一種簡單的方法可以從xcmsRaw對象中獲取保留時間列表/數組？

示例代碼：

xraw <- xcmsRaw(cdfFile) 

因此，例如，從中獲取信息：

xraw@env$intensity 

或

xraw@env$mz 

Answer 1

你可以看到槽在您的xcmsRaw實例可用

> slotNames(xraw) 
[1] "env"     "tic"     "scantime"    
[4] "scanindex"    "polarity"    "acquisitionNum"  
[7] "profmethod"   "profparam"    "mzrange"    
[10] "gradient"    "msnScanindex"   "msnAcquisitionNum"  
[13] "msnPrecursorScan"  "msnLevel"    "msnRt"     
[16] "msnPrecursorMz"  "msnPrecursorIntensity" "msnPrecursorCharge" 
[19] "msnCollisionEnergy" "filepath" 

你想要的是xraw@msnRt - 它是vector的numeric。

的env插槽是存儲3個變量環境：在class?xcmsRaw類本身

> ls(xraw@env) 
[1] "intensity" "mz"  "profile" 

更多細節。

編輯：如果指定includeMSn = TRUE和輸入文件必須在mzXML或mzML，不cdf的msnRt插槽僅填充;如果從?xcmasRaw使用faahKO例子，你會看到，

xr <- xcmsRaw(cdffiles[1], includeMSn = TRUE) 
Warning message: 
In xcmsRaw(cdffiles[1], includeMSn = TRUE) : 
    Reading of MSn spectra for NetCDF not supported 

此外，xr@msnRt只會存儲的保留時間爲MSN掃描，以n > 1。請參閱xset@rt，其中xset是xcmsSet實例，用於xcms提供的原始/更正的MS1保留時間。

EDIT2：另外，具有與mzR包

> library(mzR) 
> cdffiles[1] 
[2] "/home/lgatto/R/x86_64-unknown-linux-gnu-library/2.16/faahKO/cdf/KO/ko15.CDF" 
> xx <- openMSfile(cdffiles[1]) 
> xx 
Mass Spectrometry file handle. 
Filename: /home/lgatto/R/x86_64-unknown-linux-gnu-library/2.16/faahKO/cdf/KO/ko15.CDF 
Number of scans: 1278 
> hd <- header(xx) 
> names(hd) 
[1] "seqNum"     "acquisitionNum"   
[3] "msLevel"     "peaksCount"    
[5] "totIonCurrent"   "retentionTime"   
[7] "basePeakMZ"    "basePeakIntensity"  
[9] "collisionEnergy"   "ionisationEnergy"   
[11] "highMZ"     "precursorScanNum"   
[13] "precursorMZ"    "precursorCharge"   
[15] "precursorIntensity"  "mergedScan"    
[17] "mergedResultScanNum"  "mergedResultStartScanNum" 
[19] "mergedResultEndScanNum" 
> class(hd) 
[1] "data.frame" 
> dim(hd) 
[1] 1278 19 

一去，但你會默認xcms管道之外，如果你走這條路線（雖然斯特芬·諾伊曼，從xcms，不知道mzR很好，很明顯）。

最後，如果您想最大限度地獲得來自xcms開發人員的反饋，您最好使用xcms online forum的Bioconductor mailing list。

希望這會有所幫助。

Answer 2

很好的答案，但我一直在尋找這樣的：

xraw <- xcmsRaw(cdfFile) 
dp_index <- which(duplicated(rawMat(xraw)[,1])) 
xraw_rt_dp <- rawMat(xraw)[,1] 
xrawData.rt <- xraw_rt_dp[-dp_index] 

現在：

xrawData.rt #contains the retention time. 

觀察 - >使用MZR包：

nc  <- mzR:::netCDFOpen(cdfFile) 
ncData <- mzR:::netCDFRawData(nc) 
mzR:::netCDFClose(nc) 

ncData$rt #contains the same retention time for the same input !!!