I used the Socioeconomic data provided in the sample data and ran the following script.
//dt=($Sample data/socioeconomic.jmp")
fa=dt<<Factor Analysis(
Y(
:Total Population,
:Median School Years,
:Total Employment,
:Professional Services,
:Median House Value
),
Estimation Method( "Row-wise" ),
Variance Scaling( "Correlations" ),
Fit( "ML", "SMC", 2, "Varimax" ),
SendToReport(
Dispatch(
{"Factor Analysis on Correlations with 2 Factors: Maximum Likelihood ",
"Factor Loading Plot"},
"Loading Plot for Rotated Factor",
FrameBox,
{Frame Size( 52, 36 )}
)
)
);
fa << (Fit[1] << Save Rotated Components);