主成分分析SAS實(shí)驗(yàn).ppt

1、主成分分析,公共衛(wèi)生學(xué)院信息數(shù)據(jù)處理教學(xué)實(shí)驗(yàn)室,一、主成分分析,實(shí)際工作中原始數(shù)據(jù)的變量之間常有一定的相關(guān)性。人們希望找到較少的幾個(gè)互不相關(guān)的綜合指標(biāo)，盡可能多的反映原來(lái)的信息。 主成分分析就是由原變量X1Xp中線性組合出m個(gè)(mp)互不相關(guān)、且盡量少丟失信息的新變量(主成分)，并能給各主成分所包含的信息以恰當(dāng)?shù)膶I(yè)解釋。,主成分分析和因子分析也可以用下列各種統(tǒng)計(jì)分析的中間結(jié)果矩陣進(jìn)行分析：,CORR 相關(guān)系數(shù)矩陣 SSCP 平方和、積和矩陣 CSSCP 離均差平方和、積和矩陣 COV 方差、協(xié)方差矩陣 UCOV 為平方和、積和矩陣/n UCORR 為 矩陣 FACTOR 因子矩陣 需要在數(shù)據(jù)

2、步中指定： _TYPE_= CORR;,主成分分析過(guò)程 PROC PRINCOMP; VAR X1-X4； RUN; PRINCOMP過(guò)程后選項(xiàng)： COV 用協(xié)方差矩陣計(jì)算，默認(rèn)用相關(guān)系數(shù)矩陣 PREFIX= 指定主成分變量名前綴，默認(rèn)用prin1.print2 STD 輸出變量標(biāo)準(zhǔn)化后的主成分得分 OUT= 把原始數(shù)據(jù)和主成分得分輸出到指定數(shù)據(jù)集， 原始數(shù)據(jù)為CORR或COV矩陣時(shí)OUT無(wú)效 OUTSTAT= 把運(yùn)算結(jié)果輸出到指定數(shù)據(jù)集 NOPRINT 禁止屏幕輸出任何結(jié)果,保留主成分個(gè)數(shù)的原則： 特征根(eigenvalue)1 累計(jì)貢獻(xiàn)率(cumulative proportion of

3、 eigenvalue )70% 碎石圖(Scree plot) 能有恰當(dāng)?shù)膶I(yè)解釋,練習(xí)1：主成分分析(變量單位不同） 20例肝病患者4項(xiàng)肝功能指標(biāo)： X1：轉(zhuǎn)氨酶(SGPT); X2：肝大指數(shù)(F); X3：硫酸鋅濁度(ZnT); X4：甲胎球蛋白(AFP) 試作主成分分析 程序： unit4princomp1.sas 數(shù)據(jù)： unit4princomp1.xls,練習(xí)2：主成分分析(變量單位相同） 我國(guó)27個(gè)少數(shù)民族體型資料 X1：頭長(zhǎng)； X2: 頭寬； X3: 額最小寬； X4: 面寬； X5: 下額角間寬； X6: 容貌面高； X7: 形態(tài)面高； X8: 鼻高； X9: 鼻寬； X10: 口裂寬； X11:身長(zhǎng)； X12: 肩寬； X13: 胸圍； X14:骨盆寬； X15:全頭高； (測(cè)量單位： mm) 程序：unit4princomp2.sas 數(shù)據(jù)： unit4princomp2.xls,練習(xí)3：主成分回歸分析 某研究所調(diào)查了13名兒童的資料 X1：性別(男1，女2） X2：年齡(月） X3：身高(厘米） X4：體重(公斤） X5：胸圍(厘米） Y： 心象面積(平方厘米） 試分析性別、年齡、身高、體重和