Estimate partitioned gene-level heritability

This function estimate partitioned gene-level heritability

Usage

h2gene(susie_fit, ld, annot, n_sample = 500)

Arguments

susie_fit

susie fit object which is obtained with susie or susie_suff_stat

ld

A n_snp by n_snp linkage disequilibrium matrix which should match the one used to construct susie_fit

annot

A n_snp by n_annot binary matrix denoting the

n_sample

Number of posterior samples to simulate for estimation. Default to 500.

Value

A list with the following elements:

hsq

n_sample posterior samples of heritability

ncausal

n_sample posterior samples of number of causals