Last data update: 2014.03.03

R: Profile plots for ISA biclusters
profilePlotR Documentation

Profile plots for ISA biclusters

Description

Line plots to compare biclusters to the background, i.e. the rest of the expression matrix.

Usage

profilePlot (modules, module, eset, plot = c("samples", "features", 
     "both"), norm = "default", background = TRUE, 
    col = gray(0.7), col.mod = 1, type = "l", type.mod = type,
    mean = TRUE, meancol = "green", meancol.mod = "red",
    xlabs = c("Features", "Samples"), ylab = "Expression", 
    ...) 

Arguments

modules

An ISAModules object.

module

Numeric scalar, the module to plot.

eset

An ExpressionSet or ISAExpressionSet object. If an ExpressionSet object is supplied (and the norm argument is not set to ‘raw’), then it is normalised by calling ISANormalize on it. A subset of eset is selected that corresponds to the features included in modules.

plot

Character constant, specifies what to plot. ‘sample’ plots sample scores, ‘features’ plots feature scores. If ‘both’ is given, then the plot is divided into two subplots and both scores are plotted.

norm

Character constant, specifies how to normalize the expression matrix for plotting. It can be of length one or two, the latter for the case when plots are made both for features and samples. Possible values: ‘raw’ uses the raw expression values; ‘feature’ uses featExprs to extract the expression values from the expression set object; ‘sample’ uses sampExprs; ‘default’ means ‘feature’ for sample plots and ‘sample’ for feature plots.

background

Logical scalar, whether to plot the features/samples that are not in the module.

col

Color of lines corresponding to the background features/samples.

col.mod

Color of the lines corresponding to the features/samples included in the module.

type

Type of the plot, for the background features/samples. It is passed to plot.

type.mod

Type of the plot, for the features/samples included in the module. It is passed to plot.

mean

Logical scalar, whether to plot the mean expression for each feature/sample, separately for the samples/features that are in the module and the ones that are not.

meancol

Color of the line for the mean expression values, background.

meancol.mod

Color of the line for the mean expression values, module.

xlabs

Character vector of length one or two. The labels of the horizontal axes of the plot, the second value is used if both the feature and the sample plots are drawn.

ylab

Character vector of length one. The label of the vertical axes.

...

Additional graphical arguments. They are passed to the lines function that creates the lines of the plot.

Details

plot="both" uses the mfrow graphical parameter to create the two subplots. This does not work properly if you already have subplots.

Value

None. (Well, NULL, invisibly.)

Author(s)

Gabor Csardi csardi.gabor@gmail.com

References

Bergmann S, Ihmels J, Barkai N: Iterative signature algorithm for the analysis of large-scale gene expression data Phys Rev E Stat Nonlin Soft Matter Phys. 2003 Mar;67(3 Pt 1):031902. Epub 2003 Mar 11.

See Also

The similar parallelCoordinates function in the biclust package.

Examples

data(ALLModulesSmall)
library(ALL)
data(ALL)
if (interactive()) {
  profilePlot(ALLModulesSmall, 2, ALL, plot="samples")
}

Results


R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(eisa)
Loading required package: isa2
Loading required package: Biobase
Loading required package: BiocGenerics
Loading required package: parallel

Attaching package: 'BiocGenerics'

The following objects are masked from 'package:parallel':

    clusterApply, clusterApplyLB, clusterCall, clusterEvalQ,
    clusterExport, clusterMap, parApply, parCapply, parLapply,
    parLapplyLB, parRapply, parSapply, parSapplyLB

The following objects are masked from 'package:stats':

    IQR, mad, xtabs

The following objects are masked from 'package:base':

    Filter, Find, Map, Position, Reduce, anyDuplicated, append,
    as.data.frame, cbind, colnames, do.call, duplicated, eval, evalq,
    get, grep, grepl, intersect, is.unsorted, lapply, lengths, mapply,
    match, mget, order, paste, pmax, pmax.int, pmin, pmin.int, rank,
    rbind, rownames, sapply, setdiff, sort, table, tapply, union,
    unique, unsplit

Welcome to Bioconductor

    Vignettes contain introductory material; view with
    'browseVignettes()'. To cite Bioconductor, see
    'citation("Biobase")', and for packages 'citation("pkgname")'.

Loading required package: AnnotationDbi
Loading required package: stats4
Loading required package: IRanges
Loading required package: S4Vectors

Attaching package: 'S4Vectors'

The following objects are masked from 'package:base':

    colMeans, colSums, expand.grid, rowMeans, rowSums

> png(filename="/home/ddbj/snapshot/RGM3/R_BC/result/eisa/profilePlot.Rd_%03d_medium.png", width=480, height=480)
> ### Name: profilePlot
> ### Title: Profile plots for ISA biclusters
> ### Aliases: profilePlot
> ### Keywords: cluster
> 
> ### ** Examples
> 
> data(ALLModulesSmall)
> library(ALL)
> data(ALL)
> #if (interactive()) {
>   profilePlot(ALLModulesSmall, 2, ALL, plot="samples")
> #}
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>