Generate beta hats from standardized or non-standardized direct SNP effects and LD

Usage

gen_bhat_from_b(
  b_joint,
  N,
  trait_corr = NULL,
  R_LD = NULL,
  R_LD_eigen = NULL,
  af = NULL,
  est_s = FALSE,
  input_geno_scale = c("allele", "sd"),
  input_pheno_sd = 1,
  output_geno_scale = c("allele", "sd"),
  output_pheno_sd = 1
)

Arguments

b_joint: Matrix of joint effect sizes
N: Sample size, scalar, vector, matrix, or data.frame. See ?sim_mv for more details.
trait_corr: Matrix of population trait correlation (traits by traits)
R_LD: LD pattern (optional). See ?sim_mv for more details.
R_LD_eigen: Optionally, also supply list of eigen decompositions. The program will not check that this matches R_LD.
af: Allele frequencies (optional, allowed only if R_LD is missing). See ?sim_mv for more details.
est_s: Estimate standard errors?
input_geno_scale: Genotype scale of effects in b_joint
input_pheno_sd: Phenotype sd for effects in b_joint
output_geno_scale: Output genotype scale
output_pheno_sd: Output phenotype sd

Details

This has been made an internal function. To resample summary statistics use resample_sumstats.