Compute FDR and add it to RefSpectra Table
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 | library(bibliospec)
library(prozor)
dbfile <- file.path(path.package("bibliospec"),
"extdata/peptideStd.sqlite")
#call constructor
BS <- Bibliospec(dbfile=dbfile)
BS$listTables()
BS$annotateDB(fasta=file.path(path.package("bibliospec"),"extdata/uniprot.fasta.gz" ))
(specMeta <- BS$getSpectraMetaWithAnnotation())
specMeta <- BS$getSpectraMeta()
dim(specMeta)
specMeta
fdr <- prozor::computeFDR(specMeta$score, specMeta$proteinID,larger_better=FALSE)
head(fdr)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.