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})[_\.].*"
.