plot.fmanovatrp {fdANOVA}

Plot P-values of Tests Based on Random Projections for FMANOVA Problem

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

## 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

See Also

fmanova.trp, plot.fanovatests

Examples

# 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)