sas卡方检验编程语句

sasX2检验编程语句、复习有关内容L应用：（1） 分类变量资料；（2） 推断两个或两个以上的样本率或构成比之间有无差异；（3） 检验频数分布的拟合优度。2.四格表资料的X2检验（1） 四格表的形式（2） X2检验的基本思想（3） 基本公式尸=£以也《」T自由度卩=（行数■：!）（列数①当L，5,且nN40时使用普通烂检验，可用基本公式或专用公式(ad-bc)2nx2(ad-bc)2n(a+b)(c+d)(a+c)(b+d)a、b、c、d为实际频数②四格表为2值的校正公式当1<T<5,且n*0时，需校正%2检验a.基本公式的校正： (|A-T|-0.5)2%=x亍b■专用公式的校正W /弋:弋］(Q+/?)(c+d)(Q+測+d)3、行x列表资料的%2检验公式：2 A2?xz=n〉 1InRnc丿RXC表资料的为2检验适用条件1.理论数不能小于12,理论数1VT,且小于5的格子数不超过总格子数的1/5。若条件不适合，需作如下处理：A增大样本例数B删除理论数太小的行或列C合并（性质相同）二、完全随机设计的烂检验(SAS)例1某医院用内科疗法治疗一般类型胃溃疡患者80例，治愈63例；治疗特殊类型胃溃疡患者99例，治愈31例。问内科疗法对两种类型胃溃疡的治愈率差别有无显著意义？表两种类型胃溃疡内科疗法治疗结果dataexl;inputrccount@@;cards;11631217213122683132232procfreq;tablesr*c/chisqexpectednopercentnocol;weightcount;run;

TheFREQProcedure结果TableoTrbyc结果TableoTrbyc• + +Total9485179

Statistic 54曰-似然比X2慎DFValueProbGhi・Square 直]139.9267<-0001LikelihoodRatioChi-Square141.8605<.0001ContinuityAdj・Chi-Square38.0471<-0001Mante1-Haensze1Chi-Square-■-W^-ZQ36<-0001PhiCoefficient0.4723连续校正X2值ContingencyCoefficient0.4271Cramer'sV0.4723StatisticsforTableofrbyc一一^切概率法Fisher*sExactTestCell(1,1)Frequency(F)63

Fisher*sExactTest1.275E-101-642E-101.275E-101-642E-10结果解释本例n>40且各格子的期望值均大于5,因而选用Chi-Square的%2统计量及其显著性水平，即x2=39.927,P=0.0001,拒绝H°,认为内科疗法对两种类型胃溃疡的治愈率差别有统计学意义，一般类型的治愈率高于特殊型。

例2某省三地区花生黄曲霉素B］污染率比较地区未污染污染合计污染率（%）甲6232979.3乙30144431.8丙831127.3合计44408447.6程序1:一般输入方法dataex2;inputrccount@@;cards;116122321302214318323程序2：循环输入法dataex2;dor=lto3;doc=lto2;inputcount@@;output;end;end;cards;623301483procfreq;tablesr*c/chisqexpectednopercentnocol;weightcount;run;TableoTrbyc结果实际频戮Total23—15.1913.8120.6979.31行百分比Frequency|Expected|RowPcT结果实际频戮Total23—15.1913.8120.6979.31行百分比Frequency|Expected|RowPcT30|23.048I144430|23.048I144420.95268.1831.8268.1831.825.23815.2381Total444084

Total444084StatisticsforTableofrbycStatisticDFValueProbGhi-Square:LikelihoodRatioChi-SquareMantel-HaenszelChi-SquarePhiCoefficientContingencyCoefficientStatisticDFValueProbGhi-Square:LikelihoodRatioChi-SquareMantel-HaenszelChi-SquarePhiCoefficientContingencyCoefficientCramer'sV17,906618.754714.31550.46170.41920.46170.0001<.00010.0002本例各格子期望值均大于5,选用Chi-Square的为2统计冒及其豆薯水平，即'2=17.907,P=0.0001,按a=0.05的检验水准拒受H°,认为三地花生黄曲霉素％污染率有差别。

配对设计的产检验（SAS程序）200名已确诊的血吸虫患者，治III疗前经皮试法及粪检法检查，结果如下表，问两种检查方法的结果有无差别？表2血吸虫患者两种检查方法的结果IIII■■■■■■粪检法皮试法++.112 4012 36合计15248合计124 76 200dataex3;dor=lto2;doc=lto2;inputf@@;output;end;end;cards;112401236procfreq;tablesr*c/chisqexpectednopercentnocolagree;weightf;run;此处用Agree选项实现SAS的McNemar检验,但要注意，SAS曲McNemar检验不考虑校正问题，只适合于较大样本(b+c>40),所得的结果与2_(E)2公式的计算结果相同。 0+C)Dataex4;Inputfllf!2f21f22;Iffl2+f21v40thenChisq=(abs(fl2・f21)・l)**2/(fl2+f21);ElseChisq=abs(fl2・f21)**2/(fl2+f21);P=l-probchi(chisq,l)；Cards;112401236;procprint;Run;此法适用于b+C<=40的时候，使用公式编辑法计算%2统计量及对应的P值o

TheSASSystem16:25Tuesday,NovemberTheSASSystem16:25Tuesday,NovemberTheFREQProcedureTableofrbycFrequencyExpectedRowPct12Total111294.2473.684057.7626.3215221229.7625.003618.2475.0048Total12476200

StatisticsforTableofrbycStatisticDFValueProbChi-Square136.6991(.0001LikelihoodRatioChi-Square136.4359(.0001ContinuityAdj.Chi-Square134.6618(.0001Mantel-Haensze1Chi-Square136.5156(.0001PhiCoefficient0.4284ContingencyCoefficient0.3938Cramer!sV0.4284!1!1!此部分结果是普通四格表%2检验的结果,不适于配对乂2检验使用。

r 21!I960年Cohen等提出用Kappa值作为评价判断的一致性程度的指标。当两个诊断完全一致时，Kappa值为1。当观测一致率大于期望一致率时，Kappa值为正数，且Kappa值越大，说明一致性越好。当观察一致率小于期望一致率时，Kappa值为负数，这种情况一般来说比较少见。根据边缘概率的计算，Kappa值的范围值应在・：L~：L之间。Kappa>0.75两者一致性较好；0.75>KappaN0.4两者一敏:性一般；r 21!SimpleKappaCoefficient0.40590.06560.27740.53430.40590.06560.27740.5343955?LowerConfLimit95%UpperConfLimitSamp1eSize=200本例b+c>40,故选用Mcnemar检验的%2统计量及其显著水平，即x2=