hdl: estimate a genetic covariance matrix using High Definition...

In MichelNivard/GenomicSEM: Structural equation modeling based on GWAS summary statistics

estimate a genetic covariance matrix using High Definition Likelihood (HDL) estimation in R

Description

Function to run HDL (https://github.com/zhenin/HDL) to compute the genetic covariance between a series of traits based on genome wide summary statistics obtained from GWAS. The results generate by this function are necessary and sufficient to facilitate the fit of structural equation models (or other multivariate models) to the genetic covariance matrix. HDL is more powerfull than LDSC but if the LD structure in the reference file mismatches the GWAS LD structure, LDSC seems to perfrom better, expescially for estiamtes of heritability. For medium samples (N > 50.000) with moderate SNP-h2 (snp h2 > 0.07) where the LD structure isnt similar we would recomend ldsc, especially for GWAS. If you have small GWAS ( N < 25.000) the extra power HDL provides is worth the downward bias in snp h2 estimates relative to ldsc.

Usage

hdl(traits,sample.prev,population.prev,LDpath,Nref,trait.names=NULL,method, ...)

Arguments

traits	A vector of strings which point to munged files for trait you want to include in a Genomic SEM model. the HDL function works with standard munged files
sample.prev	A vector of sample prevalences for dichotomous traits and NA for continous traits
population.prev	A vector of population prevalences for dichotomous traits and NA for continous traits
LDpath	String which contains the path to the folder in which the LD matrices used in the analysis are located. Expects LD matirices formated as required by the original HDL software.
Nref	Sample size of the reference file, default is 335265
trait.names	A character vector specifying how the traits should be named in the genetic covariance (S) matrix. These variable names can subsequently be used in later steps for model specification. If no value is provided, the function will automatically name the variables using the generic from of V1-VX.
method	sting, either "piecewise" which estimates the heritability or genetic covariance locally in chunks across the genome and then sums these estimates, or "jackknife" which uses a genoem wide estiamte and uses a jackknife estimator for the variance of the parameter. defaults to "piecewise" the original HDL implementation is equal to "jackknife"

Value

The function returns a list with 3 named entries

S	estimated genetic covariance matrix
V	variance covariance matrix of the parameter estimates in S
I	matrix containing the "cross trait intercepts", or the error covariance between traits induced by overlap, in terms of subjects, between the GWASes on which the analyses are based

References

Ning, Z., Pawitan, Y. & Shen, X. High-definition likelihood inference of genetic correlations across human complex traits. Nat Genet (2020).

MichelNivard/GenomicSEM documentation built on Dec. 24, 2024, 3:23 a.m.