Skip to content

Commit 3b5ebce

Browse files
committed
add prior equation for LDpred2
1 parent d92e9c4 commit 3b5ebce

File tree

5 files changed

+40
-11
lines changed

5 files changed

+40
-11
lines changed

docs/articles/LDpred2.html

Lines changed: 26 additions & 11 deletions
Some generated files are not rendered by default. Learn more about customizing how changed files appear on GitHub.
11 Bytes
Loading
0 Bytes
Loading
489 Bytes
Loading

vignettes/LDpred2.Rmd

Lines changed: 14 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -193,6 +193,20 @@ In that case, you can estimate the overall $N^\text{eff}$ using e.g. `quantile(8
193193

194194
## Computing LDpred2 scores genome-wide
195195

196+
[LDpred2](https://bit.ly/ldpred2-paper) assumes the following model for effect sizes,
197+
198+
<div class="math">
199+
\[
200+
\beta_j \sim \left\{
201+
\begin{array}{ll}
202+
\mathcal N\left(0, \dfrac{h^2}{M p}\right) & \mbox{with probability } p,\\
203+
0 & \mbox{otherwise,}\end{array}
204+
\right.
205+
\]
206+
</div>
207+
208+
where $p$ is the proportion of causal variants (aka the polygenicity), $M$ the number of variants and $h^2$ the SNP-heritability (the heritability capturing by the set of variants used).
209+
196210
Note that you should **run LDpred2 genome-wide**.
197211
Just build the SFBM (the sparse LD matrix on disk) so that it contains selected variants for all chromosomes at once (see the for-loop below).
198212

0 commit comments

Comments
 (0)