康康到底哪里出了问题啊

俺是做代谢组学的,然后要俺把得到的色谱峰根据一个表格批量提取色谱峰!但是今天在处理数据上遇到了一些问题!就通过一个读取一个代谢物表信息表格,然后run下边的代码,跑完后得到一堆pdf图,但是今天run的时候它给报错了。

library(Spectra)
library(stringr)
library(data.table)
library(xcms)
library(openxlsx)
library(dplyr)
library(MsCoreUtils)
library(MSnbase)
dda_data <- MSnbase::readMSData(files = "D:/spectraData/20230406AmideSerum/正A正BullNeg80ACN2.mzXML", mode = "onDisk")

cwp <- CentWaveParam(snthresh = 3, 
                    noise = 1000, 
                    ppm = 30,
                    peakwidth = c(5, 30),
                    mzCenterFun = "wMean",
                    integrate = 1,
                    mzdiff = -0.001)

dda_data <- findChromPeaks(dda_data, param = cwp)  #

##
mass <- read.csv("C:/Users/Virus33/Desktop/parameter.csv")
ions <- data.frame(compound = mass$Metabolites,low = mass$Exact.Mass*(1-0.000010),high = mass$Exact.Mass*(1+0.000010))

mzs=ions[,c(2,3)]
par(mfrow = c(1,1))
chrs <- chromatogram(dda_data,mz = mzs)
plot(chrs[[1]]) 
pdf(file= "./ACN2.pdf")
for(i in 1:length(chrs)){
 plot(chrs[[i]],main = ions$compound[i])
}
dev.off()
getwd()

#但是在run到这一步的时候,出现了错误!昨晚上run这个代码的时候都还可以把几百个代谢物的色谱峰全部提取出来,然后今天就不行了

chrs <- chromatogram(dda_data,mz = mzs)

#然后报错的结果是这样的

Error in .loca1(object,.. .'list' object cannot be coerced to type 'double'

之前跑出来是这样的图片!

img

然后就再也出不来了!救命!