To validate our Retention Time (RT) prediction in this vignette file, we compare the predicted hydrophobicity value using the ssrc
method Krokhin et al. (2004) implemented in the
protViz package Panse and Grossmann (2019).
The following code snippet performs the comparison on the F255744 data. The file contains amino acid sequences representing the designed flycodes.
library(NestLink)
# load(url("http://fgcz-ms.uzh.ch/~cpanse/p1875/F255744.RData"))
# F255744 <- as.data.frame.mascot(F255744)
# now available through ExperimentHub
library(ExperimentHub)
eh <- ExperimentHub();
## snapshotDate(): 2021-10-18
load(query(eh, c("NestLink", "F255744.RData"))[[1]])
## see ?NestLink and browseVignettes('NestLink') for documentation
## loading from cache
.ssrc.mascot(F255744, scores = c(10, 20, 40, 50),
pch = 16,
col = rgb(0.1,0.1,0.1,
alpha = 0.1)
)
## [[1]]
##
## Call:
## lm(formula = xx.ssrc ~ xx$RTINSECONDS)
##
## Residuals:
## Min 1Q Median 3Q Max
## -38.954 -2.248 0.015 2.228 71.167
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -5.580e+00 2.030e-01 -27.48 <2e-16 ***
## xx$RTINSECONDS 8.849e-03 7.434e-05 119.04 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 5.884 on 12295 degrees of freedom
## Multiple R-squared: 0.5354, Adjusted R-squared: 0.5354
## F-statistic: 1.417e+04 on 1 and 12295 DF, p-value: < 2.2e-16
##
##
## [[2]]
##
## Call:
## lm(formula = xx.ssrc ~ xx$RTINSECONDS)
##
## Residuals:
## Min 1Q Median 3Q Max
## -37.387 -2.040 -0.042 1.930 46.035
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -6.976e+00 1.621e-01 -43.03 <2e-16 ***
## xx$RTINSECONDS 9.447e-03 6.018e-05 156.99 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 4.12 on 9835 degrees of freedom
## Multiple R-squared: 0.7148, Adjusted R-squared: 0.7147
## F-statistic: 2.464e+04 on 1 and 9835 DF, p-value: < 2.2e-16
##
##
## [[3]]
##
## Call:
## lm(formula = xx.ssrc ~ xx$RTINSECONDS)
##
## Residuals:
## Min 1Q Median 3Q Max
## -15.260 -1.963 -0.114 1.735 45.342
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -7.690e+00 1.784e-01 -43.11 <2e-16 ***
## xx$RTINSECONDS 9.781e-03 6.724e-05 145.46 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 3.506 on 5574 degrees of freedom
## Multiple R-squared: 0.7915, Adjusted R-squared: 0.7915
## F-statistic: 2.116e+04 on 1 and 5574 DF, p-value: < 2.2e-16
##
##
## [[4]]
##
## Call:
## lm(formula = xx.ssrc ~ xx$RTINSECONDS)
##
## Residuals:
## Min 1Q Median 3Q Max
## -9.570 -2.019 -0.142 1.754 45.200
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -7.827e+00 2.173e-01 -36.02 <2e-16 ***
## xx$RTINSECONDS 9.848e-03 8.271e-05 119.06 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 3.579 on 3650 degrees of freedom
## Multiple R-squared: 0.7952, Adjusted R-squared: 0.7952
## F-statistic: 1.418e+04 on 1 and 3650 DF, p-value: < 2.2e-16
Here is the output of the sessionInfo()
command.
## R version 4.1.1 (2021-08-10)
## Platform: x86_64-pc-linux-gnu (64-bit)
## Running under: Ubuntu 20.04.3 LTS
##
## Matrix products: default
## BLAS: /home/biocbuild/bbs-3.14-bioc/R/lib/libRblas.so
## LAPACK: /home/biocbuild/bbs-3.14-bioc/R/lib/libRlapack.so
##
## locale:
## [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
## [3] LC_TIME=en_GB LC_COLLATE=C
## [5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8
## [7] LC_PAPER=en_US.UTF-8 LC_NAME=C
## [9] LC_ADDRESS=C LC_TELEPHONE=C
## [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
##
## attached base packages:
## [1] stats4 stats graphics grDevices utils datasets methods
## [8] base
##
## other attached packages:
## [1] scales_1.1.1 ggplot2_3.3.5
## [3] NestLink_1.10.0 ShortRead_1.52.0
## [5] GenomicAlignments_1.30.0 SummarizedExperiment_1.24.0
## [7] Biobase_2.54.0 MatrixGenerics_1.6.0
## [9] matrixStats_0.61.0 Rsamtools_2.10.0
## [11] GenomicRanges_1.46.0 BiocParallel_1.28.0
## [13] protViz_0.7.0 gplots_3.1.1
## [15] Biostrings_2.62.0 GenomeInfoDb_1.30.0
## [17] XVector_0.34.0 IRanges_2.28.0
## [19] S4Vectors_0.32.0 ExperimentHub_2.2.0
## [21] AnnotationHub_3.2.0 BiocFileCache_2.2.0
## [23] dbplyr_2.1.1 BiocGenerics_0.40.0
## [25] BiocStyle_2.22.0
##
## loaded via a namespace (and not attached):
## [1] nlme_3.1-153 bitops_1.0-7
## [3] bit64_4.0.5 RColorBrewer_1.1-2
## [5] filelock_1.0.2 httr_1.4.2
## [7] tools_4.1.1 bslib_0.3.1
## [9] utf8_1.2.2 R6_2.5.1
## [11] KernSmooth_2.23-20 mgcv_1.8-38
## [13] colorspace_2.0-2 DBI_1.1.1
## [15] withr_2.4.2 tidyselect_1.1.1
## [17] bit_4.0.4 curl_4.3.2
## [19] compiler_4.1.1 DelayedArray_0.20.0
## [21] labeling_0.4.2 bookdown_0.24
## [23] sass_0.4.0 caTools_1.18.2
## [25] rappdirs_0.3.3 stringr_1.4.0
## [27] digest_0.6.28 rmarkdown_2.11
## [29] jpeg_0.1-9 pkgconfig_2.0.3
## [31] htmltools_0.5.2 highr_0.9
## [33] fastmap_1.1.0 rlang_0.4.12
## [35] RSQLite_2.2.8 shiny_1.7.1
## [37] farver_2.1.0 jquerylib_0.1.4
## [39] generics_0.1.1 hwriter_1.3.2
## [41] jsonlite_1.7.2 gtools_3.9.2
## [43] dplyr_1.0.7 RCurl_1.98-1.5
## [45] magrittr_2.0.1 GenomeInfoDbData_1.2.7
## [47] Matrix_1.3-4 munsell_0.5.0
## [49] Rcpp_1.0.7 fansi_0.5.0
## [51] lifecycle_1.0.1 stringi_1.7.5
## [53] yaml_2.2.1 zlibbioc_1.40.0
## [55] grid_4.1.1 blob_1.2.2
## [57] parallel_4.1.1 promises_1.2.0.1
## [59] crayon_1.4.1 lattice_0.20-45
## [61] splines_4.1.1 KEGGREST_1.34.0
## [63] magick_2.7.3 knitr_1.36
## [65] pillar_1.6.4 codetools_0.2-18
## [67] glue_1.4.2 BiocVersion_3.14.0
## [69] evaluate_0.14 latticeExtra_0.6-29
## [71] BiocManager_1.30.16 png_0.1-7
## [73] vctrs_0.3.8 httpuv_1.6.3
## [75] gtable_0.3.0 purrr_0.3.4
## [77] assertthat_0.2.1 cachem_1.0.6
## [79] xfun_0.27 mime_0.12
## [81] xtable_1.8-4 later_1.3.0
## [83] tibble_3.1.5 AnnotationDbi_1.56.1
## [85] memoise_2.0.0 ellipsis_0.3.2
## [87] interactiveDisplayBase_1.32.0
Krokhin, O. V., R. Craig, V. Spicer, W. Ens, K. G. Standing, R. C. Beavis, and J. A. Wilkins. 2004. βAn improved model for prediction of retention times of tryptic peptides in ion pair reversed-phase HPLC: its application to protein peptide mapping by off-line HPLC-MALDI MS.β Mol. Cell Proteomics 3 (9): 908β19.
Panse, Christian, and Jonas Grossmann. 2019. protViz: Visualizing and Analyzing Mass Spectrometry Related Data in Proteomics. Vienna, Austria: R Foundation for Statistical Computing. https://www.R-project.org.