统计分析与SAS实现上机实习题一、定量资料上机实习题要求：（1）先判断定量资料所对应的实验设计类型；（2）假定资料满足参数检验的前提条件，请选用相应设计的定量资料的方差分析，并用SAS软件实现统计计算；（3）摘录主要计算结果并合理解释，给出统计学结论和专业结论。【练习1】取4窝不同种系未成年的大白鼠，每窝3只，随机分配到三个实验组中，分别注射不同剂量雌激素，经过一定时间后处死大白鼠测子宫重量，资料见表1。问剂量和窝别的各自水平下子宫重量之间的差别有无统计学意义？若剂量间差别有统计学意义，请作两两比较。表1未成年大白鼠注射不同剂量雌激素后的子宫重量窝别子宫重量(g)剂量(μg/100g):0.20.40.8合计1106116145367242681152253701111333144426387192合计2603584801098定量资料的随机区组设计（区组因素：窝别；实验因素：剂量）【SAS程序】：DATAPGM15G;DOA=1TO4;/*A为窝别*/DOB=1TO3;/*B为雌激素剂量*/INPUTX@@;/*X为子宫重量*/OUTPUT;END;END;CARDS;106116145426811570111133426387;RUN;odshtml;/*将结果输出成网页格式，SAS9.0以后版本可用*/PROCGLMDATA=PGM15G;CLASSAB;MODELX=AB/SS3;MEANSAB;/*给出因素A、B各水平下的均值和标准差*/MEANSB/SNK;/*对因素B（即剂量）各水平下的均值进行两两比较*/RUN;ODSHTMLCLOSE;【练习2】一位工程师研究由钻头压力产生的冲力。考察了A（钻孔速度）和B（进料速度），两因素分别取2与3水平，各水平组合下均做了两次独立重复实验，资料见表2。假定资料满足参数检验的前提条件，且两因素对观测结果的影响地位平等，已知冲力越小越好，试作分析，尽可能给出较为明确的统计和专业结论。表2在钻孔速度和进料速度取不同水平的条件下冲力的测定结果钻孔冲力（单位）速度进料速度：0.0150.0300.0451252.702.452.602.782.492.722002.832.852.862.862.802.87两因素2x3析因设计【SAS程序】：DATAaaa;DOzs=125,200;DOrepeat=1TO2;/*每种试验条件下有2次独立重复试验*/dojs=0.015,0.030,0.045;INPUTcl@@;OUTPUT;END;END;END;CARDS;2.702.452.602.782.492.722.832.852.862.862.802.87;run;odshtml;PROCGLM;CLASSzsjs;MODELcl=zsjszs*js/SS3;MEANSzs*js;LSMEANSzs*js/TDIFFPDIFF;/*对zs和js各水平组合而成的试验条件进行均数进行两两比较*/RUN;odshtmlclose;二、定性资料上机实习题要求：（1）若题目中未给出表格，请列出标准的列联表，并对其命名；（2）若题目中已列出不规范的表格，先修改，然后对其命名；（3）根据分析目的或自己提出分析目的、资料的前提条件选用相应的统计分析方法，并用SAS软件实现计算；（4）将主要计算结果摘录出来，给出统计学和专业结论。【练习1】某卫生防疫站对屠宰场及肉食零售点的猪肉,检查其表层沙门氏菌带菌情况，结果如下表。试比较屠宰场与肉食零售点猪肉表层沙门氏菌的带菌率之间差别有无统计学意义?表1屠宰场及肉食零售点猪肉表层沙门氏菌抽检结果采样地点例数带菌率(%)阳性数值阴性数值屠宰场2267.14零售点82127.59横断面研究设计【SAS程序】：DATAPGM19A;DOA=1TO2;DOB=1TO2;INPUTF@@;OUTPUT;END;END;CARDS;226821;run;odshtml;PROCFREQ;WEIGHTF;TABLESA*B/CHISQ;RUN;odshtmlclose;【练习2】有人对某部门22707名雇员中，普查了HBsAg,其中3454名阳性，19253名为阴性。从1975年起，追踪了3年，发现在阳性组有40名患了肝癌，阴性组仅一名患肝癌。试选用合适的方法对资料进行全面分析。队列研究设计2x2表【SAS程序】：DATAPGM19B;DOA=1TO2;DOB=1TO2;INPUTF@@;OUTPUT;END;END;CARDS;403414119252;run;odshtml;PROCFREQ;WEIGHTF;TABLESA*B/CHISQcmh;RUN;odshtmlclose;【练习3】APOE-4等位基因与老年痴呆性的关联研究：以600名晚发及散发老年痴呆患者和400名正常对照为研究对象，分析APOE-4等位基因与老年痴呆性的关系。表4APOE-4等位基因与老年痴呆性病例对照关联研究AGT等位基因例数合计病例组对照组APOE-424060300非APOE-4360340700合计6004001000病例对照研究设计2x2表【SAS程序】：ATAPGM20;DOA=1TO2;DOB=1TO2;INPUTF@@;OUTPUT;END;END;CARDS;24060360340;run;odshtml;PROCFREQ;WEIGHTF;TABLESA*...
