sas卡方检验编程语句分析.ppt

《sas卡方检验编程语句分析.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、sasχ2检验编程语句一、复习有关内容1.应用:（1）分类变量资料;（2）推断两个或两个以上的样本率或构成比之间有无差异;（3）检验频数分布的拟合优度。2.四格表资料的χ2检验（1）四格表的形式（2）χ2检验的基本思想（3）基本公式自由度=(行数-1)(列数-1)①当Ti≥5,且n≥40时使用普通2检验,可用基本公式或专用公式a、b、c、d为实际频数②四格表2值的校正公式当1≤T<5,且n≥40时,需校正2检验a.基本公式的校正:b.专用公式的校正:3、行×列表资料的2检验公式:R×C表资料的2检验适用条件1.理论数不能小于12.理论数

2、1

3、weightcount;run;结果实际频数理论频数结果解释本例n>40且各格子的期望值均大于5,因而选用Chi-Square的2统计量及其显著性水平,即2＝39.927,P=0.0001,拒绝H0,认为内科疗法对两种类型胃溃疡的治愈率差别有统计学意义,一般类型的治愈率高于特殊型。例2某省三地区花生黄曲霉素B1污染率比较地区未污染污染合计污染率（%）甲乙63082314329441179.331.827.3丙合计44408447.6程序1:一般输入方法dataex2;inputrccount@@;cards;1161223213022143183

4、23;程序2:循环输入法dataex2;dor=1to3;doc=1to2;inputcount@@;output;end;end;cards;623301483;procfreq;tablesr*c/chisqexpectednopercentnocol;weightcount;run;结果本例各格子期望值均大于5,选用Chi-Square的2统计量及其显著水平,即2=17.907,P=0.0001,按=0.05的检验水准拒受H0,认为三地花生黄曲霉素B1污染率有差别。配对设计的2检验（SAS程序）200名已确诊的血吸虫患者,治疗前经皮试法

5、及粪检法检查,结果如下表,问两种检查方法的结果有无差别?粪检法合计+-皮试法+11240152-123648合计12476200表2血吸虫患者两种检查方法的结果dataex3;dor=1to2;doc=1to2;inputf@@;output;end;end;cards;112401236;procfreq;tablesr*c/chisqexpectednopercentnocolagree;weightf;run;此处用Agree选项实现SAS的McNemar检验,但要注意,SAS的McNemar检验不考虑校正问题,只适合于较大样本（b+c>40）

6、,所得的结果与公式的计算结果相同。Dataex4;Inputf11f12f21f22;Iff12+f21<40thenChisq=(abs(f12-f21)-1)**2/(f12+f21);Elsechisq=abs(f12-f21)**2/(f12+f21);P=1-probchi(chisq,1);Cards;112401236;procprint;Run;此法适用于b+c<=40的时候,使用公式编辑法计算2统计量及对应的P值。此部分结果是普通四格表2检验的结果,不适于配对2检验使用。1960年Cohen等提出用Kappa值作为评价判断的一

7、致性程度的指标。当两个诊断完全一致时，Kappa值为1。当观测一致率大于期望一致率时，Kappa值为正数，且Kappa值越大，说明一致性越好。当观察一致率小于期望一致率时，Kappa值为负数，这种情况一般来说比较少见。根据边缘概率的计算，Kappa值的范围值应在-1~1之间。Kappa≥0.75两者一致性较好；0.75>Kappa≥0.4两者一致性一般；Kappa<0.4两者一致性较差。本例b+c>40,故选用Mcnemar检验的2统计量及其显著水平,即2=15.0769,P＝0.0001,按=0.05的检验水准拒受H0,认为两种检查方法有差别

8、。公式法得到的结果:作业1.P1251;2.P1264;3.程序实现课本中确切概率法的例题的数据,并记录最后