rags2ridges (version 2.2.2)

plot.ptest: Plot the results of a fusion test

Description

Plot a histogram of the null distribution and the observed test statistic in a permutation type "fusion test".

Usage

# S3 method for ptest
plot(x, add.extra = TRUE, ...)
# S3 method for ptest
hist(x, add.extra = TRUE, ...)

Arguments

x

A ptest object (a list). Usually the output of fused.test.

add.extra

A logical. Add extra information to the plot.

…

Arguments passed to plot.

Value

Invisibly returns x with extra additions.

Details

plot.ptest is simply a wrapper for hist.ptest.

References

Bilgrau, A.E., Peeters, C.F.W., Eriksen, P.S., Boegsted, M., and van Wieringen, W.N. (2015). Targeted Fused Ridge Estimation of Inverse Covariance Matrices from Multiple High-Dimensional Data Classes, arXiv:1509.07982v1 [stat.ME].

See Also

fused.test, print.ptest

Examples

# NOT RUN {
ns <- c(10, 5, 23)
Ylist <- createS(ns, p = 15, topology = "banded", dataset = TRUE)

# Use the identity target matrix for each class
Tlist <- replicate(length(ns), diag(15), simplify = FALSE)

# Do the test
lam <- matrix(10, 3, 3)
diag(lam) <- 1
ft <- fused.test(Ylist, Tlist, lambda = lam, n.permutations = 500)

# The returned object can alo be plotted via
hist(ft)
# or via the alias
plot(ft)
# }