Compute number of peptides per 1000 aa — compute_num_peptides_per

Computing number of peptides per 1000 amino acids.

Usage

compute_num_peptides_per_1000aa(msnid, path_to_FASTA)

Arguments

msnid: (MSnID object) MS/MS ID data
path_to_FASTA: (character) Path to FASTA file.

Value

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

Examples

if (FALSE) {
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)
}