plot1D: Plot regions of the representative tree in 1D

Description Usage Arguments References Examples

View source: R/plot1D.R

Description

This function visualizes the regions of the representative tree of the output of the ARMT function. For each region the reject probability or the effect size is plotted.

Usage

1
2
plot1D(ans, type = "rej", group = 1, dim = 1, regions = rep(1,
  length(ans$RepresentativeTree$Levels)), legend = FALSE, main = "default")

Arguments

ans

An ARMT object.

type

What is represented at each node. The options are type = c("eff", "rej"). Default is type = "rej".

group

If type = "eff", which group effect size is used. Default is group = 1.

dim

If the data are multivariate, dim is the dimension plotted. Default is dim = 1.

regions

Binary vector indicating the regions to plot. The default is to plot all regions.

legend

Color legend for type. Default is legend = FALSE.

main

Overall title for the plot.

References

Soriano J. and Ma L. (2014). Multi-resolution two-sample comparison through the divide-merge Markov tree. Preprint. http://arxiv.org/abs/1404.3753

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
set.seed(1)
p = 1
n1 = 200
n2 = 200
mu1 = matrix( c(0,10), nrow = 2, byrow = TRUE)
mu2 = mu1; mu2[2] = mu1[2] + .01
sigma = c(1,.1)

Z1 = sample(2, n1, replace=TRUE, prob=c(0.9, 0.1))
Z2 = sample(2, n2, replace=TRUE, prob=c(0.9, 0.1))
X1 = mu1[Z1] + matrix(rnorm(n1*p), ncol=p)*sigma[Z1]
X2 = mu2[Z2] + matrix(rnorm(n2*p), ncol=p)*sigma[Z1]
X = rbind(X1, X2)
G = c(rep(1, n1), rep(2,n2))

ans = ARMT(X, G, K=10)
plot1D(ans, type = "rej")
plot1D(ans, type = "eff")

jacsor/ARMtree2 documentation built on May 17, 2017, 6:49 a.m.