*sizeanova_fixed.sps  ----------------------------------.
title "”팱ŽÒŠÔ—vˆö‚QŒ³•ªŽU•ª�Í(fixed effect)‚Ì•K—vƒTƒ“ƒvƒ‹�”".
*  effect size index f=.
* small f=.10 : medium f=.25: large f=.40  Cohen(1988).
preserve.
set printback=none.

data list free/name(a6) esize alpha powd f1g f2g.
begin data
small      .10 0.05 .80  3 5
medium     .25 0.05 .80  3 5 
large      .40 0.05 .80  3 5
end data.
string ftype(a8).
compute ixj=(f1g-1)*(f2g-1).
compute nixj=f1g*f2g.
compute f1g1=f1g-1.
compute f2g1=f2g-1.
do if $casenum = 1 . 
write //'                ‚QŒ³•ªŽU•ª�Í(fixed effect)‚Ì•K—vƒTƒ“ƒvƒ‹�”'
//
'—vˆö           Œø‰Ê—Êf  ƒ¿ power �…�€�”  dfh   dfe    1ŒQn      ‘Sn    �®�”‘Sn'.
end if.
do repeat
 dfh=f1g1 f2g1 ixj/ ft='f1' 'f2' 'f1xf2'/level=f1g f2g nixj.
compute ftype=ft.
*initialize values.
compute powa=0.0.
compute n=0.
compute f2=esize**2.
compute alpha2=alpha.
loop  if  (powa<=powd).
  compute n=n +100.
  compute dfe=n-nixj.

  compute lamda=f2*n.
  compute f_alpha=idf.f(1-alpha2,dfh,dfe).
  compute powa=1-ncdf.f(f_alpha,dfh,dfe,lamda).
end loop.

*begin loop for interval halving.
compute dir=-1.
compute newn=n.
compute intv=200.
compute powdiff=powa-powd.
loop if (powdiff>=.000001).
  compute  intv=intv*.5.
  compute  newn=newn +dir*intv*.5.
*compute new power.
  compute dfe=newn-nixj.
  compute lamda=f2*newn.
*compute power .
  compute f_alpha=idf.f(1-alpha2,dfh,dfe).
  compute powa=1-ncdf.f(f_alpha,dfh,dfe,lamda).
  compute powdiff=abs(powa-powd).
  do if (powa<powd).
    compute  dir=1.
    else .
    compute dir= -1.
  end if.
end loop.
compute minn=newn/nixj.
compute newnn=trunc(minn+0.99999)*nixj.
format esize powd powa alpha (f6.3).
format dfe(f8.2).
format dfh level (f6.0).
format newnn(f8.0).
format ftype (a8).
format newn(f10.2).
write /ftype name esize alpha powa level dfh dfe minn newn newnn.
do if (ft='f1xf2') .
    write /'  ---------------------------------------------------------------------------'.
end if.
end repeat.
exec.
restore.
*-------------------------------------------------.