《SAS中对应分析(含因子分析典型相关分析作业解答)》由会员分享，可在线阅读，更多相关《SAS中对应分析(含因子分析典型相关分析作业解答)（41页珍藏版）》请在金锄头文库上搜索。
1、因子分析作业解答,练习一、书p323，8-8题。用因子分析方法分析这组数据。（指出你选几个因子、你这样选择的原因、使用哪一种因子模型、做不做旋转、如何旋转、你得到的每个因子的含意。）,data data808;input x1-x5 ;label x1=力学(闭) x2=物理(闭)x3=代数(开) x4=分析(开) x5=统计(开);n=_n_;cards; 77 82 67 67 81 63 78 80 70 81 75 73 71 66 81 55 72 63 70 68 63 63 65 70 63 53 61 72 64 73 51 67 65 65 68 59 70 68 62 56
2、 62 60 58 62 70 64 72 60 62 45 52 64 60 63 54 55 67 59 62 44 50 50 64 55 63 65 63 58 56 37 31 55 60 57 73 60 64 56 54 40 44 69 53 53 53 42 69 61 55 45 62 46 61 57 45 31 49 62 63 62 44 61 52 62 46 49 41 61 49 64 12 58 61 63 67 49 53 49 62 47,54 49 56 47 53 54 53 46 59 44 44 56 55 61 36 18 44 50 57 81
3、 46 52 65 50 35 32 45 49 57 64 30 69 50 52 45 46 49 53 59 37 40 27 54 61 61 31 42 48 54 68 36 59 51 45 51 56 40 56 54 35 46 56 57 49 32 45 42 55 56 40 42 60 54 49 33 40 63 53 54 25 23 55 59 53 44 48 48 49 51 37 41 63 49 46 34 46 52 53 41 40 ;proc factor data=data808 p=0.8 method=prin priors=one simp
4、le;var x1-x5; run;,SAS程序:,Eigenvalues of the Correlation Matrix: Total = 5 Average = 1Eigenvalue Difference Proportion Cumulative1 2.61195302 1.53988724 0.5224 0.52242 1.07206578 0.50261981 0.2144 0.73683 0.56944597 0.13350618 0.1139 0.85074 0.43593980 0.12534437 0.0872 0.93795 0.31059542 0.0621 1.0
5、000,结果:,分析： 第2个因子的累计贡献率为73.68%，而且从第3个因子开始，特征根小于1。所以选2个因子。,proc factor data=data808 rotate=varimax n=2 ;var x1-x5; run;,改用这个SAS程序（不旋转结果的意义不明显）:,Rotated Factor PatternFactor1 Factor2x1 力学(闭) 0.05488 0.85565x2 物理(闭) 0.18811 0.77957x3 代数(开) 0.63399 0.56414x4 分析(开) 0.76698 0.35967x5 统计(开) 0.92948 -0.0632
6、9,结果:,因子的意义很明显，第1因子，是开卷，第2因子，是闭卷。,练习三：某医院3年中36个月的数据见hospital.xls，包括门诊人次、出院人数、病床利用率和周旋次数、平均住院天数、治愈或好转率、病死率、诊断符合率、抢救成功率，采用因子分析方法探讨综合评价指标。指出如何通过因子分析，得到你的指标。你的指标对于原数据的系数（详细写出来，如，score=0.1*门诊人次+0.2*出院人数+0.1*抢救成功率）。按照你的指标，36个样本数据排序的结果。（如24，12，33，5）。,data mylib.hospital;input t x1-x9;label t=月份 x1=门诊人次 x2=
7、出院人数 x3=病床利用率x4=病床周转数 x5=平均住院天数 x6=治愈好转率x7=病死率 x8=诊断符合率 x9=抢救成功率;cards; 1 4.34 389 99.06 1.23 25.46 93.15 3.56 97.51 61.66 2 3.45 271 88.28 0.85 23.55 94.31 2.44 97.94 73.33 3 4.38 385 103.97 1.21 26.54 92.53 4.02 98.48 76.79 36 3.77 627 86.47 1.24 23.22 91.17 3.4 98.98 89.8 ;,输入资料:,输入资料时，用了label。这将
8、极大地方便后面的分析结果。,proc factor data=mylib.hospital method=prin priors=one n=3 rotate=varimax corrscore out=hospital1;var x1-x9; run;,用这个SAS程序（不旋转结果的意义不明显）:,corr指出给出原始变量的协方差矩阵，以便分析之用。 score指出给出因子得分，以便后续的分析之用。,Eigenvalues of the Correlation Matrix: Total = 9 Average = 1Eigenvalue Difference Proportion Cumu
9、lative1 2.80742422 0.81629387 0.3119 0.31192 1.99113034 0.54280803 0.2212 0.53323 1.44832232 0.66324920 0.1609 0.69414 0.78507312 0.10436961 0.0872 0.78135 0.68070351 0.13944226 0.0756 0.85706 0.54126125 0.08823292 0.0601 0.91717 0.45302833 0.27852154 0.0503 0.96748 0.17450679 0.05595665 0.0194 0.98
10、689 0.11855013 0.0132 1.0000,选几个因子:,选3个因子时，累计贡献率超过60%，而且特征根大于1。,Rotated Factor PatternFactor1 Factor2 Factor3x1 门诊人次 0.55070 0.59527 -0.01216x2 出院人数 -0.61908 0.41630 0.23341x3 病床利用率 0.09218 0.81313 -0.01259x4 病床周转数 -0.35248 0.88650 0.07961x5 平均住院天数 0.76771 -0.22742 0.30437x6 治愈好转率 0.07020 -0.14379 0
11、.87742x7 病死率 0.19335 -0.24131 -0.73134x8 诊断符合率 0.73199 0.10834 -0.05929x9 抢救成功率 -0.80437 -0.04005 0.17998,结果:,分析： 因子1为“诊断符合率”，因子2为“病床利用率、周转数”，因子3为“治愈好转率”。 可取因子3或因子1，或它们的加权和，为评价指标。,如取因子3为评价指标，proc sort data=hospital1;by factor3; run; options ps=32 ls=85; proc print data=hospital1;var t factor3 factor
12、1 factor2 x1-x9; run;,排序结果：27, 15, 31, 19, 33, 22, 14, 16, 5, 17, 18, 2, 35, 4, 9, 32, 11, 24, 25, 28, 30, 12, 1, 34, 23, 13, 8, 3, 36, 7, 26, 20, 21, 29, 10, 6,典型相关分析作业解答,练习一、书p366，10-3题。用的是书p323，8-8题，表8.3的数据。用典型相关分析方法对闭卷（X1、X2）和开卷（X3、X4、X5）两组变量进行分析。（指出你选几对典型变量、你这样选择的原因、典型变量的意义、典型相关系数、贡献率、累计贡献率。）,d
13、ata data1003;input x1-x5 ;label x1=力学(闭) x2=物理(闭)x3=代数(开) x4=分析(开) x5=统计(开);n=_n_;cards; 77 82 67 67 81 63 78 80 70 81 75 73 71 66 81 55 72 63 70 68 63 63 65 70 63 53 61 72 64 73 41 63 49 46 34 46 52 53 41 40 ;,输入资料:,输入资料时，用了label。这将极大地方便后面的分析结果。,proc cancorr data=data1003 out=o1003;var x1 x2;with x
14、3 x4 x5; run;,用这个SAS程序进行分析:,Canonical Correlation AnalysisAdjusted Approximate SquaredCanonical Canonical Standard CanonicalCorrelation Correlation Error Correlation1 0.629913 0.600607 0.091989 0.3967902 0.095922 -.019641 0.151095 0.009201Test of H0: The canonical correlations in theEigenvalues of Inv(E)*H current row and all that follow are zero= CanRsq/(1-CanRsq)Likelihood ApproximateEigenvalue Difference Proportion Cumulative Ratio F Value Num DF Den DF Pr F1 0.6578 0.6485 0.9861 0.9861 0.59766007 3.82 6 78 0.00222 0.0093 0.0139 1.0000 0.99079906 0.19 2 40 0.8312,
