title "sas program for power estimation for csm" ;
data one ;
alpha=.05 ;  *significance level ;
rmsea0=.05 ; *null hyp value ;
rmseaa=.10 ; *alt hyp value ;
d=50 ; *degrees of freedom ;
n=200 ; *sample size ;
    ncp0=(n-1)*d*rmsea0**2 ;
    ncpa=(n-1)*d*rmseaa**2 ;
    if rmsea0<rmseaa then do ;
    cval=cinv(1-alpha,d,ncp0) ;
    power=1-probchi(cval,d,ncpa) ;
    end ;
    if rmsea0>rmseaa then do ;
    cval=cinv(alpha,d,ncp0) ;
    power=probchi(cval,d,ncpa) ;
    end ;
    output ;
proc print data=one ; var rmsea0 rmseaa alpha d n power ; run ;