title “power eatimation for csm”;
data one ;
alpha=.05 ; *significance level ;
rmsea0=.05 ; *null hyp value ;
rmseaa=.08 ; *alt hyp value ;
d=12 ; *degree of freedom
n=141 ; *sample size ;
ncp0=(n-1)*d*rmsea0**2 ;
ncpa=(n-1)*d*rmseaa**2 ;
if rmsea<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

是根据文献计算样本量统计效能的一个命令，求大佬们帮一下。