Skip to contents

Compute hertiability from standardized or non-standardized effects

Usage

compute_h2(
  b_joint,
  geno_scale = c("allele", "sd"),
  pheno_sd = 1,
  R_LD = NULL,
  af = NULL,
  full_mat = FALSE
)

Arguments

R_LD

LD pattern (optional). See ?sim_mv for more details.

af

Allele frequencies (optional, allowed only if R_LD is missing). See ?sim_mv for more details.

full_mat

If TRUE, return the full genetic variance-covariance matrix

b_joint_std, b_joint

matrix of standardized or non-standardized effects. Provide only one of these options.