compute_num_peptides_per_1000aa: Computing number of peptides per 1000 aa
In vladpetyuk/PlexedPiper: Pipeline for isobaric quantification.

Description Usage Arguments Value Examples

View source: R/compute_accession_coverage.R

Computing number of peptides per 1000 aa

1	compute_num_peptides_per_1000aa(msnid, path_to_FASTA)

`msnid`	(MSnID object) MS/MS ID data
`path_to_FASTA`	(numeric) Maximum acceptable FDR rate. Default is 0.01.

(MSnID object) MS/MS ID data with computed number of peptides per 1000 aa. Added column name - "peptides_per_1000aa".

path_to_MSGF_results <- system.file("extdata/global/msgf_output", package = "PlexedPiperTestData")
msnid <- read_msgf_data(path_to_MSGF_results)
path_to_FASTA <- system.file("extdata/Rattus_norvegicus_NCBI_RefSeq_2018-04-10.fasta.gz", package = "PlexedPiperTestData")
msnid <- compute_num_peptides_per_1000aa(msnid, path_to_FASTA)
hist(msnid$peptides_per_1000aa)