SAS數(shù)據(jù)分析應(yīng)用實(shí)例及相關(guān)程序

1、SAS數(shù)據(jù)分析應(yīng)用實(shí)例及相關(guān)程序正態(tài)性檢驗(yàn)及T檢驗(yàn)【例1】已知玉米單交種群105的平均穗重為300go噴藥后，隨機(jī)抽取9個(gè)果穗，其穗重分別為：308,305,311,298,315,300,321,294,320g。問噴藥后與噴藥前的果穗平均重量之間的差別是否具有統(tǒng)計(jì)學(xué)意義？2 .配對(duì)T檢驗(yàn)【例2】對(duì)血小板活化模型大鼠以ASAS行實(shí)驗(yàn)性治療，以血漿TXE2（ng/L）為指標(biāo)，其結(jié)果如表2-1，試進(jìn)行統(tǒng)計(jì)分析。表2-1大鼠血小板活化模型ASA臺(tái)療前后血漿TXE2的變化（ng/L）大鼠號(hào)血漿TXE2(ng/L)給藥前給藥后125011842226205999101761763 .秩和檢驗(yàn)【例3】探

2、討正己烷職業(yè)接觸人群生化指標(biāo)特征，用氣相色譜法檢測(cè)受檢者尿液2,5-己二酮濃度（mg/L）,為該人群的健康監(jiān)護(hù)尋找動(dòng)態(tài)觀察依據(jù)。正己烷職業(yè)接觸組（A組）為廣州市印刷行業(yè)彩印操作位作業(yè)人員64人，其均在同一個(gè)大的車間輪班工作，工作強(qiáng)度相當(dāng)；對(duì)照組（B組）選同廠其他車間工人53人。兩組人員除接觸正己烷因素不同外，生活水平、生活習(xí)慣、勞動(dòng)強(qiáng)度、吸煙、飲酒情況基本相同。問兩組間尿液中2,5-己二酮濃度（mg/L）平均含量之間的差別是否有統(tǒng)計(jì)學(xué)意義？數(shù)據(jù)如下所示。正己烷職業(yè)接觸組：2.89、1.85、2.27、2.07、1.62、1.77、2.53、2.02、2.07、2.07、1.93、3.01、1.

3、93、1.88、1.55、1.36、2.23、2.55、1.73、2.65、1.95、2.45、1.41、2.46、2.38、1.55、2.16、2.01、1.37、2.16、2.00、2.07、2.57、2.11、2.37、1.39、2.18、2.33、1.46、2.16、2.03、2.96、2.21、2.00、2.58、2.19、2.41、1.68、1.93、1.93、1.93、1.87、1.74、2.70、1.83、2.17、2.52、2.09、2.28、1.65、1.19、1.58、0.89、1.65對(duì)照組：0.27、0.36、0.26、0.16、0.49、0.58、0.16、0.45

4、、0.22、0.25、0.66、0.05、0.31、0.12、0.51、0.30、0.37、0.14、0.28、0.33、0.36、0.51、0.37、0.36、0.47、0.34、0.72、0.39、0.55、0.17、0.27、0.33、0.30、0.26、0.50、0.17、0.22、0.18、0.17、0.62、0.27、0.26、0.34、0.17、0.61、0.42、0.39、0.28、0.36、0.43、0.24、0.15、0.194 .兩獨(dú)立正態(tài)總體的檢驗(yàn)【例4】一個(gè)小麥新品種經(jīng)過6代選育，從第5代（A組）中抽出10株，株高為：66、65、66、68、62、65、63、66、6

5、8、62（cm）,又從第6代（B組）中抽出10株，株高為：64、61、57、65、65、63、62、63、64、60(cm),問株高性狀是否已經(jīng)達(dá)到穩(wěn)定?5 .單因素K(K3)水平方差分析【例5】從津豐小麥4個(gè)品系中分別隨機(jī)抽取10株，測(cè)量其株高(cm),數(shù)據(jù)如下所示，問不同品系津豐小麥的平均株高之間的差別是否具有統(tǒng)計(jì)學(xué)意義？品系0-3-1:63、65、64、65、61、68、65、65、63、64品系0-3-2:56、54、58、57、57、57、60、59、63、62品系0-3-3:61、61、67、62、62、60、67、66、63、65品系0-3-4:53、58、60、56、55、60

6、、59、61、60、596 .雙因素?zé)o重復(fù)試驗(yàn)的方差分析【例6】某醫(yī)生欲研究回心草各單體成分對(duì)試驗(yàn)性心肌缺血血流動(dòng)力學(xué)的影響，選取健康新西蘭家兔若干只，體重(2.00.3)kg,雌雄不計(jì)，將其隨機(jī)分成9組：胡椒堿高劑量組(100nmol/L)、胡椒堿中劑量組(10nmol/L)、胡椒堿低劑量組(1nmol/L)、胡椒酸甲酯高劑量組(100nmol/L)、胡椒酸甲酯中劑量組(10nmol/L)、胡椒酸甲酯低劑量組(1nmol/L)、咖啡酸甲酯高劑量組(100nmol/L)、咖啡酸甲酯中劑量組(10nmol/L)、咖啡酸甲酯低劑量組(1nmol/L)。所有家兔處死后，造缺血缺氧的離體心臟模型，給以

7、各試驗(yàn)組相應(yīng)種類及濃度的藥物進(jìn)行試驗(yàn)，記錄各組試驗(yàn)家兔血流動(dòng)力學(xué)指標(biāo)的平均值，結(jié)果見表4-2。試分析回心草的不同單體成分及給藥劑量對(duì)冠脈流量的影響之間的差別是否有統(tǒng)計(jì)學(xué)意義？表4-2回心草各成分對(duì)缺學(xué)缺氧后兔離體心臟血流動(dòng)力學(xué)的影響單體成分旭:脈流量(x,ml)劑量：高劑量中劑量低劑量胡椒堿8.408.477.73胡椒酸甲酯9.937.538.03咖啡酸甲酯10.736.479.477 .正交設(shè)計(jì)一元定量方差分析【例7】某研究者欲確定氧化葡萄糖的最優(yōu)制備條件，采用正交試驗(yàn)考察PH直、反應(yīng)溫度、攪拌速度3個(gè)試驗(yàn)因素的影響，因素水平見表4-10,以生成物的醛基含量和反應(yīng)時(shí)間為參考指標(biāo)進(jìn)行綜合評(píng)分，

8、評(píng)分越高說明氧化葡萄糖制備效率越高。試驗(yàn)設(shè)計(jì)及試驗(yàn)結(jié)果見表4-11，請(qǐng)進(jìn)行適當(dāng)?shù)慕y(tǒng)計(jì)分析。表4-10因素水平表因素水平反應(yīng)液pH值(A)攪拌速度(r.min-1)(B)反應(yīng)溫度(C)(C)13.610004224.512003836.1150025表4-11正交試驗(yàn)設(shè)計(jì)及試驗(yàn)結(jié)果試驗(yàn)號(hào)1(A)2(B)：4(C)I醛基含量(mmol/g)反應(yīng)時(shí)間(h)I分11118.90.916.9212：:?7.51.313.7313i37.05.79.64r2-36.76.86.65r2z：17.23.810.66r23iN6.55.47.67231N6.310.81.8823z：36.112.10.192

9、33i16.27.45.0注：評(píng)分=醛基含量x2-反應(yīng)時(shí)間8 .協(xié)方差分析【例8】某研究者欲研究三種飼料對(duì)動(dòng)物體重增長(zhǎng)的影響，按照某些重要非試驗(yàn)因素將36只大白鼠均分成12個(gè)配伍組，再將每個(gè)配伍組中的3只大白鼠隨機(jī)分入三個(gè)飼料組，各組進(jìn)食量與所增體重的測(cè)定結(jié)果見表4-18,試分析三種飼料對(duì)大鼠增重效果間的差別是否有統(tǒng)計(jì)學(xué)意義。表4-18三組白鼠的進(jìn)食量X(g)與所增體重丫(g)的試驗(yàn)結(jié)果配伍組進(jìn)食量(g)與所增體重(g)飼料1：XY飼料2:XY飼料3:XY1256.927.0260.332.0544.7160.32271.641.7271.147.1481.296.13210.225.0214

10、.736.7418.9114.64300.152.0300.165.0556.6134.85262.214.5269.739.0394.576.36304.448.8307.537.9426.672.87272.448.0278.951.5416.199.48248.29.5256.226.7549.9133.79242.837.0240.841.0580.5147.010342.956.5340.761.3608.3165.811356.976.0356.3102.1559.6169.812198.29.2199.28.1371.954.39 .有交互效應(yīng)的方差分析【例9】某研究者欲研究IL

11、-11藥對(duì)5.5Gy照射小鼠骨髓造血細(xì)胞周期（G0/G1期）的影響，選取45只小鼠并將其完全隨機(jī)地均分成3組，每組15只，分別在5.5Gy劑量照射前給IL-11藥、照射后給IL-11藥和照射對(duì)照（即不給IL-11藥）。每一組中的15只小鼠隨機(jī)等分成3組，分別在照后6h、12h、24h三個(gè)時(shí)間點(diǎn)上處死，測(cè)量其骨髓造血細(xì)胞周期（G0/G1期）。試驗(yàn)結(jié)果見表4-6,請(qǐng)進(jìn)行相應(yīng)的統(tǒng)計(jì)分析。賽4-6ILT1藥物對(duì)55G#照射小鼠骨般造血細(xì)胞周期期）的廚響給藥時(shí)機(jī)骨髓造血細(xì)胞周期5期單位）處死時(shí)間由1:61224照射前給弱65.8086.8395.8866一7384.6993.6766.59469皖676

12、2,82S7.9294.9363.69861090.84照射后給藥59.30S0.6995.53573486.70910560.164.7593.2660,46833093.67照射時(shí)照65.5886.5695.2362.36S44292.4164.03S74198,8467.598350961164盤9既4795.6510 .比例檢驗(yàn)【例10】為了調(diào)查某工廠產(chǎn)品生產(chǎn)合格率的情況，共對(duì)該工廠隨機(jī)抽取1000件產(chǎn)品進(jìn)行檢驗(yàn)，發(fā)現(xiàn)其中978件合格，22件不合格。已知業(yè)內(nèi)相同產(chǎn)品的生產(chǎn)合格率為98%問該工廠產(chǎn)品合格率是否低于業(yè)內(nèi)平均水平？11 .一致性檢驗(yàn)【例11】?jī)杉以u(píng)審

13、機(jī)構(gòu)（A和B）對(duì)同一組評(píng)價(jià)項(xiàng)目白評(píng)審結(jié)果如表7-1所示，問兩家評(píng)審機(jī)構(gòu)的評(píng)價(jià)結(jié)果是否一致。表77兩家評(píng)審機(jī)構(gòu)對(duì)同一組評(píng)價(jià)項(xiàng)目的評(píng)審結(jié)果機(jī)構(gòu)A評(píng)審結(jié)果項(xiàng)目數(shù)機(jī)構(gòu)B評(píng)審結(jié)果:通過木1U合計(jì)通過551267未通過81725合計(jì)63299212.獨(dú)立性檢驗(yàn)【例12】某研究隨機(jī)抽取了某大學(xué)四年級(jí)學(xué)生124人，調(diào)查大學(xué)英語(yǔ)六級(jí)通過情況，結(jié)果見表7-3,問該大學(xué)男生和女生英語(yǔ)六級(jí)通過率有無(wú)差別？表7-3某大學(xué)四年級(jí)男、女生大學(xué)英語(yǔ)六級(jí)通過情況性別人數(shù)英語(yǔ)六級(jí)通過情況：通過未通過人；111r男413273女43851合計(jì)844012413.列聯(lián)表分析【例13某研究為了比較三所大學(xué)大一新生中黨員比例，從三所大學(xué)

14、中隨機(jī)抽取部分學(xué)生進(jìn)行比較，得到表07-8數(shù)據(jù)，試分析三所大學(xué)新生中黨員的比例差別有無(wú)統(tǒng)計(jì)學(xué)意義？表7-8三所大學(xué)大一新生的政治面貌調(diào)查結(jié)果大學(xué)人數(shù)政治面貌:兄貝非黨員合計(jì)A5373126B38132170C168298合計(jì)10728739114.簡(jiǎn)單線性回歸【例14】土壤內(nèi)NaCl含量對(duì)植物的生長(zhǎng)有很大影響，NaCl含量過高，將增加組織內(nèi)無(wú)機(jī)鹽累積，抑制植物生長(zhǎng)。如下表所示，分別為每千克土壤中NaCl的含量（X）,植物單位葉面積干物重量（Y）,試進(jìn)行簡(jiǎn)單線性回歸分析。表107不同土NaCl含導(dǎo)對(duì)植物單位葉面積干糊打捕的幫響土塔N】含量X便0也E2,432LQ48IjfY/mgdin28090

15、9511513011513515.多重線性回歸【例15】有研究認(rèn)為血清中低密度脂蛋白增高是引起動(dòng)脈硬化的一個(gè)重要原因，現(xiàn)測(cè)量了30名被懷疑患有動(dòng)脈硬化的就診患者的載脂蛋白AI、載脂蛋白B、載脂蛋白E、載脂蛋白C低密度脂蛋白中的膽固醇含量，資料見表13-1，試分析四種載脂蛋白對(duì)低密度脂蛋白中膽固醇含量的影響。表13-13口，名患者載脂蛋白和低密度脂蛋白中膽固酶含蛋的測(cè)品結(jié)果編號(hào)我脂雷白A1(睡dl)蛇脂擊白B(mg41)我脂蛋白E(mg/dQ載脂蛋白c(mg/dl)低齡度脂蛋白tnigdl)11731067,014,713r21391326.417316231981126.916.7134.41

16、18US15.7188513994亂613.6US617516012420.3215-13115411,221,51-1S1581419.729.614Sg158137T418.219710132151T517.2113111621106.015.914512144】310.142.581131621377.220.71E514169129.516157151291386310.11971616614SILS33.41561718511B6.017.5156181551216.120.4154191751114.127,2144201361109.426.090211531338.516.92

17、15一110149910.8118241121238.0166127251471108,518,4137262041226J2L0126271311026.613一41302817012&424.719.0LES3013213113.829J12216.主成分回歸分析【例141】某種水泥在凝固時(shí)放5的熱曷F(卡/克)與該物質(zhì)中的下列4種化學(xué)成分有關(guān)：1：3CaO-Al2O3的成分(%)；X2：3CaO-SiO二的成分(%)：蒼：4CaOAbOyFjO&的成分()；X4：2CaOSiO3的成分(%)二作F對(duì)4個(gè)自變整X1、A、一七、X的

18、線性回H方程目數(shù)據(jù)結(jié)構(gòu)如表141所示口表i4-i某種水泥在凝固時(shí)放出的熱量y與其4種成分的記錄編號(hào)X蒞儲(chǔ)Y72666070.5129155274.331156320,04.34113L84787.6575二3395.9t11559109.27311170102,19131224472,5954,8n個(gè)93.1102147426115.911140233483,812116c912113.313103812.105.417主成分分析【例17】不同國(guó)家和地區(qū)的女子田徑記錄數(shù)據(jù)列于表29-1中。試對(duì)其進(jìn)行主成分分析,并對(duì)主成分進(jìn)行解釋，把不同國(guó)家和地區(qū)按它們?cè)诘谝恢鞒煞稚系牡梅峙判?。?9-1不同國(guó)

19、宗和地區(qū)女子徑賽運(yùn)動(dòng)記聚國(guó)加和地嘲100M(s)200M(s)400H(a)SOOMtm)150QM(m)3000M(ro)4拉松(m)阿根廷11.6122.9454.502.154.43979178.52澳大利亞11.2022.3551.081.984.139.0815237奧地利11432J,09?0.621.994”Tr9.34159.3比利時(shí)114123.0452.0C2.004.14157.85百慕大群島11.4623.0553J02.164.589.81169.98已兩11.3123.12.S02.104.49917165,75始匐12.1424.4755.002IS4.459.5

20、1131.02加拿大11.00222550.062.004.068.S1149.45智利12.0024.5254.902.054.239.37171.38中國(guó)1L9524.4154.972.0S4339JI168,46若倫比IkU6024.0053.262.114.35946165.42摩克群島12.9027,.1060.402.304W4JI.10233.22哥斯達(dá)黎加1L9624.G058.2?TT1&r*iJi4.S810.43171,SO捷克11,0321.9747,99l加4.1+九92158,B5丹麥1142235253.602.0J4.IS8.71151.75多米尼加共和國(guó)11.

21、7924.056,052244,749,8920,SS芬蘭1L1J22J950,142.034,108,921M,23法國(guó)11,1522.5951.732.004,148.98155,27東飛10.8121.714S,161.933.96B.75157.68飛11.0122.3949.751.954.038.59148.53惦國(guó)11.0022.1350.461.9E4.038.62149.72希臘11.7924.0854.932.0-4359.S7182.20危地馬拉11.8424.5456.092.2S4.S610.5421?.0S匈牙利11.4523.0651.502.014.148.98

22、156,3卬11.9524.2853.602.104J299818S,0S印度尼西亞11.8524.2455.347個(gè)個(gè)-a_4.6110.02201.28愛爾蘭11.4323.5053.242.054.118.S9149.38以色列1L4523.5754.902.104.259J7160.48於太利1L2923.0052.011.963取8.63151.82R本1L7J24.005173工094J59.20150.50肯尼亞1L7323.8852.02.004.159201S1.054國(guó)11.9624.4955.70265川蝦12.2525.7S51.201.974

23、.25g.ms179.17盧森耍12.0324.96561。2.074389.64174.68馬來(lái)西亞12.232421國(guó)5.092.19L6910.64182.1毛里求斯J1J625,0858,104,7910.90261,13墨西哥11.8)23.6253.762.044.259.59158.53荷蘭11.2522.E152.3S1.994.069.01152.43新西蘭11.5523,1351.602.024J88.76145.4B挪喊11.58233151122.034018.53145.4S巴布亞新幾內(nèi)MJ22525.0S6.962.244.8410.69233.00菲律賓11762

24、3.5454.6021946010.16200J7波蘭11322249.291953998.97160.82嘴牙L1.8L24.2254.302.094.163,S4151.20羅馬尼亞11.4423.4651.201.923.96S.53165.45新加坡12.3025.0055.OS2,124.529.941S2.77西班牙11.8023.9S53.592.054.149,02162,60瑞典11.J622.S251.792.024.128.84154,48瑞1：11.4523.3153.J12.024.078.77153,42臺(tái)灣11.2222.6252.502.104.389,6317

25、7.87泰國(guó)11.75244655.802.204.7210,28168,45印共11.9824.4456.452.154.379.38201,08美國(guó)10J921.8350,621.963.95S.50142J2前蘇聯(lián)11.0622.1949.J91.893.87S.45J51.22西臚糜W.12.7425,8558,732335.8113.04306.00資料來(lái)海：L4.4FATFSTrackandFieldStatisticsHandbookforthe-1981LosAngelesOlympics18探索性因子分析【例18】50個(gè)白人男性申請(qǐng)警察局職位的體檢數(shù)據(jù)（Gunst和Mason

26、（1980）,指標(biāo)括：1、身高（單位：厘米）（height）;2、體重（單位：千克）（weight）;3、肩寬（單位：米）（shldr）;4、骨盆寬（單位：厘米）（pelvic）;5、最小胸圍（單位：厘米）（chest）;6、腿皮褶厚度（單位：毫米）（thigh）;7、靜息脈率（pulse）;8、引體向上次數(shù)（chnup）;9、大肺活量（單位：公升）（breath）;10、踏板跑步休息5分鐘后脈率（recvr）;11、最大踏板度（speed）;12、體脂（fat）（資料來(lái)源：DallasE.Johnson.AppliedMultivariateMethodsforDataAnalysts）。其

27、數(shù)據(jù)結(jié)構(gòu)如表30-1所示。題3Al5。個(gè)白人男愷申請(qǐng)警殍局職位的體檢數(shù)據(jù)idheifbtcigllTvhldrpclviicchebichigjipuheclmupbreath小vpccdfai1179.674.2041.727.iS2.419.06415S1085.5011.912175.662.0437.529.1S4.15.58S20166JOB5.503.133166.272.96泗426.S稿122.010071671165.5016.S94173.3S5.9241.227.697.619.56442201205.5019.59sI84,S65承39826.18S.214.5SOg2

28、101205.507.746199A102.2643330.1101.222.0604毋916.0030.427IP1.584.Q442.825.49LQISO64l2721106.0013.70ISO.241.690.45.5741J1931175.503一049IS3,S95.M42.330.1100.213.5SO41991055.5020.26ID163.154.2537.224.280.57.04131571136.003.0411169.675.9239.427.292.016.56561801105.0012.8312171.67L7039127.086.225.56801931

29、055.5015.9?13180.080.(5540.82S.387.4r.5742181095.001LX614174470.4039825.983516.510461901295.009.9315ISIS91.4040.629.595.1,12.09212061395.003工631665,7439.726.486.0B.Oso61811205.506例17|73079.284L226.996.111.5r：61S41115.5011.5718179.892.0640.029.8100.915.06002059250024.2119176887.964L228.4100.820.57612

30、281474.0022.3920179.377.6641.4S1.690.19.558151989芯5.506.2921193.598.444L629.295.721.05402541105.5023.6322L78,S65,4239.327.133.016.58872061215.5010.5323179,6910443.830.JIQ0.822.010031941245.0020.6：24172.681.7240.927391522.07442011135.501839I7L569.6040427.S57.715J70101751105.50ll.LJ2616S.963.6639826.7

31、83.96.06871791195.505.1627183487.2443.22E.393.711.0ss72451153.509.60163.654.8637.526.684.015.56461461155.0011.9329184.384.6S40.329.09i.28.5642213J095.508.55SOIS1.075.7E42.S29.-90.38.?36111S13096.004.9431ISO.275,8441.428.88.113.?759192J445.5010.6232184.170.4S42.026.9813L4.08452311235.508.46331袍9gg42.

32、528.795L6.05412圖6ns6,0013.4734I70.Q70卵39.727._93415.05041780$5.5012.8119.典型相關(guān)分析【例19】想要研究1922歲漢族男性學(xué)生的身體形態(tài)學(xué)指標(biāo)與功能指標(biāo)之間的關(guān)系，調(diào)查得到的具體數(shù)據(jù)見表34-1，試用典型相關(guān)分析方法對(duì)該資料進(jìn)行分析。表34-1漢族男性學(xué)生的身體形態(tài)學(xué)指標(biāo)與功能指標(biāo)數(shù)據(jù)形態(tài)學(xué)指標(biāo)功能指標(biāo)褊號(hào)身痛(cm)坐高(cm)體1恤)陷W(cm)府寬(cm)骨筵筐(cm)脈搏tiw”而?？炜s壓(umiHg)舒張壓(Proto-)(nmiHg)舒張壓(Late)(nmiHg)肺活最(ml)1173.2893.6260.1

33、085.7238.?727,51753117.474.661.84508172.0&92.SJ60.3S7.39iS6227.8276.7120177.166244b9317J.4692.7859,7485.5938.8327.465.8121.875.265.443984170.0S92.2558.04859238.332-.2976.1115.173.&613406S51706192.3659.678746383827.1472.9J9477.567.143396171.6992.8559.448745381927.1072.7116274.559.31393717J.4692,9358-

34、087.0638.58273676.51廠,975.068.3438917160932g59.75S8.03i8.6&27.2275.2115174.1632450691716092.266。上087,38.7?維6374.7317,478.368J439510171.1692.6258.7287.1138.192.18732113.272.55k04462IJ170.04921756.9588.083S2427.C5.8115.976.965.641S112170.27919456.0084.52371626.8176.4113,674.J65.642321370.6192.?057.34S

35、5.6138.522-J676.4116-74J61.243051417U9924458.92537S8.8326.4774.91B.174.061.24276J5171.8392.7956.S585.3538.5827.0378.7112.472.961.44067161713692.5358.39810938,23：-,0473,9118.473.06L34421J17J.2492.6157石983.9839,042.077116.374.25LB4284181704992.0357.?6S7.1B3S.542.5772J114871.05542S9191694391.6757.22S3.8738.4126.6076.7117.5