SAS使用学习笔记(对应分析).doc

SAS使用学习笔记（对应分析）1 对应分析是不仅研究变量之间的关系、还要研究样品之间的关系。它通过在同一个直角坐标系内同时表达出变量与样品两者之间的相互关系。2 对应分析例子下面是某研究者收集到的资料，试分析各种基因频率与民族之间的关系。各民族下面的小数是44种基因出现的频率。基因型（JY）藏族(Z) 尼泊尔(N) 印度(Y) 汉族(H)。程序：DATA b; INPUT jy $ 1-3 z 6-11 n 14-19 y 22-27 h 30-35; cards;A1 0.0308 0.01800.11900.0149A2 0.3333 0.10700.14800.3492A3 0.0204 0.01900.10100.0176A9 0.3037 0.27900.15600.1414A100.0409 0.01800.03900.0313A110.1354 0.42200.12600.2977A280.0000 0.01800.08300.0094A300.0413 0.00000.00000.0217A310.0518 0.03700.02200.0121A320.0000 0.01900.03900.0013A330.0000 0.06700.08300.0608B5 0.2828 0.11800.13400.0825B7 0.0000 0.01900.08000.0244B8 0.0102 0.01180.04500.0094B120.0102 0.03700.06600.0121B130.0102 0.07700.00600.0650B140.0000 0.00000.00600.0013B150.1923 0.25400.09600.1092B180.0050 0.02800.02200.0000B270.1067 0.00000.02600.0204B350.0626 0.05700.14800.0342B370.0102 0.01800.00900.0067B380.04650.0470 0.00300.0015B390.01020.0000 0.00900.0176B460.01020.00900.00000.1813B480.05720.15000.00300.0108B500.01020.01800.03700.0000B530.00500.0000 0.00600.0000B540.01530.00000.00000.0176B550.05720.02800.02600.0217B560.01020.00900.00600.0040B570.00500.01800.03900.0341B580.00000.06700.03300.0139B600.06260.02800.02200.0723B610.08990.00000.08300.1080B700.00500.00000.00800.0000C1 0.08990.03700.02300.1716C2 0.02040.00000.07300.0397C3 0.17980.10700.08300.3269C4 0.16510.07700.13400.0495C5 0.00000.00900.01600.0054C6 0.02560.24500.04500.0081C7 0.17120.21800.11900.1152C8 0.00500.00000.00400.0027;run;PROC CORRESP data=b OUTC=ccc;VAR z n y h;LABEL z=藏族 n=尼泊尔 y=印度 h=汉族;ID jy;RUN;DATA ccc;SET ccc;X=dim1;Y=dim2;XSYS =2;YSYS =2;TEXT =jy;SIZE =2;LABEL X=Dimension 1 Y=Dimension 2;keep X Y TEXT XSYS YSYS SIZE;RUN;PROC GPLOT DATA=ccc; SYMBOL1 V=#; AXIS1 LENGTH=5 IN ORDER=-1.3 TO 1.3 BY 0.2; AXIS2 LENGTH=5 IN ORDER=-1.3 TO 1.3 BY 0.2; PLOT Y*X=1 / ANNOTATE=ccc FRAME HAXIS=AXIS1 VAXIS=AXIS2 HREF=0 VREF=0;RUN;输出： The CORRESP Procedure Inertia and Chi-Square Decomposition Singular Principal Chi- Cumulative Value Inertia Square Percent Percent 8 16 24 32 40 -+-+-+-+-+- 0.42302 0.17895 1.83072 41.61 41.61 * 0.39266 0.15418 1.57736 35.85 77.46 * 0.31137 0.09695 0.99184 22.54 100.00 * Total 0.43007 4.39992 100.00 Degrees of Freedom = 129 SAS 系统 2008年05月04日 星期日 下午02时39分56秒 The CORRESP Procedure Row Coordinates Dim1 Dim2 A1 0.5878 0.8300 A2 -0.4233 0.1209 A3 0.5708 0.7845 A9 0.0741 -0.1562 A10 -0.0466 0.2818 A11 0.0191 -0.3738 A28 0.8269 0.8844 A30 -0.7954 0.0516 A31 0.1126 -0.1366 A32 0.9506 0.5206 A33 0.2931 0.1666 B5 -0.0085 0.0766 B7 0.5508 0.7576 B8 0.5282 0.6647 B12 0.6453 0.3736 B13 -0.1439 -0.5749 B14 0.6310 1.2240 B15 0.1193 -0.3016 B18 0.8250 -0.1198 B27 -0.3381 0.2782 B35 0.4306 0.4728 B37 0.2494 -0.2481 B38 0.1990 -0.6104 B39 -0.5396 0.4166 B46 -1.2499 0.0376 B48 0.4093 -0.9132 B50 0.7714 0.4510 B53 0.3664 0.8090 B54 -0.9706 0.0684 B55 -0.0427 0.0207 B56 0.1278 -0.1051 B57 0.0754 0.3742 B58 0.6482 -0.3801 B60 -0.4505 0.0137 B61 -0.3745 0.4819 B70 0.4748 0.9100 C1 -0.7041 0.0091 C2 0.0968 0.8397 C3 -0.5193 0.0178 C4 0.1452 0.2307 C5 0.5711 0.3791 C6 0.7305 -0.8444 C7 0.1169 -0.1805 C8 -0.1612 0.5353 SAS 系统 2008年05月04日 星期日 下午02时39分56秒 The CORRESP Procedure Summary Statistics for the Row Points Quality Mass Inertia A1 0.9966 0.0179 0.0431 A2 0.9583 0.0916 0.0431 A3 0.9744 0.0154 0.0347 A9 0.3838 0.0860 0.0156 A10 0.8715 0.0126 0.0027 A11 0.6324 0.0959 0.0494 A28 0.9068 0.0108 0.0406 A30 0.5695 0.0062 0.0160 A31 0.1562 0.0120 0.0056 A32 0.9073 0.0058 0.0174 A33 0.2820 0.0206 0.0193 B5 0.0274 0.0603 0.0304 B7 0.7996 0.0121 0.0308 B8 0.9666 0.0075 0.0129 B12 0.9043 0.0122 0.0175 B13 0.5787 0.0155 0.0218 B14 0.9031 0.0007 0.0035 B15 0.8864 0.0637 0.0176 B18 0.9645 0.0054 0.0090 B27 0.1832 0.0150 0.0364 B35 0.9998 0.0295 0.0281 B37 0.9998 0.0043 0.0012 B38 0.5599 0.0096 0.0164 B39 0.9636 0.0036 0.0040 B46 0.6708 0.0196 0.1062 B48 0.9847 0.0216 0.0511 B50 0.9990 0.0064 0.0118 B53 0.7205 0.0011 0.0027 B54 0.9637 0.0032 0.0073 B55 0.0150 0.0130 0.0045 B56 0.3142 0.0029 0.0006 B57 0.3852 0.0094 0.0083 B58 0.7578 0.0111 0.0193 B60 0.9902 0.0181 0.0086 B61 0.9992 0.0275 0.0238 B70 0.8576 0.0013 0.0036 C1 0.9447 0.0314 0.0383 C2 0.9242 0.0130 0.0234 C3 0.9069 0.0681 0.0471 C4 0.3959 0.0416 0.0182 C5 0.6915 0.0030 0.0047 C6 0.9669 0.0316 0.0949 C7 0.9160 0.0609 0.0072 C8 0.7360 0.0011 0.0011 SAS 系统 2008年05月04日 星期日 下午02时39分56秒 The CORRESP Procedure Partial Contributions to Inertia for the Row Points Dim1 Dim2 A1 0.0345 0.0798 A2 0.0918 0.0087 A3 0.0281 0.0616 A9 0.0026 0.0136 A10 0.0002 0.0065 A11 0.0002 0.0869 A28 0.0412 0.0547 A30 0.0218 0.0001 A31 0.0009 0.0015 A32 0.0293 0.0102 A33 0.0099 0.0037 B5 0.0000 0.0023 B7 0.0204 0.0449 B8 0.0116 0.0214 B12 0.0285 0.0111 B13 0.0018 0.0332 B14 0.0016 0.0069 B15 0.0051 0.0376 B18 0.0204 0.0005 B27 0.0096 0.0075 B35 0.0306 0.0428 B37 0.0015 0.0017 B38 0.0021 0.0231 B39 0.0059 0.0040 B46 0.1711 0.0002 B48 0.0202 0.1168 B50 0.0212 0.0084 B53 0.0008 0.0046 B54 0.0169 0.0001 B55 0.0001 0.0000 B56 0.0003 0.0002 B57 0.0003 0.0085 B58 0.0261 0.0104 B60 0.0205 0.0000 B61 0.0215 0.0414 B70 0.0016 0.0068 C1 0.0871 0.0000 C2 0.0007 0.0595 C3 0.1026 0.0001 C4 0.0049 0.0144 C5 0.0054 0.0028 C6 0.0944 0.1463 C7 0.0047 0.0129 C8 0.0002 0.0021 S