sum_stats: Simulate summary statistics

Description Usage Arguments

View source: R/sum_stats.R

Description

Simulate summary statistics

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
sum_stats(
  snps,
  evd_list,
  n_copies = 30,
  n1,
  n2,
  h1,
  h2,
  neffect1,
  neffect2,
  gamma,
  eta,
  q,
  tau,
  omega,
  cores,
  ld_prune_pval_thresh = 0.001,
  r2_thresh = 0.1
)

Arguments

snps

SNP info data frame

evd_list

List of eigen value decompositions for LD regions

n_copies

How many times to replicate chromosome 19

n1

Sample size in study 1

n2

sample size in study 2

h1

Expected heritability of trait 1

h2

Expected heritability of trait 2

neffect1

Average number of trait 1 effect snps

neffect2

Average number of trait 2 effect snps. Note that this is in addition to SNPs that effect Y through M.

gamma

gamma

eta

eta

q

q

tau

tau = h1*(gamma^2)/h2 Proportion of Y heritability explained by M through causal mechanism

omega

omega = h1*(eta)^2/h2 omega*q^2 is proportion of Y heritability explained by M through confounder mechanism


jean997/causeSims documentation built on Sept. 4, 2020, 4:29 p.m.