SAS练习题及程序答案复习过程

1、Good is good, but better carries it.精益求精，善益求善。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各水平下的均值和标

2、准差*/MEANSB/SNK;/*对因素B（即剂量）各水平下的均值进行两两比较*/RUN;ODSHTMLCLOSE;2*3析因设计两因素完全随机统计方法2*3析因设计tiff=f的开方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;PROCGLM;CLASSzsjs;MODELcl=zsjszs*js/SS3;MEAN

3、Szs*js;LSMEANSzs*js/TDIFFPDIFF;/*对zs和js各水平组合而成的试验条件进行均数进行两两比较*/RUN;ODSHTMLCLOSE;练习一：2*2横断面研究列链表方法：卡方矫正卡方FISHERDATAPGM19A;DOA=1TO2;DOB=1TO2;INPUTF;OUTPUT;END;END;CARDS;226821;run;PROCFREQ;WEIGHTF;TABLESA*B/CHISQ;RUN;样本大小=57练习二：对裂列连表结果变量换和不换三部曲1横断面研究P0.05RDATAPGM19B;DOA=1TO2;DOB=1TO2;INPUTF;OUTPUT;END

4、;END;CARDS;403414119252;run;odshtml;PROCFREQ;WEIGHTF;TABLESA*B/CHISQcmh;RUN;odshtmlclose;样本大小=57练习三：病例对照2*2病例组中有何没有那个基因是正常的3.8倍，则有可能导致痴呆要做前瞻性研究用对裂DATAPGM20;DOA=1TO2;DOB=1TO2;INPUTF;OUTPUT;END;END;CARDS;24060360340;run;odshtml;PROCFREQ;WEIGHTF;TABLESA*B/CHISQcmh;RUN;odshtmlclose;总样本大小=1000练习四：配对设计隐含金

5、标准2*2MC卡方检验34和0在总体上（B+C40用矫正卡方）是否相等则可得甲培养基优于乙培养基一般都用矫正因卡方为近似计算DATAPGM19F;INPUTbc;chi=(ABS(b-c)-1)*2/(b+c);p=1-PROBCHI(chi,1);求概率1减掉从左侧积分到卡方的值chi=ROUND(chi,0.001);IFp0.0001THENp=ROUND(p,0.0001);FILEPRINT;PUT（打印在输出床口）#210Chisq30Pvalue（#表示行）#410chi30p;CARDS;340;run;odshtmlclose;练习五：双向有序R*C列连表用KPAdataaa

6、a;doa=1to3;dob=1to3;inputf;output;end;end;cards;582314278917;run;odshtml;*简单kappa检验;procfreqdata=aaa;weightf;（频数）tablesa*b;testkappa;run;*加权kappa检验;procfreq;weightf;tablesa*b;testwtkap;run;odshtmlclose;SAS系统FREQ过程频数百分比行百分比列百分比a*b表ab合计12315839.4692.0686.5721.363.173.7732.044.7611.116342.86210.682.001

7、.494228.5784.0079.2574.7614.0025.935034.01385.4423.5311.9496.1226.4716.981711.5650.0062.963423.13合计6745.585336.052718.37147100.00a*b表的统计量对称性检验统计量(S)2.8561自由度3PrS0.4144对称性检验指总体上主对角线的上三角数相加是否与下三角三个数相加对称性检验与KPA检验是否一致是否一个可以代替另一个检验Pe理论观察一致率独立假设性基础上计算的相互独立简单Kappa系数Kappa0.6809渐近标准误差0.050095%置信下限0.583095%置信

8、上限0.7788H0检验:Kappa=0总体的H0下的渐近标准误差0.0597Z11.4112单侧PrZ|Z|S0.4144简单Kappa系数Kappa0.6809渐近标准误差0.050095%置信下限0.583095%置信上限0.7788加权的Kappa系数加权的Kappa0.6614渐近标准误差0.056095%置信下限0.551695%置信上限0.7711H0检验:加权的Kappa=0H0下的渐近标准误差0.0646Z10.2406单侧PrZ|Z|.0001对加权的KPA检验与简单的（利用对角线上的数据分析）加权还要利用对角线以外的数据分析样本大小=147练习六：双向无序R*C列连表用卡

9、方理论频数小于5没有超过五分之一，一般用卡方实在不行用FISHER检验超过用KPA两种血型都是按小中大排列相互不影响独立的接受H0不一致行与列变量相互不影响DATAPGM20A;DOA=1TO4;DOB=1TO3;INPUTF;OUTPUT;END;END;CARDS;431490902388410800495587950137179325;run;odshtml;PROCFREQ;WEIGHTF;TABLESA*B/CHISQ;*exact;RUN;odshtmlclose;样本大小=6094练习七：单向有序R*C秩和检验*方法1;（单因素非参数HO三个药物疗效相同H1不完全相等）DATAP

10、GM20C;DOA=1TO4;DOB=1TO3;INPUTF;OUTPUT;END;END;CARDS;15414991531504552224;run;odshtml;PROCNPAR1WAYWILCOXON;FREQF;CLASSB;VARA;RUN;*方法2;（FIQCHIM）procfreqdata=PGM20C;weightf;tablesb*a/cmhscores=rank;run;odshtmlclose;总样本大小=270练习八：双向有序属性不同R*C4种目的4种方法SPEARMAN秩相关分析DATAPGM20E;DOA=1TO3;DOB=1TO3;INPUTF;OUTPUT;

11、END;END;CARDS;215131148671011284463132;run;odshtml;PROCCORRSPEARMAN;VARAB;FREQF;RUN;odshtmlclose;统计分析与SAS实现第1次上机实习题一、定量资料上机实习题要求：先判断定量资料所对应的实验设计类型；假定资料满足参数检验的前提条件，请选用相应设计的定量资料的方差分析，并用SAS软件实现统计计算；摘录主要计算结果并合理解释，给出统计学结论和专业结论。【练习1】取4窝不同种系未成年的大白鼠，每窝3只，随机分配到三个实验组中，分别注射不同剂量雌激素，经过一定时间后处死大白鼠测子宫重量，资料见表1。问剂量和窝

12、别的各自水平下子宫重量之间的差别有无统计学意义？若剂量间差别有统计学意义，请作两两比较。表1未成年大白鼠注射不同剂量雌激素后的子宫重量窝别子宫重量(g)剂量(g/100g):0.20.40.8合计1106116145367242681152253701111333144426387192合计2603584801098【SAS程序】：HYPERLINKmhtml:file:/C:Documents%20and%20Settingsdell桌面第5次课实验课2.1转学生20121215中科院研究生院统计分析与SAS实现第1次上机实习题.mht!file:/G:中科院研究生统计学讲课2011年有关资

13、料2011年春季统计分析与SAS实现教学材料第5次课实验实验课2.1转学生程序1.sas程序1【练习2】一位工程师研究由钻头压力产生的冲力。考察了A（钻孔速度）和B（进料速度），两因素分别取2与3水平，各水平组合下均做了两次独立重复实验，资料见表2。假定资料满足参数检验的前提条件，且两因素对观测结果的影响地位平等，已知冲力越小越好，试作分析，尽可能给出较为明确的统计和专业结论。表2在钻孔速度和进料速度取不同水平的条件下冲力的测定结果钻孔冲力（单位）速度进料速度：0.0150.0300.0451252.702.452.602.782.492.722002.832.852.862.862.802.

14、87【SAS程序】：HYPERLINKmhtml:file:/C:Documents%20and%20Settingsdell桌面第5次课实验课2.1转学生20121215中科院研究生院统计分析与SAS实现第1次上机实习题.mht!file:/G:中科院研究生统计学讲课2011年有关资料2011年春季统计分析与SAS实现教学材料第5次课实验实验课2.1转学生程序2.sas程序2.二、定性资料上机实习题要求：（1）若题目中未给出表格，请列出标准的列联表，并对其命名；（2）若题目中已列出不规范的表格，先修改，然后对其命名；（3）根据分析目的或自己提出分析目的、资料的前提条件选用相应的统计分析方法，

15、并用SAS软件实现计算；（4）将主要计算结果摘录出来，给出统计学和专业结论。【练习1】某卫生防疫站对屠宰场及肉食零售点的猪肉,检查其表层沙门氏菌带菌情况，结果如下表。试比较屠宰场与肉食零售点猪肉表层沙门氏菌的带菌率之间差别有无统计学意义?表1屠宰场及肉食零售点猪肉表层沙门氏菌抽检结果采样地点检查数阳性数值带菌率(%)屠宰场2827.14零售点29827.59【SAS程序】：HYPERLINKmhtml:file:/C:Documents%20and%20Settingsdell桌面第5次课实验课2.1转学生20121215中科院研究生院统计分析与SAS实现第1次上机实习题.mht!file:/

16、G:中科院研究生统计学讲课2011年有关资料2011年春季统计分析与SAS实现教学材料第5次课实验实验课2.1转学生练习1.sas练习1【练习2】有人对某部门22707名雇员中，普查了HBsAg,其中3454名阳性，19253名为阴性。从1975年起，追踪了3年，发现在阳性组有40名患了肝癌，阴性组仅一名患肝癌。试选用合适的方法对资料进行全面分析。【SAS程序】：HYPERLINKmhtml:file:/C:Documents%20and%20Settingsdell桌面第5次课实验课2.1转学生20121215中科院研究生院统计分析与SAS实现第1次上机实习题.mht!file:/G:中科院

17、研究生统计学讲课2011年有关资料2011年春季统计分析与SAS实现教学材料第5次课实验实验课2.1转学生练习2.sas练习2【练习3】APOE-4等位基因与老年痴呆性的关联研究：以600名晚发及散发老年痴呆患者和400名正常对照为研究对象，分析APOE-4等位基因与老年痴呆性的关系。表4APOE-4等位基因与老年痴呆性病例对照关联研究AGT等位基因例数合计病例组对照组APOE-424060300非APOE-4360340700合计6004001000【SAS程序】：HYPERLINKmhtml:file:/C:Documents%20and%20Settingsdell桌面第5次课实验课2.

18、1转学生20121215中科院研究生院统计分析与SAS实现第1次上机实习题.mht!file:/G:中科院研究生统计学讲课2011年有关资料2011年春季统计分析与SAS实现教学材料第5次课实验实验课2.1转学生练习3.sas练习3【练习4】请分析下表资料。已从专业上认定培养的阳性结果就是“真阳性”，而不会出现假阳性。表3两种培养基对同一批痰液标本同时培养的结果甲培养基培养结果例数乙培养基结果：+-合计+363470-0135135合计36169205配对设计22列联表资料总体率差异性检验统计量的计算公式若b+c40时若b+c40时【SAS程序】：HYPERLINKmhtml:file:/C:

19、Documents%20and%20Settingsdell桌面第5次课实验课2.1转学生20121215中科院研究生院统计分析与SAS实现第1次上机实习题.mht!file:/G:中科院研究生统计学讲课2011年有关资料2011年春季统计分析与SAS实现教学材料第5次课实验实验课2.1转学生练习4.sas练习4【练习5】请分析下表资料。表6两法检查室壁收缩运动的符合情况对比法测冠心病人数定的结果核素法正常减弱异常正常5823减弱1427异常8917合计675327【SAS程序】：HYPERLINKmhtml:file:/C:Documents%20and%20Settingsdell桌面第5

20、次课实验课2.1转学生20121215中科院研究生院统计分析与SAS实现第1次上机实习题.mht!file:/G:中科院研究生统计学讲课2011年有关资料2011年春季统计分析与SAS实现教学材料第5次课实验实验课2.1转学生练习5.sas练习5简单kappa检验和加权kappa检验这两种方法都是用来检验两种评价方法是否具有一致性的方法。其主要的区别是两种方法计算的公式不一样，更具体地说是对两个变量的打分不一样，简单kappa检验主要是利用对角线上的信息，加权kappa检验除了利用对角线上的数据外，还将对角线外的数据进行加权打分，将对角线外的信息也充分利用。所以在选择方法时应根据专业知识，如果两个变量取值的界线比较明确，如“+”“+”“+”等，这时可以选用简单的kappa检验，如果两个变量的取值不十分明确，人为的因素较多时，就可选用加权kappa检验。【练习6】请分析下表资料。