title "Power analysis: Lord's Vocabulary Data";
data lord(type=cov);
   input _type_ $ _name_ $  x1 x2 y1 y2;
datalines;
n    .  649     649     649     649
cov  x1 86.3937   .       .       .
cov  x2 57.7751 86.2632   .       .
cov  y1 56.8651 59.3177 97.2850   .
cov  y2 58.8986 59.6683 73.8201 97.8192
mean .   0       0       0       0
;
title2 "Lord's data: H1- X1 and X2 parallel, Y1 and Y2 parallel, rho=1";
proc calis data=lord cov summary outram=ram1;
   lineqs  x1  = betax F1  + e1,
           x2  = betax F1  + e2,
           y1  = betay F2  + e3,
           y2  = betay F2  + e4;
   std  F1 F2 = 1 1,
        e1 e2 e3 e4 = vex vex vey vey;
   cov  F1 F2 = 1;
run;

*-- Perform power analysis for the RMSEA statistics in this model;
title 'Retrospective power analysis';
%csmpower(data=ram1);

*--; title  'Prospective power analysis';
%csmpower(df=6, rmseaa=%str(.08 to .12 by .02),
          plot=%str(power*n =rmseaa));