Correct the association

Correct the association among SNP in SNP-set.

Usage

correct_association(results = data.frame(), GRCh = "38", output_path)

Arguments

results

data.frame A data frame containing the results of the association.

GRCh

character The genome reference version, either "19"/"37" or "38".

output_path

character The path where to save the results.

Value

data.frame A data frame containing associations p-values corrected using the Benjamini-Hochberg and Benjamini-Yekutieli methods.

Examples

tmp <- tempdir()
results <- data.frame(
  CHR = c(1, 1),
  BP = c(69091, 367658),
  p_SNPSET = c(0.01, 0.2),
  SCORE_MS = c(3.2, 1.1)
)
corrected <- correct_association(results = results, GRCh = "38", output_path = tmp)