mutualInfo: Mutual Information

Description Usage Arguments Details Value Author(s) References See Also Examples

Description

Calculate mutual information via binning

Usage

1
2

Arguments

x

an n by p matrix or ExpressionSet; if x is an ExpressionSet, then the function uses its 'exprs' slot.

...

arguments passed to mutualInfo and MIdist:

  • nbinnumber of bins to calculate discrete probabilities; default is 10.

  • diagif TRUE, then the diagonal of the distance matrix will be displayed; default is FALSE.

  • upperif TRUE, then the upper triangle of the distance matrix will be displayed; default is FALSE.

  • samplefor ExpressionSet methods, if TRUE, then distances are computed between samples, otherwise, between genes.

Details

For mutualInfo each row of x is divided into nbin groups and then the mutual information is computed, treating the data as if they were discrete.

For MIdist we use the transformation proposed by Joe (1989), delta* = (1 - exp(-2 delta))^.5 where delta is the mutual information. The MIdist is then 1-delta*. Joe argues that this measure is then similar to Kendall's tau, tau.dist.

Value

An object of class dist which contains the pairwise distances.

Author(s)

Robert Gentleman

References

H. Joe, Relative Entropy Measures of Multivariate Dependence, JASA, 1989, 157-164.

See Also

dist, KLdist.matrix, cor.dist, KLD.matrix

Examples

1
2
 x <- matrix(rnorm(100), nrow = 5)
 mutualInfo(x, nbin = 3)

Example output

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, 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

Welcome to Bioconductor

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

Loading required package: KernSmooth
KernSmooth 2.23 loaded
Copyright M. P. Wand 1997-2009
            1           2           3           4
2 0.220725846                                    
3 0.104410144 0.066556708                        
4 0.073336461 0.009320618 0.301075824            
5 0.139190106 0.038944099 0.077037020 0.049424411

bioDist documentation built on Nov. 8, 2020, 5:14 p.m.