-
Notifications
You must be signed in to change notification settings - Fork 9
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
different p-values #37
Comments
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Hi hope you are doing well,
I have tested Raremetal with two scenario's, I have a multi-sample called and qced vcf from 2 datasets. these datasets were merged and multisample called and qced together. I ran gene burden test on this final multisample called vcf using EPACTS and Raremetal. I got similar results from both tools.
then I split the two datasets, and ran firstly Raremetalworker on both datasets separately and then did aggregation test using summary statistics of the two datasets using Raremetal. and with the second approach I lost all my significant hits that I got from EPACTS and Raremetal.
I checked top hits of single.variants.score files from Combined.singlevar.score.txt and compared it with the single variant score files of separate datasets, everything is same AC, callrate, MAFs, but somehow huge difference in the pvalues. in the final burden test. I am just curious why is it so. just for your reference please see below.
result from combined.singlevar.score.txt
result from study1.singlevar.score.txt
result from study2.singlevar.score.txt
result from study1.study2.meta.score.txt
As you can see for this particular variant allele count, allele freq, N_ref, N_Het, N_Alt are all same as combined.singlevar.score.txt if we add up the columns of study1 and study2 but there is huge difference in the pvlaues. just wondering what could be the reason for that? I have noticed similar behavior with all the top variants from combined analysis and ended up losing all the significant hits. sorry for my long post, I am just curious what is going on here.
looking forward to hearing from you and thanks in advance.
Kind Regards,
Haider
The text was updated successfully, but these errors were encountered: