## ----echo = FALSE,hide=TRUE, message=FALSE, warning=FALSE---------------------
library(ELMER)
library(DT)
library(dplyr)
dir.create("result",showWarnings = FALSE)
library(BiocStyle)

## ----message=FALSE------------------------------------------------------------
# get distal probes that are 2kb away from TSS on chromosome 1
distal.probes <- get.feature.probe(
  genome = "hg19", 
  met.platform = "450K", 
  rm.chr = paste0("chr",c(2:22,"X","Y"))
)

## ----eval=TRUE, message=FALSE-------------------------------------------------
library(MultiAssayExperiment)
library(ELMER.data)
data(LUSC_RNA_refined,package = "ELMER.data")
data(LUSC_meth_refined,package = "ELMER.data")
GeneExp[1:5,1:5]
Meth[1:5,1:5]
mae <- createMAE(
  exp = GeneExp, 
  met = Meth,
  save = TRUE,
  linearize.exp = TRUE,
  save.filename = "mae.rda",
  filter.probes = distal.probes,
  met.platform = "450K",
  genome = "hg19",
  TCGA = TRUE
)
as.data.frame(colData(mae)[1:5,])  %>% datatable(options = list(scrollX = TRUE))
as.data.frame(sampleMap(mae)[1:5,])  %>% datatable(options = list(scrollX = TRUE))
as.data.frame(assay(getMet(mae)[1:5,1:5]))  %>% datatable(options = list(scrollX = TRUE))
as.data.frame(assay(getMet(mae)[1:5,1:5])) %>% datatable(options = list(scrollX = TRUE))

## ----eval=FALSE, message=FALSE------------------------------------------------
#  library(ELMER)
#  # example input
#  met <- matrix(rep(0,15),ncol = 5)
#  colnames(met) <- c(
#    "Sample1",
#    "Sample2",
#    "Sample3",
#    "Sample4",
#    "Sample5"
#  )
#  rownames(met) <- c("cg26928153","cg16269199","cg13869341")
#  
#  exp <- matrix(rep(0,15),ncol = 5)
#  colnames(exp) <- c(
#    "Sample1",
#    "Sample2",
#    "Sample3",
#    "Sample4",
#    "Sample5"
#  )
#  rownames(exp) <- c("ENSG00000073282","ENSG00000078900","ENSG00000141510")
#  
#  
#  assay <- c(
#    rep("DNA methylation", ncol(met)),
#    rep("Gene expression", ncol(exp))
#  )
#  primary <- c(colnames(met),colnames(exp))
#  colname <- c(colnames(met),colnames(exp))
#  sampleMap <- data.frame(assay,primary,colname)
#  
#  distal.probes <- get.feature.probe(
#    genome = "hg19",
#    met.platform = "EPIC"
#  )
#  
#  colData <- data.frame(sample = colnames(met))
#  rownames(colData) <- colnames(met)
#  
#  mae <- createMAE(
#    exp = exp,
#    met = met,
#    save = TRUE,
#    filter.probes = distal.probes,
#    colData = colData,
#    sampleMap = sampleMap,
#    linearize.exp = TRUE,
#    save.filename = "mae.rda",
#    met.platform = "EPIC",
#    genome = "hg19",
#    TCGA = FALSE
#  )