README.md

In sahilseth/tcga.archiver: TCGA Archive Creater (HMS - MDA) WGS sequencing center

tcga.archiver

Creating metadata for TCGA lowpass data

Installing this package

require(devtools)
install_github(repo="tcga.archiver", username="sahilseth")

Reading in example params

require(tcga.archiver)
dcc_params = file.path(path.package("tcga.archiver"), "files/DCC.params")
opt_data <- read.delim(dcc_params, sep = "\t", header = TRUE, as.is = TRUE)

Make a single SDRF row based on example data

file_type = "tsv" ## or vcf for breakdancer
disease = 'HNSC'
platform = "IlluminaHiSeq_DNASeqC"
center = "hms.harvard.edu"


## these are the TSV files we upload:
tsvs = file.path(path.package("tcga.archiver"), 
  "files/TCGA-CN-5356-01A-01D-1431___120412_SN208_0283_AD0V1GACXX---TCGA-CN-5356-10A-01D-1431___120412_SN208_0284_BD0T51ACXX.tsv")

## and example sample
vec <- c(files = tsvs, 
  NORMALSAMPLEID = "TCGA-CN-5356-10A-01D-1431-02", TUMORSAMPLEID = "TCGA-CN-5356-01A-01D-1431-02", 
  NORMALANALYSISUUID = getCGHubStatus("TCGA-CN-5356-01A-01D-1431-02", 
    by = "legacy_sample_id", get = "both")$uuid,
  TUMORANALYSISUUID = getCGHubStatus("TCGA-CN-5356-10A-01D-1431-02", 
    by = "legacy_sample_id", get = "both")$uuid)

out <- get_sdrf_row(vec = vec, disease = disease, 
  center = center, platform = platform,
  opt_data = opt_data, file_type = file_type,
  dat_level = 3, dat_batch = 0, dat_rev = 0)

sahilseth/tcga.archiver documentation built on May 29, 2019, 12:57 p.m.