FAQ/power/rmPow - CBU statistics Wiki

Upload page content

You can upload content for the page named below. If you change the page name, you can also upload content for another page. If the page name is empty, we derive the page name from the file name.

File to load page content from
Page name
Comment
Type the missing letters from: Lodon is th capial of nglnd

Revision 7 as of 2006-10-17 14:32:15

location: FAQ / power / rmPow
  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)" .