SAS练习题及程序答案.docx
1随机取组有无重复试验的两种本题是无重复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;132*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;MEANSzs*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;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练习四：配对设计隐含金标准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);IFp>0.0001THENp=ROUND(p,0.0001);FILEPRINT;PUT（打印在输出床口）#210Chisq30Pvalue（#表示行）#410chi30p;CARDS;340;run;odshtmlclose;练习五：双向有序R*C列连表用KPAdataaaa;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.862310.682.001.4985.4423.5311.944228.5784.0079.2596.1226.4716.9874.7614.0025.931711.5650.0062.965034.013423.13合计6745.585336.052718.37147100.00a*b表的统计量对称性检验统计量(S)自由度Pr>S2.856130.4144对称性检验指总体上主对角线的上三角数相加是否与下三角三个数相加对称性检验与KPA检验是否一致是否一个可以代替另一个检验Pe理论观察一致率独立假设性基础上计算的相互独立简单Kappa系数Kappa渐近标准误差0.68090.0500简单Kappa系数95%置信下限95%置信上限0.58300.7788H0检验:Kappa=0总体的H0下的渐近标准误差Z单侧Pr>Z双侧Pr>|Z|0.059711.4112<.0001<.0001总体的KPA是否为0KPA大于0两种方法的一致性有统计学意义小于0不一致性有统计学意义加权的Kappa系数加权的Kappa渐近标准误差95%置信下限95%置信上限0.66140.05600.55160.7711置信区间不包括0拒绝H0按此计算结果可以用一种取代另一种方法但要看专业要求达到多少才可以观测一致率达到多少才可以代替样本大小=147FREQ过程频数百分比行百分比列百分比a*b表ab合计123158236339.461.362.0442.8692.063.174.7686.573.7711.1121427500.6828.574.7634.012.0084.0014.001.4979.2525.9338917345.446.1211.5623.1323.5326.4750.0011.9416.9862.96合计67532714745.5836.0518.37100.00a*b表的统计量对称性检验统计量(S)2.8561自由度3Pr>S0.4144简单Kappa系数Kappa0.6809渐近标准误差0.050095%置信下限0.583095%置信上限0.7788加权的Kappa系数加权的Kappa渐近标准误差95%置信下限95%置信上限0.66140.05600.55160.7711H0检验:加权的Kappa=0H0下的渐近标准误差Z单侧Pr>Z双侧Pr>|Z|0.064610.2406<.0001<.0001对加权的KPA检验与简单的（利用对角线上的数据分析）加权还要利用对角线以外的数据分析样本大小=147练习六：双向无序R*C列连表用卡方理论频数小于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不完全相等）DATAPGM20C;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;END;END;CARDS;215131148671011284463132;run;odshtml;PROCCORRSPEARMAN;VARAB;FREQF;RUN;odshtmlclose;统计分析与SAS实现第1次上机实习题一、定量资料上机实习题要求：（1）先判断定量资料所对应的实验设计类型；（2）假定资料满足参数检验的前提条件，请选用相应设计的定量资料的方差分析，并用SAS软件实现统计计算；（3）摘录主要计算结果并合理解释，给出统计学结论和专业结论。请【练习1】取4窝不同种系未成年的大白鼠，每窝3只，随机分配到三个实验组中，分别注射不同剂量雌激素，经过一定时间后处死大白鼠测子宫重量，资料见表1。问剂量和窝别的各自水平下子宫重量之间的差别有无统计学意义？若剂量间差别有统计学意义，作两两比较。表1未成年大白鼠注射不同剂量雌激素后的子宫重量子宫重量(g)窝别1234合计剂量(g/100g):0.21064270422600.411668111633580.814511513387480合计3672253141921098【SAS程序】：程序1【练习2】一位工程师研究由钻头压力产生的冲力。考察了A（钻孔速度）和B（进料速度），两因素分别取2与3水平，各水平组合下均做了两次独立重复实验，资料见表2。假定资料满足参数检验的前提条件，且两因素对观测结果的影响地位平等，已知冲力越小越好，试作分析，尽可能给出较为明确的统计和专业结论。AGT等位基因例数表2在钻孔速度和进料速度取不同水平的条件下冲力的测定结果钻孔冲力（单位）速度进料速度：0.0150.0300.0451252.702.452.602.782.492.722002.832.852.862.862.802.87【SAS程序】：程序2.二、定性资料上机实习题要求：（1）若题目中未给出表格，请列出标准的列联表，并对其命名；（2）若题目中已列出不规范的表格，先修改，然后对其命名；（3）根据分析目的或自己提出分析目的、资料的前提条件选用相应的统计分析方法，并用SAS软件实现计算；（4）将主要计算结果摘录出来，给出统计学和专业结论。【练习1】某卫生防疫站对屠宰场及肉食零售点的猪肉,检查其表层沙门氏菌带菌情况，结果如下表。试比较屠宰场与肉食零售点猪肉表层沙门氏菌的带菌率之间差别有无统计学意义?表1屠宰场及肉食零售点猪肉表层沙门氏菌抽检结果采样地点检查数阳性数值带菌率(%)屠宰场2827.14零售点29827.59【SAS程序】：练习1【练习2】有人对某部门22707名雇员中，普查了HBsAg,其中3454名阳性，19253名为阴性。从1975年起，追踪了3年，发现在阳性组有40名患了肝癌，阴性组仅一名患肝癌。试选用合适的方法对资料进行全面分析。【SAS程序】：练习2【练习3】APOE-4等位基因与老年痴呆性的关联研究：以600名晚发及散发老年痴呆患者和400名正常对照为研究对象，分析APOE-4等位基因与老年痴呆性的关系。表4APOE-4等位基因与老年痴呆性病例对照关联研究病例组对照组合计APOE-4非APOE-424036060340300700合计6004001000a/(a+c)b/(b+d)odda/cc/(a+c)d/(b+d)oddb/dodd(比数)=odd(比数)=OR=1=22112【SAS程序】：练习3【练习4】请分析下表资料。已从专业上认定培养的阳性结果就是“真阳性”，而不会出现假阳性。表3两种培养基对同一批痰液标本同时培养的结果甲培养基例数培养结果+-合计乙培养基结果：+36036-34135169合计70135205配对设计22列联表资料总体率差异性检验统计量的计算公式若b+c40时c=n=1若b+c40时2(b-c)2b+cc=n=1(b-c-1)22b+c【SAS程序】：练习4【练习5】请分析下表资料。表6两法检查室壁收缩运动的符合情况对比法测冠心病人数定的结果核素法正常减弱异常正常5823减弱1427异常8917合计675327【SAS程序】：练习5不超过1/5的格子理论频数<5时，此类资料应采用检验处理。+简单kappa检验和加权kappa检验这两种方法都是用来检验两种评价方法是否具有一致性的方法。其主要的区别是两种方法计算的公式不一样，更具体地说是对两个变量的打分不一样，简单kappa检验主要是利用对角线上的信息，加权kappa检验除了利用对角线上的数据外，还将对角线外的数据进行加权打分，将对角线外的信息也充分利用。所以在选择方法时应根据专业知识，如果两个变量取值的界线比较明确，如“”“+”“+”等，这时可以选用简单的kappa检验，如果两个变量的取值不十分明确，人为的因素较多时，就可选用加权kappa检验。【练习6】请分析下表资料。表5某地6094人按2种血型系统划分的结果ABO人数血型MN血型:MNMNO431490902A388410800B495587950AB137179325合计145116662977【SAS程序】：练习6分析方法选择2其他情形，可采用Fisher精确检验处理。【练习7】请分析下表资料。表73种药物疗效的观察结果疗人数效药物ABC治愈1541显效49915好转315045无效52224合计1008585【SAS程序】：练习7【练习8】请分析下表资料。表8眼晶状体混浊度与年龄之关系晶状体混眼数浊程度年龄203040215131148+67101128+4463132合计326295408
