dist_perm: Burst distribution permutations

In meaRtools: Micro-Electro Array (MEA) Analysis

Description

Perform two statistical tests to quantify difference between two burst probability distributions using burst probability distribution data. Performed test are the Maximum Distance between cumulative distributions and Earth Movers Distance between the original probability distributions.

Usage

dist_perm(datafile,np,type,kotype)

Arguments

datafile	A _distributions.csv input file. Format as the output of calc_burst_distributions
np	Number of permutations to perform
type	Name of first genotype
kotype	Name of second genotype

Value

A list containing results of two statistical tests for the input probability distributions data.

data.EMD	Original value of EMD distance
data.EMD	Original value of maximum distance
perm.EMD	A permuted p.value of the EMD distance
perm_p	A permuted p.value of the maximum distance
outp	Maximum distances between genotypes for all permutations performed
out_emd	Maximum Earth Movers Distance between genotypes for all permutations performed
data.wt	Cumulative probabilities of the first genotype
data.ko	Cumulative probabilities of the second genotype
data.wt.Original	Probabilities distribution of the first genotype
data.ko.Original	Probabilities distribution of the second genotype

References

See https://redmine.igm.cumc.columbia.edu/projects/mea/wiki for further details

Examples

    result <- dist_perm(distributionFilePath,10000,"WT","KO")

meaRtools documentation built on May 1, 2019, 7:32 p.m.