hommelFast: Calculates adjusted p-values of Hommel's method efficiently.

Description Usage Arguments Value Author(s) Examples

View source: R/shortcuts.R

Description

Calculates adjusted p-values of Hommel's method efficiently.

Usage

1
 hommelFast (pvalues, simes = TRUE) 

Arguments

pvalues

A vector of p-values.

simes

If TRUE, a Simes' test is used, if FALSE Hommel's test is used.

Value

Returns a hommel object.

Author(s)

Rosa Meijer: r.j.meijer@lumc.nl

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
  
  #Generate a vector of pvalues
  set.seed(1)
  n <- 1000
  pvalues <- c(runif(0.50*n,0,0.001), runif(0.50*n,0,1))
  
  #Create an hommel object in which the adjusted p-values are stored, based on a Simes'test
  #(or Hommel's test, by choosing simes = FALSE):
  hom <- hommelFast(pvalues, simes = TRUE)
  
  #Retrieve the first 10 adjusted p-values by using the \code{p.adjust} method 
  # from the hommel package. Note that they are not sorted
  hommel::p.adjust(hom)[1:10]

cherry documentation built on May 7, 2021, 5:06 p.m.