get.phat.quadratic: Get the pvalue for the quadratic statistic

In pipeGS: Permutation p-Value Estimation for Gene Set Tests

Description

Get the pvalue for the quadratic statistic

Usage

get.phat.quadratic(x, y, method = c("is", "nested", "mc", "exact"),
  N.mc = 10^3, N.is = 10^3, N.level = 1000, q = 0.2, B = 5 *
  length(x), mc.cores = 1, do.sd = FALSE)

Arguments

x	binary vector of treatment assignment
y	gene expression measurement matrix
method	method to estimate pvalue for the linear statistic
N.mc	number of mc samples is method == "mc"
N.is	number of mc samples is method == "is"
N.level	number of samples in each level if method == "nested"
q	progression quantile if method == "nested"
B	number of burn-ins to use if method == "nested"
mc.cores	number of cores to use in mcapply
do.sd	logical to indicate whether calculating sd or not

Value

a list containing phat from different methods

Examples

x <- c(rep(0, 4), rep(1, 4))
y <- matrix(c(rnorm(4*5, 0, 1), rnorm(4*5, 2 ,1)), nrow = 8)
get.phat.quadratic(x, y, method = "is")
get.phat.quadratic(x, y, method = "nested")
get.phat.quadratic(x, y, method = "mc")
get.phat.quadratic(x, y, method = "exact")
get.phat.quadratic(x, y, method = "is", do.sd = TRUE)
get.phat.quadratic(x, y, method = "nested", do.sd = TRUE)
get.phat.quadratic(x, y, method = "mc", do.sd = TRUE)
get.phat.quadratic(x, y, method = "exact", do.sd = TRUE)

pipeGS documentation built on May 2, 2019, 5:58 a.m.