SAS练习题及程序答案.doc

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 / S

2、NK; /*对因素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.70 2.45 2.602.78 2.49 2.722.83 2.85 2.862.86 2.80 2.87;run;PROC GLM;CLASS zs js

3、;MODEL 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 RDATA PGM19B;DO A=1 TO 2;DO B=1 TO 2;INPUT F ;OUTPUT;END;END;CARDS;40 34141 19252;run;ods html;PROC 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;

5、END; END;CARDS;24060360340;run;ods html;PROC FREQ;WEIGHT F;TABLES A*B / CHISQ cmh;RUN;ods html close;总样本大小 = 1000练习四：配对设计 隐含金标准2*2 MC卡方 检验 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);I

6、F p0.0001 THEN p=ROUND(p,0.0001);FILE PRINT;PUT（打印在输出床口） #2 10 Chisq 30 P value（#表示行）#4 10 chi 30 p;CARDS;34 0;run;ods html close;练习五：双向有序R*C列连表 用KPA data aaa;do a=1 to 3;do b=1 to 3;input f ;output;end; end; cards;58 2 3 1 42 7 8 9 17 ;run;ods html;*简单kappa 检验;proc freq data=aaa; weight f;（频数）table

7、s a*b;test kappa;run;*加权kappa 检验;proc freq; weight f;tables a*b;test wtkap;run;ods html close;SAS 系统FREQ 过程频数 百分比 行百分比 列百分比 a * b 表ba1 2 3 合计1 5839.4692.0621.363.1732.044.766342.8686.57 3.77 11.11 2 10.682.001.494228.5784.0079.2574.7614.0025.935034.013 85.4423.5311.9496.1226.4716.981711.5650.0062.96

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

9、Z| S 0.4144简单 Kappa 系数Kappa 0.6809渐近标准误差 0.050095% 置信下限 0.583095% 置信上限 0.7788加权的 Kappa 系数加权的 Kappa 0.6614渐近标准误差 0.056095% 置信下限 0.551695% 置信上限 0.7711H0 检验: 加权的 Kappa = 0H0 下的渐近标准误差 0.0646Z 10.2406单侧 Pr Z |Z| .0001对加权的 KPA 检验 与简单的（利用对角线上的数据分析）加权还要利用对角线以外的数据分析 样本大小 = 147练习六：双向无序R*C 列连表 用卡方理论频数小于5没有超过五分

10、之一，一般用卡方 实在不行用FISHER 检验 超过用KPA 两种血型都是按小中大排列 相互不影响 独立的 接受H0 不一致 行与列变量相互不影响DATA PGM20A;DO A=1 TO 4;DO B=1 TO 3;INPUT F ;OUTPUT;END;END;CARDS;431 490 902388 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;（单因素非参

11、数 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 html;PROC NPAR1WAY WILCOXON;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

12、*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;PROC CORR SPEARMAN;VAR A B; FREQ F;RUN;ods html close;统计分析与 SAS 实现第 1 次上机实习题一、定量资料上机实习题要求：（1） 先判断定量资料所对应的实验设计类型；（2） 假定资料满足参数检验的前提条件，请选用相应设计的定量资料的方差分析，并用SAS 软件实现统计

13、计算；（3） 摘录主要计算结果并合理解释，给出统计学结论和专业结论。【练习1】取4窝不同种系未成年的大白鼠，每窝3只，随机分配到三个实验组中，分别注射不同剂量雌激素，经过一定时间后处死大白鼠测子宫重量，资料见表1。问剂量和窝别的各自水平下子宫重量之间的差别有无统计学意义？若剂量间差别有统计学意义，请作两两比较。表 1 未 成 年 大 白 鼠 注 射 不 同 剂 量 雌 激 素 后 的 子 宫 重 量子宫重量(g)窝别剂量(g/100g):0.2 0.4 0.8 合计1 106 116 145 3672 42 68 115 2253 70 111 133 3144 42 63 87 192合计

14、260 358 480 1098【SAS 程序】：程序 1【练习 2】一位工程师研究由钻头压力产生的冲力。考察了 A（钻孔速度）和 B（进料速度） ，两因素分别取 2 与 3 水平，各水平组合下均做了两次独立重复实验，资料见表 2。假定资料满足参数检验的前提条件，且两因素对观测结果的影响地位平等，已知冲力越小越好，试作分析，尽可能给出较为明确的统计和专业结论。表 2 在钻孔速度和进料速度取不同水平的条件下冲力的测定结果钻 孔 冲力（单位） 速 度 进料速度：0.015 0.030 0.045125 2.70 2.45 2.602.78 2.49 2.72200 2.83 2.85 2.862.

15、86 2.80 2.87【SAS 程序】： 程序 2.二、定性资料上机实习题要求：（1）若题目中未给出表格，请列出标准的列联表，并对其命名；（2）若题目中已列出不规范的表格，先修改，然后对其命名；（3）根据分析目的或自己提出分析目的、资料的前提条件选用相应的统计分析方法，并用SAS 软件实现计算；（4）将主要计算结果摘录出来，给出统计学和专业结论。【练习 1】某卫生防疫站对屠宰场及肉食零售点的猪肉,检查其表层沙门氏菌带菌情况，结果如下表。试比较屠宰场与肉食零售点猪肉表层沙门氏菌的带菌率之间差别有无统计学意义?表 1 屠宰场及肉食零售点猪肉表层沙门氏菌抽检结果采样地点 检查数 阳性数值 带菌率(%)屠宰场 28 2 7.14零售点 29 8 27.59【SAS 程序】： 练习 1【练习 2】有人对某部门 22707 名雇员中，普查了 HBsAg,其中 3454 名阳性，19253 名为阴性。从 1975 年起，追踪了 3 年，发现在阳性组有 40 名患了肝癌，阴性组仅一名患肝癌。试选用合适的方法对资料进行全面分析。