2019精选医学SAS练习题及程序答案

1、1 .随机取组有无重复试验的两种本题是无重复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 /

2、 SNK;/*对因素B （即剂量）各水平下的均值进行两两比较*/RUN;ODS HTML CLOSE;2 . 2*3析因设计两因素完全随机统计方法2*3析因设计tiff =f的开方 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.702.452.602.782.492.722.832.852.862.862.802.87JrunPROC GLM; CLASS zs js;MODEL

3、 cl=zs js zs*js /SS3 ;MEANS zs*js;LSMEANS zs*js / TDIFF PDIFF ; /*对zs和js各水平组合而成的试验条件进行均数 进行两两比较*/ RUN;ODS HTML CLOSE;练习一：2*2横断面研究列链表 方法：卡方 矫正卡方FISHER DATA PGM19A;DO A= 1 TO 2;DO B= 1 TO 2;INPUT F ;OUTPUT;END;END;CARDS;2 268 21 ;run ;PROC FREQ;WEIGHT F;TABLES A*B / CHISQ ;RUN;样本大小=57 练习二：对裂列连表结果变量换和不

4、换三部曲1横断面研究P0.05 R DATA PGM19B;DO A= 1 TO 2;DO B= 1 TO 2;INPUT F ;OUTPUT;END;END;CARDS;40 3414 1 19252 ;run ;ods htmlPROC FREQ;WEIGHT F;TABLES A*B / CHISQ cmh ;RUN;ods html close ;样本大小=57练习三：病例对照 2*2病例组中有何没有那个基因是正常的3.8倍，则有可能导致痴呆 要做前瞻性研究 用对裂DATA PGM20;DO A= 1 TO 2;DO B= 1 TO 2;INPUT F ;OUTPUT;END;END;

5、CARDS;240 60360 340;run ;ods html ;PROC FREQ;WEIGHT F;TABLES A*B / CHISQ cmh;RUN;ods html close ;总样本大小 =1000练习四：配对设计 隐含金标准2*2MC卡方检3& 34和0在总体上(B+C40用矫正卡方)是否相等则可得甲培养基优于乙培养基一般都用矫正因卡方为近似计算DATA PGM19F;INPUT b c;chi=(ABS(b-c)- 1)* 2/(b+c);p= 1-PROBCHI(chi, 1);求概率1减掉从左侧积分到卡方的值chi=ROUND(chi, 0.001 );IF

6、p> 0.0001 THEN p=ROUND(p, 0.0001 );FILE PRINT ;PUT (打印在输出床口)# 2 10 'Chisq' 30 'P value'(#表示行)#4 10 chi 30 p;CARDS;34 0run ;ods html close练习五：双向有序 R*C列连表 用KPA data aaa;do a= 1 to 3;do b= 1 to 3;input f ;output ;end ;end ;cards ;582314278917;run ;ods html ;*简单kappa检马叙;proc freq data

7、 =aaa; weight f;（频数） tables a*b;test kappa ;run ;*加权kappa检马级;proc freq ;weight f;tables a*b;test wtkap ;run ;ods html closeIsas系统FREQ过程a * b表ab合计12315839.4621.3632.046342.86频数 百分比 行百分比 列百分比92.0686.573.173.774.7611.1121427500.6828.574.7634.012.0084.0014.001.4979.2525.9338917345.446.1211.5623.1323.532

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

9、单侧Pr> Z<.0001双侧Pr>|Z|<.0001总体的KPA是否为0 KPA大于0两种方法的一致性有统计学意义小于0不一致性有统计学意义加权的Kappa系数加权的Kappa0.6614渐近标准误差0.056095%置信卜限0.551695%置信上限0.7711置信区间不包括 0拒绝H0按此计算结果可以用一种取代另一种方法但要看专业要求达到多少才可以 观测一致率达到多少才可以代替样本大小 =147FREQ过程a * b表ab合计123158236339.461.362.0442.8692.063.174.7686.573.7711.1121427500.6828.5

10、74.7634.012.0084.0014.001.4979.2525.9338917345.446.1211.5623.1323.5326.4750.0011.9416.9862.96合计6745.585336.052718.37147100.00a * b表的统计量对称性检验统计量(S)2.8561自由度3Pr > S0.4144简单Kappa系数Kappa0.6809渐近标准误差0.050095%置信卜限0.583095%置信上限0.7788加权的Kappa系数加权的Kappa0.6614渐近标准误差0.056095%置信卜限0.551695%置信上限0.7711H0检验：加权的

11、Kappa = 0H0下的渐近标准误差0.0646Z10.2406单侧Pr> Z<.0001双侧Pr>|Z|<.0001对加权的KPA检验 与简单的（利用对角线上的数据分析）加权还要利用对角线以外的数据分析样本大小 =147练习六：双向无序R*C列连表 用卡方理论频数小于5没有超过五分之一，一般用卡方 实在不行用FISHER检验超过用KPA两种血型都是按小中大排列相互不影响 独立的 接受H0不一致行与列变量相互不影响DATA PGM20A;DO A= 1 TO 4;DO B= 1 TO 3;INPUT F ;OUTPUT;END;END;CARDS;431 490 90

12、2388 410 800495 587 950137 179 325 ;run ;ods html ;PROC FREQ;WEIGHT F;TABLES A*B / CHISQ ; *exact;RUN;ods html close ;样本大小=6094练习七：单向有序 R*C秩和检验*方法1;（单因素非参数 HO三个药物疗效相同H1不完全相等）DATA PGM20C;DO A= 1 TO 4;DO B= 1 TO 3;INPUT F ;OUTPUT;END;END;CARDS;15 4 149 9 1531 50 455 22 24;run ;ods htmlPROC NPAR1WAY WI

13、LCOXON;FREQ F;CLASS B;VAR A;RUN;* 方法 2; (FIQ CHIM )proc freq data =PGM20C; |weight f;tables b*a/ cmh scores =rank;run ;ods html close ;总样本大小 =270练习八： 双向有序 属性不同 R*C 4种目的4种方法SPEARMAN秩相关分析DATA PGM20E;DO A= 1 TO 3;DO B= 1 TO 3;INPUT F ;OUTPUT;END;END;CARDS;215 131 14867 101 12844 63 132;run ;ods html ;P

14、ROC CORR SPEARMAN;VAR A B;FREQ F;RUN;ods html close ;统计分析与SAS实现第1次上机实习题一、定量资料上机实习题要求：(1) 先判断定量资料所对应的实验设计类型；(2) 假定资料满足参数检验的前提条件，请选用相应设计的定量资料的方差分析，并用SASa件实现统计计算；(3) 摘录主要计算结果并合理解释，给出统计学结论和专业结论。【练习1】取4窝不同种系未成年的大白鼠，每窝3只，随机分配到三个实验组中，分别注射不同剂量雌激素，经过一定时间后处死大白鼠测子宫重量，资料见表1。问剂量和窝别的各自水平下子宫重量之间的差别有无统计学意义？若剂量间差别有统

15、计学意义，请作两两比较。表1未成年大白鼠注射不同剂量雌激素后的子宫重量窝别(g)剂量(g g/100g):0.20.40.81106116145367242681152253701111333144426387192合计2603584801098【SAS程序】：程序1【练习2】一位工程师研究由钻头压力产生的冲力。考察了 A (钻孔速度)和 B (进料 速度)，两因素分别取 2与3水平，各水平组合下均做了两次独立重复实验，资料见表2。假定资料满足参数检验的前提条件， 且两因素对观测结果的影响地位平等， 已知冲力越小越 好，试作分析，尽可能给出较为明确的统计和专业结论。表2在钻孔速度和进料速度取不

16、同水平的条件下冲力的测定结果孔速度进料速度:0.0150.0300.0451252.702.452.602.782.492.722002.832.852.862.862.802.87【SAS程序】：程序2.二、定性资料上机实习题要求：(1)若题目中未给出表格，请列出标准的列联表，并对其命名；(2)若题目中已列出不规范的表格，先修改，然后对其命名；(3)根据分析目的或自己提出分析目的、资料的前提条件选用相应的统计分析方法，并用 SAS软件实现计算；(4)将主要计算结果摘录出来，给出统计学和专业结论。【练习1】某卫生防疫站对屠宰场及肉食零售点的猪肉，检查其表层沙门氏菌带菌情况,结果如下表。试比较屠

17、宰场与肉食零售点猪肉表层沙门氏菌的带菌率之间差别有无统计学意 义？表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等位基因与老年痴呆性的关系。表4 AP

18、OE-4等位基因与老年痴呆性病例对照关联研究AGT等位基因合计病例组对照组APOE-424060300非 APOE-43603407006004001000皿比数1）=祟Odd2（比数2）=" OR爱式【SAS程序】：练习3【练习4】请分析下表资料。已从专业上认定培养的阳性结果就是真阳性”，而不会出现假阳性。表3两种培养基对同一批痰液标本同时培养的结果甲培养基例数培养结果乙培养基结果：+-合计+363470-0135135合计36169205配对设计2X 2列联表资料总体率差异性检验统计量的计算公式若b+c>40时2 二(b - c) b c若b+cV 40时2 二(|b - c - 1 ) 2 b c【SAS程序】：练习4【练习5】请分析下表资料。表6两法检查室壁收缩运动的符合情况对比法测冠心病人数定的结果核素年云带麻 异常常弱常 十合 正减异58 1824293717675327【SAS程序】：