DepLogo (version 1.0)

getDeps: Compute dependencies between positions

Description

Computes the dependencies (as measures by mutual information) between all positions (columns) of discrete data. Specifically, it returns for each pair of positions (i,j) the mutual information I(X_i,X_j) multiplied by the number N of sequences (rows), which may also be used for testing the statistical significance of mutual information values, as for large N, 2*N*I(X_i,X_j) is approximately chi squared.

Usage

getDeps(data, ...)

Arguments

data

the data for computing mutual information. Either a DLData object or a data.frame. In the latter case, the symbols of the alphabet must be provided as a second parameter.

...

the symbols of the alphabet as character vector, only if data is a data.frame

Value

A matrix of the mutual information values, where the diagonal is fixed to zero.

Examples

# NOT RUN {
data <- DLData(c("ACGT", "ATTA"))
deps <- getDeps(data)


# }

Run the code above in your browser using DataCamp Workspace