Description Usage Arguments Details Constructor Coercion Accessors Author(s) See Also Examples
The ProteinGroup class is a container for identified peptides and proteins, and groups them to distinguish proteins with specific peptides.
1 2 3 4 | ProteinGroup(from,template=NULL,proteinInfo=data.frame())
protein.ac(x, protein.g)
protein.g(x, pattern, variables=c("AC","name"), ...)
|
from |
|
template |
'template' ProteinGroup object for grouping. |
x |
ProteinGroup object |
protein |
character string |
proteinInfo |
data.frame for proteinInfo slot |
protein.g |
character string, denoting a 'protein group'. |
pattern |
character string, see |
variables |
|
... |
Passed on to |
The ProteinGroup class stores spectrum to peptide to protein mapping.
The proteins are grouped by their evidence, i. e. peptides:
Peptides with changes only from Leucin to Isoleucin are considered the same, as they cannot be distinguished by MS.
Proteins which are detected with the same peptides are grouped together to a 'indistinguishable protein'- normally these are splice variants.
Proteins with specific peptides are 'reporters'.
Proteins with no specific peptides are grouped under these 'reporters.
This information is stored in six slots:
spectra.n.peptides
a named 'character' vector, names being spectrum identifier and values are peptides.
peptide.n.proteins
a 'data.frame' containing the number of proteins the peptides could derive from.
peptide.n.protein
a character 'matrix' linking peptides to proteins.
indistinguishable.proteins
a 'matrix' contain.
ProteinGroup(tbl.prot.pep,template=NULL)
:
Creates a ProteinGroup object.
tbl.prot.pep
A 'data.frame' with three columns: 1. Protein, 2. Peptide, 3. Spectrum.
template
Optional ProteinGroup object the grouping is based upon.
In the code snippets below, x
is a ProteinGroup object.
as(from, "ProteinGroup")
:Creates a ProteinGroup object from a data.frame.
as.data.frame(x, row.names = NULL, optional = FALSE)
:Creates a data.frame with columns protein
(character),
peptide
(character), spectrum
.
as.concise.data.frame(from)
:Creates a 'concise' data.frame with one spectrum per row, and protein ACs combined
In the following code snippets, x
is a ProteinGroup object.
spectrumToPeptide(x)
:Gets spectrum to peptide assignment.
peptideInfo(x)
:Peptide information such as protein start position.
peptideSpecificity(x)
:Gets a 'data.frame' containing the peptide specificity: they can be reporter-specific, group-specific, or non-specific.
peptideNProtein(x)
:Gets peptide to protein assignment.
indistinguishableProteins(x)
:Gets the proteins which cannot be distinguished based on peptide evidence.
proteinGroupTable
:Gets the protein grouping, listing reporters and group members.
peptides(x,protein=NULL,specificity=c("reporter-specific",
"group-specific","unspecific"),columns="peptide",set=union)
:Gets all peptides detected, or just those for a protein with the
defined specificity. columns
might define multiple columns of
peptideSpecificity(x)
. set=union returns the union of
peptides of all proteins defined, set=intersect returns the intersection.
Florian P. Breitwieser
IBSpectra
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 | tbl <- data.frame(spectrum=1:14,peptide=c(rep(letters[1:3],4),"a","x"),
modif=":",start.pos=1,
protein=c(rep(c("A","B"),each=6),"C","D"))
pg <- ProteinGroup(tbl)
pg
proteinGroupTable(pg)
data(ibspiked_set1)
pg <- proteinGroup(ibspiked_set1)
ceru.proteins <- protein.g(pg,"CERU")
## all ceru peptides
peptides(pg,ceru.proteins)
## peptides shared by all ceru proteins
peptides(pg,ceru.proteins, set=intersect)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.