{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,3,3]],"date-time":"2026-03-03T02:41:13Z","timestamp":1772505673434,"version":"3.50.1"},"reference-count":57,"publisher":"Oxford University Press (OUP)","issue":"4","license":[{"start":{"date-parts":[[2023,3,21]],"date-time":"2023-03-21T00:00:00Z","timestamp":1679356800000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"DOI":"10.13039\/501100000038","name":"Natural Sciences and Engineering Research Council of Canada","doi-asserted-by":"publisher","id":[{"id":"10.13039\/501100000038","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2023,4,3]]},"abstract":"<jats:title>Abstract<\/jats:title><jats:sec><jats:title>Motivation<\/jats:title><jats:p>Accurate power and sample size estimation is crucial to the design and analysis of genetic association studies. When analyzing a binary trait via logistic regression, important covariates such as age and sex are typically included in the model. However, their effects are rarely properly considered in power or sample size computation during study planning. Unlike when analyzing a continuous trait, the power of association testing between a binary trait and a genetic variant depends, explicitly, on covariate effects, even under the assumption of gene\u2013environment independence. Earlier work recognizes this hidden factor but the implemented methods are not flexible. We thus propose and implement a generalized method for estimating power and sample size for (discovery or replication) association studies of binary traits that (i) accommodates different types of nongenetic covariates E, (ii) deals with different types of G\u2013E relationships, and (iii) is computationally efficient.<\/jats:p><\/jats:sec><jats:sec><jats:title>Results<\/jats:title><jats:p>Extensive simulation studies show that the proposed method is accurate and computationally efficient for both prospective and retrospective sampling designs with various covariate structures. A proof-of-principle application focused on the understudied African sample in the UK Biobank data. Results show that, in contrast to studying the continuous blood pressure trait, when analyzing the binary hypertension trait ignoring covariate effects of age and sex leads to overestimated power and underestimated replication sample size.<\/jats:p><\/jats:sec><jats:sec><jats:title>Availability and implementation<\/jats:title><jats:p>The simulated datasets can be found on the online web-page of this manuscript, and the UK Biobank application data can be accessed at https:\/\/www.ukbiobank.ac.uk. The R package SPCompute that implements the proposed method is available at CRAN. The genome-wide association studies are carried out using the software PLINK 2.0 [Purcell et al. (Plink: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 2007;81:559\u201375.)].<\/jats:p><\/jats:sec>","DOI":"10.1093\/bioinformatics\/btad139","type":"journal-article","created":{"date-parts":[[2023,3,21]],"date-time":"2023-03-21T23:16:26Z","timestamp":1679440586000},"source":"Crossref","is-referenced-by-count":3,"title":["The hidden factor: accounting for covariate effects in power and sample size computation for a binary trait"],"prefix":"10.1093","volume":"39","author":[{"ORCID":"https:\/\/orcid.org\/0000-0001-7082-2323","authenticated-orcid":false,"given":"Ziang","family":"Zhang","sequence":"first","affiliation":[{"name":"Department of Statistical Science, Faculty of Arts and Science, University of Toronto , Toronto, ON M5G 1Z5, Canada"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-5640-937X","authenticated-orcid":false,"given":"Lei","family":"Sun","sequence":"additional","affiliation":[{"name":"Department of Statistical Science, Faculty of Arts and Science, University of Toronto , Toronto, ON M5G 1Z5, Canada"},{"name":"Division of Biostatistics, Dalla Lana School of Public Health, University of Toronto , Toronto, ON M5T 3M7, Canada"}],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"286","published-online":{"date-parts":[[2023,3,21]]},"reference":[{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"203","DOI":"10.1038\/s41586-018-0579-z","article-title":"The UK Biobank resource with deep phenotyping and genomic data","volume":"562","author":"Bycroft","year":"2018","journal-title":"Nature"},{"key":"2023040323221677200_","author":"Castelloe","year":"2018"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"3385","DOI":"10.1002\/sim.2771","article-title":"Sample size determination for logistic regression revisited","volume":"26","author":"Demidenko","year":"2007","journal-title":"Stat Med"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"36","DOI":"10.1002\/sim.2980","article-title":"Sample size and optimal design for logistic regression with binary interaction","volume":"27","author":"Demidenko","year":"2008","journal-title":"Stat Med"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"302","DOI":"10.1214\/13-STS456","article-title":"Pooled association tests for rare genetic variants: a review and some new results","volume":"29","author":"Derkach","year":"2014","journal-title":"Stat Sci"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"227","DOI":"10.1002\/gepi.20297","article-title":"Estimation of significance thresholds for Genomewide Association scans","volume":"32","author":"Dudbridge","year":"2008","journal-title":"Genet Epidemiol"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"431","DOI":"10.1093\/biomet\/71.3.431","article-title":"Biased estimates of treatment effect in randomized experiments with nonlinear regressions and omitted covariates","volume":"71","author":"Gail","year":"1984","journal-title":"Biometrika"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"478","DOI":"10.1093\/aje\/155.5.478","article-title":"Sample size requirements for association studies of gene-gene interaction","volume":"155","author":"Gauderman","year":"2002","journal-title":"Am J Epidemiol"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"35","DOI":"10.1002\/sim.973","article-title":"Sample size requirements for matched case-control studies of gene\u2013environment interaction","volume":"21","author":"Gauderman","year":"2002","journal-title":"Stat Med"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"E5272","DOI":"10.1073\/pnas.1419064111","article-title":"Measuring missing heritability: inferring the contribution of common variants","volume":"111","author":"Golan","year":"2014","journal-title":"Proc Natl Acad Sci USA"},{"key":"2023040323221677200_","first-page":"100","article-title":"Algorithm as 136: a k-means clustering algorithm","volume":"28","author":"Hartigan","year":"1979","journal-title":"J R Stat Soc Ser C (Appl Stat)"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"e1000008","DOI":"10.1371\/journal.pgen.1000008","article-title":"Data and theory point to mainly additive genetic variance for complex traits","volume":"4","author":"Hill","year":"2008","journal-title":"PLoS Genet"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"786","DOI":"10.1016\/j.ajhg.2021.02.014","article-title":"Estimation of non-additive genetic variance in human complex traits from a large sample of unrelated individuals","volume":"108","author":"Hivert","year":"2021","journal-title":"Am J Hum Genet"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"117","DOI":"10.5808\/GI.2012.10.2.117","article-title":"Sample size and statistical power calculation in genetic association studies","volume":"10","author":"Hong","year":"2012","journal-title":"Genomics Inform"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"795","DOI":"10.1002\/sim.4780080704","article-title":"Sample size tables for logistic regression","volume":"8","author":"Hsieh","year":"1989","journal-title":"Stat Med"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"1623","DOI":"10.1002\/(SICI)1097-0258(19980730)17:14<1623::AID-SIM871>3.0.CO;2-S","article-title":"A simple method of sample size calculation for linear and logistic regression","volume":"17","author":"Hsieh","year":"1998","journal-title":"Stat Med"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"332","DOI":"10.2174\/1389202918666170307161450","article-title":"Assessing the heritability of complex traits in humans: methodological challenges and opportunities","volume":"18","author":"Mayhew","year":"2017","journal-title":"Curr Genomics"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1017\/S0954579412000855","article-title":"Gene\u2013environment correlation in developmental psychopathology","volume":"25","author":"Knafo","year":"2013","journal-title":"Dev Psychopathol"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1186\/1746-4811-9-29","article-title":"The advantages and limitations of trait analysis with GWAS: a review","volume":"9","author":"Korte","year":"2013","journal-title":"Plant Methods"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"2528","DOI":"10.3758\/s13428-021-01546-0","article-title":"Estimating power in (generalized) linear mixed models: an open introduction and tutorial in R","volume":"53","author":"Kumle","year":"2021","journal-title":"Behav Res Methods"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"969","DOI":"10.1038\/s41588-020-0676-4","article-title":"Dynamic incorporation of multiple in silico functional annotations empowers rare variant association analysis of large whole-genome sequencing studies at scale","volume":"52","author":"Li","year":"2020","journal-title":"Nat Genet"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"802","DOI":"10.1016\/j.ajhg.2019.03.002","article-title":"Dynamic scan procedure for detecting rare-variant association regions in whole-genome sequencing studies","volume":"104","author":"Li","year":"2019","journal-title":"Am J Hum Genet"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"1599","DOI":"10.1038\/s41592-022-01640-x","article-title":"A framework for detecting noncoding rare-variant associations of large-scale whole-genome sequencing studies","volume":"19","author":"Li","year":"2022","journal-title":"Nat Methods"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1186\/s13073-020-00818-2","article-title":"Statistical power in covid-19 case-control host genomic study design","volume":"12","author":"Lin","year":"2020","journal-title":"Genome Med"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"1632","DOI":"10.1002\/sim.2617","article-title":"A practical approach to computing power for generalized linear models with nominal, count, or ordinal responses","volume":"26","author":"Lyles","year":"2007","journal-title":"Stat Med"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"e1608","DOI":"10.1002\/mpr.1608","article-title":"A tutorial on conducting genome-wide association studies: quality control and statistical analysis","volume":"27","author":"Marees","year":"2018","journal-title":"Int J Methods Psychiatr Res"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"249","DOI":"10.1375\/twin.11.3.249","article-title":"A century of Hardy\u2013Weinberg equilibrium","volume":"11","author":"Mayo","year":"2008","journal-title":"Twin Res Hum Genet"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","DOI":"10.1201\/9780203753736","volume-title":"Generalized Linear Models","author":"McCullagh","year":"2019"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1186\/s12916-019-1364-z","article-title":"GWAS and enrichment analyses of non-alcoholic fatty liver disease identify new trait-associated genes and pathways across emerge network","volume":"17","author":"Namjou","year":"2019","journal-title":"BMC Med"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"97","DOI":"10.1002\/sim.3728","article-title":"A modified approach to estimating sample size for simple logistic regression with one continuous covariate","volume":"29","author":"Novikov","year":"2010","journal-title":"Stat Med"},{"key":"2023040323221677200_","first-page":"778","author":"O\u2019Brien","year":"1986"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"539","DOI":"10.1177\/1745691616646366","article-title":"What should researchers expect when they replicate studies? a statistical view of replicability in psychological science","volume":"11","author":"Patil","year":"2016","journal-title":"Perspect Psychol Sci"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"848","DOI":"10.1038\/ng.2346","article-title":"Including known covariates can reduce power to detect genetic effects in case-control studies","volume":"44","author":"Pirinen","year":"2012","journal-title":"Nat Genet"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"309","DOI":"10.1037\/0033-2909.84.2.309","article-title":"Genotype-environment interaction and correlation in the analysis of human behavior","volume":"84","author":"Plomin","year":"1977","journal-title":"Psychol Bull"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"559","DOI":"10.1086\/519795","article-title":"Plink: a tool set for whole-genome association and population-based linkage analyses","volume":"81","author":"Purcell","year":"2007","journal-title":"Am J Hum Genet"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","DOI":"10.1002\/9780470316436","volume-title":"Linear Statistical Inference and Its Applications","author":"Rao","year":"1973"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"169","DOI":"10.2147\/TACG.S122250","article-title":"Analysis of binary responses with outcome-specific misclassification probability in genome-wide association studies","volume":"9","author":"Rekaya","year":"2016","journal-title":"Appl Clin Genet"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"227","DOI":"10.2307\/1403444","article-title":"Some surprising results about covariate adjustment in logistic regression models","volume":"59","author":"Robinson","year":"1991","journal-title":"Int Stat Rev"},{"key":"2023040323221677200_","first-page":"424","article-title":"How people make their own environments: a theory of genotype \u2192 environment effects","volume":"54","author":"Scarr","year":"1983","journal-title":"Child Dev"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"79","DOI":"10.2307\/2531897","article-title":"Power\/sample size calculations for generalized linear models","volume":"44","author":"Self","year":"1988","journal-title":"Biometrics"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"31","DOI":"10.2307\/2532736","article-title":"Power calculations for likelihood ratio tests in generalized linear models","volume":"48","author":"Self","year":"1992","journal-title":"Biometrics"},{"key":"2023040323221677200_","volume-title":"Approximation Theorems of Mathematical Statistics","author":"Serfling","year":"2009"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"1192","DOI":"10.1111\/j.0006-341X.2000.01192.x","article-title":"On power and sample size calculations for likelihood ratio tests in generalized linear models","volume":"56","author":"Shieh","year":"2000","journal-title":"Biometrics"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"4696","DOI":"10.1002\/sim.5879","article-title":"Ignoring the matching variables in cohort studies\u2013when is it valid and why?","volume":"32","author":"Sj\u00f6lander","year":"2013","journal-title":"Stat Med"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"e1001779","DOI":"10.1371\/journal.pmed.1001779","article-title":"UK Biobank: an open access resource for identifying the causes of a wide range of complex diseases of middle and old age","volume":"12","author":"Sudlow","year":"2015","journal-title":"PLoS Med"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"352","DOI":"10.1002\/gepi.20068","article-title":"Reduction of selection bias in genomewide studies by resampling","volume":"28","author":"Sun","year":"2005","journal-title":"Genet Epidemiol"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"545","DOI":"10.1007\/s00439-011-0948-2","article-title":"Br-squared: a practical solution to the winner\u2019s curse in genome-wide scans","volume":"129","author":"Sun","year":"2011","journal-title":"Hum Genet"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"229","DOI":"10.1038\/s41588-017-0009-4","article-title":"Multi-trait analysis of genome-wide association summary statistics using MTAG","volume":"50","author":"Turley","year":"2018","journal-title":"Nat Genet"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"e1010231","DOI":"10.1371\/journal.pgen.1010231","article-title":"Major sex differences in allele frequencies for x chromosomal variants in both the 1000 genomes project and gnomAD","volume":"18","author":"Wang","year":"2022","journal-title":"PLoS Genet"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"89","DOI":"10.1016\/j.ajhg.2018.06.002","article-title":"Estimating SNP-based heritability and genetic correlation in case-control studies directly and with summary statistics","volume":"103","author":"Weissbrod","year":"2018","journal-title":"Am J Hum Genet"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"27","DOI":"10.1080\/01621459.1981.10477597","article-title":"Sample size for logistic regression with small response probability","volume":"76","author":"Whittemore","year":"1981","journal-title":"J Am Stat Assoc"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"565","DOI":"10.1038\/ng.608","article-title":"Common SNPs explain a large proportion of the heritability for human height","volume":"42","author":"Yang","year":"2010","journal-title":"Nat Genet"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"1304","DOI":"10.1038\/ng.3941","article-title":"Concepts, estimation and interpretation of SNP-based heritability","volume":"49","author":"Yang","year":"2017","journal-title":"Nat Genet"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"3700","DOI":"10.1002\/sim.8688","article-title":"Genetic association studies with bivariate mixed responses subject to measurement error and misclassification","volume":"39","author":"Zhang","year":"2020","journal-title":"Stat Med"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"78","DOI":"10.1002\/gepi.20437","article-title":"Correcting \u201cwinner\u2019s curse\u201d in odds ratios from genomewide association findings for major complex human diseases","volume":"34","author":"Zhong","year":"2010","journal-title":"Genet Epidemiol"},{"key":"2023040323221677200_","first-page":"1","article-title":"Causal associations between risk factors and common diseases inferred from GWAS summary data","volume":"9","author":"Zhu","year":"2018","journal-title":"Nat Commun"},{"key":"2023040323221677200_","doi-asserted-by":"crossref","first-page":"605","DOI":"10.1086\/512821","article-title":"Overcoming the winner\u2019s curse: estimating penetrance parameters from case-control data","volume":"80","author":"Z\u00f6llner","year":"2007","journal-title":"Am J Hum Genet"}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/bioinformatics\/advance-article-pdf\/doi\/10.1093\/bioinformatics\/btad139\/49583758\/btad139.pdf","content-type":"application\/pdf","content-version":"am","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/39\/4\/btad139\/49736179\/btad139.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/39\/4\/btad139\/49736179\/btad139.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2023,12,9]],"date-time":"2023-12-09T09:01:57Z","timestamp":1702112517000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/doi\/10.1093\/bioinformatics\/btad139\/7082519"}},"subtitle":[],"editor":[{"given":"Russell","family":"Schwartz","sequence":"additional","affiliation":[],"role":[{"role":"editor","vocabulary":"crossref"}]}],"short-title":[],"issued":{"date-parts":[[2023,3,21]]},"references-count":57,"journal-issue":{"issue":"4","published-print":{"date-parts":[[2023,4,3]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/btad139","relation":{},"ISSN":["1367-4811"],"issn-type":[{"value":"1367-4811","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2023,4,1]]},"published":{"date-parts":[[2023,3,21]]},"article-number":"btad139"}}