如何在SAS中实现R×C列联表的两两比较.pdf

《如何在SAS中实现R×C列联表的两两比较.pdf》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、A14-如何在SAS中实现R×C列联表的两两比较内容来自网络，侵删在分析R×C列联表时，在卡方检验有统计意义的情况下常常需要做进一步的多重比较。可以采用的方法为1）卡方分割（具体见本人另外一篇文章《R×C行列表卡方值分割的概念及运用》）将原表卡方值分割成独立的子卡方值，分割后的子卡方值和对应的自由度相加会和原表的卡方值和自由度相等。2）或者采用彼此之间非独立的两两比较。但是两两比较的卡方值和对应的自由度相加不会等于原表卡方值和自由度，所以此类比较不能称为卡方分割法。尽管后者更为灵活但需要调整多重比较的次数以避免增加第一类错误。本文将具体讲解如何在SAS中

2、实现R×C列联表的两两比较。1．研究数据-血型和疾病类型假定某医学课题想研究血型(O,A,B)和疾病类型(PepticUlcer,GastricCancer,Control)是否相互关联,具体临床数据见表1。表1：血型和疾病类型行列表DiseaseTypeBloodTypePepticUlcer(I)-1GastricCancer(II)-2Control(III)-3TotalO-1983(14.024)383(4.9139)2892(1.4159)4258A-2679(9.0743)416(4.5484)2625(0.679)3720B-3134(4

3、.6663)84(0.2695)570(0.9519)788Totals17968836087N=876612．初步分析我们先用SAS/PROCFREQ和PROCCORRESP先对数据做初步的分析，来判断行列变量之间的关系。2CELLCHI2选项是计算按公式2-1计算每个单元格在总体值的组成。2(EO)ijijC（2-1）i,jEijprocfreqdata=paper14;weightcount;tabler*c/chisqcellchi2nopercent;run;proccorrespdata=paper14;weightcount;tabl

4、esr,c;run;表2–卡方检验结果TableofrbycrcFrequencyCellChi-SquareRowPctColPct123Total19833832892425814.0244.91391.415923.098.9967.9254.7343.3747.512679416262537209.07434.54840.67918.2511.1870.5637.8147.1143.123134845707884.66630.26950.951917.0110.6672.347.469.519.36Total179688360878766Stat

5、isticsforTableofrbycStatisticDFValueProb------------------------------------------------------Chi-Square440.5434<.0001LikelihoodRatioChi-Square440.6401<.0001Mantel-HaenszelChi-Square121.0035<.0001SampleSize=87662从表２中可以看出第１行（血型Ｏ）的卡方检验分值在所有行中所占比重最高（(14.024+4.9139+1.4159)/40.5434≈50

6、％）;第１列（PepticUlcer(I)）的卡方检验分值在所有列中所占比重最高（(14.024+9.0743+4.6663)/40.5434≈68.5%）。这表明第1行和第1列与其他行列相比有显著不同。同时，在图1中行1和列1很明显地远离其他行列。这说明2了表1中的行列表行变量（血型）和列变量（疾病类型）并非相互独立（Pearson检验的P值<0.0001）。那么如何运用统计检验来说明血型O和PepticUlcer(I)与其他组存在有统计意义的区别？下面讨论如何运用两两比较来回答这个问题。3.在SASPROCFREQ中实现两两比较具体做法是：1)从R

7、×C行列表中行变量中任取两不同的行和原行列表的列变量组成新的2×C子行列R表，共有次。分别计算各子表的2值。232)从R×C行列表中列变量中任取两不同的列和原行列表的行变量组成新的R×2子行列C表，共有次。分别计算各子表的2值。23)最后作根据比较次数做Bonferroniadjustment。a）如果只做行变量之间的比较，共RCRC次;b)如果只做列变量之间的比较，共次;c)如果都做比较，则+2222次；表1-血型和疾病类型数据的两两比较的SAS程序如下

8、：procfreqdata=paper14;whererin(1,2);weig