# plot1D: Plot regions of the representative tree in 1D In jacsor/MRS: Multi-Resolution Scanning for Cross-Sample Differences

 plot1D R Documentation

## Plot regions of the representative tree in 1D

### Description

This function visualizes the regions of the representative tree of the output of the `mrs` function. For each region the posterior probability of difference (PMAP) or the effect size is plotted.

### Usage

```plot1D(ans, type = "prob", group = 1, dim = 1, regions = rep(1,
length(ans\$RepresentativeTree\$Levels)), legend = FALSE, main = "default",
abs = TRUE)
```

### Arguments

 `ans` An `mrs` object. `type` What is represented at each node. The options are `type = c("eff", "prob")`. Default is `type = "prob"`. `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. `abs` If `TRUE`, plot the absolute value of the effect size. Only used when `type = "eff"`.

### References

Soriano J. and Ma L. (2016). Probabilistic multi-resolution scanning for two-sample differences. Journal of the Royal Statistical Society: Series B (Statistical Methodology). http://onlinelibrary.wiley.com/doi/10.1111/rssb.12180/abstract

Ma L. and Soriano J. (2016). Analysis of distributional variation through multi-scale Beta-Binomial modeling. arXiv. http://arxiv.org/abs/1604.01443

### Examples

```set.seed(1)
p = 1
n1 = 200
n2 = 200
mu1 = matrix( c(0,10), nrow = 2, byrow = TRUE)
mu2 = mu1; mu2 = mu1 + .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 = mrs(X, G, K=10)
plot1D(ans, type = "prob")
plot1D(ans, type = "eff")
```

jacsor/MRS documentation built on Oct. 12, 2022, 8:33 p.m.