Caclulate effect estimates and SE from Z scores
1
2
Entering edit mode
5.3 years ago
muralinmars ▴ 100

Hi,

I am using a post-gwas analysis tools to test genetic overlap between different traits based on summary level data, and it requires the effect size(beta) and SE.

So for summary stats that has the Z score measures, beta and SE are estimated by using the following formula that was published in a large scale study on Alzheimers.

SE~=sqrt(VP/(ES×2pq))
    Beta=SE×Z

where:

  • VP: phenotypic variance approximated to 1;
  • ES: effective sample size
  • pq: allele frequencies (p allele, q allele)
  • Z: Z scores

The summary stats that I currently have access to does not have allele frequency information.

head of the summary stats:

> Chr     Pos     MarkerName      EffectAllele    OtherAllele     Z    P.value TotalSampleSize 
> 1       744055  rs3131000               a           g         -1.60      0.10  23546

Is there a way to calculate the Betas and SE without contacting the study, with just relying on the information that is currently publicly available.

gwas meta-analysis R • 7.8k views
ADD COMMENT
0
Entering edit mode

Hi, Would it be possible to send the link of the Alzheimers paper you found the equation in please?

ADD REPLY
0
Entering edit mode

I am somewhat pressed for time but it's not 100% clear which data-points you have. Please take a look here and see if you can work back from the Z-score calculation: A: SNP dataset and Z Score

Do you have the allele tallies?

ADD REPLY
0
Entering edit mode
2.4 years ago

You can use this equation:

Beta = z / sqrt(2p(1− p)(n + z^2)) and

SE =1 / sqrt(2p(1− p)(n + z^2))

Where p is the frequency of the imputed SNP, you could use out reference panel to calculate p. For reference please go to

https://images.nature.com/full/nature-assets/ng/journal/v48/n5/extref/ng.3538-S1.pdf

ADD COMMENT

Login before adding your answer.

Traffic: 1332 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6