SAS数据分析应用实例及相关程序

1、SAS数据分析应用实例及相关程序正态性检验及T检验【例1】已知玉米单交种群105的平均穗重为300go喷药后，随机抽取9个果穗，其穗重分别为：308,305,311,298,315,300,321,294,320g。问喷药后与喷药前的果穗平均重量之间的差别是否具有统计学意义？2 .配对T检验【例2】对血小板活化模型大鼠以ASAS行实验性治疗，以血浆TXE2（ng/L）为指标，其结果如表2-1，试进行统计分析。表2-1大鼠血小板活化模型ASA台疗前后血浆TXE2的变化（ng/L）大鼠号血浆TXE2(ng/L)给药前给药后125011842226205999101761763 .秩和检验【例3】探

2、讨正己烷职业接触人群生化指标特征，用气相色谱法检测受检者尿液2,5-己二酮浓度（mg/L）,为该人群的健康监护寻找动态观察依据。正己烷职业接触组（A组）为广州市印刷行业彩印操作位作业人员64人，其均在同一个大的车间轮班工作，工作强度相当；对照组（B组）选同厂其他车间工人53人。两组人员除接触正己烷因素不同外，生活水平、生活习惯、劳动强度、吸烟、饮酒情况基本相同。问两组间尿液中2,5-己二酮浓度（mg/L）平均含量之间的差别是否有统计学意义？数据如下所示。正己烷职业接触组：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对照组：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 .两独立正态总体的检验【例4】一个小麦新品种经过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),问株高性状是否已经达到稳定?5 .单因素K(K3)水平方差分析【例5】从津丰小麦4个品系中分别随机抽取10株，测量其株高(cm),数据如下所示，问不同品系津丰小麦的平均株高之间的差别是否具有统计学意义？品系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 .双因素无重复试验的方差分析【例6】某医生欲研究回心草各单体成分对试验性心肌缺血血流动力学的影响，选取健康新西兰家兔若干只，体重(2.00.3)kg,雌雄不计，将其随机分成9组：胡椒碱高剂量组(100nmol/L)、胡椒碱中剂量组(10nmol/L)、胡椒碱低剂量组(1nmol/L)、胡椒酸甲酯高剂量组(100nmol/L)、胡椒酸甲酯中剂量组(10nmol/L)、胡椒酸甲酯低剂量组(1nmol/L)、咖啡酸甲酯高剂量组(100nmol/L)、咖啡酸甲酯中剂量组(10nmol/L)、咖啡酸甲酯低剂量组(1nmol/L)。所有家兔处死后，造缺血缺氧的离体心脏模型，给以

7、各试验组相应种类及浓度的药物进行试验，记录各组试验家兔血流动力学指标的平均值，结果见表4-2。试分析回心草的不同单体成分及给药剂量对冠脉流量的影响之间的差别是否有统计学意义？表4-2回心草各成分对缺学缺氧后兔离体心脏血流动力学的影响单体成分旭:脉流量(x,ml)剂量：高剂量中剂量低剂量胡椒碱8.408.477.73胡椒酸甲酯9.937.538.03咖啡酸甲酯10.736.479.477 .正交设计一元定量方差分析【例7】某研究者欲确定氧化葡萄糖的最优制备条件，采用正交试验考察PH直、反应温度、搅拌速度3个试验因素的影响，因素水平见表4-10,以生成物的醛基含量和反应时间为参考指标进行综合评分，

8、评分越高说明氧化葡萄糖制备效率越高。试验设计及试验结果见表4-11，请进行适当的统计分析。表4-10因素水平表因素水平反应液pH值(A)搅拌速度(r.min-1)(B)反应温度(C)(C)13.610004224.512003836.1150025表4-11正交试验设计及试验结果试验号1(A)2(B)：4(C)I醛基含量(mmol/g)反应时间(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注：评分=醛基含量x2-反应时间8 .协方差分析【例8】某研究者欲研究三种饲料对动物体重增长的影响，按照某些重要非试验因素将36只大白鼠均分成12个配伍组，再将每个配伍组中的3只大白鼠随机分入三个饲料组，各组进食量与所增体重的测定结果见表4-18,试分析三种饲料对大鼠增重效果间的差别是否有统计学意义。表4-18三组白鼠的进食量X(g)与所增体重丫(g)的试验结果配伍组进食量(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 .有交互效应的方差分析【例9】某研究者欲研究IL

11、-11药对5.5Gy照射小鼠骨髓造血细胞周期（G0/G1期）的影响，选取45只小鼠并将其完全随机地均分成3组，每组15只，分别在5.5Gy剂量照射前给IL-11药、照射后给IL-11药和照射对照（即不给IL-11药）。每一组中的15只小鼠随机等分成3组，分别在照后6h、12h、24h三个时间点上处死，测量其骨髓造血细胞周期（G0/G1期）。试验结果见表4-6,请进行相应的统计分析。赛4-6ILT1药物对55G#照射小鼠骨般造血细胞周期期）的厨响给药时机骨髓造血细胞周期5期单位）处死时间由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照射时照65.5886.5695.2362.36S44292.4164.03S74198,8467.598350961164盘9既4795.6510 .比例检验【例10】为了调查某工厂产品生产合格率的情况，共对该工厂随机抽取1000件产品进行检验，发现其中978件合格，22件不合格。已知业内相同产品的生产合格率为98%问该工厂产品合格率是否低于业内平均水平？11 .一致性检验【例11】两家评审

13、机构（A和B）对同一组评价项目白评审结果如表7-1所示，问两家评审机构的评价结果是否一致。表77两家评审机构对同一组评价项目的评审结果机构A评审结果项目数机构B评审结果:通过木1U合计通过551267未通过81725合计63299212.独立性检验【例12】某研究随机抽取了某大学四年级学生124人，调查大学英语六级通过情况，结果见表7-3,问该大学男生和女生英语六级通过率有无差别？表7-3某大学四年级男、女生大学英语六级通过情况性别人数英语六级通过情况：通过未通过人；111r男413273女43851合计844012413.列联表分析【例13某研究为了比较三所大学大一新生中党员比例，从三所大学

14、中随机抽取部分学生进行比较，得到表07-8数据，试分析三所大学新生中党员的比例差别有无统计学意义？表7-8三所大学大一新生的政治面貌调查结果大学人数政治面貌:兄贝非党员合计A5373126B38132170C168298合计10728739114.简单线性回归【例14】土壤内NaCl含量对植物的生长有很大影响，NaCl含量过高，将增加组织内无机盐累积，抑制植物生长。如下表所示，分别为每千克土壤中NaCl的含量（X）,植物单位叶面积干物重量（Y）,试进行简单线性回归分析。表107不同土NaCl含导对植物单位叶面积干糊打捕的帮响土塔N】含量X便0也E2,432LQ48IjfY/mgdin28090

15、9511513011513515.多重线性回归【例15】有研究认为血清中低密度脂蛋白增高是引起动脉硬化的一个重要原因，现测量了30名被怀疑患有动脉硬化的就诊患者的载脂蛋白AI、载脂蛋白B、载脂蛋白E、载脂蛋白C低密度脂蛋白中的胆固醇含量，资料见表13-1，试分析四种载脂蛋白对低密度脂蛋白中胆固醇含量的影响。表13-13口，名患者载脂蛋白和低密度脂蛋白中胆固酶含蛋的测品结果编号我脂雷白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】某种水泥在凝固时放5的热曷F(卡/克)与该物质中的下列4种化学成分有关：1：3CaO-Al2O3的成分(%)；X2：3CaO-SiO二的成分(%)：苍：4CaOAbOyFjO&的成分()；X4：2CaOSiO3的成分(%)二作F对4个自变整X1、A、一七、X的

18、线性回H方程目数据结构如表141所示口表i4-i某种水泥在凝固时放出的热量y与其4种成分的记录编号X莅储Y72666070.5129155274.331156320,04.34113L84787.6575二3395.9t11559109.27311170102,19131224472,5954,8n个93.1102147426115.911140233483,812116c912113.313103812.105.417主成分分析【例17】不同国家和地区的女子田径记录数据列于表29-1中。试对其进行主成分分析,并对主成分进行解释，把不同国家和地区按它们在第一主成分上的得分排序。泰29-1不同国

19、宗和地区女子径赛运动记聚国加和地嘲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比利时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中国1L9524.4154.972.0S4339JI168,46若伦比IkU6024.0053.262.114.35946165.42摩克群岛12.9027,.1060.402.304W4JI.10233.22哥斯达黎加1L9624.G058.2?TT1&r*iJi4.S810.43171,SO捷克11,0321.9747,99l加4.1+九92158,B5丹麦1142235253.602.0J4.IS8.71151.75多米尼加共和国11.

21、7924.056,052244,749,8920,SS芬兰1L1J22J950,142.034,108,921M,23法国11,1522.5951.732.004,148.98155,27东飞10.8121.714S,161.933.96B.75157.68飞11.0122.3949.751.954.038.59148.53惦国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个个-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国11.9624.4955.70265川虾12.2525.7S51.201.974

23、.25g.ms179.17卢森耍12.0324.96561。2.074389.64174.68马来西亚12.232421国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巴布亚新几内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台湾11.2222.6252.502.104.389,6317

25、7.87泰国11.75244655.802.204.7210,28168,45印共11.9824.4456.452.154.379.38201,08美国10J921.8350,621.963.95S.50142J2前苏联11.0622.1949.J91.893.87S.45J51.22西胪糜W.12.7425,8558,732335.8113.04306.00资料来海：L4.4FATFSTrackandFieldStatisticsHandbookforthe-1981LosAngelesOlympics18探索性因子分析【例18】50个白人男性申请警察局职位的体检数据（Gunst和Mason

26、（1980）,指标括：1、身高（单位：厘米）（height）;2、体重（单位：千克）（weight）;3、肩宽（单位：米）（shldr）;4、骨盆宽（单位：厘米）（pelvic）;5、最小胸围（单位：厘米）（chest）;6、腿皮褶厚度（单位：毫米）（thigh）;7、静息脉率（pulse）;8、引体向上次数（chnup）;9、大肺活量（单位：公升）（breath）;10、踏板跑步休息5分钟后脉率（recvr）;11、最大踏板度（speed）;12、体脂（fat）（资料来源：DallasE.Johnson.AppliedMultivariateMethodsforDataAnalysts）。其

27、数据结构如表30-1所示。题3Al5。个白人男恺申请警殍局职位的体检数据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.典型相关分析【例19】想要研究1922岁汉族男性学生的身体形态学指标与功能指标之间的关系，调查得到的具体数据见表34-1，试用典型相关分析方法对该资料进行分析。表34-1汉族男性学生的身体形态学指标与功能指标数据形态学指标功能指标褊号身痛(cm)坐高(cm)体1恤)陷W(cm)府宽(cm)骨筵筐(cm)脉搏tiw”而。快缩压(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