plot.ptest: Plot the results of a fusion test

Description Usage Arguments Details Value Author(s) References See Also Examples

View source: R/rags2ridgesFused.R

Description

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

Usage

1
2
3
4
5
## S3 method for class 'ptest'
plot(x, add.extra = TRUE, ...)

## S3 method for class '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.

Details

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

Value

Invisibly returns x with extra additions.

Author(s)

Anders Ellern Bilgrau, Carel F.W. Peeters <[email protected]>, Wessel N. van Wieringen

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

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
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)

CFWP/rags2ridges documentation built on Sept. 23, 2017, 6:38 a.m.