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

corrected_effect_se is unusually large #5

Open
desert0eagle opened this issue Dec 20, 2022 · 2 comments
Open

corrected_effect_se is unusually large #5

desert0eagle opened this issue Dec 20, 2022 · 2 comments

Comments

@desert0eagle
Copy link

desert0eagle commented Dec 20, 2022

Hi,

I`m new in this field and not good at the mathematics.

Recently I`m doing two sample MR analysis between metabolome sumstats by using MRlap. I found most of my results have extremely large corrected_effect_se value (corrected_effect_p is nearly 1) while the observed_effect_se seems usual (observed_effect_p is very low and significant) and p_difference is nearly 1 (as below). May I ask how to interpret the results and which value (corrected or observed) should I use in this case.
For the GWAS sumstats, both of them are analyzed by about 10000 samples, and sample overlap rate is 100%, number of IVs is 14 to 22. I also attached my R code below. This is a common issue and I also found the similar case by using other sumstats (n=100,000, sample overlap rate 10%).
スクリーンショット 2022-12-20 17 07 32
スクリーンショット 2022-12-20 17 23 40

@desert0eagle
Copy link
Author

@n-mounier HI, could you please response the above issue, Im afraid its a common issue, if that this model may not be useful anymore.

@n-mounier
Copy link
Owner

Hello,
Sorry for not being able to reply before.
It would have been helpful to see the log, not just the function call and the results. As described in the paper, the heritability and the polygenicity are very important parameters to estimate the corrected effect, and the variance of the corrected causal effects is estimated through a sampling strategy (the number of simulations needed to estimate the variance would help detect abnormal situations), which could fail if the heritability is too low or if the trait is not polygenic enough. Could this be the case with metabolomics data?

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