第8章 R×C表资料分析.pptVIP

中医药统计学;第一节一般四格表?2检验;第r行第c列位置上理论数的计算公式Trc＝nrnc/n ;H0：两种药物的疗效相同, 有效率345/376＝91.76% ;df＝ (2－1)×(2－1)＝1 ;解 H0：两组疗效相同；H1：两疗效不相同 ;按自由度df＝1查?2界值表，?20.01(1) ＝6.63，因?2＞ ?20.01(1) ，P＜0.01，按α＝0.01水准拒绝H0，认为两药治愈率不同 ;【SPSS操作】以组别、疗效及例数为变量名，建立数据文件L8-1.sav;data L27_20; do r=1 to 2; do c=1 to 2; input f@@; output; end; end; cards; 271 5 74 26 ; proc freq; weight f; tables r*c/chisq expected norow nocol nopct; run; ;匿陈过晒樱堰熙庚埂涌俐颓杂剖睦锅撤塌奉练荔失馒密究解嗜盅纵玲莱驮第8章 R×C表资料分析第8章 R×C表资料分析;【例8.2】甲组用活血温经汤，乙组用通塞脉1号治疗，问两药的疗效有无差别 ;查?2界值表，得P＞0.05，不拒绝H0，两种疗法效果无统计学差异，不能认为两法疗效不同 ;⑵ 用专有公式计算?2值，要先计算最小理论数。自由度df＝1时，u20.05/2＝?20.05(1)（即1.962 ＝3.84），u20.01/2＝?20.01(1)（即2.57582 ＝6.63），故n较大时的两样本率比较，既可用u检验，又可用?2检验 ;当①总样本含量n40时；或者② 用其它假设检验方法所得概率接近检验水准时；或者③ 四格表中有实际频数为0或有理论数T1时，都无法用?2检验，常用以超几何分布为前提的确切概率法（exact test） ;组别;【SPSS操作】四格表数据合为一列输入为变量x的值，变量m、n以值1、2分别代表数据所在的行或列。选择Data→Weight Cases命令，把x选入Frequency Variable框进行加权。 ;由n＝38＜40，宜用Fishers Exact Test。双侧确切概率P＝0.003，单侧确切概率P＝0.001，故拒绝H0，认为该药能够改善周围血管闭塞性病变患者的皮肤微循环状况 ;8.3 分层四格表分析;【DPS操作】分层四格表，输入数据以行属性表示分层（中心、疗效）、列属性表示不分层（分组）;选定数据块，选择分类数据统计菜单“分层四格表”命令，可得 ;由校正卡方＝3.4508，概率P＝0.0632＞0.05，表明去除中心混杂作用后，疗效与分组无关。在一致性检验中，概率P＝0.8462＞0.05，表明不同层间的OR值相同，去除中心混杂作用后，正常的危险度大约为异常的1.7136倍 ;双向无序R×C表，即行列皆无序的表。R×C表资料，在没有小于5的理论频数且总例数很大时，用?2检验法，可以推断样本率或构成比间的差异，也可推断行与列的属性间有无关系 ;【例8.5】某医院将用三种治疗方法治疗199例消化性溃疡住院病人，资料如表9-10。为避免中医不同证型对疗效比较的影响，试分析三种疗法治疗的病人按中医分型构成比有无差别 ;解 H0：三组的中医分型构成比相同；H1：三组的中医分型构成比不同 ;1. ?2检验要求理论数T不太小 一般认为R×C表中不宜有20%以上格子的T＜5，不能有T＜1。理论数太小时的处理办法：⑴最好增加样本例数以增大理论数。⑵删去理论数太小的行和列。⑶将太小理论数所在行或列的实际数与性质相近的邻行邻列合并，使重新计算的理论数增大。须知，后两法可能会损失信息，也会损害样本的随机性，不同的合并方式有可能影响推断结论，故不宜作为常规方法。⑷不满足?2检验条件的R×C表资料，最好计算确切概率 ;2. R×C表的确切概率法步骤 ⑴确定一个统计量，例如?2值 ⑵列出行合计与列合计固定的全部可能的R×C表，对于每一个表计算统计量和P值； ⑶凡统计量大于原始表的统计量值的那些R×C表的P值之和，即为确切概率P值，一般用统计软件完成 SPSS，Analyze→Dsecriptive statistic→Crosstabs，击Statistics按钮，指定Chi-square及Exact:Monte Carlo（确切概率：蒙特卡洛法） ;3. ?2检验拒绝H0，只能说明有差异，不能说明优劣 ;把第1，2列对换，计算结果完全一样，这时虽可得出治愈与药物剂量有关的结论，但不能推论治愈率随药物剂量增加而升高 ;4. 实验分组变量无序，指标变量有序时，宜用秩和或检验Ridit分析。但当实验分组变量有序，指标分组变量无序时，仍可用双向无序R×C表分析方法 ;在卡方检验有统计