医学统计学实习课：STATA软件的应用3)

1、医学统计学实习课STATA软件及其应用III课件下载账号：stata_密码：stata123456课件：stata3（4次）分类资料的统计分析分类变量的统计描述分类资料的统计推断总体率和事件数的可信区间估计两个率的假设检验卡方检验四格表卡方检验RC表卡方检验配对四格表卡方检验确切概率法分类变量的统计描述tabulate 变量1 变量2 , cell column row nofreqtab 分组变量 ， sum(数值变量)tab 分组变量1 分组变量2 ， sum(数值变量)cell /* 打印每个格子的频数占总频数的百分比column /* 打印每个格子的频数占相应列合计的百分比row /*

2、 打印每个格子的频数占相应行合计的百分比nofreq /* 不打印频数例题（ex1） 有三组(group)患者，男女(sex)若干人，sex=1表示男性，sex=0表示女性。测得其血红蛋白浓度(x,%) x group sex 1. 3.9 1 0 2. 4.2 1 0 3. 3.7 1 0 4. 4 1 0 5. 4.4 1 0 例题tab grouptab sextab sex grouptab group,sum(x)tab sex,sum(x)tab group sex,sum(x)tab groupgroup | Freq. Percent Cum.-+- 1 | 12 40.00

3、40.00 2 | 11 36.67 76.67 3 | 7 23.33 100.00-+- Total | 30 100.00tab sex sex | Freq. Percent Cum.-+- 0 | 15 50.00 50.00 1 | 15 50.00 100.00-+- Total | 30 100.00tab group sex | sex group | 0 1 | Total-+-+- 1 | 8 4 | 12 2 | 4 7 | 11 3 | 3 4 | 7 -+-+- Total | 15 15 | 30tab group,sum(x) | Summary of x gr

4、oup | Mean Std. Dev. Freq.-+- 1 | 3.8166667 .93889033 12 2 | 4.8727273 .52932203 11 3 | 3.7714286 .58513326 7-+- Total | 4.1933333 .88236879 30tab sex,sum(x) | Summary of x sex | Mean Std. Dev. Freq.-+- 0 | 4.1266667 .8224238 15 1 | 4.26 .9627342 15-+- Total | 4.1933333 .88236879 30tab group sex,sum

5、(x)Means, Standard Deviations and Frequencies of x | sex group | 0 1 | Total-+-+- 1 | 3.8125 3.825 | 3.8166667 | .90307009 1.1528949 | .93889033 | 8 4 | 12-+-+- 2 | 4.85 4.8857143 | 4.8727273 | .3872984 .6256425 | .52932203 | 4 7 | 11-+-+- 3 | 4 3.6 | 3.7714286 | .43588989 .68313003 | .58513326 | 3

6、4 | 7-+-+- Total | 4.1266667 4.26 | 4.1933333 | .8224238 .9627342 | .88236879 | 15 15 | 30率和事件计数的可信区间仅有原始资料，统计量未知时 ci 变量名统计量已知时cii 样本含量 样本率（发生数），选项level(#) /* 指定可信度，缺失时为 95(%)例题矫正25名学生的近视眼，3人有效，求总体率的95的可信区间；cii 25 3 样本量 发生数cii 25 3/25 样本量 发生率某医师治疗10名患者，治愈8人，求治愈率的90可信区间cii 10 8, level(90)四格表的卡方检验tabu

7、late var1 var2 fw=频数变量，选择项tabi a bc dchi2 /* (Pearson) 2检验lrchi2 /* 似然比2检验gamma /*Goodman-Kruskal的系数 taub /*Kendall的相关系数bV /* Cramer的列联系数Vall /* 同时给出以上结果exact /* Fisher的确切概率cell /* 打印每个格子的频数占总频数的百分比 column /* 打印每个格子的频数占相应列合计的百分比 row /* 打印每个格子的频数占相应行合计的百分比 expected /* 计算期望频数卡方检验方法的选择卡方检验：N40，T 5(stat

8、a有)校正的卡方检验： N40，1T 5确切概率法：N40或者T 1，无须关心理论频数(stata有)两药治疗贫血有效率的比较组别有效人数无效人数合计试验组431053对照组401656合计8326109Tmin=T12=53*26/109=12.64数据录入格式：rcf1143121021402216例题(ex2)tab r c fw=f,row all | c r | 1 2 | Total-+-+- 1 | 43 10 | 53 | 81.13 18.87 | 100.00 -+-+- 2 | 40 16 | 56 | 71.43 28.57 | 100.00 -+-+- Total |

9、 83 26 | 109 | 76.15 23.85 | 100.00 Pearson chi2(1) = 1.4115 Pr = 0.235 likelihood-ratio chi2(1) = 1.4231 Pr = 0.233 Cramers V = 0.1138 gamma = 0.2647 ASE = 0.213 Kendalls tau-b = 0.1138 ASE = 0.094tabi 43 1040 16,row chi2 exact | col row | 1 2 | Total-+-+- 1 | 43 10 | 53 | 81.13 18.87 | 100.00 -+-+

10、- 2 | 40 16 | 56 | 71.43 28.57 | 100.00 -+-+- Total | 83 26 | 109 | 76.15 23.85 | 100.00 Pearson chi2(1) = 1.4115 Pr = 0.235 Fishers exact = 0.267 1-sided Fishers exact = 0.168mcci #a #b #c #d配对四个表的应用条件：1，b+c40时，使用该公式：2，20b+c40时，使用该公式：3，b+c20，使用确切概率法：配对四格表甲种培养法+-乙种培养法 +3624-10135配对四格表 mcci 36 24 10

11、135有205份咽喉涂抹样本，把每份标本依同样的条件分别接种于甲、乙两种白喉杆菌培养基上，观察白喉杆菌生长的情况，结果如下，问两种培养基的阳性率有无差别？R*C表的卡方检验tabulate var1 var2 fw=频数变量 ，选择项tabi多个率的比较 多组构成比的比较 多个率的比较分 组 阳性例数 阴性例数 合 计阳性率(%)a=11881620492.3a=21023 3330.3a=34933338212.8合 计24737261939.9 tabi 188 1610 2349 333 , row chi2 lrchi2多组构成比的比较 A B O AB 合计 (xx=1) (xx=2

12、) (xx=3) (xx=4) 傣 族 (mz=1) 112 150 205 40 507佤 族 (mz=2) 200 112 135 73 520土家族(mz=3) 362 219 310 69 960合 计 674 481 650 182 1987tabi 112 150 205 40200 112 135 73362 219 310 69, row chi2 lrchi2 无效有效合计甲法617乙法112031合计172138tabi 6 111 20,row all exact tab r c fw=freq ,row all exact确切概率法练习1某医生用甲乙两法检测同一组患者,具体结果如下,试比较这两种方法的阳性率是否相同.甲法 乙法+-合计+351853-533285合计8850138练习2比较男性和女性暴露率是否相同。暴露未暴露合计男性13518女性21728合计341246知识点考察形式例：某医院用A，B两种仪器测了16名健康男青年的血红蛋白含量（g/L），检测结果见下表。问两种血红蛋白测定仪器的检测结果是否有差别。1234567891011121314151611312512613015014513510512813510013011011512015514015013