    ## <script type="text/javascript">
    ## document.addEventListener("DOMContentLoaded", function() {
    ##   document.querySelector("h1").style.marginTop = "0";
    ## });
    ## </script>
    ## <script type="text/javascript">
    ## document.addEventListener("DOMContentLoaded", function() {
    ##   var links = document.links;
    ##   for (var i = 0, linksLength = links.length; i < linksLength; i++)
    ##     if (links[i].hostname != window.location.hostname)
    ##       links[i].target = '_blank';
    ## });
    ## </script>

# Pbase example data

Laurent Gatto

The original data is a 10 fmol Peptide Retention Time Calibration Mixture spiked into 50 ng HeLa background acquired on a Thermo Orbitrap Q Exactive instrument. A restricted set of high scoring human proteins from the UniProt release 2014_06 were searched using the MSGF+ engine.

## The fasta database

    library("Biostrings")
    fafile <- system.file("extdata/HUMAN_2014_06-09prots.fasta", package = "Pbase")
    fa <- readAAStringSet(fafile)
    fa

    ##   A AAStringSet instance of length 9
    ##     width seq                                          names
    ## [1]   375 MDDDIAALVVDNGSGMCKAGF...WISKQEYDESGPSIVHRKCF sp|P60709|ACTB_HU...
    ## [2]   664 METPSQRRATRSGAQASSTPL...SYLLGNSSPRTQSPQNCSIM sp|P02545|LMNA_HU...
    ## [3]   417 MSLSNKLTLDKLDVKGKRVVM...ASLELLEGKVLPGVDALSNI sp|P00558|PGK1_HU...
    ## [4]  2602 MPVTEKDLAEDAPWKKIQQNT...LAVKWGEEHIPGSPFHVTVP sp|O75369|FLNB_HU...
    ## [5]   911 MVDYHAANQSYQYGPSSAGNG...VPGALDYKSFSTALYGESDL sp|O43707|ACTN4_H...
    ## [6]  2624 MFRRARLSVKPNVRPGVGARG...ATTVSEYFFNDIFIEVDETE sp|A6H8Y1|BDP1_HU...
    ## [7]  3374 MSPESGHSRIFEATAGPNKPE...TLSKDSLSNGVPSGRQAEFS sp|A4UGR9|XIRP2_H...
    ## [8]   364 MPYQYPALTPEQKKELSDIAH...PSGQAGAAASESLFVSNHAY sp|P04075|ALDOA_H...
    ## [9]   418 MARRKPEGSSFNMTHLSMAMA...PSGQAGAAASESLFVSNHAY sp|P04075-2|ALDOA...

## The PSM data

    library("mzID")
    idfile <- system.file("extdata/Thermo_Hela_PRTC_1_selected.mzid", package = "Pbase")
    id <- flatten(mzID(idfile))

    ## reading Thermo_Hela_PRTC_1_selected.mzid... DONE!

    dim(id)

    ## [1] 137  41

    head(id)

    ##     spectrumid scan number(s)
    ## 1    index=173          12256
    ## 1.1  index=173          12256
    ## 2    index=163          11860
    ## 2.1  index=163          11860
    ## 3    index=200          13408
    ## 3.1  index=200          13408
    ##                                                                                spectrum title
    ## 1   msLevel 2; retentionTime 2094.56706; scanNum 12256; precMz 1137.06665029649; precCharge 2
    ## 1.1 msLevel 2; retentionTime 2094.56706; scanNum 12256; precMz 1137.06665029649; precCharge 2
    ## 2   msLevel 2; retentionTime 2039.84424; scanNum 11860; precMz 1136.57450195803; precCharge 2
    ## 2.1 msLevel 2; retentionTime 2039.84424; scanNum 11860; precMz 1136.57450195803; precCharge 2
    ## 3   msLevel 2; retentionTime 2258.27868; scanNum 13408; precMz 703.038108542133; precCharge 3
    ## 3.1 msLevel 2; retentionTime 2258.27868; scanNum 13408; precMz 703.038108542133; precCharge 3
    ##     acquisitionnum passthreshold rank calculatedmasstocharge
    ## 1              173          TRUE    1               1136.574
    ## 1.1            173          TRUE    1               1136.574
    ## 2              163          TRUE    1               1136.574
    ## 2.1            163          TRUE    1               1136.574
    ## 3              200          TRUE    1                703.037
    ## 3.1            200          TRUE    1                703.037
    ##     experimentalmasstocharge chargestate ms-gf:denovoscore
    ## 1                  1137.0667           2               132
    ## 1.1                1137.0667           2               132
    ## 2                  1136.5745           2               230
    ## 2.1                1136.5745           2               230
    ## 3                   703.0381           3               145
    ## 3.1                 703.0381           3               145
    ##                   ms-gf:evalue ms-gf:rawscore
    ## 1   0.000000000000000002597097            118
    ## 1.1 0.000000000000000002597097            118
    ## 2   0.000000000000000049426636            186
    ## 2.1 0.000000000000000049426636            186
    ## 3   0.000000000408042900000000             98
    ## 3.1 0.000000000408042900000000             98
    ##                   ms-gf:specevalue assumeddissociationmethod
    ## 1   0.0000000000000000000002276758                       CID
    ## 1.1 0.0000000000000000000002276758                       CID
    ## 2   0.0000000000000000000043330090                       CID
    ## 2.1 0.0000000000000000000043330090                       CID
    ## 3   0.0000000000000357806840000000                       CID
    ## 3.1 0.0000000000000357806840000000                       CID
    ##     ctermioncurrentratio explainedioncurrentratio isotopeerror
    ## 1              0.5061673                0.5888804            1
    ## 1.1            0.5061673                0.5888804            1
    ## 2             0.48202217               0.57532895            0
    ## 2.1           0.48202217               0.57532895            0
    ## 3             0.32245573                0.5909792            0
    ## 3.1           0.32245573                0.5909792            0
    ##     meanerrorall meanerrortop7 meanrelerrorall meanrelerrortop7
    ## 1      17.757027     4.3787823     -15.7646885        -3.255431
    ## 1.1    17.757027     4.3787823     -15.7646885        -3.255431
    ## 2      18.029726     5.0083427       -16.59046        -3.333969
    ## 2.1    18.029726     5.0083427       -16.59046        -3.333969
    ## 3      22.256987      8.589616       -17.53643       -1.7168152
    ## 3.1    22.256987      8.589616       -17.53643       -1.7168152
    ##     ms2ioncurrent ntermioncurrentratio nummatchedmainions stdeverrorall
    ## 1        526864.0           0.08271315                 25     23.135044
    ## 1.1      526864.0           0.08271315                 25     23.135044
    ## 2       3342169.2           0.09330681                 30     28.647552
    ## 2.1     3342169.2           0.09330681                 30     28.647552
    ## 3       5008038.0            0.2685235                 22      38.71648
    ## 3.1     5008038.0            0.2685235                 22      38.71648
    ##     stdeverrortop7 stdevrelerrorall stdevrelerrortop7 isdecoy post pre end
    ## 1         3.582861         24.53603          4.627396   FALSE    C   K 134
    ## 1.1       3.582861         24.53603          4.627396   FALSE    C   K 188
    ## 2        2.9929223        29.504406          4.788082   FALSE    C   K 134
    ## 2.1      2.9929223        29.504406          4.788082   FALSE    C   K 188
    ## 3        2.4928854        41.070827          8.777729   FALSE    Y   K 173
    ## 3.1      2.4928854        41.070827          8.777729   FALSE    Y   K 227
    ##     start               accession length
    ## 1     112   sp|P04075|ALDOA_HUMAN    364
    ## 1.1   166 sp|P04075-2|ALDOA_HUMAN    418
    ## 2     112   sp|P04075|ALDOA_HUMAN    364
    ## 2.1   166 sp|P04075-2|ALDOA_HUMAN    418
    ## 3     154   sp|P04075|ALDOA_HUMAN    364
    ## 3.1   208 sp|P04075-2|ALDOA_HUMAN    418
    ##                                                                description
    ## 1      Fructose-bisphosphate aldolase A OS=Homo sapiens GN=ALDOA PE=1 SV=2
    ## 1.1 Isoform 2 of Fructose-bisphosphate aldolase A OS=Homo sapiens GN=ALDOA
    ## 2      Fructose-bisphosphate aldolase A OS=Homo sapiens GN=ALDOA PE=1 SV=2
    ## 2.1 Isoform 2 of Fructose-bisphosphate aldolase A OS=Homo sapiens GN=ALDOA
    ## 3      Fructose-bisphosphate aldolase A OS=Homo sapiens GN=ALDOA PE=1 SV=2
    ## 3.1 Isoform 2 of Fructose-bisphosphate aldolase A OS=Homo sapiens GN=ALDOA
    ##                      pepseq modified modification
    ## 1   GVVPLAGTNGETTTQGLDGLSER    FALSE         <NA>
    ## 1.1 GVVPLAGTNGETTTQGLDGLSER    FALSE         <NA>
    ## 2   GVVPLAGTNGETTTQGLDGLSER    FALSE         <NA>
    ## 2.1 GVVPLAGTNGETTTQGLDGLSER    FALSE         <NA>
    ## 3      IGEHTPSALAIMENANVLAR    FALSE         <NA>
    ## 3.1    IGEHTPSALAIMENANVLAR    FALSE         <NA>
    ##                      spectrumFile                databaseFile
    ## 1   Thermo_Hela_PRTC_selected.mgf HUMAN_2014_06-09prots.fasta
    ## 1.1 Thermo_Hela_PRTC_selected.mgf HUMAN_2014_06-09prots.fasta
    ## 2   Thermo_Hela_PRTC_selected.mgf HUMAN_2014_06-09prots.fasta
    ## 2.1 Thermo_Hela_PRTC_selected.mgf HUMAN_2014_06-09prots.fasta
    ## 3   Thermo_Hela_PRTC_selected.mgf HUMAN_2014_06-09prots.fasta
    ## 3.1 Thermo_Hela_PRTC_selected.mgf HUMAN_2014_06-09prots.fasta

## The Proteins object

    library("Pbase")
    p <- Proteins(fafile)
    p <- addIdentificationData(p, idfile)

    ## Reading 1 identification files:
    ##   1. /private/tmp/RtmpXTFOYc/Rinst10c40378b8d32/Pbase/extdata/Thermo_Hela_PRTC_1_selected.mzid
    ## done.

    p

    ## S4 class type     : Proteins
    ## Class version     : 0.1
    ## Created           : Mon Oct 13 21:06:50 2014
    ## Number of Proteins: 9
    ## Sequences:
    ##   [1] A4UGR9 [2] A6H8Y1 ... [8] P04075-2 [9] P60709
    ## Sequence features:
    ##   [1] DB [2] AccessionNumber ... [11] Filename [12] npeps
    ## Peptide features:
    ##   [1] DB [2] AccessionNumber ... [25] DatabaseSeq [26] filenames

A `Proteins` object is composed of a set of protein sequences accessible with the `aa` accessor as well as an optional set of peptides features that are mapped as coordinates along the proteins, available with `pranges`. The actual peptide sequences can be extraced with `pfeatures`.

    aa(p)

    ##   A AAStringSet instance of length 9
    ##     width seq                                          names
    ## [1]  3374 MSPESGHSRIFEATAGPNKPE...TLSKDSLSNGVPSGRQAEFS A4UGR9
    ## [2]  2624 MFRRARLSVKPNVRPGVGARG...ATTVSEYFFNDIFIEVDETE A6H8Y1
    ## [3]   911 MVDYHAANQSYQYGPSSAGNG...VPGALDYKSFSTALYGESDL O43707
    ## [4]  2602 MPVTEKDLAEDAPWKKIQQNT...LAVKWGEEHIPGSPFHVTVP O75369
    ## [5]   417 MSLSNKLTLDKLDVKGKRVVM...ASLELLEGKVLPGVDALSNI P00558
    ## [6]   664 METPSQRRATRSGAQASSTPL...SYLLGNSSPRTQSPQNCSIM P02545
    ## [7]   364 MPYQYPALTPEQKKELSDIAH...PSGQAGAAASESLFVSNHAY P04075
    ## [8]   418 MARRKPEGSSFNMTHLSMAMA...PSGQAGAAASESLFVSNHAY P04075-2
    ## [9]   375 MDDDIAALVVDNGSGMCKAGF...WISKQEYDESGPSIVHRKCF P60709

    pranges(p)

    ## IRangesList of length 9
    ## $A4UGR9
    ## IRanges of length 36
    ##      start  end width  names
    ## [1]   2743 2760    18 A4UGR9
    ## [2]    307  318    12 A4UGR9
    ## [3]   1858 1870    13 A4UGR9
    ## [4]   1699 1708    10 A4UGR9
    ## [5]   2622 2637    16 A4UGR9
    ## ...    ...  ...   ...    ...
    ## [32]    20   31    12 A4UGR9
    ## [33]  1712 1729    18 A4UGR9
    ## [34]    48   61    14 A4UGR9
    ## [35]  2082 2094    13 A4UGR9
    ## [36]  2743 2756    14 A4UGR9
    ##
    ## ...
    ## <8 more elements>

    pfeatures(p)

    ## AAStringSetList of length 9
    ## [["A4UGR9"]] QEITQNKSFFSSVKESQR LPVPKDVYSKQR ... QEITQNKSFFSSVK
    ## [["A6H8Y1"]] EDAEQVALEVDLNQKKRR ETGKENFR ... ARLSVKPNVRPGVGARGSTASNPQRGR
    ## [["O43707"]] QQRKTFTAWCNSHLR CQKICDQWDALGSLTHSR ... VGWEQLLTTIAR
    ## [["O75369"]] DLDIIDNYDYSHTVK PFDLVIPFAVRK ... VQAQGPGLKEAFTNK
    ## [["P00558"]] ELNYFAKALESPER DLMSKAEK QIVWNGPVGVFEWEAFAR FHVEEEGKGKDASGNK GTKALMDEVVK
    ## [["P02545"]] METPSQRRATR DTSRRLLAEKEREMAEMR ... RATRSGAQASSTPLSPTR
    ## [["P04075"]] GVVPLAGTNGETTTQGLDGLSER ... YTPSGQAGAAASESLFVSNHAY
    ## [["P04075-2"]] GVVPLAGTNGETTTQGLDGLSER ... YTPSGQAGAAASESLFVSNHAY
    ## [["P60709"]] DLTDYLMKILTER

A Proteins instance is further described by general `metadata`. Protein sequence and peptide features annotations can be accessed with `ametadata` and `pmetadata` (or `acols` and `pcols`) respectively.

    metadata(p)

    ## $created
    ## [1] "Mon Oct 13 21:06:50 2014"

    head(acols(p))

    ## DataFrame with 6 rows and 12 columns
    ##      DB AccessionNumber   EntryName IsoformName
    ##   <Rle>     <character> <character>       <Rle>
    ## 1    sp          A4UGR9 XIRP2_HUMAN          NA
    ## 2    sp          A6H8Y1  BDP1_HUMAN          NA
    ## 3    sp          O43707 ACTN4_HUMAN          NA
    ## 4    sp          O75369  FLNB_HUMAN          NA
    ## 5    sp          P00558  PGK1_HUMAN          NA
    ## 6    sp          P02545  LMNA_HUMAN          NA
    ##                                         ProteinName OrganismName GeneName
    ##                                         <character>        <Rle>    <Rle>
    ## 1     Xin actin-binding repeat-containing protein 2 Homo sapiens    XIRP2
    ## 2 Transcription factor TFIIIB component B'' homolog Homo sapiens     BDP1
    ## 3                                   Alpha-actinin-4 Homo sapiens    ACTN4
    ## 4                                         Filamin-B Homo sapiens     FLNB
    ## 5                         Phosphoglycerate kinase 1 Homo sapiens     PGK1
    ## 6                                      Prelamin-A/C Homo sapiens     LMNA
    ##            ProteinExistence SequenceVersion Comment
    ##                       <Rle>           <Rle>   <Rle>
    ## 1 Evidence at protein level               2      NA
    ## 2 Evidence at protein level               3      NA
    ## 3 Evidence at protein level               2      NA
    ## 4 Evidence at protein level               2      NA
    ## 5 Evidence at protein level               3      NA
    ## 6 Evidence at protein level               1      NA
    ##                                                                               Filename
    ##                                                                                  <Rle>
    ## 1 /private/tmp/RtmpXTFOYc/Rinst10c40378b8d32/Pbase/extdata/HUMAN_2014_06-09prots.fasta
    ## 2 /private/tmp/RtmpXTFOYc/Rinst10c40378b8d32/Pbase/extdata/HUMAN_2014_06-09prots.fasta
    ## 3 /private/tmp/RtmpXTFOYc/Rinst10c40378b8d32/Pbase/extdata/HUMAN_2014_06-09prots.fasta
    ## 4 /private/tmp/RtmpXTFOYc/Rinst10c40378b8d32/Pbase/extdata/HUMAN_2014_06-09prots.fasta
    ## 5 /private/tmp/RtmpXTFOYc/Rinst10c40378b8d32/Pbase/extdata/HUMAN_2014_06-09prots.fasta
    ## 6 /private/tmp/RtmpXTFOYc/Rinst10c40378b8d32/Pbase/extdata/HUMAN_2014_06-09prots.fasta
    ##       npeps
    ##   <integer>
    ## 1        36
    ## 2        23
    ## 3         6
    ## 4        13
    ## 5         5
    ## 6        12

    head(pcols(p))

    ## SplitDataFrameList of length 6
    ## $A4UGR9
    ## DataFrame with 36 rows and 26 columns
    ##        DB AccessionNumber   EntryName IsoformName
    ##     <Rle>     <character> <character>       <Rle>
    ## 1      sp          A4UGR9 XIRP2_HUMAN          NA
    ## 2      sp          A4UGR9 XIRP2_HUMAN          NA
    ## 3      sp          A4UGR9 XIRP2_HUMAN          NA
    ## 4      sp          A4UGR9 XIRP2_HUMAN          NA
    ## 5      sp          A4UGR9 XIRP2_HUMAN          NA
    ## ...   ...             ...         ...         ...
    ## 32     sp          A4UGR9 XIRP2_HUMAN          NA
    ## 33     sp          A4UGR9 XIRP2_HUMAN          NA
    ## 34     sp          A4UGR9 XIRP2_HUMAN          NA
    ## 35     sp          A4UGR9 XIRP2_HUMAN          NA
    ## 36     sp          A4UGR9 XIRP2_HUMAN          NA
    ##                                                             ProteinName
    ##                                                             <character>
    ## 1   sp|A4UGR9|XIRP2_HUMAN Xin actin-binding repeat-containing protein 2
    ## 2   sp|A4UGR9|XIRP2_HUMAN Xin actin-binding repeat-containing protein 2
    ## 3   sp|A4UGR9|XIRP2_HUMAN Xin actin-binding repeat-containing protein 2
    ## 4   sp|A4UGR9|XIRP2_HUMAN Xin actin-binding repeat-containing protein 2
    ## 5   sp|A4UGR9|XIRP2_HUMAN Xin actin-binding repeat-containing protein 2
    ## ...                                                                 ...
    ## 32  sp|A4UGR9|XIRP2_HUMAN Xin actin-binding repeat-containing protein 2
    ## 33  sp|A4UGR9|XIRP2_HUMAN Xin actin-binding repeat-containing protein 2
    ## 34  sp|A4UGR9|XIRP2_HUMAN Xin actin-binding repeat-containing protein 2
    ## 35  sp|A4UGR9|XIRP2_HUMAN Xin actin-binding repeat-containing protein 2
    ## 36  sp|A4UGR9|XIRP2_HUMAN Xin actin-binding repeat-containing protein 2
    ##     OrganismName GeneName          ProteinExistence SequenceVersion
    ##            <Rle>    <Rle>                     <Rle>           <Rle>
    ## 1   Homo sapiens    XIRP2 Evidence at protein level               2
    ## 2   Homo sapiens    XIRP2 Evidence at protein level               2
    ## 3   Homo sapiens    XIRP2 Evidence at protein level               2
    ## 4   Homo sapiens    XIRP2 Evidence at protein level               2
    ## 5   Homo sapiens    XIRP2 Evidence at protein level               2
    ## ...          ...      ...                       ...             ...
    ## 32  Homo sapiens    XIRP2 Evidence at protein level               2
    ## 33  Homo sapiens    XIRP2 Evidence at protein level               2
    ## 34  Homo sapiens    XIRP2 Evidence at protein level               2
    ## 35  Homo sapiens    XIRP2 Evidence at protein level               2
    ## 36  Homo sapiens    XIRP2 Evidence at protein level               2
    ##     Comment spectrumID chargeState      rank passThreshold
    ##       <Rle>   <factor>   <integer> <integer>     <logical>
    ## 1        NA  index=124           3         1          TRUE
    ## 2        NA   index=28           2         1          TRUE
    ## 3        NA   index=20           2         1          TRUE
    ## 4        NA  index=187           2         1          TRUE
    ## 5        NA  index=211           3         1          TRUE
    ## ...     ...        ...         ...       ...           ...
    ## 32       NA   index=99           2         1          TRUE
    ## 33       NA    index=9           2         1          TRUE
    ## 34       NA  index=122           2         1          TRUE
    ## 35       NA   index=87           2         1          TRUE
    ## 36       NA   index=77           2         1          TRUE
    ##     experimentalMassToCharge calculatedMassToCharge           sequence
    ##                    <numeric>              <numeric>           <factor>
    ## 1                   715.0305               715.0308 QEITQNKSFFSSVKESQR
    ## 2                   715.9177               715.4117       LPVPKDVYSKQR
    ## 3                   786.9066               786.9081      EQNNDALEKSLRR
    ## 4                   629.8380               629.3386         SLKESSHRWK
    ## 5                   645.3429               645.3511   LKMVPRKQREFSGSDR
    ## ...                      ...                    ...                ...
    ## 32                  619.2888               618.7782       PESGFAEDSAAR
    ## 33                 1014.0198              1013.5117 QPDAIPGDIEKAIECLEK
    ## 34                  821.4005               820.8909     MARYQAAVSRGDCR
    ## 35                  720.3445               720.3527      TNTSTGLKMAMER
    ## 36                  821.9231               821.9254     QEITQNKSFFSSVK
    ##        modNum   isDecoy     post      pre     start       end
    ##     <integer> <logical> <factor> <factor> <integer> <integer>
    ## 1           0     FALSE        D        K      2743      2760
    ## 2           0     FALSE        N        R       307       318
    ## 3           0     FALSE        L        R      1858      1870
    ## 4           0     FALSE        E        K      1699      1708
    ## 5           0     FALSE        G        K      2622      2637
    ## ...       ...       ...      ...      ...       ...       ...
    ## 32          0     FALSE        G        K        20        31
    ## 33          1     FALSE        A        K      1712      1729
    ## 34          1     FALSE        S        R        48        61
    ## 35          0     FALSE        S        K      2082      2094
    ## 36          0     FALSE        E        K      2743      2756
    ##            DatabaseAccess DatabaseSeq
    ##                  <factor>    <factor>
    ## 1   sp|A4UGR9|XIRP2_HUMAN
    ## 2   sp|A4UGR9|XIRP2_HUMAN
    ## 3   sp|A4UGR9|XIRP2_HUMAN
    ## 4   sp|A4UGR9|XIRP2_HUMAN
    ## 5   sp|A4UGR9|XIRP2_HUMAN
    ## ...                   ...         ...
    ## 32  sp|A4UGR9|XIRP2_HUMAN
    ## 33  sp|A4UGR9|XIRP2_HUMAN
    ## 34  sp|A4UGR9|XIRP2_HUMAN
    ## 35  sp|A4UGR9|XIRP2_HUMAN
    ## 36  sp|A4UGR9|XIRP2_HUMAN
    ##                                                                                     filenames
    ##                                                                                         <Rle>
    ## 1   /private/tmp/RtmpXTFOYc/Rinst10c40378b8d32/Pbase/extdata/Thermo_Hela_PRTC_1_selected.mzid
    ## 2   /private/tmp/RtmpXTFOYc/Rinst10c40378b8d32/Pbase/extdata/Thermo_Hela_PRTC_1_selected.mzid
    ## 3   /private/tmp/RtmpXTFOYc/Rinst10c40378b8d32/Pbase/extdata/Thermo_Hela_PRTC_1_selected.mzid
    ## 4   /private/tmp/RtmpXTFOYc/Rinst10c40378b8d32/Pbase/extdata/Thermo_Hela_PRTC_1_selected.mzid
    ## 5   /private/tmp/RtmpXTFOYc/Rinst10c40378b8d32/Pbase/extdata/Thermo_Hela_PRTC_1_selected.mzid
    ## ...                                                                                       ...
    ## 32  /private/tmp/RtmpXTFOYc/Rinst10c40378b8d32/Pbase/extdata/Thermo_Hela_PRTC_1_selected.mzid
    ## 33  /private/tmp/RtmpXTFOYc/Rinst10c40378b8d32/Pbase/extdata/Thermo_Hela_PRTC_1_selected.mzid
    ## 34  /private/tmp/RtmpXTFOYc/Rinst10c40378b8d32/Pbase/extdata/Thermo_Hela_PRTC_1_selected.mzid
    ## 35  /private/tmp/RtmpXTFOYc/Rinst10c40378b8d32/Pbase/extdata/Thermo_Hela_PRTC_1_selected.mzid
    ## 36  /private/tmp/RtmpXTFOYc/Rinst10c40378b8d32/Pbase/extdata/Thermo_Hela_PRTC_1_selected.mzid
    ##
    ## ...
    ## <5 more elements>

Specific proteins can be extracted by index of name using `[` and proteins and their peptide features can be plotted with the default plot method.

    seqnames(p)

    ## [1] "A4UGR9"   "A6H8Y1"   "O43707"   "O75369"   "P00558"   "P02545"
    ## [7] "P04075"   "P04075-2" "P60709"

    plot(p[c(1,9)])

More details can be found in `?Proteins`. The object generated above is also directly available as `data(p)`.

## Session information

    sessionInfo()

    ## R version 3.1.1 Patched (2014-09-25 r66681)
    ## Platform: x86_64-apple-darwin10.8.0 (64-bit)
    ##
    ## locale:
    ## [1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
    ##
    ## attached base packages:
    ##  [1] stats4    parallel  grid      stats     graphics  grDevices utils
    ##  [8] datasets  methods   base
    ##
    ## other attached packages:
    ##  [1] mzID_1.4.0
    ##  [2] TxDb.Hsapiens.UCSC.hg19.knownGene_3.0.0
    ##  [3] GenomicFeatures_1.18.0
    ##  [4] AnnotationDbi_1.28.0
    ##  [5] Biobase_2.26.0
    ##  [6] ggplot2_1.0.0
    ##  [7] BSgenome.Hsapiens.UCSC.hg19_1.4.0
    ##  [8] BSgenome_1.34.0
    ##  [9] rtracklayer_1.26.0
    ## [10] biomaRt_2.22.0
    ## [11] Biostrings_2.34.0
    ## [12] XVector_0.6.0
    ## [13] Pbase_0.4.0
    ## [14] Rcpp_0.11.3
    ## [15] Gviz_1.10.0
    ## [16] GenomicRanges_1.18.0
    ## [17] GenomeInfoDb_1.2.0
    ## [18] IRanges_2.0.0
    ## [19] S4Vectors_0.4.0
    ## [20] BiocGenerics_0.12.0
    ## [21] BiocStyle_1.4.0
    ##
    ## loaded via a namespace (and not attached):
    ##  [1] acepack_1.3-3.3          affy_1.44.0
    ##  [3] affyio_1.34.0            base64enc_0.1-2
    ##  [5] BatchJobs_1.4            BBmisc_1.7
    ##  [7] BiocInstaller_1.16.0     BiocParallel_1.0.0
    ##  [9] biovizBase_1.14.0        bitops_1.0-6
    ## [11] brew_1.0-6               checkmate_1.4
    ## [13] chron_2.3-45             cleaver_1.4.0
    ## [15] cluster_1.15.3           codetools_0.2-9
    ## [17] colorspace_1.2-4         data.table_1.9.4
    ## [19] DBI_0.3.1                dichromat_2.0-0
    ## [21] digest_0.6.4             doParallel_1.0.8
    ## [23] evaluate_0.5.5           fail_1.2
    ## [25] foreach_1.4.2            foreign_0.8-61
    ## [27] formatR_1.0              Formula_1.1-2
    ## [29] GenomicAlignments_1.2.0  gtable_0.1.2
    ## [31] Hmisc_3.14-5             impute_1.40.0
    ## [33] iterators_1.0.7          knitr_1.7
    ## [35] labeling_0.3             lattice_0.20-29
    ## [37] latticeExtra_0.6-26      limma_3.22.0
    ## [39] MALDIquant_1.11          markdown_0.7.4
    ## [41] MASS_7.3-35              matrixStats_0.10.0
    ## [43] mime_0.2                 MSnbase_1.14.0
    ## [45] munsell_0.4.2            mzR_2.0.0
    ## [47] nnet_7.3-8               pcaMethods_1.56.0
    ## [49] plyr_1.8.1               preprocessCore_1.28.0
    ## [51] proto_0.3-10             Pviz_1.0.0
    ## [53] R.methodsS3_1.6.1        RColorBrewer_1.0-5
    ## [55] RCurl_1.95-4.3           reshape2_1.4
    ## [57] rpart_4.1-8              Rsamtools_1.18.0
    ## [59] RSQLite_0.11.4           scales_0.2.4
    ## [61] sendmailR_1.2-1          splines_3.1.1
    ## [63] stringr_0.6.2            survival_2.37-7
    ## [65] tools_3.1.1              VariantAnnotation_1.12.0
    ## [67] vsn_3.34.0               XML_3.98-1.1
    ## [69] zlibbioc_1.12.0