中科院研究生院统计分析与SAS实现第1次上机实习题
一、定量资料上机实习题
要求：
（1）先判断定量资料所对应的实验设计类型；
（2）假定资料满足参数检验的前提条件，请选用相应设计的定量资料的方差分析，并用SAS软件实现统计计算；
（3）摘录主要计算结果并合理解释，给出统计学结论和专业结论。
【练习1】取4窝不同种系未成年的大白鼠，每窝3只，随机分配到三个实验组中，分别注射不同剂量雌激素，经过一定时间后处死大白鼠测子宫重量，资料见表1。问剂量和窝别的各自水平下子宫重量之间的差别有无统计学意义？若剂量间差别有统计学意义，请作两两比较。
表1 未成年大白鼠注射不同剂量雌激素后的子宫重量
子宫重量(g)
窝别
剂量(μg/100g):0.2 0.4 0.8 合计
1 106 116 145 367
2 42 68 115 225
3 70 111 133 314
4 42 63 87 192
合计260 358 480 1098
定量资料的随机区组设计（区组因素：窝别；实验因素：剂量）【SAS程序】：
DATA PGM15G;
DO A=1TO4; /*A为窝别*/
DO B=1TO3; /*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.045
125 2.70 2.45 2.60
2.78 2.49 2.72
200 2.83 2.85 2.86
2.86 2.80 2.87
两因素2 x 3析因设计
【SAS程序】：
DATA aaa;
DO zs=125,200;
DO repeat=1TO2; /*每种试验条件下有2次独立重复试验*/
do js=0.015,0.030,0.045;
INPUT cl @@;
OUTPUT;
END;
END;
END;
CARDS;
2.70 2.45 2.60
2.78 2.49 2.72
2.83 2.85 2.86
2.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 屠宰场及肉食零售点猪肉表层沙门氏菌抽检结果
采样地点带菌率(%)
阳性数值阴性数值
屠宰场 2 26 7.14
零售点8 21 27.59
横断面研究设计
【SAS程序】：
DATA PGM19A;
DO A=1TO2;
DO B=1TO2;
INPUT F @@;
OUTPUT;
END;
END;
CARDS;
2 26
8 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=1TO2;
DO B=1TO2;
INPUT F @@;
OUTPUT;
END;
END;
CARDS;
40 3414
1 19252
run;
ods html;
PROC FREQ;
WEIGHT F;
TABLES A*B / CHISQ cmh;
RUN;
