topD3M: Subsetting Top Significant Distribution Pairs and Output...

Description Usage Arguments Details Value Author(s) Examples

Description

Subsetting Top Significant Distribution Pairs and Output Graphics.

Usage

1
2
topD3M(obj, id = seq_along(obj[[1]]), ntop = 10, qval = TRUE,
  plot.it = FALSE)

Arguments

obj

D3M object resulted from d3m.

id

identifier of each row. Default is sequential integer number starting from 1.

ntop

the number of top distribution pairs to be investigated. Default is 10.

qval

if TRUE, ordering the result by false discovery rate.

plot.it

whether output the grphics based on Q-Q plot by JPEG format. Defarult is plot.it = FALSE.

Details

this function extract a subset of the most significant distribution pairs based on p-values. The graphical representation is based on Q-Q plots, which represents the shape difference being tested by d3m function.

Value

top.cases top significant distributions in case group.

top.control top significant distributions in control group.

pval p-values related to top significant distribution pairs.

Q-Q plots of distribution pairs.

Author(s)

Yusuke Matsui & Teppei Shimamura

Examples

1
2
3
4
5
6
library(D3M)
nrep <- 12
cases <- Map(rbeta,rep(30,nrep),rep(1,nrep),rep(5,nrep)); cases <- do.call("rbind",cases)
control <- Map(rbeta,rep(30,nrep),rep(1,nrep),rep(5,nrep)); control <- do.call("rbind",control)
obj <- d3m(cases,control,paranum = 101, q = 2, bsn = 1000)
topD3M(obj,ntop = 10,qval = TRUE, plot.it = FALSE)

ymatts/D3M documentation built on May 4, 2019, 5:30 p.m.