Last data update: 2014.03.03

R: Clustered Image Maps (CIMs) (heat maps)
plotcimR Documentation

Clustered Image Maps (CIMs) (heat maps)

Description

This function generates color-coded Clustered Image Maps (CIMs) representing the pairwise correlation between a subset of predictors and the outcomes.

Usage

  plotcim(x, select.variable, labelX = NULL)

Arguments

x

an object of class ESS.

select.variable

name of the selected predictors.

labelX

names of the predictors used when running R2GUESS. This argument is optional if a MAP.file has been provided.

Author(s)

Benoit Liquet, b.liquet@uq.edu.au,
Marc Chadeau-Hyam m.chadeau@imperial.ac.uk,
Leonardo Bottolo l.bottolo@imperial.ac.uk,
Gianluca Campanella g.campanella11@imperial.ac.uk

Examples

modelY_Hopx <- example.as.ESS.object()
MPPI.Hopx <-
plotMPPI(modelY_Hopx,threshold.model=20,threshold.variable=0.25,Figure=FALSE)
plotcim(modelY_Hopx,select.variable=MPPI.Hopx$var.TOP.MPI)

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(R2GUESS)
Loading required package: fields
Loading required package: spam
Loading required package: grid
Spam version 1.3-0 (2015-10-24) is loaded.
Type 'help( Spam)' or 'demo( spam)' for a short introduction 
and overview of this package.
Help for individual functions is also obtained by adding the
suffix '.spam' to the function name, e.g. 'help( chol.spam)'.

Attaching package: 'spam'

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

    backsolve, forwardsolve

Loading required package: maps

 # maps v3.1: updated 'world': all lakes moved to separate new #
 # 'lakes' database. Type '?world' or 'news(package="maps")'.  #


Loading required package: MCMCpack
Loading required package: coda
Loading required package: MASS
##
## Markov Chain Monte Carlo Package (MCMCpack)
## Copyright (C) 2003-2016 Andrew D. Martin, Kevin M. Quinn, and Jong Hee Park
##
## Support provided by the U.S. National Science Foundation
## (Grants SES-0350646 and SES-0350613)
##
Loading required package: mixOmics
Loading required package: lattice
Loading required package: ggplot2

Loaded mixOmics 6.0.0

Visit http://www.mixOmics.org for more details about our methods.
Any bug reports or comments? Notify us at mixomics at math.univ-toulouse.fr or https://bitbucket.org/klecao/package-mixomics/issues

Thank you for using mixOmics!

Attaching package: 'mixOmics'

The following object is masked from 'package:maps':

    map

Loading required package: mvtnorm
Loading required package: snowfall
Loading required package: snow
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/R2GUESS/plotcim.Rd_%03d_medium.png", width=480, height=480)
> ### Name: plotcim
> ### Title: Clustered Image Maps (CIMs) (heat maps)
> ### Aliases: plotcim
> 
> ### ** Examples
> 
> modelY_Hopx <- example.as.ESS.object()
The run is ok 
You can now analyse the results 
> MPPI.Hopx <-
+ plotMPPI(modelY_Hopx,threshold.model=20,threshold.variable=0.25,Figure=FALSE)
> plotcim(modelY_Hopx,select.variable=MPPI.Hopx$var.TOP.MPI)
Read 22330 items
Read 116 items
$mat.sim
           [,1]        [,2]        [,3]        [,4]        [,5]        [,6]
[1,]  0.2490801  0.02647336 -0.10415155 -0.55957534  0.22833102  0.06747719
[2,]  0.1787624 -0.04225032 -0.07122495 -0.17350509  0.53617164 -0.52150815
[3,] -0.3751483  0.46440155  0.02093906 -0.06100608 -0.03048752  0.14129624
[4,]  0.1615677 -0.08140259 -0.25226534 -0.35580078  0.22630945 -0.04989970
          [,7]
[1,] 0.4785538
[2,] 0.2688677
[3,] 0.9002820
[4,] 0.4827163

$labelY
[1] "ADR"    "Fat"    "Heart"  "Kidney"

$labelX
[1] "D2Cebr204s17" "D3Mit16"      "D6Cebrp40s27" "D10Rat33"     "Dcp1"        
[6] "D11Mit4"      "D14Mit3"     

> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>