fdANOVA (version 0.1.2)

plot.fanovatests: Plot P-values of Tests Based on Random Projections for FANOVA Problem

Description

The p-values of the tests based on random projections for the one-way analysis of variance for (univariate) functional data against the number of projections are plotted.

Usage

# S3 method for fanovatests
plot(x, y, …)

Arguments

x

an "fanovatests" object. More precisely, a result of the function fanova.tests for the standard tests based on random projections.

y

an "fanovatests" object. More precisely, a result of the function fanova.tests for the permutation tests based on random projections.

...

additional arguments not used.

See Also

fanova.tests, plot.fmanovatrp

Examples

# NOT RUN {
# Some of the examples may run some time.

# gait data (the first feature)
library(fda)
gait.data.frame <- as.data.frame(gait)
x.gait <- as.matrix(gait.data.frame[, 1:39])

# vector of group labels
group.label.gait <- rep(1:3, each = 13)
# }
# NOT RUN {
set.seed(123)
fanova4 <- fanova.tests(x.gait, group.label.gait, test = "TRP",
                        parallel = TRUE, nslaves = 2,
                        params = list(paramTRP = list(k = c(10, 20, 30), B.TRP = 1000)))
set.seed(123)
fanova5 <- fanova.tests(x.gait, group.label.gait, test = "TRP",
                        parallel = TRUE, nslaves = 2,
                        params = list(paramTRP = list(k = c(10, 20, 30),
                                                      permutation = TRUE, B.TRP = 1000)))
plot(x = fanova4)
plot(y = fanova5)
plot(x = fanova4, y = fanova5)

set.seed(123)
fanova6 <- fanova.tests(x.gait, group.label.gait, test = "TRP",
                        parallel = TRUE, nslaves = 2,
                        params = list(paramTRP = list(k = c(10, 20, 30), projection = "BM",
                                                      B.TRP = 1000)))
set.seed(123)
fanova7 <- fanova.tests(x.gait, group.label.gait, test = "TRP",
                        parallel = TRUE, nslaves = 2,
                        params = list(paramTRP = list(k = c(10, 20, 30), projection = "BM",
                                                      permutation = TRUE, B.TRP = 1000)))
plot(x = fanova6)
plot(y = fanova7)
plot(x = fanova6, y = fanova7)
# }