Simulation of molecular regulatory networks with graphical models
|
|
- Brent Banks
- 6 years ago
- Views:
Transcription
1 Simulation of molecular regulatory networks with graphical models Inma Tur 1 inma.tur@upf.edu Alberto Roverato 2 alberto.roverato@unibo.it Robert Castelo 1 robert.castelo@upf.edu 1 Universitat Pompeu Fabra, Barcelona, Spain 2 Università di Bologna, Bologna, Italy user! UCLM, Albacete, Spain - July Inma Tur, Alberto Roverato, Robert Castelo Simulation of molecular regulatory networks with graphical models 1 / 27
2 Motivation - Genomics data High-throughput genomics technologies produce high-dimensional (p n) multivariate data sets of continuous and discrete random variables. Gene expression data Network of molecular regulatory interactions Genetical genomics data Inma Tur, Alberto Roverato, Robert Castelo Simulation of molecular regulatory networks with graphical models 2 / 27
3 Motivation - Graphical Markov Models (GMM) Gaussian GMMs Homogeneous Mixed GMMs X V N p (µ, Σ) 1 X V N p (µ(i), Σ(i)) with Σ(i) Σ I1 2 3 Y1 Y2 > library(qpgraph) > set.seed(12345) > gmm <- ruggmm(dregulargraphparam()) > round(solve(gmm$sigma), digits=1) Y3 > library(qpgraph) > set.seed(12345) > gmm <- rhmgmm(dregularmarkedgraphparam()) > round(solve(gmm$sigma), digits=1) Y1 Y2 Y3 Y Y Y > gmm$mean() Y1 Y2 Y Inma Tur, Alberto Roverato, Robert Castelo Simulation of molecular regulatory networks with graphical models 3 / 27
4 Motivation - Graphical Markov Model (GMM) Testing GMM estimation procedures on simulated data is a fundamental step to verify properties such as correctness or asymptotic behavior. The R/Bioconductor package qpgraph implements algorithms to simulate Gaussian GMMs, homogeneous mixed GMMs and data from them. Simulating Gaussian GMMs requires simulating covariance matrices whose inverse matches a zero-pattern defined by missing edges in a graph. Simulating homogeneous mixed GMMs requires simulating conditional covariance matrices whose inverse matches a zero-pattern defined by a graph, and conditional mean vectors satisfying additive effects. The rest of this talk is based on the vignette entitled Simulating molecular regulatory networks using qpgraph from the qpgraph package. Inma Tur, Alberto Roverato, Robert Castelo Simulation of molecular regulatory networks with graphical models 4 / 27
5 Outline 1 Simulation of graphs for GMMs 2 Simulation of undirected Gaussian Graphical Markov Models 3 Simulation of Homogeneous Mixed Graphical Markov Models 4 Simulation of eqtl models of experimental crosses 5 Conclusions Inma Tur, Alberto Roverato, Robert Castelo Simulation of molecular regulatory networks with graphical models 5 / 27
6 Outline 1 Simulation of graphs for GMMs 2 Simulation of undirected Gaussian Graphical Markov Models 3 Simulation of Homogeneous Mixed Graphical Markov Models 4 Simulation of eqtl models of experimental crosses 5 Conclusions Inma Tur, Alberto Roverato, Robert Castelo Simulation of molecular regulatory networks with graphical models 6 / 27
7 Simulation of graphs for GMMs Consider undirected labeled graphs G = (V, E) where V = {1,..., p} is the vertex set indexing a vector of random variables (r.v. s) X V = {X 1,..., X p } that belong to some probability distribution P. E (V V ) is the edge set, such that (i, j ) E X i X j X V \{X i, X j }, holds in P (pairwise Markov property w.r.t G). Let {A, B, S} V, S separates A from B in G (A G B S) if every path between A and B intersects S. G should be such that A G B S X A X B X S, holds in P (global Markov property w.r.t G). In the context of GMMs we want to simulate graphs in which we can control separation and sparseness. Inma Tur, Alberto Roverato, Robert Castelo Simulation of molecular regulatory networks with graphical models 7 / 27
8 Simulation of graphs for GMMs qpgraph simulates undirected graphs according to 1 the type of graph pure single type of vertices marked two subsets of vertices (associated to discrete and continuous random variables) 2 the model to simulate the random graph Erdös-Rényi edges occur with equal probability, or graphs are chosen uniformly at random with given number of vertices & edges. d-regular vertices have a constant degree d (Harary, 1969). It follows that graph density is a linear function of d: D = d/(p 1). This is implemented in the function rgraphbam(n, param,...) which returns objects of the class graphbam defined in the graph package. Inma Tur, Alberto Roverato, Robert Castelo Simulation of molecular regulatory networks with graphical models 8 / 27
9 Simulation of graphs for GMMs Input parameters in rgraphbam() are defined by S4 classes graphparam p labels markedgraphparam pi Ilabels py Ylabels ergraphparam m [# edges] prob [Pr(edge)] dregulargraphparam d [vertex degree] exclude ermarkedgraphparam dregularmarkedgraphparam Inma Tur, Alberto Roverato, Robert Castelo Simulation of molecular regulatory networks with graphical models 9 / 27
10 Simulation of graphs for GMMs d-regular graphs are simulated with the Steger and Wormald (1999) algorithm > library(qpgraph) > set.seed(1234) > g <- rgraphbam(dregularmarkedgraphparam(pi=2, py=10, d=3)) > plot(g, lwd=3) I2 Y1 I1 Y2 Y4 Y3 Y10 Y6 Y5 Y7 Y8 Y9 plot() is overloaded to use the functionality from the Rgraphviz package. Inma Tur, Alberto Roverato, Robert Castelo Simulation of molecular regulatory networks with graphical models 10 / 27
11 Outline 1 Simulation of graphs for GMMs 2 Simulation of undirected Gaussian Graphical Markov Models 3 Simulation of Homogeneous Mixed Graphical Markov Models 4 Simulation of eqtl models of experimental crosses 5 Conclusions Inma Tur, Alberto Roverato, Robert Castelo Simulation of molecular regulatory networks with graphical models 11 / 27
12 Simulation of undirected Gaussian GMMs Undirected Gaussian GMMs are multivariate models on continuous r.v. s X V = {X 1,..., X p } determined by an undirected graph G = (V, E) with V = {1,..., p} and E V V such that X V N p (µ, Σ) where {Σ 1 } ij = 0 for i j and (i, j ) G Therefore, to simulate an undirected Gaussian GMM we need to build a matrix Σ such that 1 Σ is positive definite (Σ S + ), 2 the off-diagonal cells of the scaled Σ corresponding to the present edges in G match a given marginal correlation ρ, 3 the zero pattern of Σ 1 matches the missing edges in G. This is not straightforward since setting directly off-diagonal cells to zero in some initial Γ S + will not typically lead to a positive definite matrix. Inma Tur, Alberto Roverato, Robert Castelo Simulation of molecular regulatory networks with graphical models 12 / 27
13 Simulation of undirected Gaussian GMMs Let Γ G be an incomplete matrix with elements {γ ij } for i = j or (i, j ) G. 1 γ 11 γ 12 γ Γ G = γ 21 γ 22 γ 24 γ 31 γ 33 γ 34 γ 42 γ 43 γ 44 4 Γ is a positive completion of Γ G if Γ S + and {Γ 1 } ij =0 for i j, (i, j ) G. Draw Γ G from a Wishart distributionw p (Λ, p); Λ= R, =diag({ 1/p} p ) and R = {R ij } p p where R ij = 1 for i = j and R ij = ρ for i j. It is required that Λ S + and this happens if and only if 1/(p 1) < ρ < 1. Finally, to obtain Σ Γ from Γ G, qpgraph uses the regression algorithm by Hastie, Tibshirani and Friedman (2009, pg. 634) as matrix completion algorithm. See functions qprndwishart(), qphtf() and qpg2sigma() for further details. Inma Tur, Alberto Roverato, Robert Castelo Simulation of molecular regulatory networks with graphical models 13 / 27
14 Simulation of undirected Gaussian GMMs ruggmm() simulates undirected Gaussian GMMs taking as input either a graphparam or a graphbam object. It returns an S4 object of class UGgmm. > set.seed(12345) > gmm <- ruggmm(n=1, g=dregulargraphparam(p=4, d=2), rho=0.75) > class(gmm) [1] "UGgmm" attr(,"package") [1] "qpgraph" > names(gmm) ## it behaves pretty much like a 'list' object [1] "X" "p" "g" "mean" "sigma" > round(solve(gmm$sigma), digits=1) > gmm2 <- ruggmm(n=1, g=gmm$g) ## fix 'g' simulate covariance only > round(solve(gmm2$sigma), digits=1) Inma Tur, Alberto Roverato, Robert Castelo Simulation of molecular regulatory networks with graphical models 14 / 27 4
15 Simulation of undirected Gaussian GMMs rmvnorm() from the mvtnorm package is overloaded in qpgraph to simulate multivariate normal observations from the undirected GMM: > set.seed(123) > X <- rmvnorm(n=10000, gmm) > head(x, n=3) [1,] [2,] [3,] > round(solve(cov(x)), digits=1) > mean(cor(x)[upper.tri(cor(x)) & as(gmm$g, "matrix")]) [1] As the number of present edges grow, their mean marginal correlation approaches the given one. Inma Tur, Alberto Roverato, Robert Castelo Simulation of molecular regulatory networks with graphical models 15 / 27
16 Outline 1 Simulation of graphs for GMMs 2 Simulation of undirected Gaussian Graphical Markov Models 3 Simulation of Homogeneous Mixed Graphical Markov Models 4 Simulation of eqtl models of experimental crosses 5 Conclusions Inma Tur, Alberto Roverato, Robert Castelo Simulation of molecular regulatory networks with graphical models 16 / 27
17 Simulation of Homogeneous Mixed GMMs Simplifying assumptions (in the context of genetical genomics data): 1 Discrete genotypes affect gene expression and not the other way around. 2 Joint distribution of X is a conditional Gaussian distribution X V N py (µ(i), Σ(i)) with i I. 3 Genotype alleles affect only mean expression levels of genes and not the correlations between them, i.e., Σ(i) Σ is constant throughout i I. 4 Currently, discrete r.v. s are simulated as being marginally independent between them. 5 Currently, every continuous r.v. cannot depend on more than one discrete r.v., where µ(i) = Σ h(i) and h(i) = {h 1 (i),..., h p (i)} contain the mixed linear parameters that define the additive effects. Y1 I1 Y3 Y2 Inma Tur, Alberto Roverato, Robert Castelo Simulation of molecular regulatory networks with graphical models 17 / 27
18 Simulation of Homogeneous Mixed GMMs rhmgmm() simulates homogeneous mixed GMMs taking as input either a markedgraphparam or graphbam object. It returns an S4 object of class HMgmm. Here we simulate one involving 2 50 joint discrete levels: > set.seed(12345) > gmm <- rhmgmm(n=1, dregularmarkedgraphparam(pi=50, py=100, d=2), a=2) > names(gmm) [1] "X" "I" "Y" "p" "pi" "py" "g" "mean" "sigma" [10] "a" "eta2" > print(object.size(gmm), units="kb") 90.1 Kb rcmvnorm() simulates data from a homogeneous mixed GMM using rmvnorm() for the continuous observations (mvtnorm package): > set.seed(123) > X <- rcmvnorm(n=1000, gmm) > Ymix <- gmm$y[which(rowsums(as(gmm$g, "matrix")[gmm$y, gmm$i]) == 1)] > Imix <- gmm$i[apply(as(gmm$g, "matrix")[ymix, gmm$i] == 1, 1, which)] > smean <- sapply(1:length(ymix), function(i, Y, I) tapply(y[, i], I[, i], mean), + X[, Ymix], X[, Imix]) > mean(abs(smean[1, ] - smean[2, ])) ## sample additive effects approach the given one [1] Inma Tur, Alberto Roverato, Robert Castelo Simulation of molecular regulatory networks with graphical models 18 / 27
19 Outline 1 Simulation of graphs for GMMs 2 Simulation of undirected Gaussian Graphical Markov Models 3 Simulation of Homogeneous Mixed Graphical Markov Models 4 Simulation of eqtl models of experimental crosses 5 Conclusions Inma Tur, Alberto Roverato, Robert Castelo Simulation of molecular regulatory networks with graphical models 19 / 27
20 Simulation of eqtl models of experimental crosses S4 class eqtlcross: a genetic map (map class in qtl package) a homogeneous mixed GMM (HMgmm) reqtlcross simulates n eqtlcross models: > reqtlcross(n=2, network=eqtlcrossparam()) [[1]] eqtl backcross model with 20 markers, 20 genes, 20 eqtl and 20 gene-gene expression associations Genetic map [[2]] eqtl backcross model with 20 markers, 20 genes, 20 eqtl and 20 gene-gene expression associations. Location (cm) X Chromosome Inma Tur, Alberto Roverato, Robert Castelo Simulation of molecular regulatory networks with graphical models 20 / 27
21 Simulation of eqtl models of experimental crosses eqtl model with cis-associations only > data(map10, package="qtl") ## load example genetic map from the 'qtl' package > (eqtl <- reqtlcross(eqtlcrossparam(map=map10, genes=100))) eqtl backcross model with 187 markers, 100 genes, 100 eqtl and 100 gene-gene expression associations. > eqtl$model Homogeneous mixed graphical Markov model with 100 discrete and 100 continuous r.v., and 200 edges. > plot(eqtl, main="eqtl model with cis-associations only", cex.lab=1.5, cex.main=1.5) eqtl model with cis associations only Gene location X eqtl location Inma Tur, Alberto Roverato, Robert Castelo Simulation of molecular regulatory networks with graphical models 21 / 27
22 Simulation of eqtl models of experimental crosses eqtl model with cis and trans-associations > set.seed(123) > (eqtl <- reqtlcross(eqtlcrossparam(map=map10, genes=100, cis=0.5, trans=c(5,5)), a=5)) eqtl backcross model with 187 markers, 100 genes, 60 eqtl and 100 gene-gene expression associations. > head(ciseqtl(eqtl), n=3) chrom location QTL gene a QTL1 g QTL2 g QTL3 g6 5 > head(transeqtl(eqtl), n=10) chrom location QTL gene a QTL17 g QTL17 g QTL17 g QTL17 g QTL17 g QTL19 g QTL19 g QTL19 g QTL19 g QTL19 g79 5 > plot(eqtl, cex.lab=1.5, cex.main=1.5, + main="eqtl model with trans-eqtl") Gene location X eqtl model with trans eqtl eqtl location Inma Tur, Alberto Roverato, Robert Castelo Simulation of molecular regulatory networks with graphical models 22 / 27
23 Simulation of eqtl models of experimental crosses sim.cross() from package qtl is overloaded to simulate genotypes from experimental crosses and place eqtl using an input homogeneous mixed GMM: > set.seed(123) > cross <- sim.cross(map10, eqtl) > eqtlgene <- "g5" ## take gene 'g5' with a cis-eqtl > connectedgenes <- graph::edges(eqtl$model$g)[[eqtlgene]] > (connectedgenes <- connectedgenes[connectedgenes %in% eqtl$model$y]) [1] "g15" "g22" > out.mr <- qtl::scanone(cross, method="mr", pheno.col=c(eqtlgene, connectedgenes)) 15 g5 QTL2 LOD scores 10 5 g15 g22 g5 0 QTL2 g Map position (cm) g15 Inma Tur, Alberto Roverato, Robert Castelo Simulation of molecular regulatory networks with graphical models 23 / 27
24 Outline 1 Simulation of graphs for GMMs 2 Simulation of undirected Gaussian Graphical Markov Models 3 Simulation of Homogeneous Mixed Graphical Markov Models 4 Simulation of eqtl models of experimental crosses 5 Conclusions Inma Tur, Alberto Roverato, Robert Castelo Simulation of molecular regulatory networks with graphical models 24 / 27
25 Conclusions The R/Bioconductor package qpgraph enables simulating Gaussian GMMs, homogeneous mixed GMMs and data from them. The implemented algorithms can simulate GMMs with given marginal correlations and additive effects. This provides flexibility when simulating data to verify correctness or asymptotic behavior of GMM estimation algorithms. The use of S4 classes hides algorithmic complexity to the non-expert user and provides a seamless integration with other packages such as graph, Rgraphviz, mvtnorm and qtl. By overloading the function sim.cross() from qtl, qpgraph allows one to easily simulate genetical genomics data from experimental crosses. Inma Tur, Alberto Roverato, Robert Castelo Simulation of molecular regulatory networks with graphical models 25 / 27
26 Session information > tolatex(sessioninfo()) R version ( ), x86_64-unknown-linux-gnu Locale: LC_CTYPE=en_US.UTF8, LC_NUMERIC=C, LC_TIME=en_US.UTF8, LC_COLLATE=en_US.UTF8, LC_MONETARY=en_US.UTF8, LC_MESSAGES=en_US.UTF8, LC_PAPER=C, LC_NAME=C, LC_ADDRESS=C, LC_TELEPHONE=C, LC_MEASUREMENT=en_US.UTF8, LC_IDENTIFICATION=C Base packages: base, datasets, graphics, grdevices, methods, stats, utils Other packages: colorout 1.0-0, qpgraph , setwidth 1.0-3, vimcom Loaded via a namespace (and not attached): annotate , AnnotationDbi , Biobase , BiocGenerics 0.6.0, DBI 0.2-7, genefilter , GGBase , graph , grid 3.0.0, IRanges , lattice , Matrix , mvtnorm , parallel 3.0.0, qtl , Rgraphviz 2.4.1, RSQLite , snpstats , splines 3.0.0, stats , survival , tools 3.0.0, XML , xtable Inma Tur, Alberto Roverato, Robert Castelo Simulation of molecular regulatory networks with graphical models 26 / 27
27 Acknowledgements and advertisements Funding: FPI predoctoral fellowship [BES ] to I. Tur Spanish MINECO project grant [TIN ] to R. Castelo The qpgraph package is available at Follow news and bugfixes about qpgraph in Twitter Thanks for listening! Inma Tur, Alberto Roverato, Robert Castelo Simulation of molecular regulatory networks with graphical models 27 / 27
How to Use pkgdeptools
How to Use pkgdeptools Seth Falcon February 10, 2018 1 Introduction The pkgdeptools package provides tools for computing and analyzing dependency relationships among R packages. With it, you can build
More informationAnalysis of two-way cell-based assays
Analysis of two-way cell-based assays Lígia Brás, Michael Boutros and Wolfgang Huber April 16, 2015 Contents 1 Introduction 1 2 Assembling the data 2 2.1 Reading the raw intensity files..................
More informationHow to use CNTools. Overview. Algorithms. Jianhua Zhang. April 14, 2011
How to use CNTools Jianhua Zhang April 14, 2011 Overview Studies have shown that genomic alterations measured as DNA copy number variations invariably occur across chromosomal regions that span over several
More informationHow To Use GOstats and Category to do Hypergeometric testing with unsupported model organisms
How To Use GOstats and Category to do Hypergeometric testing with unsupported model organisms M. Carlson October 30, 2017 1 Introduction This vignette is meant as an extension of what already exists in
More informationsurvsnp: Power and Sample Size Calculations for SNP Association Studies with Censored Time to Event Outcomes
survsnp: Power and Sample Size Calculations for SNP Association Studies with Censored Time to Event Outcomes Kouros Owzar Zhiguo Li Nancy Cox Sin-Ho Jung Chanhee Yi June 29, 2016 1 Introduction This vignette
More informationCount outlier detection using Cook s distance
Count outlier detection using Cook s distance Michael Love August 9, 2014 1 Run DE analysis with and without outlier removal The following vignette produces the Supplemental Figure of the effect of replacing
More informationHow to use cghmcr. October 30, 2017
How to use cghmcr Jianhua Zhang Bin Feng October 30, 2017 1 Overview Copy number data (arraycgh or SNP) can be used to identify genomic regions (Regions Of Interest or ROI) showing gains or losses that
More informationThe qp Package. December 21, 2006
The qp Package December 21, 2006 Type Package Title q-order partial correlation graph search algorithm Version 0.2-1 Date 2006-12-18 Author Robert Castelo , Alberto Roverato
More informationShrinkage of logarithmic fold changes
Shrinkage of logarithmic fold changes Michael Love August 9, 2014 1 Comparing the posterior distribution for two genes First, we run a DE analysis on the Bottomly et al. dataset, once with shrunken LFCs
More informationAnalysis of screens with enhancer and suppressor controls
Analysis of screens with enhancer and suppressor controls Lígia Brás, Michael Boutros and Wolfgang Huber April 16, 2015 Contents 1 Introduction 1 2 Assembling the data 2 2.1 Reading the raw intensity files..................
More informationMIGSA: Getting pbcmc datasets
MIGSA: Getting pbcmc datasets Juan C Rodriguez Universidad Católica de Córdoba Universidad Nacional de Córdoba Cristóbal Fresno Instituto Nacional de Medicina Genómica Andrea S Llera Fundación Instituto
More informationBioconductor: Annotation Package Overview
Bioconductor: Annotation Package Overview April 30, 2018 1 Overview In its current state the basic purpose of annotate is to supply interface routines that support user actions that rely on the different
More informationThe analysis of rtpcr data
The analysis of rtpcr data Jitao David Zhang, Markus Ruschhaupt October 30, 2017 With the help of this document, an analysis of rtpcr data can be performed. For this, the user has to specify several parameters
More informationPreprocessing and Genotyping Illumina Arrays for Copy Number Analysis
Preprocessing and Genotyping Illumina Arrays for Copy Number Analysis Rob Scharpf September 18, 2012 Abstract This vignette illustrates the steps required prior to copy number analysis for Infinium platforms.
More informationHowTo: Querying online Data
HowTo: Querying online Data Jeff Gentry and Robert Gentleman November 12, 2017 1 Overview This article demonstrates how you can make use of the tools that have been provided for on-line querying of data
More informationAnalysis of Pairwise Interaction Screens
Analysis of Pairwise Interaction Screens Bernd Fischer April 11, 2014 Contents 1 Introduction 1 2 Installation of the RNAinteract package 1 3 Creating an RNAinteract object 2 4 Single Perturbation effects
More informationPackage qpgraph. January 22, 2018
Package qpgraph January 22, 2018 Title Estimation of genetic and molecular regulatory networks from high-throughput genomics data Version 2.12.0 Estimate gene and eqtl networks from high-throughput expression
More informationImage Analysis with beadarray
Mike Smith October 30, 2017 Introduction From version 2.0 beadarray provides more flexibility in the processing of array images and the extraction of bead intensities than its predecessor. In the past
More informationIntroduction to the Codelink package
Introduction to the Codelink package Diego Diez October 30, 2018 1 Introduction This package implements methods to facilitate the preprocessing and analysis of Codelink microarrays. Codelink is a microarray
More informationcrlmm to downstream data analysis
crlmm to downstream data analysis VJ Carey, B Carvalho March, 2012 1 Running CRLMM on a nontrivial set of CEL files To use the crlmm algorithm, the user must load the crlmm package, as described below:
More informationAn Introduction to the methylumi package
An Introduction to the methylumi package Sean Davis and Sven Bilke October 13, 2014 1 Introduction Gene expression patterns are very important in understanding any biologic system. The regulation of gene
More informationDesign Issues in Matrix package Development
Design Issues in Matrix package Development Martin Maechler and Douglas Bates R Core Development Team maechler@stat.math.ethz.ch, bates@r-project.org Spring 2008 (typeset on November 16, 2017) Abstract
More informationIntroduction: microarray quality assessment with arrayqualitymetrics
Introduction: microarray quality assessment with arrayqualitymetrics Audrey Kauffmann, Wolfgang Huber April 4, 2013 Contents 1 Basic use 3 1.1 Affymetrix data - before preprocessing.......................
More informationDesign Issues in Matrix package Development
Design Issues in Matrix package Development Martin Maechler and Douglas Bates R Core Development Team maechler@stat.math.ethz.ch, bates@r-project.org Spring 2008 (typeset on October 26, 2018) Abstract
More informationCardinal design and development
Kylie A. Bemis October 30, 2017 Contents 1 Introduction.............................. 2 2 Design overview........................... 2 3 iset: high-throughput imaging experiments............ 3 3.1 SImageSet:
More informationsegmentseq: methods for detecting methylation loci and differential methylation
segmentseq: methods for detecting methylation loci and differential methylation Thomas J. Hardcastle October 13, 2015 1 Introduction This vignette introduces analysis methods for data from high-throughput
More informationComparing Least Squares Calculations
Comparing Least Squares Calculations Douglas Bates R Development Core Team Douglas.Bates@R-project.org November 16, 2017 Abstract Many statistics methods require one or more least squares problems to be
More informationTriform: peak finding in ChIP-Seq enrichment profiles for transcription factors
Triform: peak finding in ChIP-Seq enrichment profiles for transcription factors Karl Kornacker * and Tony Håndstad October 30, 2018 A guide for using the Triform algorithm to predict transcription factor
More informationsegmentseq: methods for detecting methylation loci and differential methylation
segmentseq: methods for detecting methylation loci and differential methylation Thomas J. Hardcastle October 30, 2018 1 Introduction This vignette introduces analysis methods for data from high-throughput
More informationSimBindProfiles: Similar Binding Profiles, identifies common and unique regions in array genome tiling array data
SimBindProfiles: Similar Binding Profiles, identifies common and unique regions in array genome tiling array data Bettina Fischer October 30, 2018 Contents 1 Introduction 1 2 Reading data and normalisation
More informationAnalyse RT PCR data with ddct
Analyse RT PCR data with ddct Jitao David Zhang, Rudolf Biczok and Markus Ruschhaupt October 30, 2017 Abstract Quantitative real time PCR (qrt PCR or RT PCR for short) is a laboratory technique based on
More informationAdvanced analysis using bayseq; generic distribution definitions
Advanced analysis using bayseq; generic distribution definitions Thomas J Hardcastle October 30, 2017 1 Generic Prior Distributions bayseq now offers complete user-specification of underlying distributions
More informationThe rtracklayer package
The rtracklayer package Michael Lawrence January 22, 2018 Contents 1 Introduction 2 2 Gene expression and microrna target sites 2 2.1 Creating a target site track..................... 2 2.1.1 Constructing
More informationWorking with aligned nucleotides (WORK- IN-PROGRESS!)
Working with aligned nucleotides (WORK- IN-PROGRESS!) Hervé Pagès Last modified: January 2014; Compiled: November 17, 2017 Contents 1 Introduction.............................. 1 2 Load the aligned reads
More informationThe bumphunter user s guide
The bumphunter user s guide Kasper Daniel Hansen khansen@jhsph.edu Martin Aryee aryee.martin@mgh.harvard.edu Rafael A. Irizarry rafa@jhu.edu Modified: November 23, 2012. Compiled: April 30, 2018 Introduction
More informationCreating a New Annotation Package using SQLForge
Creating a New Annotation Package using SQLForge Marc Carlson, Herve Pages, Nianhua Li November 19, 2013 1 Introduction The AnnotationForge package provides a series of functions that can be used to build
More informationRDAVIDWebService: a versatile R interface to DAVID
RDAVIDWebService: a versatile R interface to DAVID Cristóbal Fresno 1,2 and Elmer A. Fernández 1,2 1 Bio-science Data Mining Group, Catholic University of Córdoba, Córdoba, Argentina. 2 CONICET, Argentina
More informationSIBER User Manual. Pan Tong and Kevin R Coombes. May 27, Introduction 1
SIBER User Manual Pan Tong and Kevin R Coombes May 27, 2015 Contents 1 Introduction 1 2 Using SIBER 1 2.1 A Quick Example........................................... 1 2.2 Dealing With RNAseq Normalization................................
More informationmocluster: Integrative clustering using multiple
mocluster: Integrative clustering using multiple omics data Chen Meng Modified: June 03, 2015. Compiled: October 30, 2018. Contents 1 mocluster overview......................... 1 2 Run mocluster............................
More informationGetting started with flowstats
Getting started with flowstats F. Hahne, N. Gopalakrishnan October, 7 Abstract flowstats is a collection of algorithms for the statistical analysis of flow cytometry data. So far, the focus is on automated
More informationINSPEcT - INference of Synthesis, Processing and degradation rates in Time-course analysis
INSPEcT - INference of Synthesis, Processing and degradation rates in Time-course analysis Stefano de Pretis April 7, 20 Contents 1 Introduction 1 2 Quantification of Exon and Intron features 3 3 Time-course
More informationgenefu: a package for breast cancer gene expression analysis
genefu: a package for breast cancer gene expression analysis Benjamin Haibe-Kains 1, Markus Schröder 2, Christos Sotiriou 3, Gianluca Bontempi 4, and John Quackenbush 5,6 1 Bioinformatics and Computational
More informationEvaluation of VST algorithm in lumi package
Evaluation of VST algorithm in lumi package Pan Du 1, Simon Lin 1, Wolfgang Huber 2, Warrren A. Kibbe 1 December 22, 2010 1 Robert H. Lurie Comprehensive Cancer Center Northwestern University, Chicago,
More informationVisualisation, transformations and arithmetic operations for grouped genomic intervals
## Warning: replacing previous import ggplot2::position by BiocGenerics::Position when loading soggi Visualisation, transformations and arithmetic operations for grouped genomic intervals Thomas Carroll
More informationCreating a New Annotation Package using SQLForge
Creating a New Annotation Package using SQLForge Marc Carlson, HervÃľ PagÃĺs, Nianhua Li April 30, 2018 1 Introduction The AnnotationForge package provides a series of functions that can be used to build
More informationriboseqr Introduction Getting Data Workflow Example Thomas J. Hardcastle, Betty Y.W. Chung October 30, 2017
riboseqr Thomas J. Hardcastle, Betty Y.W. Chung October 30, 2017 Introduction Ribosome profiling extracts those parts of a coding sequence currently bound by a ribosome (and thus, are likely to be undergoing
More informationGxE.scan. October 30, 2018
GxE.scan October 30, 2018 Overview GxE.scan can process a GWAS scan using the snp.logistic, additive.test, snp.score or snp.matched functions, whereas snp.scan.logistic only calls snp.logistic. GxE.scan
More informationUsing crlmm for copy number estimation and genotype calling with Illumina platforms
Using crlmm for copy number estimation and genotype calling with Illumina platforms Rob Scharpf November, Abstract This vignette illustrates the steps necessary for obtaining marker-level estimates of
More informationAn Overview of the S4Vectors package
Patrick Aboyoun, Michael Lawrence, Hervé Pagès Edited: February 2018; Compiled: June 7, 2018 Contents 1 Introduction.............................. 1 2 Vector-like and list-like objects...................
More informationThe Command Pattern in R
The Command Pattern in R Michael Lawrence September 2, 2012 Contents 1 Introduction 1 2 Example pipelines 2 3 sessioninfo 8 1 Introduction Command pattern is a design pattern used in object-oriented programming,
More informationMultipleAlignment Objects
MultipleAlignment Objects Marc Carlson Bioconductor Core Team Fred Hutchinson Cancer Research Center Seattle, WA April 30, 2018 Contents 1 Introduction 1 2 Creation and masking 1 3 Analytic utilities 7
More informationMMDiff2: Statistical Testing for ChIP-Seq Data Sets
MMDiff2: Statistical Testing for ChIP-Seq Data Sets Gabriele Schwikert 1 and David Kuo 2 [1em] 1 The Informatics Forum, University of Edinburgh and The Wellcome
More informationOverview of ensemblvep
Overview of ensemblvep Lori Shepherd and Valerie Obenchain October 30, 2017 Contents 1 Introduction 1 2 Results as R objects 1 3 Write results to a file 4 4 Configuring runtime options 5 5 sessioninfo()
More informationAn Introduction to ShortRead
Martin Morgan Modified: 21 October, 2013. Compiled: April 30, 2018 > library("shortread") The ShortRead package provides functionality for working with FASTQ files from high throughput sequence analysis.
More informationCreating IGV HTML reports with tracktables
Thomas Carroll 1 [1em] 1 Bioinformatics Core, MRC Clinical Sciences Centre; thomas.carroll (at)imperial.ac.uk June 13, 2018 Contents 1 The tracktables package.................... 1 2 Creating IGV sessions
More informationAn Introduction to the methylumi package
An Introduction to the methylumi package Sean Davis and Sven Bilke April 30, 2018 1 Introduction Gene expression patterns are very important in understanding any biologic system. The regulation of gene
More informationSparse Model Matrices
Sparse Model Matrices Martin Maechler R Core Development Team maechler@r-project.org July 2007, 2008 (typeset on April 6, 2018) Introduction Model matrices in the very widely used (generalized) linear
More informationDupChecker: a bioconductor package for checking high-throughput genomic data redundancy in meta-analysis
DupChecker: a bioconductor package for checking high-throughput genomic data redundancy in meta-analysis Quanhu Sheng, Yu Shyr, Xi Chen Center for Quantitative Sciences, Vanderbilt University, Nashville,
More informationOverview of ensemblvep
Overview of ensemblvep Lori Shepherd and Valerie Obenchain October 30, 2018 Contents 1 Introduction 1 2 Results as R objects 1 3 Write results to a file 4 4 Configuring runtime options 5 5 sessioninfo()
More informationSeminar III: R/Bioconductor
Leonardo Collado Torres lcollado@lcg.unam.mx Bachelor in Genomic Sciences www.lcg.unam.mx/~lcollado/ August - December, 2009 1 / 50 Class outline Public Data Intro biomart GEOquery ArrayExpress annotate
More informationStep-by-Step Guide to Advanced Genetic Analysis
Step-by-Step Guide to Advanced Genetic Analysis Page 1 Introduction In the previous document, 1 we covered the standard genetic analyses available in JMP Genomics. Here, we cover the more advanced options
More informationmatter: Rapid prototyping with data on disk
matter: Rapid prototyping with data on disk Kylie A. Bemis March 17, 2017 Contents 1 Introduction 1 2 Installation 1 3 Basic use and data manipulation 1 4 Linear regression for on-disk datasets 5 5 Principal
More informationClassification by Nearest Shrunken Centroids and Support Vector Machines
Classification by Nearest Shrunken Centroids and Support Vector Machines Florian Markowetz florian.markowetz@molgen.mpg.de Max Planck Institute for Molecular Genetics, Computational Diagnostics Group,
More informationRaman Spectra of Chondrocytes in Cartilage: hyperspec s chondro data set
Raman Spectra of Chondrocytes in Cartilage: hyperspec s chondro data set Claudia Beleites CENMAT and DI3, University of Trieste Spectroscopy Imaging, IPHT Jena e.v. February 13,
More informationCQN (Conditional Quantile Normalization)
CQN (Conditional Quantile Normalization) Kasper Daniel Hansen khansen@jhsph.edu Zhijin Wu zhijin_wu@brown.edu Modified: August 8, 2012. Compiled: April 30, 2018 Introduction This package contains the CQN
More informationMLSeq package: Machine Learning Interface to RNA-Seq Data
MLSeq package: Machine Learning Interface to RNA-Seq Data Gokmen Zararsiz 1, Dincer Goksuluk 2, Selcuk Korkmaz 2, Vahap Eldem 3, Izzet Parug Duru 4, Turgay Unver 5, Ahmet Ozturk 5 1 Erciyes University,
More informationComputer Vision Group Prof. Daniel Cremers. 4. Probabilistic Graphical Models Directed Models
Prof. Daniel Cremers 4. Probabilistic Graphical Models Directed Models The Bayes Filter (Rep.) (Bayes) (Markov) (Tot. prob.) (Markov) (Markov) 2 Graphical Representation (Rep.) We can describe the overall
More informationIntroduction to R (BaRC Hot Topics)
Introduction to R (BaRC Hot Topics) George Bell September 30, 2011 This document accompanies the slides from BaRC s Introduction to R and shows the use of some simple commands. See the accompanying slides
More informationIntroduction: microarray quality assessment with arrayqualitymetrics
: microarray quality assessment with arrayqualitymetrics Audrey Kauffmann, Wolfgang Huber October 30, 2017 Contents 1 Basic use............................... 3 1.1 Affymetrix data - before preprocessing.............
More informationSeminar III: R/Bioconductor
Leonardo Collado Torres lcollado@lcg.unam.mx Bachelor in Genomic Sciences www.lcg.unam.mx/~lcollado/ August - December, 2009 1 / 25 Class outline Working with HTS data: a simulated case study Intro R for
More informationNacoStringQCPro. Dorothee Nickles, Thomas Sandmann, Robert Ziman, Richard Bourgon. Modified: April 10, Compiled: April 24, 2017.
NacoStringQCPro Dorothee Nickles, Thomas Sandmann, Robert Ziman, Richard Bourgon Modified: April 10, 2017. Compiled: April 24, 2017 Contents 1 Scope 1 2 Quick start 1 3 RccSet loaded with NanoString data
More informationBiobase development and the new eset
Biobase development and the new eset Martin T. Morgan * 7 August, 2006 Revised 4 September, 2006 featuredata slot. Revised 20 April 2007 minor wording changes; verbose and other arguments passed through
More informationA complete analysis of peptide microarray binding data using the pepstat framework
A complete analysis of peptide microarray binding data using the pepstat framework Greg Imholte, Renan Sauteraud, Mike Jiang and Raphael Gottardo April 30, 2018 This document present a full analysis, from
More informationUsing the SRAdb Package to Query the Sequence Read Archive
Using the SRAdb Package to Query the Sequence Read Archive Jack Zhu and Sean Davis Genetics Branch, Center for Cancer Research, National Cancer Institute, National Institutes of Health January 24, 2013
More informationIntroduction to GenomicFiles
Valerie Obenchain, Michael Love, Martin Morgan Last modified: October 2014; Compiled: October 30, 2018 Contents 1 Introduction.............................. 1 2 Quick Start..............................
More informationRAPIDR. Kitty Lo. November 20, Intended use of RAPIDR 1. 2 Create binned counts file from BAMs Masking... 1
RAPIDR Kitty Lo November 20, 2014 Contents 1 Intended use of RAPIDR 1 2 Create binned counts file from BAMs 1 2.1 Masking.................................................... 1 3 Build the reference 2 3.1
More informationLearning DAGs from observational data
Learning DAGs from observational data General overview Introduction DAGs and conditional independence DAGs and causal effects Learning DAGs from observational data IDA algorithm Further problems 2 What
More informationBasic4Cseq: an R/Bioconductor package for the analysis of 4C-seq data
Basic4Cseq: an R/Bioconductor package for the analysis of 4C-seq data Carolin Walter April 16, 2015 Contents 1 Introduction 1 1.1 Loading the package...................................... 2 1.2 Provided
More informationBasic4Cseq: an R/Bioconductor package for the analysis of 4C-seq data
Basic4Cseq: an R/Bioconductor package for the analysis of 4C-seq data Carolin Walter October 30, 2017 Contents 1 Introduction 1 1.1 Loading the package...................................... 2 1.2 Provided
More informationhighscreen: High-Throughput Screening for Plate Based Assays
highscreen: High-Throughput Screening for Plate Based Assays Ivo D. Shterev Cliburn Chan Gregory D. Sempowski August 22, 2018 1 Introduction This vignette describes the use of the R extension package highscreen
More informationThe fastclime Package for Linear Programming and Large-Scale Precision Matrix Estimation in R
Journal of Machine Learning Research (2013) Submitted ; Published The fastclime Package for Linear Programming and Large-Scale Precision Matrix Estimation in R Haotian Pang Han Liu Robert Vanderbei Princeton
More informationMissing Data Analysis for the Employee Dataset
Missing Data Analysis for the Employee Dataset 67% of the observations have missing values! Modeling Setup For our analysis goals we would like to do: Y X N (X, 2 I) and then interpret the coefficients
More informationUsing metama for differential gene expression analysis from multiple studies
Using metama for differential gene expression analysis from multiple studies Guillemette Marot and Rémi Bruyère Modified: January 28, 2015. Compiled: January 28, 2015 Abstract This vignette illustrates
More informationINSPEcT - INference of Synthesis, Processing and degradation rates in Time-course analysis
INSPEcT - INference of Synthesis, Processing and degradation rates in Time-course analysis Stefano de Pretis October 30, 2017 Contents 1 Introduction.............................. 1 2 Quantification of
More informationA spatio-temporal model for extreme precipitation simulated by a climate model.
A spatio-temporal model for extreme precipitation simulated by a climate model. Jonathan Jalbert Joint work with Anne-Catherine Favre, Claude Bélisle and Jean-François Angers STATMOS Workshop: Climate
More informationbayseq: Empirical Bayesian analysis of patterns of differential expression in count data
bayseq: Empirical Bayesian analysis of patterns of differential expression in count data Thomas J. Hardcastle October 30, 2017 1 Introduction This vignette is intended to give a rapid introduction to the
More informationPROPER: PROspective Power Evaluation for RNAseq
PROPER: PROspective Power Evaluation for RNAseq Hao Wu [1em]Department of Biostatistics and Bioinformatics Emory University Atlanta, GA 303022 [1em] hao.wu@emory.edu October 30, 2017 Contents 1 Introduction..............................
More informationAn Introduction to Bioconductor s ExpressionSet Class
An Introduction to Bioconductor s ExpressionSet Class Seth Falcon, Martin Morgan, and Robert Gentleman 6 October, 2006; revised 9 February, 2007 1 Introduction Biobase is part of the Bioconductor project,
More informationThe 2-core of a Non-homogeneous Hypergraph
July 16, 2012 k-cores A Hypergraph G on vertex set V is a collection E of subsets of V. E is the set of hyperedges. For ordinary graphs, e = 2 for all e E. The k-core of a (hyper)graph is the maximal subgraph
More informationUsing the dorng package
Using the dorng package dorng package Version 1.6 Renaud Gaujoux March 5, 2014 Contents Introduction............. 1 1 The %dorng% operator...... 3 1.1 How it works......... 3 1.2 Seeding computations.....
More informationNetwork Traffic Measurements and Analysis
DEIB - Politecnico di Milano Fall, 2017 Sources Hastie, Tibshirani, Friedman: The Elements of Statistical Learning James, Witten, Hastie, Tibshirani: An Introduction to Statistical Learning Andrew Ng:
More informationAn Introduction to GenomeInfoDb
Martin Morgan, Hervé Pagès, Marc Carlson, Sonali Arora Modified: 17 January, 2014. Compiled: January 21, 2018 Contents 1 Introduction.............................. 2 2 Functionality for all existing organisms...............
More informationBig Data Analytics. Special Topics for Computer Science CSE CSE Feb 11
Big Data Analytics Special Topics for Computer Science CSE 4095-001 CSE 5095-005 Feb 11 Fei Wang Associate Professor Department of Computer Science and Engineering fei_wang@uconn.edu Clustering II Spectral
More informationComputer Vision Group Prof. Daniel Cremers. 4. Probabilistic Graphical Models Directed Models
Prof. Daniel Cremers 4. Probabilistic Graphical Models Directed Models The Bayes Filter (Rep.) (Bayes) (Markov) (Tot. prob.) (Markov) (Markov) 2 Graphical Representation (Rep.) We can describe the overall
More informationNested Loop Cross Validation for Classification using nlcv
Nested Loop Cross Validation for Classification using nlcv Willem Talloen February 15, 2013 1 Introduction Microarrays may provide clinicians with valuable biomarkers for disease status or treatment sensitivity.
More informationNotes on QTL Cartographer
Notes on QTL Cartographer Introduction QTL Cartographer is a suite of programs for mapping quantitative trait loci (QTLs) onto a genetic linkage map. The programs use linear regression, interval mapping
More informationThe fastliquidassociation Package
The fastliquidassociation Package Tina Gunderson October 30, 2018 1 Introduction The fastliquidassociation package is intended to provide more efficient tools for genome-wide application of liquid association.
More informationClustering K-means. Machine Learning CSEP546 Carlos Guestrin University of Washington February 18, Carlos Guestrin
Clustering K-means Machine Learning CSEP546 Carlos Guestrin University of Washington February 18, 2014 Carlos Guestrin 2005-2014 1 Clustering images Set of Images [Goldberger et al.] Carlos Guestrin 2005-2014
More informationAn introduction to gaucho
An introduction to gaucho Alex Murison Alexander.Murison@icr.ac.uk Christopher P Wardell Christopher.Wardell@icr.ac.uk October 30, 2017 This vignette serves as an introduction to the R package gaucho.
More informationThe Genominator User Guide
The Genominator User Guide James Bullard Kasper D. Hansen Modified: April 18, 2010, Compiled: April 3, 2013 > ## -- For Debugging / Timing > ## require(rsqlite) > ## source("../../r/genominator.r") > ##
More information