matrixSNP: Correlation between methylation states

Description Usage Arguments Details Value Author(s) Examples

Description

Correlation between methylation states for each CpG position and each sequence under study

Usage

1
matrixSNP(methData,correlation)

Arguments

methData

List; contains information on the pairwise alignments, and methylated CpG motifs.

correlation

calculation of corraltion on Matrix, if TRUE: calculation of spearman correlation matrix. if FALSE: computes a matrix with absolut numberes of methylation over all CpG position given methylation in a certain CpG sites.

Details

The function enable the user to explore cooccurrence between non neighbored CpG sites, which can be made by calculating all pairwise cooccurrences, due to correlation over all methylation CpG sites. Furthermore it is possible to display the computed distance cooccurrences in a plot containing all CpG sites and their correlation with other CpG sites. The correlation values are color coded (gray levels) and the color coding bar is given beside the graph. The red numbers in the diagonal give the genomic position of each displayed CpG site.

Value

Returns a correlation matrix displaying dependencies values for all CpG positions

Author(s)

Arie Zackay <arie.zackay@mail.huji.ac.il>, Christine Steinhoff <steinhof@molgen.mpg.de>

Examples

1
2
3
## using methData
data(methData)
matrixSNP(methData)

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, basename, cbind, colMeans, colSums, colnames,
    dirname, 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
Loading required package: grid
Loading required package: sqldf
Loading required package: RSQLite
Warning message:
no DISPLAY variable so Tk is not available 
      [,1]  [,2] [,3]  [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13]
 [1,] 1.00  0.63   NA  0.63   NA 1.00   NA   NA 1.00  0.63  0.25  1.00  1.00
 [2,] 0.63  1.00   NA  1.00   NA 0.63   NA   NA 0.63  1.00 -0.32  0.63  0.63
 [3,]   NA    NA   NA    NA   NA   NA   NA   NA   NA    NA    NA    NA    NA
 [4,] 0.63  1.00   NA  1.00   NA 0.63   NA   NA 0.63  1.00 -0.32  0.63  0.63
 [5,]   NA    NA   NA    NA   NA   NA   NA   NA   NA    NA    NA    NA    NA
 [6,] 1.00  0.63   NA  0.63   NA 1.00   NA   NA 1.00  0.63  0.25  1.00  1.00
 [7,]   NA    NA   NA    NA   NA   NA   NA   NA   NA    NA    NA    NA    NA
 [8,]   NA    NA   NA    NA   NA   NA   NA   NA   NA    NA    NA    NA    NA
 [9,] 1.00  0.63   NA  0.63   NA 1.00   NA   NA 1.00  0.63  0.25  1.00  1.00
[10,] 0.63  1.00   NA  1.00   NA 0.63   NA   NA 0.63  1.00 -0.32  0.63  0.63
[11,] 0.25 -0.32   NA -0.32   NA 0.25   NA   NA 0.25 -0.32  1.00  0.25  0.25
[12,] 1.00  0.63   NA  0.63   NA 1.00   NA   NA 1.00  0.63  0.25  1.00  1.00
[13,] 1.00  0.63   NA  0.63   NA 1.00   NA   NA 1.00  0.63  0.25  1.00  1.00
[14,] 1.00  0.63   NA  0.63   NA 1.00   NA   NA 1.00  0.63  0.25  1.00  1.00
[15,] 0.71  0.45   NA  0.45   NA 0.71   NA   NA 0.71  0.45  0.71  0.71  0.71
[16,] 0.71  0.45   NA  0.45   NA 0.71   NA   NA 0.71  0.45  0.71  0.71  0.71
[17,] 0.71  0.45   NA  0.45   NA 0.71   NA   NA 0.71  0.45  0.71  0.71  0.71
[18,] 0.71  0.45   NA  0.45   NA 0.71   NA   NA 0.71  0.45  0.71  0.71  0.71
[19,] 0.32  0.20   NA  0.20   NA 0.32   NA   NA 0.32  0.20  0.32  0.32  0.32
[20,] 1.00  0.63   NA  0.63   NA 1.00   NA   NA 1.00  0.63  0.25  1.00  1.00
[21,] 0.63  1.00   NA  1.00   NA 0.63   NA   NA 0.63  1.00 -0.32  0.63  0.63
[22,] 0.71  0.45   NA  0.45   NA 0.71   NA   NA 0.71  0.45  0.71  0.71  0.71
[23,] 0.25 -0.32   NA -0.32   NA 0.25   NA   NA 0.25 -0.32  1.00  0.25  0.25
[24,] 0.25 -0.32   NA -0.32   NA 0.25   NA   NA 0.25 -0.32  1.00  0.25  0.25
[25,] 0.63 -0.20   NA -0.20   NA 0.63   NA   NA 0.63 -0.20  0.63  0.63  0.63
[26,]   NA    NA   NA    NA   NA   NA   NA   NA   NA    NA    NA    NA    NA
[27,] 1.00  0.63   NA  0.63   NA 1.00   NA   NA 1.00  0.63  0.25  1.00  1.00
[28,] 0.71  0.45   NA  0.45   NA 0.71   NA   NA 0.71  0.45  0.71  0.71  0.71
[29,] 0.63 -0.20   NA -0.20   NA 0.63   NA   NA 0.63 -0.20  0.63  0.63  0.63
[30,] 0.25  0.63   NA  0.63   NA 0.25   NA   NA 0.25  0.63  0.25  0.25  0.25
[31,] 0.63  1.00   NA  1.00   NA 0.63   NA   NA 0.63  1.00 -0.32  0.63  0.63
[32,]   NA    NA   NA    NA   NA   NA   NA   NA   NA    NA    NA    NA    NA
[33,] 0.63 -0.20   NA -0.20   NA 0.63   NA   NA 0.63 -0.20  0.63  0.63  0.63
[34,] 1.00  0.63   NA  0.63   NA 1.00   NA   NA 1.00  0.63  0.25  1.00  1.00
[35,] 0.63 -0.20   NA -0.20   NA 0.63   NA   NA 0.63 -0.20  0.63  0.63  0.63
      [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25]
 [1,]  1.00  0.71  0.71  0.71  0.71  0.32  1.00  0.63  0.71  0.25  0.25  0.63
 [2,]  0.63  0.45  0.45  0.45  0.45  0.20  0.63  1.00  0.45 -0.32 -0.32 -0.20
 [3,]    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA
 [4,]  0.63  0.45  0.45  0.45  0.45  0.20  0.63  1.00  0.45 -0.32 -0.32 -0.20
 [5,]    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA
 [6,]  1.00  0.71  0.71  0.71  0.71  0.32  1.00  0.63  0.71  0.25  0.25  0.63
 [7,]    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA
 [8,]    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA
 [9,]  1.00  0.71  0.71  0.71  0.71  0.32  1.00  0.63  0.71  0.25  0.25  0.63
[10,]  0.63  0.45  0.45  0.45  0.45  0.20  0.63  1.00  0.45 -0.32 -0.32 -0.20
[11,]  0.25  0.71  0.71  0.71  0.71  0.32  0.25 -0.32  0.71  1.00  1.00  0.63
[12,]  1.00  0.71  0.71  0.71  0.71  0.32  1.00  0.63  0.71  0.25  0.25  0.63
[13,]  1.00  0.71  0.71  0.71  0.71  0.32  1.00  0.63  0.71  0.25  0.25  0.63
[14,]  1.00  0.71  0.71  0.71  0.71  0.32  1.00  0.63  0.71  0.25  0.25  0.63
[15,]  0.71  1.00  1.00  1.00  1.00  0.45  0.71  0.45  1.00  0.71  0.71  0.45
[16,]  0.71  1.00  1.00  1.00  1.00  0.45  0.71  0.45  1.00  0.71  0.71  0.45
[17,]  0.71  1.00  1.00  1.00  1.00  0.45  0.71  0.45  1.00  0.71  0.71  0.45
[18,]  0.71  1.00  1.00  1.00  1.00  0.45  0.71  0.45  1.00  0.71  0.71  0.45
[19,]  0.32  0.45  0.45  0.45  0.45  1.00  0.32  0.20  0.45  0.32  0.32  0.20
[20,]  1.00  0.71  0.71  0.71  0.71  0.32  1.00  0.63  0.71  0.25  0.25  0.63
[21,]  0.63  0.45  0.45  0.45  0.45  0.20  0.63  1.00  0.45 -0.32 -0.32 -0.20
[22,]  0.71  1.00  1.00  1.00  1.00  0.45  0.71  0.45  1.00  0.71  0.71  0.45
[23,]  0.25  0.71  0.71  0.71  0.71  0.32  0.25 -0.32  0.71  1.00  1.00  0.63
[24,]  0.25  0.71  0.71  0.71  0.71  0.32  0.25 -0.32  0.71  1.00  1.00  0.63
[25,]  0.63  0.45  0.45  0.45  0.45  0.20  0.63 -0.20  0.45  0.63  0.63  1.00
[26,]    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA
[27,]  1.00  0.71  0.71  0.71  0.71  0.32  1.00  0.63  0.71  0.25  0.25  0.63
[28,]  0.71  1.00  1.00  1.00  1.00  0.45  0.71  0.45  1.00  0.71  0.71  0.45
[29,]  0.63  0.45  0.45  0.45  0.45  0.20  0.63 -0.20  0.45  0.63  0.63  1.00
[30,]  0.25  0.71  0.71  0.71  0.71  0.32  0.25  0.63  0.71  0.25  0.25 -0.32
[31,]  0.63  0.45  0.45  0.45  0.45  0.20  0.63  1.00  0.45 -0.32 -0.32 -0.20
[32,]    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA
[33,]  0.63  0.45  0.45  0.45  0.45  0.20  0.63 -0.20  0.45  0.63  0.63  1.00
[34,]  1.00  0.71  0.71  0.71  0.71  0.32  1.00  0.63  0.71  0.25  0.25  0.63
[35,]  0.63  0.45  0.45  0.45  0.45  0.20  0.63 -0.20  0.45  0.63  0.63  1.00
      [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35]
 [1,]    NA  1.00  0.71  0.63  0.25  0.63    NA  0.63  1.00  0.63
 [2,]    NA  0.63  0.45 -0.20  0.63  1.00    NA -0.20  0.63 -0.20
 [3,]    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA
 [4,]    NA  0.63  0.45 -0.20  0.63  1.00    NA -0.20  0.63 -0.20
 [5,]    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA
 [6,]    NA  1.00  0.71  0.63  0.25  0.63    NA  0.63  1.00  0.63
 [7,]    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA
 [8,]    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA
 [9,]    NA  1.00  0.71  0.63  0.25  0.63    NA  0.63  1.00  0.63
[10,]    NA  0.63  0.45 -0.20  0.63  1.00    NA -0.20  0.63 -0.20
[11,]    NA  0.25  0.71  0.63  0.25 -0.32    NA  0.63  0.25  0.63
[12,]    NA  1.00  0.71  0.63  0.25  0.63    NA  0.63  1.00  0.63
[13,]    NA  1.00  0.71  0.63  0.25  0.63    NA  0.63  1.00  0.63
[14,]    NA  1.00  0.71  0.63  0.25  0.63    NA  0.63  1.00  0.63
[15,]    NA  0.71  1.00  0.45  0.71  0.45    NA  0.45  0.71  0.45
[16,]    NA  0.71  1.00  0.45  0.71  0.45    NA  0.45  0.71  0.45
[17,]    NA  0.71  1.00  0.45  0.71  0.45    NA  0.45  0.71  0.45
[18,]    NA  0.71  1.00  0.45  0.71  0.45    NA  0.45  0.71  0.45
[19,]    NA  0.32  0.45  0.20  0.32  0.20    NA  0.20  0.32  0.20
[20,]    NA  1.00  0.71  0.63  0.25  0.63    NA  0.63  1.00  0.63
[21,]    NA  0.63  0.45 -0.20  0.63  1.00    NA -0.20  0.63 -0.20
[22,]    NA  0.71  1.00  0.45  0.71  0.45    NA  0.45  0.71  0.45
[23,]    NA  0.25  0.71  0.63  0.25 -0.32    NA  0.63  0.25  0.63
[24,]    NA  0.25  0.71  0.63  0.25 -0.32    NA  0.63  0.25  0.63
[25,]    NA  0.63  0.45  1.00 -0.32 -0.20    NA  1.00  0.63  1.00
[26,]    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA
[27,]    NA  1.00  0.71  0.63  0.25  0.63    NA  0.63  1.00  0.63
[28,]    NA  0.71  1.00  0.45  0.71  0.45    NA  0.45  0.71  0.45
[29,]    NA  0.63  0.45  1.00 -0.32 -0.20    NA  1.00  0.63  1.00
[30,]    NA  0.25  0.71 -0.32  1.00  0.63    NA -0.32  0.25 -0.32
[31,]    NA  0.63  0.45 -0.20  0.63  1.00    NA -0.20  0.63 -0.20
[32,]    NA    NA    NA    NA    NA    NA    NA    NA    NA    NA
[33,]    NA  0.63  0.45  1.00 -0.32 -0.20    NA  1.00  0.63  1.00
[34,]    NA  1.00  0.71  0.63  0.25  0.63    NA  0.63  1.00  0.63
[35,]    NA  0.63  0.45  1.00 -0.32 -0.20    NA  1.00  0.63  1.00

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