SAS中典型相關(guān)分析可以先上本章,再上對(duì)應(yīng)分析ppt課件

1、用SAS進(jìn)行典型相關(guān)分析書p358，例10.2.1，表10.1，礦區(qū)下部礦Ptx1）、Pdx2與Cuy1）、Niy2的數(shù)據(jù)。典型相關(guān)分析研究共生規(guī)律。nx1x2y1y21 0.140.3 0.03 0.1420.20.5 0.14 0.223 0.06 0.11 0.03 0.024 0.07 0.11 0.04 0.135 0.12 0.22 0.06 0.126 0.52 0.87 0.190.27 0.23 0.47 0.140.18 1.19 0.38 0.09 0.119 0.37 0.66 0.14 0.1510 0.360.6 0.12 0.1411 0.42 0.77 0.17

2、0.112 0.35 0.850.3 0.19130.5 0.87 0.23 0.2214 0.56 1.15 0.29 0.28nx1x2y1y2150.430.90.130.22160.470.970.260.22170.490.790.210.2180.470.770.510.22190.40.880.330.19200.661.30.210.3210.631.30.450.28220.521.430.310.23230.440.870.170.25240.030.070.050.08250.20.280.040.08260.040.10.110.07270.170.280.150.09

3、輸入資料:n=_n_;這一行指出，對(duì)數(shù)據(jù)集加入一個(gè)新變量一列），變量名為n，其值等于樣本的實(shí)際編號(hào)。data d1021; input x1 x2 y1 y2; n=_n_; cards;0.14 0.30 0.03 0.140.20 0.50 0.14 0.220.06 0.11 0.03 0.020.07 0.11 0.04 0.130.12 0.22 0.06 0.120.52 0.87 0.19 0.200.23 0.47 0.14 0.101.19 0.38 0.09 0.110.37 0.66 0.14 0.150.36 0.60 0.12 0.140.42 0.77 0.17 0

4、.100.35 0.85 0.30 0.190.50 0.87 0.23 0.220.56 1.15 0.29 0.280.43 0.90 0.13 0.220.47 0.97 0.26 0.22/*未完未完*/*續(xù)前*/0.49 0.79 0.21 0.200.47 0.77 0.51 0.220.40 0.88 0.33 0.190.66 1.30 0.21 0.300.63 1.30 0.45 0.280.52 1.43 0.31 0.230.44 0.87 0.17 0.250.03 0.07 0.05 0.080.20 0.28 0.04 0.080.04 0.10 0.11 0.0

5、70.17 0.28 0.15 0.09;典型相關(guān)分析在SAS中用cancorr過(guò)程：proc cancorr data=d1021 simple corr vprefix=v wprefix=w out=o1021 ; var x1 x2; with y1 y2;run;proc plot data=o1021; plot w1*v1 $ n=* / vref=0 href=0;run;cancorr過(guò)程中的simple指出，分析各原始變量的均值和標(biāo)準(zhǔn)差并輸出。cancorr過(guò)程中的corr指出，計(jì)算原始變量間的協(xié)方差矩陣，并輸出。vprefix=v （這是默認(rèn)值，故可省指出，第一組典型變量

6、對(duì)應(yīng)第一組原始變量的名稱叫作，v1、v2、。wprefix=w這是默認(rèn)值，故可省指出，第二組典型變量對(duì)應(yīng)第二組原始變量的名稱叫作，w1、w2、。var x1 x2;with y1 y2;表示第一組原始變量是x1、x2。第二組變量是y1、y2。plot過(guò)程把各樣本以v1、w1為坐標(biāo)畫出來(lái)。href=0表示在橫坐標(biāo)v1=0處畫一條垂線，vref=0表示在縱坐標(biāo)w1=0處畫一條垂線。$n=*表示每個(gè)點(diǎn)在圖上用*表示，并且在*后顯示該樣本點(diǎn)的n變量的值。這是各原始變量的均值、標(biāo)準(zhǔn)差。這是cancorr過(guò)程中的simple選項(xiàng)所要求輸出的。輸出的一部分： Means and Standard Devia

7、tions Standard Variable Mean Deviation x1 0.371852 0.249014 x2 0.659259 0.396435 y1 0.181481 0.123870 y2 0.168519 0.073626這是各原始變量間的協(xié)方差矩陣。這是cancorr過(guò)程中的corr選項(xiàng)所要求輸出的。輸出的一部分： Correlations Among the Original Variables Correlations Among the VAR Variables x1 x2 x1 1.0000 0.6209 x2 0.6209 1.0000 Correlatio

8、ns Among the WITH Variables y1 y2 y1 1.0000 0.6818 y2 0.6818 1.0000 Correlations Between the VAR Variables and the WITH Variables y1 y2 x1 0.4477 0.5273 x2 0.7471 0.8691我們看到，第一典型相關(guān)系數(shù)為0.894618，貢獻(xiàn)率為100%，以99.99%的顯著性非0。而第二典型相關(guān)系數(shù)的貢獻(xiàn)率為0.00%，只以3.67%的顯著性非0。典型相關(guān)系數(shù)，及其檢驗(yàn)： Canonical Correlation Analysis Adjuste

9、d Approximate Squared Canonical Canonical Standard Canonical Correlation Correlation Error Correlation 1 0.894618 0.889469 0.039156 0.800341 2 0.009496 . 0.196098 0.000090 Test of H0: The canonical correlations in the Eigenvalues of Inv(E)*H current row and all that follow are zero = CanRsq/(1-CanRs

10、q) Likelihood Approximate Eigenvalue Difference Proportion Cumulative Ratio F Value Num DF Den DF Pr F 1 4.0085 4.0084 1.0000 1.0000 0.20194137 14.24 4 46 F 1 1.7247 1.6828 0.9734 0.9734 0.35039053 2.05 9 34.223 0.0635 2 0.0419 0.0366 0.0237 0.9970 0.95472266 0.18 4 30 0.9491 3 0.0053 0.0030 1.0000

11、0.99473355 0.08 1 16 0.7748PHYS1=-0.031404688 weight + 0.4932416756 waist -0.008199315 pulseEXER1=-0.066113986 chins -0.016846231 situps + 0.0139715689 jumps典型變量關(guān)于原始變量的系數(shù)： Raw Canonical Coefficients for the 生理指標(biāo) PHYS1 PHYS2 PHYS3 weight 體重 -0.031404688 -0.076319506 -0.007735047 waist 腰圍 0.4932416756

12、 0.3687229894 0.1580336471 pulse 脈博 -0.008199315 -0.032051994 0.1457322421 Raw Canonical Coefficients for the 訓(xùn)練指標(biāo) EXER1 EXER2 EXER3 chins 單杠 -0.066113986 -0.071041211 -0.245275347 situps 仰臥起坐 -0.016846231 0.0019737454 0.0197676373 jumps 跳高 0.0139715689 0.0207141063 -0.008167472典型變量關(guān)于標(biāo)準(zhǔn)化的原始變量的系數(shù)：PHY

13、S1近似地代表waist與weight的加權(quán)差；EXER1主要代表situps仰臥起坐）。PHYS1=-0.7754 weight* + 1.5793 waist* -0.0591 pulse*EXER1=-0.3495 chins* -1.0540 situps* + 0.7164 jumps* Standardized Canonical Coefficients for the 生理指標(biāo) PHYS1 PHYS2 PHYS3 weight 體重 -0.7754 -1.8844 -0.1910 waist 腰圍 1.5793 1.1806 0.5060 pulse 脈博 -0.0591 -0

14、.2311 1.0508 Standardized Canonical Coefficients for the 訓(xùn)練指標(biāo) EXER1 EXER2 EXER3 chins 單杠 -0.3495 -0.3755 -1.2966 situps 仰臥起坐 -1.0540 0.1235 1.2368 jumps 跳高 0.7164 1.0622 -0.4188典型結(jié)構(gòu)一，典型變量與同組的原始變量的相關(guān)系數(shù)： Correlations Between the 生理指標(biāo) and Their Canonical Variables PHYS1 PHYS2 PHYS3 weight 體重 0.6206 -0.

15、7724 -0.1350 waist 腰圍 0.9254 -0.3777 -0.0310 pulse 脈博 -0.3328 0.0415 0.9421 Correlations Between the 訓(xùn)練指標(biāo) and Their Canonical Variables EXER1 EXER2 EXER3 chins 單杠 -0.7276 0.2370 -0.6438 situps 仰臥起坐 -0.8177 0.5730 0.0544 jumps 跳高 -0.1622 0.9586 -0.2339典型結(jié)構(gòu)二，典型變量與異組的原始變量的相關(guān)系數(shù)： Correlations Between the

16、 生理指標(biāo) and the Canonical Variables of the 訓(xùn)練指標(biāo) EXER1 EXER2 EXER3 weight 體重 0.4938 -0.1549 -0.0098 waist 腰圍 0.7363 -0.0757 -0.0022 pulse 脈博 -0.2648 0.0083 0.0684 Correlations Between the 訓(xùn)練指標(biāo) and the Canonical Variables of the 生理指標(biāo) PHYS1 PHYS2 PHYS3 chins 單杠 -0.5789 0.0475 -0.0467 situps 仰臥起坐 -0.6506

17、0.1149 0.0040 jumps 跳高 -0.1290 0.1923 -0.0170 Standardized Variance of the 生理指標(biāo) Explained by Their Own The Opposite Canonical Variables Canonical Variables Canonical Variable Cumulative Canonical Cumulative Number Proportion Proportion R-Square Proportion Proportion 1 0.4508 0.4508 0.6330 0.2854 0.2

18、854 2 0.2470 0.6978 0.0402 0.0099 0.2953 3 0.3022 1.0000 0.0053 0.0016 0.2969 Standardized Variance of the 訓(xùn)練指標(biāo) Explained by Their Own The Opposite Canonical Variables Canonical Variables Canonical Variable Cumulative Canonical Cumulative Number Proportion Proportion R-Square Proportion Proportion 1

19、 0.4081 0.4081 0.6330 0.2584 0.2584 2 0.4345 0.8426 0.0402 0.0175 0.2758 3 0.1574 1.0000 0.0053 0.0008 0.2767典型冗余分析，原始變量被同組、異組的典型變量所解釋的比例：練習(xí)一、書p366，10-3題。用的是書p323，8-8題，表8.3的數(shù)據(jù)。用典型相關(guān)分析方法對(duì)閉卷X1、X2和開卷X3、X4、X5兩組變量進(jìn)行分析。（指出你選幾對(duì)典型變量、你這樣選擇的原因、典型變量的意義、典型系數(shù)、貢獻(xiàn)率、累計(jì)貢獻(xiàn)率。）x1x2x3x4x5778267678163788070817573716681557263706863636570635361726473516765656859706862566260586270647260624552646063545567596244505064556365635856373155605773x1x2x3x4x5606456544044695353534269615545624661574531496263624461526246494161496412586163674953496247544956475354534659444456556136184450578146526550353245495764x1x2x3x