quantFile <- dir(system.file(package = "MSnbase", dir = "extdata"),
full.name = TRUE, pattern = "mzXML$")
identFile <- dir(system.file(package = "MSnbase", dir = "extdata"),
full.name = TRUE, pattern = "dummyiTRAQ.mzid")
msexp <- readMSData(quantFile)
msexp <- addIdentificationData(msexp, identFile)
fData(msexp)$pepseq
length(msexp)
## using default fcol
msexp2 <- removeNoId(msexp)
length(msexp2)
fData(msexp2)$pepseq
## using keep
print(fvarLabels(msexp))
(k <- fData(msexp)$'ms-gf:evalue' > 75)
k[is.na(k)] <- FALSE
k
msexp3 <- removeNoId(msexp, keep = k)
length(msexp3)
fData(msexp3)$pepseq
Run the code above in your browser using DataLab