sas卡方检验编程语句【行业内容】
sas 2 检验编程语句,1,一、复习有关内容,1.应用: (1)分类变量资料; (2)推断两个或两个以上的样本率或构成比之间有无差异; (3)检验频数分布的拟合优度。,2,2. 四格表资料的2 检验,(1)四格表的形式 (2)2 检验的基本思想 (3) 基本公式 自由度 =(行数-1)(列数-1),3, 当Ti5,且n 40时 使用普通2 检验,可用基本 公式或专用公式,a、b、c、d为实际频数,4,四格表2 值的校正公式,当1T 5,且n40时,需校正2 检验 a .基本公式的校正: b.专用公式的校正:,5,3、行列表资料的2 检验,公式:,6,RC表资料的2检验适用条件,1.理论数不能小于1 2.理论数1T,且小于5的格子数不超过总格子数的1/5。 若条件不适合,需作如下处理: A 增大样本例数 B 删除理论数太小的行或列 C 合并(性质相同),7,二、完全随机设计的2检验(SAS),例1 某医院用内科疗法治疗一般类型胃溃疡患者80例,治愈63例;治疗特殊类型胃溃疡患者99例,治愈31例。 问内科疗法对两种类型胃溃疡的治愈率差别有无显著意义?,8,表 两种类型胃溃疡内科疗法治疗结果,9,data ex1; input r c count ; cards; 1 1 63 1 2 17 2 1 31 2 2 68 3 1 322 3 2 ; proc freq; tables r*c/chisq expected nopercent nocol; weight count; run;,10,结果,11,12,结果解释,本例n40且各格子的期望值均大于5,因而选用Chi-Square的2统计量及其显著性水平,即239.927,P=0.0001,拒绝H0,认为内科疗法对两种类型胃溃疡的治愈率差别有统计学意义,一般类型的治愈率高于特殊型。,13,例2 某省三地区花生黄曲霉素B1污染率比较,地区,未污染,污染,合计,污染率(%),甲 乙,6 30 8,23 14 3,29 44 11,79.3 31.8 27.3,丙,合计,44,40,84,47.6,14,程序1:一般输入方法,data ex2; input r c count; cards; 1 1 6 1 2 23 2 1 30 2 2 14 3 1 8 3 2 3 ;,15,程序2:循环输入法,data ex2; do r=1 to 3; do c=1 to 2; input count; output;end;end; cards; 6 23 30 14 8 3 ;,16,proc freq; tables r*c/chisq expected nopercent nocol; weight count; run;,17,结果,18,本例各格子期望值均大于5,选用Chi-Square的2统计量及其显著水平,即2=17.907,P=0.0001,按=0.05的检验水准拒受H0,认为三地花生黄曲霉素B1污染率有差别。,19,配对设计的2检验(SAS程序),200名已确诊的血吸虫患者,治疗前经皮试法及粪检法检查,结果如下表,问两种检查方法的结果有无差别?,20,表2 血吸虫患者两种检查方法的结果,21,data ex3; do r=1 to 2; do c=1 to 2; input f; output;end;end; cards; 112 40 12 36 ;,22,proc freq; tables r*c/chisq expected nopercent nocol agree; weight f; run;,此处用Agree选项实现SAS的McNemar检验,但要注意,SAS的McNemar检验不考虑 校正问题,只适合于较大样本(b+c40),所得的结果与 公式的计算结果相同。,23,Data ex4; Input f11 f12 f21 f22; If f12+f2140 then Chisq=(abs(f12-f21)-1)*2/(f12+f21); Else chisq=abs(f12-f21)*2/(f12+f21); P=1-probchi(chisq,1); Cards; 112 40 12 36 ;proc print; Run;,此法适用于b+c=40的时候, 使用公式编辑法计算2统计量及对应的P值。,24,25,此部分结果是普通四格表2检验的结果, 不适于配对2检验使用。,26,27,本例b+c40,故选用Mcnemar检验的2统计量及其显著水平,即2=15.0769,P0.0001,按=0.05的检验水准拒受H0,认为两种检查方法有差别。,28,公式法得到的结果:,29,作业 1. P125 1; 2. P126 4; 3. 程序实现课本中确切概率法的例题的数据,并记录最后的结果。,30,