library(dplyr)
# Read in SILAC passage 0 peptide data
# 100% light cells (R0K0)
pepg <- read.delim(
here::here("data-raw/Molm_13_P0_PeptideGroups.txt")
)
# Remove extraneous periods from colnames
colnames(pepg) <- remove_dots(colnames(pepg))
# Set seed for reproducible sampling
set.seed(2021)
# Get all peptides for a random sample of 100 Master proteins
pep_silac_p0 <- pepg %>%
filter(Master.Protein.Accessions %in% sample(unique(Master.Protein.Accessions), 100))
# Output .rda file
usethis::use_data(pep_silac_p0, overwrite = TRUE)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.