Diff for "FAQ/power/rmPow" - CBU statistics Wiki
location: Diff for "FAQ/power/rmPow"
Differences between revisions 15 and 17 (spanning 2 versions)
Revision 15 as of 2007-01-15 10:04:39
Size: 1618
Editor: PeterWatson
Comment:
Revision 17 as of 2007-01-15 17:11:35
Size: 1733
Editor: PeterWatson
Comment:
Deletions are marked like this. Additions are marked like this.
Line 5: Line 5:
 * num(erator) is product over factors comprising term of interest of (numer of levels of factor -1) For B between subject factors with levels $$b_{i}$$, i=1, ..., B and W with subject factors with levels $$w_{i}$$, j=1, ..., W
Line 7: Line 7:
 * d1 is sum over the factors of (number of levels of '''between''' subject factor - 1)  * num(erator) = $$ \prod_{\mbox{factors}} $$ (number of levels of factor -1) in term of interest
Line 9: Line 9:
 * d2 is product over the factors of (number of levels of '''within''' subject factor - 1)  *
d1 = $$\sum_{i}^{B} (b_{i} - 1) $$ if B > 0 in anova
   = 0 otherwise
Line 11: Line 13:
 * d2 = $$ \prod_{j} (w_{j} - 1) $$ if W > 0 in term of interest
      = 1 otherwise
  
Line 13: Line 18:
Power can also be computed using a [attachment:ao.xls spreadsheet.] Power can also be computed using a [attachment:aov.xls spreadsheet.]
  • alpha is likelihood of making a type I error (usually = 0.05)
  • etasq is partial eta-squared/100 so, for example, 5.9% = 0.059

For B between subject factors with levels $$b_{i}$$, i=1, ..., B and W with subject factors with levels $$w_{i}$$, j=1, ..., W

  • num(erator) = $$ \prod_{\mbox{factors}} $$ (number of levels of factor -1) in term of interest

d1 = $$\sum_{i}^{B} (b_{i} - 1) $$ if B > 0 in anova

  • = 0 otherwise
  • d2 = $$ \prod_{j} (w_{j} - 1) $$ if W > 0 in term of interest

    • = 1 otherwise
  • ntot is the total sample size

Power can also be computed using a [attachment:aov.xls spreadsheet.]

[ COPY AND PASTE THE BOXED BELOW SYNTAX BELOW INTO A SPSS SYNTAX WINDOW AND RUN; ADJUST INPUT DATA AS REQUIRED]

DATA LIST free
/alpha etasq num d1 d2 ntot. 
BEGIN DATA.
.05 0.02 2 3 2 58 
.05 0.12 2 3 2 58
END DATA.
set errors=none. 
matrix.
get m /variables=alpha etasq num d1 d2 ntot /missing=omit.
compute alpha=make(1,1,0).
compute etasq=make(1,1,0).
compute num=make(1,1,0).
compute d1=make(1,1,0).
compute d2=make(1,1,0).
compute ntot=make(1,1,0).
compute alpha=m(:,1).
compute etasq=m(:,2).
compute num=m(:,3).
compute d1=m(:,4).
compute d2=m(:,5).
compute ntot=m(:,6).
end matrix.
compute denom=(ntot-d1-1)*d2.
COMPUTE power = 1 - NCDF.F(IDF.F(1-ALPHA,num,denom),num,denom,denom*etasq/(1-etasq)).
EXE.
formats alpha (f5.2) num (f5.2) denom (f5.2) etasq (f5.2) power (f5.2).
variable labels ntot 'Total Sample' /alpha 'Alpha' /num 'Numerator F' /denom 'Denominator F' /etasq 'Partial Eta-Sq' /power 'Power'.
report format=list automatic align(center)
  /variables=ntot alpha num denom etasq power 
  /title "ANOVA power (any anova)" .

None: FAQ/power/rmPow (last edited 2013-03-08 10:18:01 by localhost)