The functions detailed inside the piece of code below (in a Gist) has been useful for me when I had to calculate many possible scenarios of statistical power and sample size. The formulae were taken from the article of Samuels et al., AJHG 2006, and the script showed even useful for making a variety of comparative plots. This is intended for estimating power/ sample size in association studies, involving mitochondrial DNA haplogroups (which are categories whose frequencies depend on each other), on a Chi-square test basis.

Continue reading

Author's picture

aurora-mareviv

Anesthesiologist, MD, postdoc. Utter Rstats geek

Universidade de Santiago de Compostela

Spain