heatMapMeth: HeatMap diagram over methylation data

In methVisual: Methods for visualization and statistics on DNA methylation data

Description

HeatMap diagram of methylation classes over CpG position and sequence level

Usage

heatMapMeth(methData,file)

Arguments

methData	List; contains information on the pairwise alignments, and methylated CpG motifs
file	optionally, quoted character string for specification of path and file name for saving the result. By default, the result file's format is .pdf. If argument is omitted, screen output is provided only.

Details

Clustering is a prominent method for visualizing and studying groups of similar features, and is also widely used for the analysis of microarrays. In the case of analyzing DNA methylation datasets the matrix to be explored is a I*J binary matrix, were I are the clone sequences and J are the CpG positions. Every index in this matrix has the value 1 for methylated or 0 for non methylated CpG sites. Providing a hierarchical clustering option based on the quality checked methylation data can be useful in finding clusters in the explored data in two dimensions, of the methylation state of CpG sites (J) and for distribution of methylation state over explored clone sequences (I). Very importantly, one has to keep in mind, that this method does not take into account the genomic ordering of CpG sites. Clustering methods can be also useful in quality control check. Observing clones from the same PCR product in different clusters can point to bad quality clone sequences. The clustering method which was used in methVisual R package is the heatmap() function from stat package which is a color image with two dendrograms added to the sides of the columns and rows. Since the data is binary, the distances that are calculated among the columns and among the rows are computed with a binary distance function which is the asymmetric binary function. This method assumes that non zero elements are ON and zero elements are OFF. The distance is the proportion of bits in which only one is ON amongst those in which at least one is ON. Based on aligned sequences under study a heatmap is created that displays two way clustering of methylation status of all sequences and all aligned CpG positions.

Value

Heat-Map image is displayed and optionally saved as postscript in given path and name.

Author(s)

Arie Zackay <[e-mail]>,Christine Steinhoff <[e-mail]>

Examples

## using methData, file is the path to R home directory.
## In order to save heatMapMeth.pdf, make sure that you have writing 
## permission under R.home() directory. If you do not have permission
## choose your own path.  
#dir.create(file.path(R.home(component="home"),"/BiqAnalyzer"))
BiqAnalyzer_path <- file.path(tempdir(), "BiqAnalyzer")
dir.create(BiqAnalyzer_path) 
data(methData)
heatMapMeth(methData,file=file.path(BiqAnalyzer_path, "heatMapMeth.pdf"))

Example output

Loading required package: Biostrings
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, sd, var, xtabs

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

    Filter, Find, Map, Position, Reduce, anyDuplicated, append,
    as.data.frame, cbind, colMeans, colSums, 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, rowMeans, rowSums, rownames, sapply,
    setdiff, sort, table, tapply, union, unique, unsplit, which,
    which.max, which.min

Loading required package: S4Vectors
Loading required package: stats4

Attaching package: 'S4Vectors'

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

    expand.grid

Loading required package: IRanges
Loading required package: XVector

Attaching package: 'Biostrings'

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

    strsplit

Loading required package: plotrix
Loading required package: gsubfn
Loading required package: proto
Could not load tcltk.  Will use slower R code instead.
Loading required package: grid
Loading required package: sqldf
Loading required package: RSQLite
  On_Plot_Nr     Experiment
1          1 QC_seq_B.fasta
2          2 QC_seq_C.fasta
3          3 QC_seq_E.fasta
4          4 QC_seq_G.fasta
5          5 QC_seq_H.fasta
6          6 QC_seq_I.fasta

methVisual documentation built on April 28, 2020, 7:08 p.m.