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

SAS数据分析应用实例及相关程序 正态性检验及 T 检验 【例 1】已知玉米单交种群105的平均穗重为 300g。喷药后，随机抽取9个果穗，其穗重分别 为： 308，305，311，298，315，300，321，294， 320g。问喷药后与喷药前的果穗平均重量之 间的差别是否具有统计学意义？ 2.配对 T 检验 【例 2】对血小板活化模型大鼠以ASA 进行实验性治疗，以血浆TXB2（ng/L ）为指标，其结 果如表 2-1，试进行统计分析。 表2-1 大鼠血小板活化模型ASA 治疗前后血浆 TXB2的变化（ ng/L ） 大 鼠号 血浆TXB2 （ng/L） 给 药前 给 药后 1 250 184 2 226 205 10 176 176 3. 秩和检验 【例 3】探讨正己烷职业接触人群生化指标特征，用气相色谱法检测受检者尿液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.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、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.19 4. 两独立正态总体的检验 【例 4】一个小麦新品种经过6代选育，从第 5代（A组）中抽出 10株，株高为： 66、65、66、 68、 62、65、63、 66、68、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、59、 61、60、59 6. 双因素无重复试验的方差分析 【例 6】某医生欲研究回心草各单体成分对试验性心肌缺血血流动力学的影响，选取健康新 西兰家兔若干只，体重（2.0 ± 0.3 ） kg，雌雄不计，将其随机分成9组：胡椒碱高剂量组 （100nmol/L ）、胡椒碱中剂量组（10nmol/L）、胡椒碱低剂量组（1nmol/L）、胡椒酸甲酯高 剂量组（ 100nmol/L）、胡椒酸甲酯中剂量组（10nmol/L）、胡椒酸甲酯低剂量组（1nmol/L）、 咖啡酸甲酯高剂量组（100nmol/L）、咖啡酸甲酯中剂量组（10nmol/L）、咖啡酸甲酯低剂量组 （1nmol/L）。所有家兔处死后，造缺血缺氧的离体心脏模型，给以各试验组相应种类及浓度的 药物进行试验，记录各组试验家兔血流动力学指标的平均值，结果见表4-2 。试分析回心草的不 同单体成分及给药剂量对冠脉流量的影响之间的差别是否有统计学意义？ 表4-2 回心草各成分对缺学缺氧后兔离体心脏血流动力学的影响 单 体 成 分 冠脉流量（x，ml） 剂 量： 高 剂量 中剂 量 低 剂量 胡椒碱8.40 8.47 7.73 胡椒酸甲酯9.93 7.53 8.03 咖啡酸甲酯10.7 3 6.47 9.47 7. 正交设计一元定量方差分析 【例 7】某研究者欲确定氧化葡萄糖的最优制备条件，采用正交试验考察PH 值、反应温度、 搅拌速度 3个试验因素的影响，因素水平见表 4-10 ，以生成物的醛基含量和反应时间为参考指标 进行综合评分，评分越高说明氧化葡萄糖制备效率越高。试验设计及试验结果见表4-11，请进 行适当的统计分析。 表4-10 因素水平表 因 素水平 反应液 pH 值(A) 搅 拌 速 度 ( r min-1) (B) 反 应 温 度 ( ) (C) 1 3.6 1000 42 2 4.5 1200 38 3 6.1 1500 25 表4-11 正交试验设计及试验结果 试验号1 (A) 2 (B) 4 (C) 醛 基含 量 (mmol/g) 反应时 间 (h) 评 分 1 1118.9 0.9 1 6.9 2 1227.5 1.3 1 3.7 3 1337.0 5.7 9 .6 4 2136.7 6.8 6 .6 5 2217.2 3.8 1 0.6 6 2326.5 5.4 7 .6 7 3126.3 10.8 1 .8 8 3236.1 12.1 0 .1 9 3316.2 7.4 5 .0 注：评分 =醛基含量× 2-反应时间 8. 协方差分析 【例 8】某研究者欲研究三种饲料对动物体重增长的影响，按照某些重要非试验因素将36 只大白鼠均分成12个配伍组， 再将每个配伍组中的3只大白鼠随机分入三个饲料组，各组进食量 与所增体重的测定结果见表4-18，试分析三种饲料对大鼠增重效果间的差别是否有统计学意义。 表4-18 三组白鼠的进食量X(g) 与所增体重 Y(g) 的试验结果 配伍 组 进食量 (g)与所增体重 (g) 饲料 1: XY 饲 料 2: X Y 饲 料 3: X Y 1 256.9 27.0 260.3 32.0 544.7 160.3 2 271.6 41.7 271.1 47.1 481.2 96.1 3 210.2 25.0 214.7 36.7 418.9 114.6 4 300.1 52.0 300.1 65.0 556.6 134.8 5 262.2 14.5 269.7 39.0 394.5 76.3 6 304.4 48.8 307.5 37.9 426.6 72.8 7 272.4 48.0 278.9 51.5 416.1 99.4 8 248.2 9.5 256.2 26.7 549.9 133.7 9 242.8 37.0 240.8 41.0 580.5 147.0 10 342.9 56.5 340.7 61.3 608.3 165.8 11 356.9 76.0 356.3 102.1 559.6 169.8 12 198.2 9.2 199.2 8.1 371.9 54.3 9. 有交互效应的方差分析 【例 9】某研究者欲研究IL-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， 请进行相应的统计分析。 10 比例检验 【例 10】为了调查某工厂产品生产合格率的情况，共对该工厂随机抽取1000件产品进行检 验，发现其中978件合格， 22件不合格。已知业内相同产品的生产合格率为98% ，问该工厂产品 合格率是否低于业内平均水平？ 11. 一致性检验 【例 11】两家评审机构(A和B)对同一组评价项目的评审结果如表7-1所示， 问两家评审机构 的评价结果是否一致。 12. 独立性检验 【例 12】某研究随机抽取了某大学四年级学生124人，调查大学英语六级通过情况，结果见 表7-3 ，问该大学男生和女生英语六级通过率有无差别？ 13. 列联表分析 【例 13】某研究为了比较三所大学大一新生中党员比例，从三所大学中随机抽取部分学生 进行比较，得到表07-8数据，试分析三所大学新生中党员的比例差别有无统计学意义？ 14简单线性回归 【例 14】土壤内 NaCl含量对植物的生长有很大影响，NaCl含量过高，将增加组织内无机 盐累积，抑制植物生长。如下表所示，分别为每千克土壤中NaCl的含量（ X），植物单位叶面 积干物重量（ Y），试进行简单线性回归分析。 15. 多重线性回归 【例 15】有研究认为血清中低密度脂蛋白增高是引起动脉硬化的一个重要原因，现测量了 30名被怀疑患有动脉硬化的就诊患者的载脂蛋白A、载脂蛋白 B、载脂蛋白 E、载脂蛋白 C、低 密度脂蛋白中的胆固醇含量，资料见表13-1 ，试分析四种载脂蛋白对低密度脂蛋白中胆固醇含 量的影响。 16. 主成分回归分析 17 主成分分析 【例 17】 不同国家和地区的女子田径记录数据列于表29-1 中。试对其进行主成分分析， 并对主成分进行解释，把不同国家和地区按它们在第一主成分上的得分排序。 18 探索性因子分析 【例 18】 50个白人男性申请警察局职位的体检数据(Gunst 和 Mason(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)（资料来源： Dallas E.Johnson. Applied Multivariate Methods for Data Analysts）。 其数据结构如表30-1所示。 19. 典型相关分析 【例 19】想要研究 1922岁汉族男性学生的身体形态学指标与功能指标之间的关系，调 查得到的具体数据见表34-1，试用典型相关分析方法对该资料进行分析。 20 变量聚类 【例 20】为研究人脑老化的严重成度，有人测定了不同年龄的60名正常男性 10项有关指 标，数据见表 35-1 。各变量的含义如下：AGE 为年龄、 TJ为图片记忆、SG 为数字广度记忆、 ?TS为图形顺序记忆、XX 为心算位数、XS 为心算时间、 CK 为规定时间内穿孔数、BJ为步距、 JJ为步行时双下肢夹角、BS 步速。试对这些指标作变量聚类分析。 21 样品聚类 【例 21】 有一项对美国39座城市空气污染情况的调查数据，在这39个城市中，对每个 城市检测了二氧化硫（SO2）、平均气温（ temperature）、20人以上的手工企业（factories）、 人口数（ population）、平均风速（windspeed）、平均降雨量（rain）、每年平均降雨天数 （rainydays）7项指标，其中二氧化硫（SO2）是结果变量，其余的六个变量为原因变量，即 对二氧化硫（ SO2）有影响的因素。数据见表36-1，试对这 39个城市的空气污染程度按原因 变量进行分类。 22. 判别分析 【例 22】 从心电图的 5个不同指标中对健康人（C1）、硬化症患者（C 2）和冠心 病患者（ C 3）的数据 ( 见 cl.dbf). 试进行判别分析。 23 时间序列分析指数平滑分析 【例 23】 某药品公司自1990年以来生产的某种抗生素的出厂数量时间序列如表20-1所 示，试用指数平滑法预测20062010年该药的出厂数量，并计算模型参数、预测值及其置信 区间。 24. 时间序列 ARIMA 模型 【例 24】为了对某区级医院体检中心1970-2005年的收入序列（连续性变量）进行预测 （数据详见表 20-2），选用 ARIMA 过程进行。通过ARIMA 模型预测的三个主要阶段，通过 识别（自相关、逆自相关及偏自相关函数，并最终根据AIC 最小的原则选择合适的模型）、 建模和预测，完成对该医院收入序列的预测结果。