sim_pants: Simulate 'pants' to estimate size & power

Description Usage Arguments See Also

View source: R/sim_pants.R

Description

Simulate pants to estimate size (effect=0) & power (effect>0) with score_fcn=identity vs. H1: "greater". Does not support type=="mediation".

Usage

1
2
3
sim_pants(Gmat, phenotype, type = c("contrasts", "correlation"),
  effect.v = c(0, 0.2), alpha = 0.05, nsim = 10^3, nperm = 10^3,
  ncores = 1, seed = 1, verbose = TRUE, ker = NULL)

Arguments

Gmat

Binary feature (e.g. gene) by pathway inclusion matrix, indicating which features are in which pathways.

phenotype

Vector of sample characteristics (correlation: numeric; contrasts: character). Should be same length as ncol(object).

type

Type of ezlimma analysis per feauture; must be one of"contrasts" (limma_contrasts), "correlation" (limma_cor), or "mediation" (hitman). You can specify just the initial letter.

effect.v

Numeric vector of log fold-changes or percent of phentotypes to add.

alpha

Alpha level.

nsim

Number of simulations.

nperm

Number of sample permutations to evaluate significance of pathways.

ncores

Integer. If > 1, number of cores to use for parallel computing. You can detect how many are available for your system using detectCores.

seed

Random seed.

verbose

Logical; should the number of simulations be printed every 100 simulations?

ker

Laplacian kernel matrix representing the interaction network.

See Also

sim_pants_mediation


jdreyf/PANTS documentation built on July 18, 2019, 10:12 a.m.