plot.fmanovatrp: Plot P-values of Tests Based on Random Projections for... In fdANOVA: Analysis of Variance for Univariate and Multivariate Functional Data

Description

The p-values of the tests based on random projections for multivariate analysis of variance for functional data against the number of projections are plotted.

Usage

 ```1 2``` ```## S3 method for class 'fmanovatrp' plot(x, y, withoutRoy = FALSE, ...) ```

Arguments

 `x` an "`fmanovatrp`" object. More precisely, a result of the function `fmanova.trp` for the standard tests based on random projections. `y` an "`fmanovatrp`" object. More precisely, a result of the function `fmanova.trp` for the permutation tests based on random projections. `withoutRoy` a logical indicating whether to plot the p-values of the Rp test. `...` additional arguments not used.

Note

We can use only one of the arguments `x` and `y`, or both simultaneously.

Author(s)

Tomasz Gorecki, Lukasz Smaga

`fmanova.trp`, `plot.fanovatests`
 ``` 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33``` ```# Some of the examples may run some time. # gait data (both features) library(fda) gait.data.frame <- as.data.frame(gait) x.gait <- vector("list", 2) x.gait[[1]] <- as.matrix(gait.data.frame[, 1:39]) x.gait[[2]] <- as.matrix(gait.data.frame[, 40:78]) # vector of group labels group.label.gait <- rep(1:3, each = 13) # the tests based on random projections with the Gaussian white noise generated for projections set.seed(123) fmanova1 <- fmanova.trp(x.gait, group.label.gait, k = c(1, 5, 10, 15, 20)) fmanova2 <- fmanova.trp(x.gait, group.label.gait, k = c(1, 5, 10, 15, 20), permutation = TRUE, B = 1000, parallel = TRUE, nslaves = 2) plot(x = fmanova1) plot(x = fmanova1, withoutRoy = TRUE) plot(y = fmanova2) plot(x = fmanova1, y = fmanova2) plot(x = fmanova1, y = fmanova2, withoutRoy = TRUE) # the tests based on random projections with the Brownian motion generated for projections set.seed(123) fmanova3 <- fmanova.trp(x.gait, group.label.gait, k = c(1, 5, 10, 15, 20), projection = "BM") fmanova4 <- fmanova.trp(x.gait, group.label.gait, k = c(1, 5, 10, 15, 20), projection = "BM", permutation = TRUE, B = 1000, parallel = TRUE, nslaves = 2) plot(x = fmanova3) plot(x = fmanova3, withoutRoy = TRUE) plot(y = fmanova4) plot(x = fmanova3, y = fmanova4) plot(x = fmanova3, y = fmanova4, withoutRoy = TRUE) ```