Diff for "FAQ/power/rmPow" - CBU statistics Wiki
location: Diff for "FAQ/power/rmPow"
Differences between revisions 5 and 7 (spanning 2 versions)
Revision 5 as of 2006-10-05 15:18:40
Size: 199
Editor: PeterWatson
Comment:
Revision 7 as of 2006-10-17 14:32:15
Size: 1460
Editor: PeterWatson
Comment:
Deletions are marked like this. Additions are marked like this.
Line 1: Line 1:
 1. [:FAQ/power/rmPow/bpow: Sample Size for anova term involving only between subject factors]  1. [:FAQ/power/rmPowN/bN: Sample Size for anova term involving only between subject factors]
Line 4: Line 4:

 * Num(erator) is (nb(1)-1)...(nb(i)-1) for i between subject factors with categories nb(1) to nb(i)

 * Denom(inator) is (Ntot- 1- (sum of between subject factor levels) + number of between subject factors) for a total sample of size Ntot

 * Rsq is r-squared/100 so 5.9% = 0.059

{{{

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)" .
}}}
  1. [:FAQ/power/rmPowN/bN: Sample Size for anova term involving only between subject factors]
  2. [:FAQ/power/rmPowN/wN: Sample Size for anova term involving at least one within subject factor]
  3. Num(erator) is (nb(1)-1)...(nb(i)-1) for i between subject factors with categories nb(1) to nb(i)
  4. Denom(inator) is (Ntot- 1- (sum of between subject factor levels) + number of between subject factors) for a total sample of size Ntot
  5. Rsq is r-squared/100 so 5.9% = 0.059

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)