October 2019
Intermediate to advanced
316 pages
9h 45m
English
Proteomics data can be viewed in a genome browser using the following steps:
library(MSnID)library(EnsDb.Hsapiens.v86)library(rtracklayer)
msnid <- MSnID() msnid <- read_mzIDs(msnid, file.path(getwd(), "datasets", "ch6", "HeLa_180123_m43_r2_CAM.mzid.gz"))
real_hits <- msnid@psms[! msnid@psms$isDecoy, ]
required_info <- real_hits[, c("spectrumID", "pepSeq", "accession", "start", "end")]
uniprot_ids <- unlist(lapply(strsplit(required_info$accession, "\\|"), function(x){x[2]}) ) uniprot_ids <- uniprot_ids[!is.na(uniprot_ids)] ...Read now
Unlock full access