典型相关分析的SAS实现由刀豆文库小编整理，希望给你工作、学习、生活带来方便，猜你可能喜欢“sas典型相关分析”。
典型相关分析的SAS实现
data ex20_1(type=corr);input _name_$3.x1 x2 y1-y4;_type_='corr';cards;x1 1.0 0.8491 0.5106 0.2497 0.5285 0.3019 x2 0.8491 1.0 0.8062 0.5438 0.7887 0.6064 y1 0.5106 0.8062 1.0 0.7833 0.9284 0.8364 y2 0.2497 0.5438 0.7833 1.0 0.6457 0.9051 y3 0.5285 0.7887 0.9284 0.6457 1.0 0.7097 y4 0.3019 0.6064 0.8364 0.9051 0.7097 1.0;proc cancorr edf=102;var x1 x2;with y1-y4;run;
data ex20_2(type=corr);input _name_$ x1-x9 y1-y5;_type_='corr';cards;x1 1-0.027 0.052 0.037 0.003 0.009-0.030 0.014 0.020 0.034 0.071 0.054-0.022 0.012 x2-0.027 1-0.065-0.057-0.016 0.014 0.053-0.080-0.067-0.023 0.019 0.046-0.034 0.012 x3 0.052-0.065 1 0.772 0.018 0.042-0.094 0.279 0.163 0.221 0.093-0.030 0.040-0.118 x4 0.037-0.057 0.772 1 0.003 0.044-0.080 0.098 0.040 0.224 0.116-0.013-0.013-0.096 x5 0.003-0.016 0.018 0.003 1-0.017-0.030 0.019 0.025-0.044-0.021 0.002-0.019 0.002 x6 0.009 0.014 0.042 0.044-0.017 1-0.020 0.009-0.018 0.032 0.047 0.009-0.037-0.015 x7-0.030 0.053-0.094-0.080-0.030-0.020 1-0.057 0.013 0.009-0.032 0.089 0.027 0.014 x8 0.014-0.080 0.279 0.098 0.019 0.009-0.057 1 0.447-0.018-0.036-0.031 0.156-0.041 x9 0.020-0.067 0.163 0.040 0.025-0.018 0.013 0.447 1-0.063-0.017-0.027 0.117-0.019 y1 0.034-0.023 0.221 0.224-0.044 0.032 0.009-0.018-0.063 1 0.156-0.070-0.070-0.116
y2 0.071 0.019 0.093 0.116-0.021 0.047-0.032-0.036-0.017 0.156 1-0.039-0.046-0.052 y3 0.054 0.046-0.030-0.013 0.002 0.009 0.089-0.031-0.027-0.070-0.039 1 0.001 0.013 y4-0.022-0.034 0.040-0.013-0.019-0.037 0.027 0.156 0.117-0.070-0.046 0.001 1-0.019 y5 0.012 0.012-0.118-0.096 0.002-0.015 0.014-0.041-0.019-0.116-0.052 0.013-0.019 1;proc cancorr data=ex20_2 edf=1494 red;var x1-x9;with y1-y5;run;
data ex21_1;input group num;do i=1 to num;input value@@;output;end;cards;1 25 6.5 13.5 12.8 6.2 13.9 14.7 9.5 9.0 6.9 16.8 13.3 10.8 12.2 14.9 13.7 12.8 5.3 11.8 12.4 7.6 13.3 11.9 11.2 12.3 12.7 0 20 8.5 6.4 4.6 1.7 9.7 5.3 4.9 5.7 3.8 6.5 6.3 5.4 3.3 4.7 8.6 6.3 5.9 4.8 4.5 5.2;proc logistic descending;model group=value/scale=none outroc=roc1;run;proc print;run;proc gplot;plot _SENSIT_*_1MSPEC_;run;
