Skip to content

Instantly share code, notes, and snippets.

@hbeale

hbeale/plot TCGA RINs.R

Created March 28, 2018 00:34
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save hbeale/c3a65a57f4f30573f86c3213e63463ae to your computer and use it in GitHub Desktop.
Save hbeale/c3a65a57f4f30573f86c3213e63463ae to your computer and use it in GitHub Desktop.
Download ZIP
plot TCGA RIN scores for some samples
Raw
plot TCGA RINs.R
# to reproduce
# The biospecimen_analyte_combined file can be downloaded from Synapse.org under the id: `syn7499603`
library(tidyverse)
biospecimen_values=read_tsv("~/Documents/Dropbox/ucsc/operations/registry/TCGA biospecimen_analyted_combined.tsv")
biospecimen_values %>%
mutate(
rinNumeric=as.numeric(rinvalue),
a260_a280_ratioNumeric=as.numeric(a260_a280_ratio),
source=gsub("^([A-Z]*)[-_].*$", "\\1", bcr_sample_barcode)
) %>%
ggplot() + geom_histogram(aes(rinNumeric), binwidth=0.1) + ggtitle("RIN scores in some TCGA samples") + scale_x_continuous(breaks=0:10)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment