BTestSpearman: Bootstrap Hypothesis Test on Spearman Correlation...
In roahd: Robust Analysis of High Dimensional Data

Description Usage Arguments Details Value See Also Examples

View source: R/bootstrap_spearman_inference.R

This function performs a bootstrap test that checks whether the Spearman correlation structures (e.g. matrices) of two populations of compatible multivariate functional data are equal or not.

BTestSpearman(
  mfD1,
  mfD2,
  bootstrap_iterations = 1000,
  ordering = "MEI",
  normtype = "f",
  verbose = FALSE
)

`mfD1`	is the first functional dataset, specified in form of `mfData` object; it must be compatible with `mfD2`.
`mfD2`	is the second functional dataset, specified in form of `mfData` object; it must be compatible with `mfD1`.
`bootstrap_iterations`	is the number of bootstrap iterations to be performed.
`ordering`	is the kind of ordering to be used in the computation of Spearman's correlation coefficient (default is `MEI`).
`normtype`	is the norm to be used when comparing the Spearman correlation matrices of the two functional datasets (default is Frobenius, allowed values are the same as for parameter `type` in the base function `norm`).
`verbose`	a boolean flag specifying whether to print the progress of bootstrap iterations or not (default is FALSE).

Given a first multivariate functional population, X_1^(i), …, X_n^(i) with i=1, …, L, defined on the grid I = t_1, …, t_P, and a second multivariate functional population, Y_1^(i), …, Y_m^(i) with i=1, …, L, defined on the same grid I, the function performs a bootstrap test to check the hypothesis:

H_0: R_X = R_Y

H_1: R_X \neq R_Y,

where R_X and R_Y denote the L x L matrices of Spearman correlation coefficients of the two populations.

The two functional samples must have the same number of components and must be defined over the same discrete interval t_1, …, t_P.
The test is performed through a bootstrap argument, so a number of bootstrap iterations must be specified as well. A high value for this parameter may result in slow performances of the test (you may consider setting verbose to TRUE to get hints on the process).

The function returns the estimates of the test's p-value and statistics.

BCIntervalSpearman, BCIntervalSpearmanMultivariate, mfData

set.seed(1)

N  <- 200
P <- 100
L <- 2

grid <- seq(0, 1, length.out = P)

# Creating an exponential covariance function to simulate Gaussian data
Cov <- exp_cov_function(grid, alpha = 0.3, beta = 0.4)

# Simulating two populations of bivariate functional data
#
# The first population has very high correlation between first and second component
centerline_1 <- matrix(
  data = rep(sin(2 * pi * grid)),
  nrow = L,
  ncol = P,
  byrow = TRUE
)
values1 <- generate_gauss_mfdata(
  N = N,
  L = L,
  correlations = 0.9,
  centerline = centerline_1,
  listCov = list(Cov, Cov)
)
mfD1 <- mfData(grid, values1)

# Pointwise estimate
cor_spearman(mfD1)

# The second population has zero correlation between first and second component
centerline_2 <- matrix(
  data = rep(cos(2 * pi * grid)),
  nrow = L,
  ncol = P,
  byrow = TRUE
)
values2 <- generate_gauss_mfdata(
  N = N,
  L = L,
  correlations = 0,
  centerline = centerline_2,
  listCov = list(Cov, Cov)
)
mfD2 <- mfData(grid, values2)

# Pointwise estimate
cor_spearman(mfD2)

# Applying the test

BTestSpearman(mfD1, mfD2)