inst/doc/bioqc-signedGenesets.R

## ----setup, include=FALSE------------------------------------------------
library(knitr)
opts_chunk$set(out.extra='style="display:block; margin: auto"', fig.align="center")

## ----lib, warning=FALSE, message=FALSE, results="hide"-------------------
library(BioQC)

## ----gmt-----------------------------------------------------------------
gmtFile <- system.file("extdata/test.gmt", package="BioQC")
## print the file context
cat(readLines(gmtFile), sep="\n")
## read the file into SignedGenesets
genesets <- readSignedGmt(gmtFile)
print(genesets)

## ----gmtPos--------------------------------------------------------------
genesets <- readSignedGmt(gmtFile, nomatch="pos")
print(genesets)

## ----data----------------------------------------------------------------
set.seed(1887)
testN <- 100L
testSampleCount <- 3L
testGenes <- c("AKT1", "AKT2", "ERBB2", "ERBB3", "EGFR","TSC1", "TSC2", "GATA2", "GATA4", "GATA1", "GATA3")
testRows <- c(testGenes, paste("Gene", (length(testGenes)+1):testN, sep=""))
testMatrix <- matrix(rnorm(testN*testSampleCount, sd=0.1), nrow=testN, dimnames=list(testRows, NULL))
testMatrix[1:2,] <- testMatrix[1:2,]+10
testMatrix[6:7,] <- testMatrix[6:7,]-10
testMatrix[3:4,] <- testMatrix[3:4,]-5
testMatrix[5,] <- testMatrix[5,]+5
testEset <- new("ExpressionSet", exprs=testMatrix)

## ----index---------------------------------------------------------------
testIndex <- matchGenes(genesets, testEset, col=NULL)
print(testIndex)

## ----runWmwGreater-------------------------------------------------------
wmwResult.greater <- wmwTest(testEset, testIndex, valType="p.greater")
print(wmwResult.greater)

## ----runWmwLess----------------------------------------------------------
wmwResult.less <- wmwTest(testEset, testIndex, valType="p.less")
print(wmwResult.less)

## ----runWmwTwoSided------------------------------------------------------
wmwResult.two.sided <- wmwTest(testEset, testIndex, valType="p.two.sided")
print(wmwResult.two.sided)

## ----runWmwQ-------------------------------------------------------------
wmwResult.Q <- wmwTest(testEset, testIndex, valType="Q")
print(wmwResult.Q)

## ----session-------------------------------------------------------------
sessionInfo()

Try the BioQC package in your browser

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

BioQC documentation built on Nov. 1, 2018, 2:12 a.m.