Skip to contents

Re-build the LD reference panel with given SNPs.

Usage

sHDL.rebuild.ref(
  LD.path,
  gwas.snps,
  LD.path.new,
  nthreads = 1,
  lam.cut = NULL,
  pattern = ".*chr(\\d{1,2})\\.(\\d{1,2})[_\\.].*"
)

Arguments

LD.path

Path to the .rda file where the Eigen decomposition of LD matrix is stored.

gwas.snps

A vector contains SNP IDs from the GWAS summary statistics.

LD.path.new

The new path to the re-built LD reference.

nthreads

Number of threads to use, default nthreads = 1.

lam.cut

Eigenvalue cutoff for LD matrices, default lam.cut = NULL, which means the same cutoff in the original LD reference panel.

pattern

Chromosome and picece pattern of LD files, default is ".*chr(\d{1,2})\.(\d{1,2})[_\.].*".

Value

The new path to the re-built LD reference.