Benchmark function for a new type of non-functional dependence

Share:

Description

This function is used to benchmark test for independence on a new type of non-function dependence called patchwork copula

Usage

1
benchmark.patchwork.copula(fun, args, cvals, n = 320, nsim = 500, bins = 20)

Arguments

fun

function or character naming a function. A function should have two vectors of coordinates as first two arguments

args

list of additional arguments to the functions fun. If a function does not need any arguments use an empty list.

cvals

target mutual information values vector of concentration factors, these represent mutual information values (see generate.patchwork.copula)

n

numeric, size of the data sets to generate (default 320 points)

nsim

numeric, how many replicate simulations to run under the null model and H1, default 500

bins

decimal, number of bins of the bins*bins grid, (see generate.patchwork.copula)

Value

This function returns a list data structure that can be further processed with the functions of this package, calculate.power,generate.roc

Author(s)

Sebastian D├╝mcke duemcke@mpipz.mpg.de

See Also

calculate.power,generate.roc

Examples

1
2
3
4
mycor = function(...) cor(...)^2
copula.vals = benchmark.patchwork.copula(mycor,list(),c(.3,1,10))
drop(calculate.power(copula.vals,.95))
roc.plot(generate.roc(copula.vals))