<?xml version="1.0" encoding="utf-8"?><!DOCTYPE article  PUBLIC '-//OASIS//DTD DocBook XML V4.4//EN'  'http://www.docbook.org/xml/4.4/docbookx.dtd'><article><articleinfo><title>FAQ/power/rmPow</title><revhistory><revision><revnumber>56</revnumber><date>2013-03-08 10:18:01</date><authorinitials>localhost</authorinitials><revremark>converted to 1.6 markup</revremark></revision><revision><revnumber>55</revnumber><date>2009-01-28 10:32:51</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>54</revnumber><date>2009-01-28 10:31:49</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>53</revnumber><date>2009-01-28 10:30:44</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>52</revnumber><date>2009-01-28 10:30:06</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>51</revnumber><date>2009-01-28 10:25:06</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>50</revnumber><date>2009-01-28 10:23:11</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>49</revnumber><date>2009-01-28 09:54:19</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>48</revnumber><date>2009-01-27 16:31:57</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>47</revnumber><date>2009-01-27 16:30:48</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>46</revnumber><date>2008-09-25 13:35:49</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>45</revnumber><date>2008-09-25 10:46:07</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>44</revnumber><date>2008-09-25 10:45:47</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>43</revnumber><date>2008-09-25 10:42:49</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>42</revnumber><date>2008-09-25 10:40:37</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>41</revnumber><date>2008-09-25 10:39:51</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>40</revnumber><date>2008-09-25 10:36:59</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>39</revnumber><date>2008-09-25 10:34:20</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>38</revnumber><date>2008-09-25 09:09:29</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>37</revnumber><date>2008-09-25 09:07:55</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>36</revnumber><date>2008-09-24 16:21:48</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>35</revnumber><date>2008-09-24 16:21:31</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>34</revnumber><date>2008-09-24 16:14:33</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>33</revnumber><date>2008-09-24 16:09:10</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>32</revnumber><date>2008-09-24 16:08:48</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>31</revnumber><date>2007-09-05 12:54:26</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>30</revnumber><date>2007-09-05 12:54:06</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>29</revnumber><date>2007-09-05 12:48:17</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>28</revnumber><date>2007-09-05 09:17:36</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>27</revnumber><date>2007-09-04 15:31:55</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>26</revnumber><date>2007-09-04 15:31:18</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>25</revnumber><date>2007-09-04 15:22:21</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>24</revnumber><date>2007-03-20 15:12:50</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>23</revnumber><date>2007-03-06 15:06:33</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>22</revnumber><date>2007-03-06 15:00:59</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>21</revnumber><date>2007-03-06 15:00:09</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>20</revnumber><date>2007-02-19 17:21:30</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>19</revnumber><date>2007-02-19 17:21:04</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>18</revnumber><date>2007-02-19 17:20:51</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>17</revnumber><date>2007-01-15 17:11:35</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>16</revnumber><date>2007-01-15 10:05:05</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>15</revnumber><date>2007-01-15 10:04:39</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>14</revnumber><date>2006-12-08 10:31:40</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>13</revnumber><date>2006-12-08 10:31:24</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>12</revnumber><date>2006-12-06 13:11:58</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>11</revnumber><date>2006-12-06 13:11:33</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>10</revnumber><date>2006-11-30 16:00:42</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>9</revnumber><date>2006-10-17 14:53:55</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>8</revnumber><date>2006-10-17 14:43:04</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>7</revnumber><date>2006-10-17 14:32:15</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>6</revnumber><date>2006-10-05 15:19:04</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>5</revnumber><date>2006-10-05 15:18:40</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>4</revnumber><date>2006-10-05 15:18:30</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>3</revnumber><date>2006-10-05 15:11:52</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>2</revnumber><date>2006-10-05 15:09:58</date><authorinitials>PeterWatson</authorinitials></revision><revision><revnumber>1</revnumber><date>2006-10-03 11:14:59</date><authorinitials>PeterWatson</authorinitials></revision></revhistory></articleinfo><para><emphasis role="underline">Spreadsheet and SPSS macro inputs</emphasis> </para><itemizedlist><listitem><para>alpha is the likelihood of making a type I error (usually = 0.05) </para></listitem><listitem><para>etasq is partial $$\eta^text{2}$$/100 so, for example, 5.9% = 0.059 </para></listitem></itemizedlist><para>Partial $$\eta^text{2}$$ =  $$ \frac{\mbox{SS(effect)}}{\mbox{SS(effect) + SS(its error)}}$$ where SS represents the sum of squares associated with a particular term in the anova. </para><para>In other words, partial $$\eta^text{2}$$ represents the proportion of variance in outcome predicted by the effect after adjusting for other terms in the anova.  Click <ulink url="https://imaging.mrc-cbu.cam.ac.uk/statswiki/FAQ/power/rmPow/statswiki/FAQ/power/rmPow?action=AttachFile&amp;do=get&amp;target=etasqrp.pdf">here</ulink> for further details on partial $$\eta^text{2}$$ and <ulink url="https://imaging.mrc-cbu.cam.ac.uk/statswiki/FAQ/power/rmPow/statswiki/FAQ/power/rmPow?action=AttachFile&amp;do=get&amp;target=etasq.pdf">here.</ulink> </para><para>If the Sums of Squares are not available you can<ulink url="https://imaging.mrc-cbu.cam.ac.uk/statswiki/FAQ/power/rmPow/statswiki/FAQ/power/rsqform#">construct eta-squared</ulink>. </para><itemizedlist><listitem><para>num(erator) = df of term of interest = the product of the (number of levels of each factor -1) in term of interest </para></listitem><listitem><para>sum (B-1) = sum of dfs involving <emphasis role="strong">only</emphasis> between subject factors in the anova or zero otherwise. df = Product of number of levels minus 1 of each between subject factor in term of interest. e.g. For a three way interaction involving three between subject factors, abc, we sum the dfs of the six lower order combinations: ab, ac and bc, a, b and c to that of abc </para></listitem><listitem><para>Prod (W-1) = df of within subject effect if in term of interest or 1 otherwise. df = Product of number of levels minus 1 of each within subject factor in term of interest </para></listitem><listitem><para>Prod W = The product of all the ithin subject levels in term of interest </para></listitem><listitem><para>corr is the average correlation between levels of the repeated measures (=0 if no within subjects factors) </para></listitem><listitem><para>Total sample size </para></listitem></itemizedlist><para><ulink url="https://imaging.mrc-cbu.cam.ac.uk/statswiki/FAQ/power/rmPow/statswiki/FAQ/power/powexampleN#">Example input</ulink> </para><para>Power can be computed using an EXCEL <ulink url="https://imaging.mrc-cbu.cam.ac.uk/statswiki/FAQ/power/rmPow/statswiki/FAQ/power/rmPow?action=AttachFile&amp;do=get&amp;target=aovp.xls">spreadsheet</ulink> or the SPSS syntax below. Power analysis software using Winer (1991, pp 136-138) for balanced anovas may be downloaded from <ulink url="http://www.soton.ac.uk/~cpd/anovas/datasets/">here</ulink> with details of how to compute inputs <ulink url="https://imaging.mrc-cbu.cam.ac.uk/statswiki/FAQ/power/rmPow/statswiki/FAQ/Doncaster#">here.</ulink> </para><para>[ COPY AND PASTE THE SYNTAX IN THE BOX BELOW INTO A SPSS SYNTAX WINDOW AND RUN; ADJUST INPUT DATA AS REQUIRED] </para><screen><![CDATA[DATA LIST free
/alpha num bsum wdf corr ntot rsq. 
BEGIN DATA. 
.05 2 1 2 0.0 60 0.0588
.05 2 1 2 0.3 67 0.0588
END DATA.
set errors=none. 
matrix.
get m /variables=alpha num bsum wdf corr ntot rsq  /missing=omit.
compute alpha=make(1,1,0).
compute num=make(1,1,0).
compute bsum=make(1,1,0).
compute wdf=make(1,1,0).
compute corr=make(1,1,0).
compute ntot=make(1,1,0).
compute rsq=make(1,1,0).
compute alpha=m(:,1).
compute num=m(:,2).
compute bsum=m(:,3).
compute wdf=m(:,4).
compute corr=m(:,5).
compute ntot=m(:,6).
compute rsq=(m:,7).  
end matrix.
compute denom = (ntot-1-bsum)*wdf.
COMPUTE power = 1 - NCDF.F(IDF.F(1-ALPHA,num,denom),num,denom,((NTOT-1-bsum)*wdf*RSQ/(1-RSQ))/(1-corr)).
EXE.
formats ntot (f7.0) alpha (f5.2) num (f5.2) denom (f5.2) corr
 (f5.2)rsq (f5.2) power (f5.2).
variable labels ntot 'Total Sample Size' /alpha 'Alpha' /num 'Numerator F' /denom 'Denominator F' /corr 'Correlation' /rsq 'R-squared' /power 'Power'.
report format=list automatic align(center)
  /variables=ntot alpha num denom corr rsq power 
  /title "ANOVA power, between subjects factor possibly in a mixed design for given total sample size" .]]></screen><para><emphasis role="underline">Reference</emphasis> </para><para>Doncaster CP and Davey AJH (2007) Analysis of covariance. How to choose and construct models for the life sciences. CUP:Cambridge. </para><para>Winer BJ, Brown DR and Michels KM (1991) Statistical principles in experimental design, 3rd edition. <ulink url="https://imaging.mrc-cbu.cam.ac.uk/statswiki/FAQ/power/rmPow/statswiki/McGraw#">McGraw</ulink>-Hill:New York. </para></article>