| createUCSCtrack {sesame} | R Documentation |
Turn beta values into a UCSC browser track
createUCSCtrack(betas, output = NULL, platform = "HM450", refversion = "hg38")
betas |
a named numeric vector |
output |
output file name |
platform |
HM450, EPIC etc. |
refversion |
hg38, hg19 etc. |
when output is null, return a data.frame, otherwise NULL
betas.tissue <- sesameDataGet('HM450.1.TCGA.PAAD')$betas
## add output to create an actual file
df <- createUCSCtrack(betas.tissue)
## to convert to bigBed
## sort -k1,1 -k2,2n output.bed >output_sorted.bed
## bedToBigBed output_sorted.bed hg38.chrom output.bb