SEM

data illfl (type=corr);
input _type_ $1-4 _name_ $6-13
exercise 15-20 hardy 22-27
fitness 29-34 stress 36-41
illness 43-48;
cards;
n             373    373    373    373    373
mean          40.90  0.00   67.10  4.80   716.7
std           66.50  3.80   18.40  6.70   624.8
corr exercise 1.00    
corr hardy    -0.03  1.00
corr fitness  0.39   0.07   1.00
corr stress    -0.05  -0.23  -0.13  1.00
corr illness  -0.08  -0.16  -0.29  0.34   1.00
;;;;

/* Confirmatory Factor Analysis */

proc calis data=illfl corr;
lineqs
hardy = p1 F1 + e1,
stress = p2 F1 + e2,
illness = p3 F1 + e3,
fitness = p4 F2 + e4,
exercise= p5 F2 + e5;
std
e1-e5 = vare1-vare5,
F1 = 1,
F2 = 1;
cov
F1 F2 = covf1f2;
VAR exercise fitness hardy stress illness;
run;

/* Exploratory Factor Analysis */

proc factor data=illfl method=ml scree priors=smc;
  var hardy stress illness fitness exercise;
run;

proc factor data=illfl method=ml rotate=v n=2 reorder plot priors=smc;
  var hardy stress illness fitness exercise;
run;