Nothing
## ----style, echo = FALSE, results = 'asis'------------------------------------
BiocStyle::markdown()
## ---- message=FALSE-----------------------------------------------------------
library(Chicago)
library(PCHiCdata)
## -----------------------------------------------------------------------------
dataPath <- system.file("extdata", package="PCHiCdata")
testDesignDir <- file.path(dataPath, "hg19TestDesign")
dir(testDesignDir)
## -----------------------------------------------------------------------------
testDataPath <- file.path(dataPath, "GMchinputFiles")
dir(testDataPath)
files <- c(
file.path(testDataPath, "GM_rep1.chinput"),
file.path(testDataPath, "GM_rep2.chinput"),
file.path(testDataPath, "GM_rep3.chinput")
)
## -----------------------------------------------------------------------------
settingsFile <- file.path(system.file("extdata", package="PCHiCdata"),
"sGM12878Settings", "sGM12878.settingsFile")
## ---- message=FALSE-----------------------------------------------------------
library(Chicago)
cd <- setExperiment(designDir = testDesignDir, settingsFile = settingsFile)
## ---- message=FALSE-----------------------------------------------------------
cd <- readAndMerge(files=files, cd=cd)
## ---- eval=FALSE--------------------------------------------------------------
# cd <- chicagoPipeline(cd)
## ---- echo=FALSE, message=FALSE-----------------------------------------------
cd <- chicagoPipeline(cd)
## -----------------------------------------------------------------------------
outputDirectory <- tempdir()
exportResults(cd, file.path(outputDirectory, "vignetteOutput"))
## ----echo=FALSE---------------------------------------------------------------
a <- read.table(file.path(outputDirectory, "vignetteOutput.ibed"), header=TRUE)
head(a)
## ----echo=FALSE---------------------------------------------------------------
a <- read.table(file.path(outputDirectory, "vignetteOutput_seqmonk.txt"), header=FALSE)
head(a)
## ----echo=FALSE---------------------------------------------------------------
a <- read.table(file.path(outputDirectory, "vignetteOutput_washU_text.txt"), header=FALSE)
head(a)
## ----eval=FALSE---------------------------------------------------------------
#
# plottedBaitIDs <- plotBaits(cd, n=6)
#
## ----echo=FALSE, fig.height=10------------------------------------------------
invisible(plotBaits(cd, baits=c(415839, 404491, 425847, 417632, 409335, 414114)))
## -----------------------------------------------------------------------------
featuresFolder <- file.path(dataPath, "GMfeatures")
dir(featuresFolder)
featuresFile <- file.path(featuresFolder, "featuresGM.txt")
featuresTable <- read.delim(featuresFile, header=FALSE, as.is=TRUE)
featuresList <- as.list(featuresTable$V2)
names(featuresList) <- featuresTable$V1
featuresList
## ---- message=FALSE, fig.width=12, fig.height=7-------------------------------
no_bins <- ceiling(max(abs(intData(cd)$distSign), na.rm = TRUE)/1e4)
enrichmentResults <- peakEnrichment4Features(cd, folder=featuresFolder,
list_frag=featuresList, no_bins=no_bins, sample_number=100)
## -----------------------------------------------------------------------------
enrichmentResults
## ---- message=FALSE-----------------------------------------------------------
library(GenomicInteractions)
library(GenomicRanges)
gi <- exportToGI(cd)
## ---- message=FALSE-----------------------------------------------------------
library(AnnotationHub)
ah <- AnnotationHub()
hs <- query(ah, c("GRanges", "EncodeDCC", "Homo sapiens", "H3k4me1"))
enhancerTrack <- hs[["AH23254"]]
## -----------------------------------------------------------------------------
otherEnds <- anchorTwo(gi)
otherEnds <- renameSeqlevels(otherEnds, c("chr20","chr21"))
## -----------------------------------------------------------------------------
findOverlaps(otherEnds, enhancerTrack)
## -----------------------------------------------------------------------------
hs["AH23254"]$genome
## -----------------------------------------------------------------------------
head(intData(cd), 2)
## -----------------------------------------------------------------------------
newCd = copyCD(cd)
## -----------------------------------------------------------------------------
weightsPath <- file.path(system.file("extdata", package="Chicago"),
"weights")
dir(weightsPath)
## ---- message=FALSE-----------------------------------------------------------
weightSettings <- file.path(weightsPath, "GM12878-2reps.settings")
cd <- setExperiment(designDir = testDesignDir, settingsFile = weightSettings)
## -----------------------------------------------------------------------------
sessionInfo()
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.