seqqvalhh: Sequential generalized Hochberg and Hommel procedures based...

View source: R/seqqvalhh.r

seqqvalhhR Documentation

Sequential generalized Hochberg and Hommel procedures based on q-values

Description

Sequential generalized Hochberg and Hommel procedures based on q-values

Usage

seqqvalhh(
  pm = matrix(rep(c(0.03, 0.04, 0.01), times = 2), ncol = 3, nrow = 2),
  alpham = matrix(rep(c(0.02, 0.03, 0.05), each = 2), ncol = 3, nrow = 2),
  epsilon = 1e-10,
  precision = 10,
  method = "Hochberg"
)

Arguments

pm

Matrix of group-sequential p-values for different hypotheses (in row) at different times (in column).

alpham

Matrix of alpha spending corresponding to the p-values pm. For each row, alpha levels must be non-decreasing.

epsilon

Numeric scalar indicating the lower bound for alpha.

precision

Integer scalar for precision of the values, obsolete for backward compatibility.

method

Character scalar "Hochberg" or "Hommel".

Value

List with elements

  • rejected: the index set of rejected hypotheses

  • decisionsm: rejection decision for each endpoint (row) at each timepoint (column)

  • cumdecisionsm: cumulative rejection decision for each endpoint (row) at each timepoint (column);

  • alphaused: alpha levels actually used for each endpoint (row) at each timepoint (column).

Author(s)

Xiaodong Luo

Examples

pm <- matrix(rep(c(0.03, 0.04, 0.01), times = 2), ncol = 3, nrow = 2)
alpham <- matrix(rep(c(0.02, 0.03, 0.05), each = 2), ncol = 3, nrow = 2)
seqqvalhh(pm = pm, alpham = alpham, method = "Hochberg")
seqqvalhh(pm = pm, alpham = alpham, method = "Hommel")

MultiGroupSequential documentation built on Aug. 9, 2023, 5:10 p.m.