Description Usage Arguments Value Examples
peptides which do not have protein assignment drop out
1 |
file |
where to write |
... |
fasta list or single file |
writes a file.
1 2 3 4 5 6 7 8 9 10 11 12 | #example how to create a protein db with decoy sequences
library(seqinr)
#library(prozor)
file = system.file("extdata/fgcz_contaminants2021_20210929.fasta.gz",package = "prozor")
fasta = readPeptideFasta(file = file)
revfasta <- reverseSeq(fasta)
decoyDB <- c(fasta,revfasta)
stopifnot(length(decoyDB) == 2 * length(fasta))
## Not run:
writeFasta(decoyDB, file="test.fasta")
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.