Calculate Z-score for motif frequencies

Calculate Z-score for motif frequencies

Usage

calculate_Z(observed = NULL, nulls = NULL)

Arguments

observed

A list of observed motif frequencies for each motif type. List elements must be named 'lambda', 'bifan', 'V', 'PPI_V', and 'delta' (not necessarily in that order).

nulls

A list of null distributions for each motif type as returned by generate_nulls.

Value

A numeric vector with the Z-score for each motif type.

Examples

# Simulating it for test purposes
null <- rnorm(1000, mean = 5, sd = 1)
nulls <- list(
    lambda = null, V = null, PPI_V = null, delta = null, bifan = null
)
observed <- list(lambda = 7, bifan = 13, delta = 9, V = 5, PPI_V = 10)
z <- calculate_Z(observed, nulls)
# Check for motif enrichment (Z > 5)
z[which(z > 5)]
#>    bifan 
#> 8.003716