def test_surrogate_analysis_fdr():
rng = np.random.RandomState(0)
data = np.load("../examples/data/eeg_32chans_10secs.npy")
ts1 = data[0, 0:512].ravel()
ts2 = data[1, 0:512].ravel()
p_val, surr_vals, surrogates, r_value = surrogate_analysis(
ts1,
ts2,
num_surr=1000,
estimator_func=None,
ts1_no_surr=False,
rng=rng)
num_ts = 2
num_pvals = num_ts * (num_ts - 1) / 2.0
num_pvals = np.int32(num_pvals)
p_vals = np.ones([num_pvals, 1]) * p_val
h, crit_p = fdr(p_vals, 0.01, 'pdep')
expected_result = np.load("data/test_ts_surrogates_fdr_h.npy")
assert h == expected_result
expected_result = np.load("data/test_ts_surrogates_fdr_crit_p.npy")
assert crit_p == expected_result
def test_fdr():
pvals = np.ones((1, 2 * (2 - 1) / 2)) * 0.0025
q = 0.01
method = "pdep"
h, crit_p = fdr(pvals, q, method)
h = h.ravel()[0]
crit_p = crit_p.ravel()[0]
assert h == True
assert crit_p == 0.0025
def test_fdr():
num_pvals = 2 * (2 - 1) / 2.0
num_pvals = np.int32(num_pvals)
pvals = np.ones((1, num_pvals)) * 0.0025
q = 0.01
method = "pdep"
h, crit_p = fdr(pvals, q, method)
h = h.ravel()[0]
crit_p = crit_p.ravel()[0]
assert h == True
assert crit_p == 0.0025
def test_surrogate_analysis2_fdr():
rng = np.random.RandomState(0)
data = np.load("../examples/data/eeg_32chans_10secs.npy")
ts1 = data[0, 0:512].ravel()
ts2 = data[1, 0:512].ravel()
p_val, surr_vals, surrogates, r_value = surrogate_analysis(
ts1, ts2, num_surr=1000, ts1_no_surr=True, rng=rng)
num_ts = 2
p_vals = np.ones([num_ts * (num_ts - 1) / 2, 1]) * p_val
h, crit_p = fdr(p_vals, 0.01, 'pdep')
expected_result = np.load("data/test_ts_surrogates2_fdr_h.npy")
assert h == expected_result
# -*- coding: utf-8 -*-
import numpy as np
np.set_printoptions(precision=3, linewidth=256)
from dyfunconn.ts import fdr, surrogate_analysis
if __name__ == "__main__":
rng = np.random.RandomState(0)
data = np.load(
"/home/makism/Github/dyfunconn/examples/data/eeg_32chans_10secs.npy")
ts1 = data[0, :].ravel()
ts2 = data[1, :].ravel()
p_val, corr_surr, surrogates, r_value = surrogate_analysis(
ts1, ts2, num_surr=1000, ts1_no_surr=True, rng=rng)
num_ts = 2
p_vals = np.ones([num_ts * (num_ts - 1) / 2, 1]) * p_val
q = 0.01
method = 'pdep'
h, crit_p = fdr(p_vals, q, method)
print("p-value: {0}, h: {1} (critical p-value: {2})".format(
p_val, h, crit_p))