{"id":482,"date":"2018-03-04T18:01:19","date_gmt":"2018-03-04T23:01:19","guid":{"rendered":"https:\/\/umaine.edu\/research-compliance\/?page_id=482"},"modified":"2025-06-04T15:28:38","modified_gmt":"2025-06-04T19:28:38","slug":"how-to-do-a-power-analysis","status":"publish","type":"page","link":"https:\/\/umaine.edu\/research-compliance\/animal-care\/how-to-do-a-power-analysis\/","title":{"rendered":"How to do a Power Analysis"},"content":{"rendered":"\n<p>To determine the number of replicates needed to detect a \u2018true\u2019 difference between sample means use the following formula (Sokal and Rohlf, 1981 Biometry: The Principles and Practice of Statistics in Biological Research. W.H. Freeman &amp; Co., New York. p 262, Box 9.13):<\/p>\n\n\n\n<p>N = 2(s\/d)<sup>2<\/sup> * {t<sub>\u03b1<\/sub> [v] + t<sub>2(1-P)<\/sub>[v]}<sup>2<\/sup><\/p>\n\n\n\n<p>where<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>N=number of replications<\/li>\n\n\n\n<li>s= true Standard deviation<\/li>\n\n\n\n<li>d=smallest true difference that is desired to be detected (note: one only need know the ratio of s\/d, not their actual values)<\/li>\n\n\n\n<li>v=degrees of freedom (df) with \u2018a\u2019 groups and \u2018n\u2019 replicates per group<\/li>\n\n\n\n<li>\u03b1=significance level<\/li>\n\n\n\n<li>P=desired probability that a significant difference will be found (power of the test)<\/li>\n<\/ul>\n\n\n\n<p>To determine how many animals you will need (N<sub>stable<\/sub>), you need some idea of what the variability is for the parameters you will be measuring. &nbsp;N<sub>stable<\/sub> is obtained through an interative method.&nbsp; To calculate N<sub>stable<\/sub>, you need the df (which is a function of N).&nbsp; Use an estimate of N (N<sub>intitial<\/sub>) to obtain an estimate of df (df<sub>initial<\/sub> ) which is then used to calculate &nbsp;N<sub>2<\/sub>. &nbsp;This new estimate of N (N<sub>2<\/sub>) is then used to calculate a new df, (df<sub>2<\/sub>), which is used, in turn, for the calculation of a new estimate of N (N<sub>3<\/sub>). This method is repeated until a \u2018stable\u2019 N (N<sub>stable<\/sub>) is found.<\/p>\n\n\n\n<h3 class=\"wp-block-heading\">Example<\/h3>\n\n\n\n<p>In the following example, we want to measure how chemicals affect CYP1A enzyme activity in fish. We look up the data we have (or the data of others) on the variability of this enzyme in fish:<\/p>\n\n\n\n<p>Control fish: 269 \u00b1 49 pmol product\/min\/mg protein, n = 3 individual fish<\/p>\n\n\n\n<p>Treated fish: 1,453 \u00b1 139 pmol product\/min\/mg protein, n =3 individual fish<\/p>\n\n\n\n<p>If unknown, the term s\/d can be replaced by CV\/D were CV is the Coefficient of Variation (in %) and D is d in %.&nbsp; The coefficient of variation for these activities varied between ~ 10 \u2013 18% (that is, 49\/269 * 100 = 18.2%; 139\/1453*100 = 9.6%); in other studies we conducted, the CV for this enzyme varied between 34 \u2013 55%. We would like to detect at least a 50% difference between means. Using an average coefficient variation of 30%, an \u03b1 level = 0.05, a desired probability of P = 0.8 and 32 degrees of freedom<sup>1<\/sup>:<\/p>\n\n\n\n<p>N<sub>2<\/sub> = 2 (30\/50)<sup>2<\/sup> * {t<sub>0.05<\/sub>[32] + t<sub>2(1-0.8)<\/sub>[32]}<sup>2<\/sup><\/p>\n\n\n\n<p>= 0.72 * {2.037 + 0.853}<sup>2<\/sup><\/p>\n\n\n\n<p>= 6<\/p>\n\n\n\n<p>Note: Since the enzyme activity could go either up or down with treatment, look up the \u2018t\u2019 values in a Two-tailed Student\u2019s t Table. The value of \u2018t\u2019 for an \u03b1 of 0.05 and 32 df (t<sub>0.05<\/sub>[32]) = 2.037; the value of \u2018 t\u2032 for a P of 0.8 and 32 df (t<sub>2(1-0.8)<\/sub>[32]) = (t<sub>0.4<\/sub>[32]) = 0.853.<\/p>\n\n\n\n<p>The second \u2018round\u2019 of calculations (using N<sub>2 <\/sub>= 6) leads to a df<sub>2<\/sub> of 40 and a new N<sub>3<\/sub> of 6, it is then stable and N<sub>stable <\/sub>= 6. So, 6 fish per group is the number required to detect a significant difference of at least 50 % between the treatments (at a 0.05 level \u03b1), with a likelihood of detecting this difference 80% of the time, if this difference truly exists (this is the power of the test, P).<\/p>\n\n\n\n<p class=\"has-small-font-size\"><sup>1<\/sup> Each fish experiment consists of 8 groups (2 sites x 2 sexes x 2 treatments = 8 groups) with 5 replicates per group (N<sub>intitial<\/sub> for the calculation of the df<sub>initial<\/sub>), yielding v = 8 * (5 \u2013 1) = 32 df.<\/p>\n\n\n\n<h3 class=\"wp-block-heading\">G*Power<\/h3>\n\n\n\n<p><a href=\"https:\/\/www.psychologie.hhu.de\/en\/arbeitsgruppen\/allgemeine-psychologie-und-arbeitspsychologie\/gpower\" data-type=\"link\" data-id=\"https:\/\/www.psychologie.hhu.de\/en\/arbeitsgruppen\/allgemeine-psychologie-und-arbeitspsychologie\/gpower\" target=\"_blank\" rel=\"noreferrer noopener\">G*Power (external link)<\/a>\u00a0is a free software to estimate the sample size and to conduct statistical power analysis.<\/p>\n\n\n\n<p><\/p>\n","protected":false},"excerpt":{"rendered":"<p>To determine the number of replicates needed to detect a \u2018true\u2019 difference between sample means use the following formula (Sokal and Rohlf, 1981 Biometry: The Principles and Practice of Statistics in Biological Research. W.H. Freeman &amp; Co., New York. p 262, Box 9.13): N = 2(s\/d)2 * {t\u03b1 [v] + t2(1-P)[v]}2 where To determine how [&hellip;]<\/p>\n","protected":false},"author":430,"featured_media":0,"parent":6660,"menu_order":0,"comment_status":"closed","ping_status":"closed","template":"templates\/page-withsidebar.php","meta":{"_seopress_robots_primary_cat":"","_seopress_titles_title":"","_seopress_titles_desc":"","_seopress_robots_index":"","_kad_blocks_custom_css":"","_kad_blocks_head_custom_js":"","_kad_blocks_body_custom_js":"","_kad_blocks_footer_custom_js":"","_monsterinsights_skip_tracking":false,"_monsterinsights_sitenote_active":false,"_monsterinsights_sitenote_note":"","_monsterinsights_sitenote_category":0,"footnotes":""},"class_list":["post-482","page","type-page","status-publish","hentry"],"taxonomy_info":[],"featured_image_src_large":false,"author_info":{"display_name":"aashe","author_link":"https:\/\/umaine.edu\/research-compliance\/author\/aashe\/"},"comment_info":"","_links":{"self":[{"href":"https:\/\/umaine.edu\/research-compliance\/wp-json\/wp\/v2\/pages\/482","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/umaine.edu\/research-compliance\/wp-json\/wp\/v2\/pages"}],"about":[{"href":"https:\/\/umaine.edu\/research-compliance\/wp-json\/wp\/v2\/types\/page"}],"author":[{"embeddable":true,"href":"https:\/\/umaine.edu\/research-compliance\/wp-json\/wp\/v2\/users\/430"}],"replies":[{"embeddable":true,"href":"https:\/\/umaine.edu\/research-compliance\/wp-json\/wp\/v2\/comments?post=482"}],"version-history":[{"count":15,"href":"https:\/\/umaine.edu\/research-compliance\/wp-json\/wp\/v2\/pages\/482\/revisions"}],"predecessor-version":[{"id":7899,"href":"https:\/\/umaine.edu\/research-compliance\/wp-json\/wp\/v2\/pages\/482\/revisions\/7899"}],"up":[{"embeddable":true,"href":"https:\/\/umaine.edu\/research-compliance\/wp-json\/wp\/v2\/pages\/6660"}],"wp:attachment":[{"href":"https:\/\/umaine.edu\/research-compliance\/wp-json\/wp\/v2\/media?parent=482"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}