Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

African LD reference #305

Closed
rameez500 opened this issue Jan 29, 2022 · 5 comments
Closed

African LD reference #305

rameez500 opened this issue Jan 29, 2022 · 5 comments

Comments

@rameez500
Copy link

I have a quick question; I am working in African ancestry for both target dataset and summary statistics.
I noticed that in the ldpred2 paper, the European LD reference is available at https://doi.org/10.6084/m9.figshare.13034123.
Is there any African LD reference for LDPRED2 that I can find?

@privefl
Copy link
Owner

privefl commented Jan 29, 2022

No, we have not provided any LD ref for African ancestry yet.
If you have more than 2000 individuals in your target data, you can always use this as reference panel.

Otherwise, if you have access to the UK Biobank and your data is from West/South Africa, you can use the "Nigeria" group of individuals that you can reproduce from https://github.com/privefl/UKBB-PGS#code-to-reproduce-ancestry-groups (from this paper).

@privefl privefl pinned this issue Jan 29, 2022
@rameez500
Copy link
Author

Hi Florian,

Thanks for the prompt reply.

I am working on African Americans ancestry samples (target dataset). This target dataset is imputed by CAAPA from the Michigan imputation server.
The total number of samples is about 4000 with 10 million variants.

I know Hapmap3 variants are strongly recommended to be filtered in LDpred2 because they offer good coverage of the genome.

I was reading one of the issue of using Using Alternate LD Reference Panels: #225

Which one of the following would you recommend me to follow:

  1. The code used to provide the LD ref in the paper: https://github.com/privefl/paper-ldpred2/blob/master/code/provide-ld-ref.R.
    OR
  2. European LD reference which contains about 1.1M variants

@privefl
Copy link
Owner

privefl commented Jan 30, 2022

If you have African sumstats, you need to use an African LD ref.
So option 1.

@rameez500
Copy link
Author

Thanks for your help.
I wanted to ask something regarding option 1.

I have two datasets African summary statistics and African target dataset; but I'm seeking to find African LD ref.

In provide-ld-ref.R line 74 , you can compute LD matrices using snp_cor() and save the result in "ld-ref/LD_chr"

If I want to calculate PRS using the LDpred2-auto, I can use example-with-provided-ldref.R with provided LD reference.

My question is:
Are the variants provided in the LD reference in example-with-provided-ldref.R line 7, generated
from provide-ld-ref.R line 74?

@privefl
Copy link
Owner

privefl commented Jan 31, 2022

Yes, this is the code I used to prepare the LD ref.
And then made this other script as an example/tutorial how to use it.

@privefl privefl closed this as completed Mar 1, 2022
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants