### R code from vignette source 'vignettes/methylumi/inst/doc/methylumi450k.Rnw' ################################################### ### code chunk number 1: methylumi450k.Rnw:23-25 ################################################### library(knitr) opts_chunk$set(tidy=FALSE,cache=TRUE,size='scriptsize') ################################################### ### code chunk number 2: setup ################################################### options('width'=50) ################################################### ### code chunk number 3: loadLibraries ################################################### suppressPackageStartupMessages(require('methylumi')) suppressPackageStartupMessages(require('IlluminaHumanMethylation450k.db')) suppressPackageStartupMessages(require('TCGAMethylation450k')) ################################################### ### code chunk number 4: loadData ################################################### ## read in 10 BRCA IDATs idatPath <- system.file('extdata/idat',package='TCGAMethylation450k') mset450k <- methylumIDAT(getBarcodes(path=idatPath), idatPath=idatPath) sampleNames(mset450k) <- paste0('TCGA', seq_along(sampleNames(mset450k))) show(mset450k) ################################################### ### code chunk number 5: controls ################################################### library(ggplot2) ## for larger datasets, the by.type argument be set to FALSE ## positional effects will manifest as a wave-like pattern p <- qc.probe.plot(mset450k, by.type=TRUE) print(p) ################################################### ### code chunk number 6: preprocess ################################################### mset450k.proc <- stripOOB(normalizeMethyLumiSet(methylumi.bgcorr(mset450k))) ################################################### ### code chunk number 7: controls2 ################################################### library(ggplot2) p2 <- qc.probe.plot(mset450k.proc, by.type=TRUE) print(p2) ################################################### ### code chunk number 8: coerceLumi ################################################### suppressPackageStartupMessages(require(lumi)) mset450k.lumi <- as(mset450k.proc, 'MethyLumiM') show(mset450k.lumi) ################################################### ### code chunk number 9: coerceBack ################################################### mset450k.andBack <- as(mset450k.lumi, 'MethyLumiSet') show(mset450k.andBack) ################################################### ### code chunk number 10: coerceMinfi ################################################### suppressPackageStartupMessages(require(FDb.InfiniumMethylation.hg19)) rgSet450k <- as(mset450k, 'RGChannelSet') show(rgSet450k) ################################################### ### code chunk number 11: coerceMinfi2 ################################################### suppressPackageStartupMessages(require(minfi)) suppressPackageStartupMessages(require(IlluminaHumanMethylation450kanno.ilmn12.hg19)) grSet450k <- mapToGenome(mset450k.andBack) sexChroms <- GRanges( seqnames=c('chrX','chrY'), IRanges(start=c(1, 1), end=c(155270560, 59373566)), strand=c('*','*') ) summary(subsetByOverlaps(grSet450k, sexChroms)) dim(subsetByOverlaps(grSet450k, sexChroms)) ################################################### ### code chunk number 12: subsetMinfi ################################################### ## perhaps more topical: suppressPackageStartupMessages(require(TxDb.Hsapiens.UCSC.hg19.knownGene)) suppressPackageStartupMessages(require(Homo.sapiens)) txdb <- TxDb.Hsapiens.UCSC.hg19.knownGene KDM6AEntrezID=org.Hs.egSYMBOL2EG[['KDM6A']] txs.KDM6A <- transcriptsBy(txdb, 'gene')[[KDM6AEntrezID]] tss.KDM6A <- unique(resize(txs.KDM6A, 1, fix='start')) ## two start sites promoters.KDM6A <- flank(tss.KDM6A, 100) ## an arbitrary distance upstream show( subsetByOverlaps(grSet450k, promoters.KDM6A) ) ## probes in this window