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