aSPC_dcor_get_pvals: An Adaptive Sum of Powered Correlation Test (aSPC) with dcor

In jasonzyx/aSPC: An Adaptive Sum of Powered Correlation Test (aSPC) to Test Global Association Between Two Groups of Variables

Description

An Adaptive Sum of Powered Correlation Test (aSPC) with dcor

Usage

aSPC_dcor_get_pvals(df1, df2, pow = pow, T0s = T0s)

Arguments

df1,	first matrix
df2,	second matrix
pow,	power integer candidates, default c(1:8, Inf)
T0s,	A B by length(pow) T statistcs matrix based on B permutations

Value

Pvalues of SPC and aSPC tests

References

Xu Z., Pan W. An adaptive and powerful test for two groups of variables with high dimension

Examples

library(mvtnorm)
sigma = diag(0.6, 10) + 0.4
n = 260 # sample size
p = 3000; q = 300;

Z = rmvnorm(n=n, mean=rep(0,10), sigma=sigma)
X = rmvnorm(n=n, mean=rep(0,p), sigma=diag(1, p))
Y = rmvnorm(n=n, mean=rep(0,q), sigma=diag(1, q))
X = as.data.frame(cbind(Z[,1:5], X))
Y = as.data.frame(cbind(Z[,6:10], Y))
dim(X)
dim(Y)
set.seed(123) # to ensure we can replicate the permutation P-value
T_nulls = aSPC_dcor_get_Tstat(X, Y, pow = c(1:8, Inf), B = 100, show_b = T)
fit = aSPC_dcor_get_pvals(X, Y, pow = c(1:8, Inf), T0s = T_nulls)
fit

jasonzyx/aSPC documentation built on May 18, 2019, 5:55 p.m.