Skip to contents

Compute multiple statistics

comp_stat.Rd

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

See also

stat_mo(), stat_med(), stat_wmw(), stat_hkr(), stat_cff()

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.9989137
#> 
#> $med
#> [1] 0.9992842
#> 
#> $wmw
#> [1] 0.998584
#> 
#> $hkr
#>            [,1]
#> T1 2588074354.4
#> T2     302154.8
#> 
#> $cff
#> [1] 511743.7
#>