## ----style, echo = FALSE, results = 'asis', message=FALSE------------------ BiocStyle::markdown() ## ----env, echo=FALSE, message=FALSE---------------------------------------- suppressPackageStartupMessages(library("AnnotationHub")) suppressPackageStartupMessages(library("ProteomicsAnnotationHubData")) suppressPackageStartupMessages(library("mzR")) ## ----ahinit---------------------------------------------------------------- library("AnnotationHub") ah <- AnnotationHub() ah ## ----provider-------------------------------------------------------------- query(ah, "PRIDE") ## ----title----------------------------------------------------------------- query(ah, "PXD000001") ## ----ah49008--------------------------------------------------------------- ah["AH49008"] ## ----rawmsdata------------------------------------------------------------- library("mzR") rw <- ah[["AH49008"]] rw ## ----msdataplot------------------------------------------------------------ plot(peaks(rw, 1), type = "h", xlab = "M/Z", ylab = "Intensity") ## ----availablepahd--------------------------------------------------------- library("ProteomicsAnnotationHubData") availableProteomicsAnnotationHubData ## ---- echo=FALSE----------------------------------------------------------- cat(readLines(dir(system.file("extdata", package = "ProteomicsAnnotationHubData"), full.names = TRUE, pattern = "PXD000001.dcf"))[10:30], sep = "\n") ## -------------------------------------------------------------------------- writePahdTemplate() ## ----PXD000001.dcf, echo=FALSE--------------------------------------------- f <- list.files(system.file("extdata", package = "ProteomicsAnnotationHubData"), full.names = TRUE, pattern = "PXD000001.dcf") cat(readLines(f)[10:30], sep = "\n") ## ----providers, echo=FALSE------------------------------------------------- tab <- do.call(rbind, ProteomicsAnnotationHubData:::ProteomicsAnnotationHubDataProviders) knitr::kable(tab, row.names = FALSE) ## ----tags, echo=FALSE------------------------------------------------------ ProteomicsAnnotationHubData:::ProteomicsAnnotationHubDataTags ## ----PXD000001.R, echo=FALSE----------------------------------------------- f <- list.files(system.file("scripts", package = "ProteomicsAnnotationHubData"), full.names = TRUE, pattern = "PXD000001.R") cat(readLines(f)[1:4], sep = "\n") ## ---- eval=FALSE----------------------------------------------------------- # setClass("mzRpwizResource", contains="AnnotationHubResource") # setMethod(".get1", "mzRpwizResource", # function(x, ...) # { # .require("mzR") # yy <- cache(.hub(x)) # mzR::openMSfile(yy, backend = "pwiz") # }) ## ---- eval=FALSE----------------------------------------------------------- # setClass("mzRidentResource", contains="AnnotationHubResource") # setMethod(".get1", "mzRidentResource", # function(x, ...) # { # .require("mzR") # yy <- cache(.hub(x)) # mzR::openIDfile(yy) # }) ## ---- eval=FALSE----------------------------------------------------------- # setClass("MSnSetResource", contains="RdaResource") # setMethod(".get1", "MSnSetResource", # function(x, ...) # { # .require("MSnbase") # callNextMethod(x, ...) # }) ## ---- eval=FALSE----------------------------------------------------------- # setClass("AAStringSetResource", contains="AnnotationHubResource") # setMethod(".get1", "AAStringSetResource", # function(x, ...) # { # .require("Biostrings") # yy <- cache(.hub(x)) # Biostrings::readAAStringSet(yy) # })