This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How do I know the SNPS which are removed when clumping are in a LD block or SNPs not existed in 1000 Genome

Hi everyone, I'm not sure how to know why my SNPs are being removed applying clumping when reading the exposure data (in MR analysis). I would like to know as if the SNPs are being deleted because they're not in 1000 Genome I can use proxies. And this way more IVs will be left for MR analysis. Does anybody take any approaches or we need to do it manually? Any help would be appreciated Paria

randomization data clumping mendelian gwas
• 2,325 views • link
updated by Kevin Blighe Kevin Blighe • written by paria • paria

Please let us know which program you are using.

• 0 views • link
written by Kevin Blighe Kevin Blighe

I use TwoSampleMR. And the code is as below:

exp_data <- read_exposure_data("ALS_sig.txt", sep="\t", snp_col= "SNP", beta_col= "b", eaf_col= "effectAlleleFreq", se_col= "StdErr", effect_allele_col= "A1", other_allele_col= "A2", pval_col= "p", samplesize_col= "SampleSize", ncase_col= "N_cases", ncontrol_col= "N_controls", clump=TRUE)
• 0 views • link
updated by Kevin Blighe Kevin Blighe • written by paria • paria

I see, I am not familiar with this package, unfortunately. Can you not derive the list of excluded SNPs via an indirect process of deduction / elimination? That is, save the list before pruning into a vector, then save the pruned list into another vector, and then determine the differences between both of these.

• 0 views • link
written by Kevin Blighe Kevin Blighe

0 answers

No answers yet.

Log in to answer this question.

More posts like this