统计分析与统计分析与 SAS 实现上机实习题实现上机实习题 一、定量资料上机实习题一、定量资料上机实习题 要求： （1）先判断定量资料所对应的实验设计类型；先判断定量资料所对应的实验设计类型； （2）假定资料满足参数检验的前提条件，请选用相应设计的定量资料的方差分析，并用假定资料满足参数检验的前提条件，请选用相应设计的定量资料的方差分析，并用 SAS 软件实现统计计算；软件实现统计计算； （3）摘录主要计算结果并合理解释，给出统计学结论和专业结论。摘录主要计算结果并合理解释，给出统计学结论和专业结论。【练习1】取4窝不同种系未成年的大白鼠，每窝3只，随机分配到三个实验组中，分 别注射不同剂量雌激素，经过一定时间后处死大白鼠测子宫重量，资料见表1。问剂量和窝 别的各自水平下子宫重量之间的差别有无统计学意义？若剂量间差别有统计学意义，请作 两两比较。表表1 未未成成年年大大白白鼠鼠注注射射不不同同剂剂量量雌雌激激素素后后的的子子宫宫重重量量子宫重量(g) 窝别 剂量(g/100g):0.20.40.8合计1106116145367242681152253701111333144426387192合计2603584801098定量定量资资料的料的随机区随机区组设计组设计（区（区组组因素：因素：窝别窝别； ；实验实验因素：因素：剂剂量）量）【SAS 程序】： DATA PGM15G;DO A=1 TO 4; /*A为窝别*/ DO B=1 TO 3; /*B为雌激素剂量*/INPUT X ; /*X为子宫重量*/OUTPUT;END; END;CARDS;106 116 145 42 68 115 70 111 133 42 63 87 ;RUN;ods html; /*将结果输出成网页格式，SAS9.0以后版本可用*/PROC GLM DATA=PGM15G;CLASS A B;MODEL X=A B / SS3;MEANS A B; /*给出因素A、B各水平下的均值和标准差*/MEANS B / SNK; /*对因素B（即剂量）各水平下的均值进行两两比较*/RUN;ODS HTML CLOSE;【练习 2】一位工程师研究由钻头压力产生的冲力。考察了 A（钻孔速度）和 B（进 料速度） ，两因素分别取 2 与 3 水平，各水平组合下均做了两次独立重复实验，资料见表 2。假定资料满足参数检验的前提条件，且两因素对观测结果的影响地位平等，已知冲力越 小越好，试作分析，尽可能给出较为明确的统计和专业结论。 表表 2 在钻孔速度和进料速度取不同水平的条件下冲力的测定结果在钻孔速度和进料速度取不同水平的条件下冲力的测定结果钻 孔 冲力（单位） 速 度 进料速度：0.015 0.030 0.045125 2.70 2.45 2.602.78 2.49 2.72200 2.83 2.85 2.862.86 2.80 2.87两因素两因素 2 x 3 析因析因设计设计【SAS 程序】： DATA aaa;DO zs=125,200;DO repeat=1 TO 2; /*每种试验条件下有2次独立重复试验*/do js=0.015,0.030,0.045;INPUT cl ;OUTPUT;END; END; END;CARDS;2.70 2.45 2.602.78 2.49 2.722.83 2.85 2.862.86 2.80 2.87;run;ods html;PROC GLM;CLASS zs js;MODEL cl=zs js zs*js / SS3;MEANS zs*js;LSMEANS zs*js / TDIFF PDIFF; /*对 zs和js各水平组合而成的试验条件进行均数进行两两比较*/ RUN; ods html close;二、定性资料上机实习题二、定性资料上机实习题 要求：要求：（1）若题目中未给出表格，请列出标准的列联表，并对其命名；）若题目中未给出表格，请列出标准的列联表，并对其命名； （2）若题目中已列出不规范的表格，先修改，然后对其命名；）若题目中已列出不规范的表格，先修改，然后对其命名； （3）根据分析目的或自己提出分析目的、资料的前提条件选用相应的统计分析方法，并用）根据分析目的或自己提出分析目的、资料的前提条件选用相应的统计分析方法，并用 SAS 软件实现计算；软件实现计算； （4）将主要计算结果摘录出来，给出统计学和专业结论。）将主要计算结果摘录出来，给出统计学和专业结论。【练习 1】某卫生防疫站对屠宰场及肉食零售点的猪肉,检查其表层沙门氏菌带菌情况， 结果如下表。试比较屠宰场与肉食零售点猪肉表层沙门氏菌的带菌率之间差别有无统计学 意义? 表表 1 屠宰场及肉食零售点猪肉表层沙门氏菌抽检结果屠宰场及肉食零售点猪肉表层沙门氏菌抽检结果 例数采样地点阳性数值阴性数值带菌率(%)屠宰场226 7.14 零售点82127.59横断面研究横断面研究设计设计【SAS 程序】： DATA PGM19A;DO A=1 TO 2;DO B=1 TO 2;INPUT F ;OUTPUT;END;END; CARDS;2 268 21;run;ods html;PROC FREQ;WEIGHT F;TABLES A*B / CHISQ;RUN;ods html close;【练习 2】有人对某部门 22707 名雇员中，普查了 HBsAg,其中 3454 名阳性，19253 名 为阴性。从 1975 年起，追踪了 3 年，发现在阳性组有 40 名患了肝癌，阴性组仅一名患肝 癌。试选用合适的方法对资料进行全面分析。 队队列研究列研究设计设计 2 x 2 表表【SAS 程序】： DATA PGM19B;DO A=1 TO 2;DO B=1 TO 2;INPUT F ;OUTPUT;END;END;CARDS;40 34141 19252;run;ods html;PROC FREQ;WEIGHT F;TABLES A*B / CHISQ cmh;RUN;ods html close;【练习 3】APOE-4 等位基因与老年痴呆性的关联研究：以 600 名晚发及散发老年痴呆 患者和 400 名正常对照为研究对象，分析 APOE-4 等位基因与老年痴呆性的关系。表表 4 APOE-4 等位基因与老年痴呆性病例对照关联研究等位基因与老年痴呆性病例对照关联研究例数例数AGT 等位基因等位基因病例组病例组对照组对照组合计合计APOE-424060300 非非 APOE-4360340700合计合计6004001000 病例病例对对照研究照研究设计设计 2 x 2 表表【SAS 程序】： ATA PGM20;DO A=1 TO 2;DO B=1 TO 2;INPUT F ;OUTPUT;END; END;CARDS;240 60360 340;run;ods html;PROC FREQ;WEIGHT F;TABLES A*B / CHISQ cmh;RUN; ods html close; caccaaodd/)(11比数 dbca oddoddOR/21 dbddbbodd/)(22比数【练习 4】请分析下表资料。已从专业上认定培养的阳性结果就是“真阳性”，而不会 出现假阳性。 表表 3 两种培养基对同一批痰液标本同时培养的结果两种培养基对同一批痰液标本同时培养的结果例数甲培养基 培养结果乙培养基结果：+-合计+363470-0135135合计36169205 隐隐含金含金标标准配准配对设计对设计 2 x 2 表表配对设计配对设计 22 列联表资料总体率差异性检验统计量的计算公式列联表资料总体率差异性检验统计量的计算公式 若若 b+c40 时时若若 b+c40 时时【SAS 程序】： DATA PGM19F;INPUT b c;chi=(ABS(b-c)-1)*2/(b+c);p=1-PROBCHI(chi,1);chi=ROUND(chi, 0.001);IF p0.0001 THEN p=ROUND(p,0.0001);FILE PRINT;PUT #2 10 Chisq 30 P value#4 10 chi 30 p;CARDS;34 0; run; （甲培养基的阳性数大于乙培养基的阳性数）【练习 5】请分析下表资料。 表 6 两法检查室壁收缩运动的符合情况 对比法测 冠心病人数 定的结果核素法正常 减弱 异常 正 常 58 2 3 减 弱 1 42 7 异 常 8 9 17 合计 67 53 27 1)(2 2cbcb1)1(2 2cbcb双向有序且属性相同双向有序且属性相同 R x C 列列联联表（表（Kappa 检验检验） ）【SAS 程序】： data aaa;do a=1 to 3;do b=1 to 3;input f ;output;end; end; cards;58 2 3 1 42 7 8 9 17 ;run;ods html;*简单kappa检验;proc freq data=aaa; weight f;tables a*b;test kappa;run;*加权kappa检验;proc freq; weight f;tables a*b;test wtkap;run;ods html close; （两种方法的一致性检测有统计学意义） 简单 kappa 检验和加权 kappa 检验这两种方法都是用来检验两种评价方法是否具有一检验两种评价方法是否具有一 致性的方法致性的方法。其主要的区别是两种方法计算的公式不一样，更具体地说是对两个变量的打 分不一样，简单 kappa 检验主要是利用对角线上的信息，加权 kappa 检验除了利用对角线 上的数据外，还将对角线外的数据进行加权打分，将对角线外的信息也充分利用。所以在 选择方法时应根据专业知识，如果两个变量取值的界线比较明确，如“+”“+”“+”等，这 时可以选用简单的 kappa 检验，如果两个变量的取值不十分明确，人为的因素较多时，就 可选用加权 kappa 检验。【练习 6】请分析下表资料。 表 5 某地 6094 人按 2 种血型系统划分的结果