Compute multiple statistics

Computation of the different statistics defined in the package. See Smida et al (2022) for more details.

Usage

comp_stat(MatX, MatY, stat = c("mo", "med"))

Arguments

MatX: numeric matrix of dimension n_point x n containing n trajectories (in columns) of size n_point (in rows).
MatY: numeric matrix of dimension n_point x m containing m trajectories (in columns) of size n_point (in rows).
stat: character string or vector of character string, name of the statistics for which the p-values will be computed, among "mo", "med", "wmw", "hkr", "cff".

Value

list of named numeric value corresponding to the statistic values listed in stat input.

Details

For HKR statistics, only the values of the two statistics, namely HKR1 and HKR2 and not the eigen values (see stat_hkr() for more details).

References

Zaineb Smida, Lionel Cucala, Ali Gannoun & Ghislain Durif (2022) A median test for functional data, Journal of Nonparametric Statistics, 34:2, 520-553, doi:10.1080/10485252.2022.2064997 , hal-03658578

Author

Zaineb Smida, Ghislain DURIF, Lionel Cucala

Examples

simu_data <- simul_data(
    n_point = 100, n_obs1 = 50, n_obs2 = 75, c_val = 10, 
    delta_shape = "constant", distrib = "normal"
)

MatX <- simu_data$mat_sample1
MatY <- simu_data$mat_sample2

res <- comp_stat(MatX, MatY, stat = c("mo", "med", "wmw", "hkr", "cff"))
res
#> $mo
#> [1] 0.9989724
#> 
#> $med
#> [1] 0.9993059
#> 
#> $wmw
#> [1] 0.9987091
#> 
#> $hkr
#>          [,1]
#> T1 18720331.9
#> T2   297319.7
#> 
#> $cff
#> [1] 505329.9
#>