inst/doc/doppelgangR.R

## ---- echo=FALSE, results='hide', eval=FALSE-------------------------------
#  library(knitr)
#  opts_knit$set(cache=TRUE)

## ---- echo=FALSE, message=FALSE--------------------------------------------
library(doppelgangR)

## ---- message = FALSE------------------------------------------------------
library(curatedOvarianData)
data(GSE32062.GPL6480_eset)
data(GSE17260_eset)

## ----testesets-------------------------------------------------------------
testesets <- list(JapaneseA=GSE32062.GPL6480_eset,
                  Yoshihara2010=GSE17260_eset)

## ----rundopp, results="hide", message=FALSE, cache=TRUE--------------------
results1 <- doppelgangR(testesets, phenoFinder.args=NULL)

## ----summarizedop, results='hide'------------------------------------------
summary(results1)

## ----plotdop, fig.cap="Doppelgängers identified on the basis of similar expression profiles.  The vertical red lines indicate samples that were flagged."----
par(mfrow=c(2,2), las=1)
plot(results1)

## --------------------------------------------------------------------------
plot(results1, plot.pair=c("JapaneseA", "JapaneseA"))

## ----rundopp2, results="hide", message=FALSE, eval=FALSE-------------------
#  results1 <- doppelgangR(testesets,
#          outlierFinder.expr.args = list(bonf.prob = 1.0, transFun = atanh,
#                                         tail = "upper"))

## --------------------------------------------------------------------------
mat <- matrix(1:4, ncol=2)
library(Biobase)
eset <- ExpressionSet(mat)
class(eset)

## ---- eval=FALSE-----------------------------------------------------------
#  results2 <- doppelgangR(testesets, BPPARAM = MulticoreParam(workers = 8))

Try the doppelgangR package in your browser

Any scripts or data that you put into this service are public.

doppelgangR documentation built on May 2, 2018, 4:47 a.m.