Power estimation and sample size determination for replication studies of genome-wide association studies

BMC Genomics

Methodology

Weichuan Yu¹  Wei Jiang¹ 

[1] Department of Electronic and Computer Engineering, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China;

关键词: Replication study;    Power;    Empirical Bayes;   

DOI  :  10.1186/s12864-015-2296-4

来源: Springer

【 摘 要 】

BackgroundReplication study is a commonly used verification method to filter out false positives in genome-wide association studies (GWAS). If an association can be confirmed in a replication study, it will have a high confidence to be true positive. To design a replication study, traditional approaches calculate power by treating replication study as another independent primary study. These approaches do not use the information given by primary study. Besides, they need to specify a minimum detectable effect size, which may be subjective. One may think to replace the minimum effect size with the observed effect sizes in the power calculation. However, this approach will make the designed replication study underpowered since we are only interested in the positive associations from the primary study and the problem of the “winner’s curse” will occur.ResultsAn Empirical Bayes (EB) based method is proposed to estimate the power of replication study for each association. The corresponding credible interval is estimated in the proposed approach. Simulation experiments show that our method is better than other plug-in based estimators in terms of overcoming the winner’s curse and providing higher estimation accuracy. The coverage probability of given credible interval is well-calibrated in the simulation experiments. Weighted average method is used to estimate the average power of all underlying true associations. This is used to determine the sample size of replication study. Sample sizes are estimated on 6 diseases from Wellcome Trust Case Control Consortium (WTCCC) using our method. They are higher than sample sizes estimated by plugging observed effect sizes in power calculation.ConclusionsOur new method can objectively determine replication study’s sample size by using information extracted from primary study. Also the winner’s curse is alleviated. Thus, it is a better choice when designing replication studies of GWAS. The R-package is available at: http://bioinformatics.ust.hk/RPower.html.

【 授权许可】

CC BY   
© Jiang and Yu. 2015

【 图 表 】

12951_2017_270_Article_IEq10.gif

Fig. 8

Fig. 3

Fig. 3

Fig. 7

12951_2015_111_Article_IEq5.gif

Fig. 1

12951_2015_155_Article_IEq30.gif

Fig. 5

Fig. 6

12951_2017_270_Article_IEq9.gif

Fig. 2

12951_2017_270_Article_IEq8.gif

Fig. 1

12951_2015_111_Article_IEq4.gif

Fig. 1

13100_2023_302_Article_IEq7.gif

Fig. 2

Fig. 4

Fig. 5

12951_2015_155_Article_IEq28.gif

42004_2023_1020_Article_IEq18.gif

Fig. 3

Fig. 7

13100_2023_302_Article_IEq2.gif

13100_2023_302_Article_IEq1.gif

12951_2017_270_Article_IEq7.gif

42004_2023_1020_Article_IEq12.gif

Fig. 2

Fig. 7

Fig. 1

Fig. 3

12951_2015_111_Article_IEq2.gif

Fig. 1

Fig. 5

Fig. 6

12951_2015_155_Article_IEq27.gif

12951_2017_315_Article_IEq2.gif

12951_2015_155_Article_IEq26.gif

12951_2015_155_Article_IEq25.gif

12951_2015_155_Article_IEq24.gif

Fig. 2

12951_2017_270_Article_IEq5.gif

Fig. 3

Fig. 1

12951_2015_111_Article_IEq1.gif

Fig. 6

Fig. 4

Fig. 1

12951_2017_315_Article_IEq1.gif

Fig. 1

12951_2015_155_Article_IEq23.gif

12951_2015_155_Article_IEq22.gif

12951_2015_155_Article_IEq21.gif

12951_2015_155_Article_IEq20.gif

12951_2015_155_Article_IEq19.gif

12951_2015_155_Article_IEq18.gif

12951_2015_155_Article_IEq16.gif

Fig. 3

12951_2015_155_Article_IEq14.gif

12951_2015_155_Article_IEq13.gif

12951_2015_155_Article_IEq12.gif

Fig. 2

12951_2015_155_Article_IEq9.gif

Fig. 2

12864_2017_3527_Article_IEq9.gif

Fig. 7

Fig. 2

Fig. 3

Fig. 5

Fig. 5

Fig. 10

Fig. 1

Fig. 9

12951_2015_155_Article_IEq7.gif

Fig. 4

Fig. 8

Fig. 4

Fig. 3

Fig. 2

Fig. 4

Fig. 1

Fig. 3

Fig. 1

40517_2023_273_Article_IEq9.gif

Fig. 2

40517_2023_273_Article_IEq6.gif

40517_2023_273_Article_IEq4.gif

Fig. 1

40517_2023_273_Article_IEq2.gif

Fig. 2

42004_2023_1025_Article_IEq7.gif

Fig. 9

Fig. 9

Fig. 4

Fig. 1

Fig. 4

Fig. 3

Fig. 4

41512_2023_158_Article_IEq26.gif

Fig. 2

Fig. 2

12951_2017_270_Article_IEq3.gif

Fig. 1

41512_2023_158_Article_IEq20.gif

Fig. 1

Fig. 8

Fig. 6

Fig. 2

Fig. 2

Fig. 6

Fig. 6

12951_2015_155_Article_IEq6.gif

Fig. 8

12951_2017_323_Article_IEq1.gif

Fig. 5

Fig. 3

41512_2023_158_Article_IEq2.gif

Fig. 7

41512_2023_158_Article_IEq1.gif

Fig. 1

Fig. 1

Fig. 7

Fig. 4

Fig. 3

Fig. 3

Fig. 6

12951_2015_155_Article_IEq4.gif

12867_2016_60_Article_IEq1.gif

Fig. 4

Fig. 2

Fig. 9