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

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 K 3 水平方差分析 例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 3 4 C 醛基含量 mmol g 反应时 间 h 评 分 1 1 1 1 1 8 9 0 9 1 6 9 2 1 2 2 2 7 5 1 3 1 3 7 3 1 3 3 3 7 0 5 7 9 6 4 2 1 2 3 6 7 6 8 6 6 5 2 2 3 1 7 2 3 8 1 0 6 6 2 3 1 2 6 5 5 4 7 6 7 3 1 3 2 6 3 10 8 1 8 8 3 2 1 3 6 1 12 1 0 1 9 3 3 2 1 6 2 7 4 5 0 注 评分 醛基含量 2 反应时间 8 协方差分析 例8 某研究者欲研究三种饲料对动物体重增长的影响 按照某些重要非试验因素将36 只大白鼠均分成12个配伍组 再将每个配伍组中的3只大白鼠随机分入三个饲料组 各组进食 量与所增体重的测定结果见表4 18 试分析三种饲料对大鼠增重效果间的差别是否有统计学意 义 表4 18 三组白鼠的进食量X g 与所增体重Y g 的试验结果 进食量 g 与所增体重 g 配伍 组 饲料1 X Y 饲料 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 想要研究19 22岁汉族男性学生的身体形态学指标与功能指标之间的关系 调 查得到的具体数据见表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个不同指标中对健康人 C 1 硬化症患者 C 2 和冠心 病患者 C 3 的数据 见cl dbf 试进行判别分析 23 时间