研究生统计学讲义第10讲中医临床实验设计与统计分析.ppt

《研究生统计学讲义第10讲中医临床实验设计与统计分析.ppt》由会员分享，可在线阅读，更多相关《研究生统计学讲义第10讲中医临床实验设计与统计分析.ppt（48页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、中医临床研究设计及统计分析中医临床研究设计及统计分析第一章第一章制定中医证候宏观辩证量表制定中医证候宏观辩证量表:统一辩证标准统一辩证标准第一节第一节现时中医证候宏观标准产生的背景现时中医证候宏观标准产生的背景1.古代医家论述、医案记载；古代医家论述、医案记载；2.现代教材描述现代教材描述;3.学术机构的制定学术机构的制定;4.政府有关部门组织专家制定政府有关部门组织专家制定;5.来源于一定范围的临床流行病学研究来源于一定范围的临床流行病学研究;就第就第2、3、4项与第项与第5项结合项结合,阐述其制定方法阐述其制定方法.第二节第二节中医证候宏观辩证量表制定方法中医证候宏观辩证量表制定方法1.利

2、用现有四诊项目组合法利用现有四诊项目组合法+临床流行病学调查证候加临床流行病学调查证候加权积分法权积分法:国际或行标规定的主要症状都具备国际或行标规定的主要症状都具备(如果无主要症如果无主要症状状,则要求任意三项症状则要求任意三项症状)+舌苔舌苔+脉象脉象.合计至少合计至少5项阳项阳性才能辩为某证候性才能辩为某证候.根据此项规定可以计算出每个证根据此项规定可以计算出每个证候候诊断的最低加权积分临界值诊断的最低加权积分临界值.优点优点:宏观辩证统一宏观辩证统一化、规范化、量化。化、规范化、量化。现场调查证候现场调查证候加权积分步骤：加权积分步骤：每项症状按严重程度或频次制定每项症状按严重程度或频

3、次制定0、2、4、6分的得分的得分量表分量表.分别对应无分别对应无、轻、中、重症状；使用文字进、轻、中、重症状；使用文字进行描述得分具体内容行描述得分具体内容.是症状和体征是症状和体征(四诊四诊)的定量方法的定量方法.舌和脉象如果不容易打分舌和脉象如果不容易打分,也可按照无也可按照无=0分分,有有=1分分.证候加权积分证候加权积分(计量资料计量资料)=每项四诊得分每项四诊得分*每项四诊每项四诊对证候的权重对证候的权重.(=同一证候中每项四诊中出现频率的百同一证候中每项四诊中出现频率的百分构成比分构成比,可根据文献或现场调查计算可根据文献或现场调查计算)证候加权积分作为证候疗效评价的定量依据证候

4、加权积分作为证候疗效评价的定量依据.若再将证候加权积分转化为证候分级值，若再将证候加权积分转化为证候分级值，则可以把则可以把证候分级为证候分级为0、1、2、3级级,分别代表某证候的无、轻、分别代表某证候的无、轻、中、重中、重.2.利用古今传统文献统计法利用古今传统文献统计法:因为四诊资料和证候在原因为四诊资料和证候在原始文献始文献(医案医案)中没有定量值中没有定量值,只能查到有或无的记载只能查到有或无的记载.因此可以统计各种证候的每项四诊资料发生的频数因此可以统计各种证候的每项四诊资料发生的频数,换算成为频率换算成为频率(条件概率条件概率),再用对数转换为诊断概率再用对数转换为诊断概率指数指数

5、,从而制定出计量诊断与鉴别诊断表来从而制定出计量诊断与鉴别诊断表来.3.只按照现代流行病学调查和统计分析只按照现代流行病学调查和统计分析特点特点:不受现有证候诊断不受现有证候诊断(宏观辩证宏观辩证)标准的限制标准的限制,完全完全按照无偏的自然人群样本调查收集到的四诊资料进行按照无偏的自然人群样本调查收集到的四诊资料进行团结分析团结分析.优点优点:不仅可以筛选出某证主要四诊资料不仅可以筛选出某证主要四诊资料(指标指标),而且而且可以得到每个指标对证候的权重可以得到每个指标对证候的权重.例如主成分及因子分例如主成分及因子分析和指标聚类分析析和指标聚类分析.问题问题:对证候的重新分类可能与现有辩证体

6、系不全吻合对证候的重新分类可能与现有辩证体系不全吻合上述第上述第1种方法可做暂行种方法可做暂行“金金”标准标准,其最大优点是其最大优点是:1.辩证统一辩证统一;2.辩证结果可以重复辩证结果可以重复,即重测信度即重测信度(retestreliability)好好.3.症状和证候可定量症状和证候可定量,均可作为症状和证候疗效评价均可作为症状和证候疗效评价的依据的依据.第三节第三节制定方法举例制定方法举例以中医内科学以中医内科学(张伯臾张伯臾)对胃痛的辩证分型和四诊对胃痛的辩证分型和四诊资料项目作为基础资料项目作为基础.再按照四诊资料的严重程度或发再按照四诊资料的严重程度或发生频次评定得分生频次评定

7、得分,拟定四诊资料表拟定四诊资料表,称为四诊资料得称为四诊资料得分分证候加权积分证候加权积分=各项四诊资料得分各项四诊资料得分*各项四诊项目权重的各项四诊项目权重的积和积和(代数和代数和).加权方法加权方法:某一证候中某一证候中,每一项四诊资料出现频率的每一项四诊资料出现频率的百分构成比百分构成比(根据临床流行病学调查数据计算根据临床流行病学调查数据计算).上述四诊资料得分为等级资料上述四诊资料得分为等级资料,证候加权积分是计量证候加权积分是计量资料资料.表表1胃痛各证候的宏观辩证量表胃痛各证候的宏观辩证量表(金标准金标准)症状症状|得分得分0分分2分分4分分6分分权重权重1.寒邪客胃寒邪客胃

8、1.1胃脘疼痛胃脘疼痛无无轻微轻微较痛能忍较痛能忍不能忍不能忍0.31.2得温痛减得温痛减从不从不偶尔偶尔经常经常总是总是0.3 遇寒痛增遇寒痛增1.3口和不渴口和不渴非常渴非常渴很渴很渴有点渴有点渴不渴不渴0.31.4喜热饮喜热饮从不从不偶尔偶尔经常经常总是总是0.11.5舌苔白舌苔白不是不是是是0.11.6脉弦紧脉弦紧不是不是是是0.12.饮食停滞饮食停滞.3.肝气犯胃肝气犯胃4.肝胃郁热肝胃郁热5.瘀血停滞瘀血停滞6.胃阴亏虚胃阴亏虚7.脾胃虚寒脾胃虚寒.参考文献参考文献:中医内科学张伯臾主编中医内科学张伯臾主编,第第1版版上海科上海科技出版社技出版社,1986年年临床流行病学王家良主编

9、临床流行病学王家良主编,第第2版版,上海科技出版上海科技出版社社,2001,437页页临床流行病学在中医药研究中的临床流行病学在中医药研究中的应用应用:赖世隆赖世隆第二章第二章中医证候宏观辩证量表中医证候宏观辩证量表的表面效度评价的表面效度评价第一节第一节Kappa检验的原理检验的原理在在诊诊断断检检验验中中,对对水水平平数数3的的配配伍伍设设计计问问题题.通通常常进进行行一一 致致 性性(intraobserver aggrement)检检 验验,可可 使使 用用McNemar检验和检验和Kappa检验检验.1.Kappa检验的意义检验的意义Kappa检验（检验（Kappatest）或称一）

10、或称一致性致性(intraobserveragreement)检验，包括检验，包括同一医同一医药工作者对同一受试对象进行两次以上观察作出判断药工作者对同一受试对象进行两次以上观察作出判断的一致性。例如：同一临床医师对同一患者的两次胃的一致性。例如：同一临床医师对同一患者的两次胃镜检查的结论、两次阅读同一镜检查的结论、两次阅读同一X光片的诊断结论，病理光片的诊断结论，病理学医师两次阅读同一张切片的诊断结论，同一医药工学医师两次阅读同一张切片的诊断结论，同一医药工作者两次测定同一药材某特征物含量的结论，都可能作者两次测定同一药材某特征物含量的结论，都可能产生不一致的情况。产生不一致的情况。两个及两

11、个以上的医药工作者两个及两个以上的医药工作者对同一对象进行观察作出判断的一致性。对同一对象进行观察作出判断的一致性。例如：两个或更多的医药工作者检查同一患者、阅读例如：两个或更多的医药工作者检查同一患者、阅读同一同一X光片、同一病理切片也可能发生判断的不一致。光片、同一病理切片也可能发生判断的不一致。因此在临床实际工作中，特别是临床试验研究中，有因此在临床实际工作中，特别是临床试验研究中，有必要对判断的一致性作出估计。例如，在利用胃镜检必要对判断的一致性作出估计。例如，在利用胃镜检查结果作慢性胃炎的诊断之前，研究者很有必要了解查结果作慢性胃炎的诊断之前，研究者很有必要了解重复检查的一致程度，若

12、一致性较差，说明诊断结果重复检查的一致程度，若一致性较差，说明诊断结果的可靠性受到重复检查结果判断不一致的影响；若一的可靠性受到重复检查结果判断不一致的影响；若一致性较好，说明重复检查结果是可信的。由于实际工致性较好，说明重复检查结果是可信的。由于实际工作中重复观察造成的判断不一致很难避免，只能尽量作中重复观察造成的判断不一致很难避免，只能尽量采取措施控制在最低限度，同时，选用适当的统计方采取措施控制在最低限度，同时，选用适当的统计方法评价其判断的一致性。法评价其判断的一致性。2Kappa检验的统计量检验的统计量为为Kappa值，简记为值，简记为K，计算，计算公式为：公式为：式中式中PA为实际

13、观察到的一致率；为实际观察到的一致率；Pe为期望一致率，简为期望一致率，简称期望率，即两次检查结果由于偶然机会所造成的一称期望率，即两次检查结果由于偶然机会所造成的一致率。致率。PA和和Pe的计算公式分别为：的计算公式分别为：可可见见，K值值实实际际为为两两个个差差值值之之比比，分分母母为为(1期期望望率率)，即即非非偶偶然然机机会会造造成成的的一一致致率率；分分子子称称为为校校正正后后的的一一致致率率，它它是是实实际际观观察察到到的的一一致致率率和和可可能能由由于于偶偶然然机机会会造造成成的的期期望望率率之之差差值值，差差值值越越大大，说说明明实实际际观观察察到到的的一一致致率率远远比比由由

14、于于机机会会造造成成的的期期望望一一致致率率高高。K值值为为经经校校正后一致率与非机会造成一致率的比值。正后一致率与非机会造成一致率的比值。配对四格表资料的期望率配对四格表资料的期望率Pe、实际观察一致率、实际观察一致率PA分别分别为：为：Pe=(a+b)(a+c)/n+(b+d)(c+d)/n/n=(a+b)(a+c)+(b+d)(c+d)/n2(10.15)PA=(a+d)/n(10.14)配对四格表资料的配对四格表资料的Kappa值可简化为式值可简化为式(10.16)：K=2(adbc)/(a+b)(a+c)+(b+d)(c+d)(10.16Kappa的取值范围在的取值范围在-1至至+1

15、之间，之间，K值愈接近值愈接近+1，表明，表明一致程度愈好，一致程度愈好，K值值+1，说明两次判断的结果完全一，说明两次判断的结果完全一致；致；K值值1，表示二者完全相反，说明两次判断的，表示二者完全相反，说明两次判断的结果完全不一致；若结果完全不一致；若K值值0，表示一致率完全由机会，表示一致率完全由机会引起，说明两次检查结果很不一致。在实际运用中，只引起，说明两次检查结果很不一致。在实际运用中，只有有K值在值在0到到+1区间时，对判断一致性才有意义。区间时，对判断一致性才有意义。K值越值越大说明二者的一致性越好。一般认为，若大说明二者的一致性越好。一般认为，若K值值0.75，可，可说明已经

16、达到相当满意的一致程度；若小于说明已经达到相当满意的一致程度；若小于0.4，则可，则可说明一致程度不够理想。说明一致程度不够理想。Landis和和Koch对四格表资料对四格表资料提出判断一致性程度的提出判断一致性程度的Kappa量判断表，见表量判断表，见表10-15。表表10-15一致性程度的一致性程度的Kappa量判断表量判断表*Kappa值值一致性程度一致性程度Kappa值值一致性程一致性程度度0弱弱0.410.60中度中度0.000.20轻轻0.610.80高度高度0.210.40尚好尚好0.811.00极强极强*Kappa值值0.4，Kappa值的标值的标准误准误=0.049，近似概率

17、，近似概率P=0.000。说明样本。说明样本Kappa值来值来自总体自总体Kappa值为值为0的概率近似的概率近似0.000，故可认为两者，故可认为两者辩证结果的一致性较好。辩证结果的一致性较好。第三节第三节定性观察定性观察(二分类变量二分类变量)的一致性评价方法的一致性评价方法如果将证候分度简化为二分类变量如果将证候分度简化为二分类变量,则资料整理成表则资料整理成表乙辩证例数乙辩证例数是是非非合计合计甲辩证例数甲辩证例数 是是a b a+b非非 c d c+d 合计合计 a+c b+d n例例10.14湖南中医学院研究的血吸虫病肝纤维化气虚血湖南中医学院研究的血吸虫病肝纤维化气虚血瘀证判别模

18、型回代检验结果见表瘀证判别模型回代检验结果见表10-17，用，用Kappa值评值评价该判别模型诊断的吻合性与一致性。价该判别模型诊断的吻合性与一致性。乙辩证例数乙辩证例数是是非非合计合计甲辩证例数甲辩证例数 是是36541非非 34043 合计合计 394584H0：两者吻合。：两者吻合。H1：两者不吻合。：两者不吻合。2(xji-xij)2/(xji+xij)=(5-3)2/(5+3)=0.5，自由度自由度df=k(k-1)/2=2(2-1)/21,查查 2界值表得界值表得P0.05K=2(adbc)/(a+b)(a+c)+(b+d)(c+d)2(3640-53)/(4145+3943)=0

19、.81，按表，按表10-15的一致性程度的的一致性程度的Kappa量判断表，可认为该判别模型诊断的吻合性极量判断表，可认为该判别模型诊断的吻合性极强。强。第三章第三章中医证候宏观辩证量表的重测信度评价中医证候宏观辩证量表的重测信度评价(R值法值法)及反应度评价及反应度评价信度又称为可靠性信度又称为可靠性(reliability):指在相同条件指在相同条件下下,对同一客观事件重复测量若干次对同一客观事件重复测量若干次,测量结果相互测量结果相互符合的程度符合的程度;说明数据的可靠性说明数据的可靠性.常用信度评价的方法有四种常用信度评价的方法有四种:1.客观信度客观信度(objectivityrel

20、iability)2.精确程度精确程度(precisionreliability)3.和谐程度和谐程度(congruencereliablity)4.重测信度重测信度(re-testreliability):用同样测量方法对同用同样测量方法对同一批对象在不同时间段进行重复测量一批对象在不同时间段进行重复测量,目的是考核测目的是考核测量工具的可靠性量工具的可靠性.下面介绍使用组内相关可靠性系数下面介绍使用组内相关可靠性系数(intra-classcorrelationcoeficientofreliability,R值值)对中医证候对中医证候宏观辩证量表进行重测信度评价的方法和反应度评价宏观辩证

21、量表进行重测信度评价的方法和反应度评价的方法的方法.第一节第一节R的设计计算和统计分析的设计计算和统计分析一、例一、例1某研究者使用中医证候宏观辩证量表对同一某研究者使用中医证候宏观辩证量表对同一批受测者批受测者(10例例)在上在上,下午两次盲法观测证候的直接下午两次盲法观测证候的直接积分积分(四诊资料总得分四诊资料总得分),资料见下表资料见下表:受试编号受试编号12345678910第第1次积分次积分545312791586第第2次积分次积分644410581878 二、计算组内相关可靠性系数二、计算组内相关可靠性系数:可使用可使用SPSS统计分析统计分析0.7R0.4,可靠性较好可靠性较好

22、,R0.7可靠性甚佳可靠性甚佳.这里这里R=0.92可靠性好可靠性好.第二节第二节量表反应度的评价量表反应度的评价量表不但要评价其真实性量表不但要评价其真实性(效度效度)和可靠度和可靠度,还要还要评价反应度评价反应度.一、反应度的概念一、反应度的概念反应度反应度(responsibilitytochange):反应内外环境和反应内外环境和受试者有一定变化受试者有一定变化(例如中医四诊资料的严重度或频度例如中医四诊资料的严重度或频度有所变化有所变化)的情况下该量表的应变能力的情况下该量表的应变能力.二、二、反应度指标反应度指标1.实施干预措施后实施干预措施后,量表两次得分的差别有无统计学量表两次

23、得分的差别有无统计学意义的假设检验意义的假设检验:既可以按照计量资料进行假设检验既可以按照计量资料进行假设检验,如如t检验检验(正态性正态性,方差齐性方差齐性)或或t检验检验(方差不齐方差不齐).也可以进行非参数检验也可以进行非参数检验.2.效应尺度统计量效应尺度统计量(effectsizestatistics):除开假设检除开假设检验外验外,还可以把实施干预措施前后还可以把实施干预措施前后,量表两次得分转量表两次得分转换为换为“效应尺度效应尺度”统计量来评价量表的反应度统计量来评价量表的反应度.效应尺度效应尺度=(干预得分干预得分-干预后得分干预后得分)/干预前得分的标准差干预前得分的标准差

24、一般认为一般认为,效应尺度小于效应尺度小于0.2,则反应度甚差则反应度甚差;0.20.50之间反应度较小之间反应度较小;0.510.80为反应度较大为反应度较大.如果如果量表在干预前后两次得分的假设检验有统计学意义的量表在干预前后两次得分的假设检验有统计学意义的差别差别,而效应尺度却小于而效应尺度却小于0.20,则显示该量表的反应度则显示该量表的反应度甚差甚差.第四章病证结合的观察性描述研究：横断面调查分析第四章病证结合的观察性描述研究：横断面调查分析观察性描述性研究主要包括横断面调查、诊断试观察性描述性研究主要包括横断面调查、诊断试验评价与筛查、生态学研究、临床医疗调查报告、病验评价与筛查、

25、生态学研究、临床医疗调查报告、病例总结分析，专家述评等，属于观察性描述性研究。例总结分析，专家述评等，属于观察性描述性研究。其科学价值可参考加拿大其科学价值可参考加拿大MCMaster大学的四级分类方大学的四级分类方法，属于第三级（较低一级）研究设计方案，论证强法，属于第三级（较低一级）研究设计方案，论证强度较弱。度较弱。第一节横断面调查的概念和特点第一节横断面调查的概念和特点横断面调查又称现况调查（横断面调查又称现况调查（cross-sectionalsurvey）或患病率（现患病，或患病率（现患病，prevalencerate）调查。是在某个）调查。是在某个时点或较短时段内（看成一个时点，

26、一般在一个月内时点或较短时段内（看成一个时点，一般在一个月内完成调查），调查一个特定人群样本中医证候的分布完成调查），调查一个特定人群样本中医证候的分布频率以及与某些因素的相关关系，可以得到证候的现频率以及与某些因素的相关关系，可以得到证候的现患率（患率（prevalencerate）和一些可疑的影响因素。）和一些可疑的影响因素。现况调查可以将证候按人群各种特征间、时间、地现况调查可以将证候按人群各种特征间、时间、地区间（合称区间（合称“三间三间”）的横断面分布频率显示出来，）的横断面分布频率显示出来，而且可以获得影响证候而且可以获得影响证候“三间三间”分布频率的相关因素分布频率的相关因素信息

27、。信息。由于证候与相关因素是在同一次调查中得到的，即由于证候与相关因素是在同一次调查中得到的，即果与因同时存在，无法判断发生时间的先后顺序，除非果与因同时存在，无法判断发生时间的先后顺序，除非一些持续不变的因素如家族史、性别、血型、基因型外，一些持续不变的因素如家族史、性别、血型、基因型外，不能推断证候与相关因素之间存在因果联系，即不能确不能推断证候与相关因素之间存在因果联系，即不能确定病因或危险因子。定病因或危险因子。尽管不能进行病因推断，但却是病因研究必须的第尽管不能进行病因推断，但却是病因研究必须的第一步，能为病或证候的病因（或危险因子）研究提供工一步，能为病或证候的病因（或危险因子）研

28、究提供工作假设的线索，为第二步观察性分析性研究（包括回顾作假设的线索，为第二步观察性分析性研究（包括回顾性病例对照研究和前瞻性的队列研究）以及第三步的干性病例对照研究和前瞻性的队列研究）以及第三步的干预研究（包括治疗性临床试验、预防性临床试验、以群预研究（包括治疗性临床试验、预防性临床试验、以群体为单位的社区试验）提供基础。体为单位的社区试验）提供基础。第二节第二节横断面调查的目的和用途横断面调查的目的和用途1.证候按人群各种特征的分布频率：年龄、性别、种证候按人群各种特征的分布频率：年龄、性别、种族、基因型、经济收入、职业、婚姻、家庭状况等一切族、基因型、经济收入、职业、婚姻、家庭状况等一切

29、可能影响证候的分布频率。可能影响证候的分布频率。一、了解证候人群人中的现况和描述证候按一、了解证候人群人中的现况和描述证候按“三间三间”的分布频率的分布频率3.证候按地区的分布频率：不同地理区域可能影响证证候按地区的分布频率：不同地理区域可能影响证候的分布频率。候的分布频率。二、了解影响证候二、了解影响证候“三间三间”分布频率的相关因素分布频率的相关因素2.证候按时间的分布频率：人群抽样调查时间的不同证候按时间的分布频率：人群抽样调查时间的不同（例如季节和年份）可能影响证候的分布频率。（例如季节和年份）可能影响证候的分布频率。1.证候按人群各种特征的分布频率：年龄、性别、种证候按人群各种特征的

30、分布频率：年龄、性别、种族、基因型、经济收入、职业、婚姻、家庭状况等一族、基因型、经济收入、职业、婚姻、家庭状况等一切可能影响证候的分布频率。切可能影响证候的分布频率。2.证候按时间的分布频率：人群抽样调查时间的不同证候按时间的分布频率：人群抽样调查时间的不同（例如季节和年份）可能影响证候的分布频率。（例如季节和年份）可能影响证候的分布频率。3.证候按地区的分布频率：不同地理区域可能影响证证候按地区的分布频率：不同地理区域可能影响证候的分布频率。候的分布频率。三、了解影响证候三、了解影响证候“三间三间”分布频率的相关因素分布频率的相关因素如果是病证结合研究，影响病的患病率的相关因素，如果是病证

31、结合研究，影响病的患病率的相关因素，也可能是证候影响因素。例如：软水质使冠心病高发，也可能是证候影响因素。例如：软水质使冠心病高发，硬水质使龋齿高发，黄曲霉素硬水质使龋齿高发，黄曲霉素B在食品中的浓度与肝癌在食品中的浓度与肝癌患病率呈正相关，吸烟与支气管肺癌，某些地区的民患病率呈正相关，吸烟与支气管肺癌，某些地区的民族某种病低发和高发，即使移民到某国几代人后仍与族某种病低发和高发，即使移民到某国几代人后仍与该国居民的患病率（或发病率）水平有显著差别。所该国居民的患病率（或发病率）水平有显著差别。所有这些相关因素都可能影响证候的分布频率，因为证有这些相关因素都可能影响证候的分布频率，因为证候与病

32、都是反映疾病和健康状况的概念。候与病都是反映疾病和健康状况的概念。四、了解证候分布频率的动态变化四、了解证候分布频率的动态变化如果能像全球疾病监测那样，在监测点同样可以做如果能像全球疾病监测那样，在监测点同样可以做到证候的监测，在同监测的相同人群样本中到证候的监测，在同监测的相同人群样本中2次以上的次以上的调查，就可以提供时间动态变化。调查，就可以提供时间动态变化。第三节第三节人群抽样调查设计及随机抽样的人群抽样调查设计及随机抽样的SAS程序程序一、目标明确一、目标明确1、主要了解某病的证候分布频率：使用医院确诊为某、主要了解某病的证候分布频率：使用医院确诊为某病的连续病例，直到满足样本含量为

33、止。医院病例的病的连续病例，直到满足样本含量为止。医院病例的缺点：不能计算病或证候的患病率（时点患病率某缺点：不能计算病或证候的患病率（时点患病率某时点新老病例数之和时点新老病例数之和/该时点平均人口数），只能计算该时点平均人口数），只能计算某病的某证候分布某病的某证候分布频率。因为医院的病例不能包含各频率。因为医院的病例不能包含各年龄组、各职业、各病期、种种病情（轻、中、重、年龄组、各职业、各病期、种种病情（轻、中、重、并发症）的全部病例，所以医院病例是该病的有偏样并发症）的全部病例，所以医院病例是该病的有偏样本。理想的人群样本来自自然人群（人群各种特征的本。理想的人群样本来自自然人群（人群

34、各种特征的比例符合所在的市或县最近一次人口普查的特征比例）。比例符合所在的市或县最近一次人口普查的特征比例）。2、了解某病及证候的患病率及其影响因子（相关因素）、了解某病及证候的患病率及其影响因子（相关因素）：来自自然人群的样本。城镇的居民委员会或农村的行：来自自然人群的样本。城镇的居民委员会或农村的行政村均为自然人群。政村均为自然人群。二、统一诊断标准二、统一诊断标准1、病的诊断标准：国际标准或国家标准。、病的诊断标准：国际标准或国家标准。2、证候的辩证标准：笔者建议制定、证候的辩证标准：笔者建议制定“统一的辩证用四统一的辩证用四诊量表诊量表”作为目前阶段统一的辩证标准，即暂行作为目前阶段统

35、一的辩证标准，即暂行“金金”标准。制定方法可参阅本书第一章。标准。制定方法可参阅本书第一章。三、调查表三、调查表应包含一般项目和研究项目。应包含一般项目和研究项目。一般项目：姓名；性别；年龄组（按世界卫生组织提一般项目：姓名；性别；年龄组（按世界卫生组织提供的供的“标准人口构成标准人口构成”构成年龄分组方法，成年人构成年龄分组方法，成年人5岁岁一个年龄组，一个年龄组，75岁以上只一个年龄组，例如岁以上只一个年龄组，例如20253075）；职业。）；职业。研研究究项项目目：凡凡是是可可疑疑的的病病或或证证的的一一切切相相关关因因素素，而而且且打算统计分析的，都不应遗漏；打算统计分析的，都不应遗漏

36、；不打算统计分析的，一项也不要列入。不打算统计分析的，一项也不要列入。各各种种项项目目都都要要写写明明计计算算机机编编码码：三三分分类类名名义义变变量量如如男男和和女女，编编码码为为1和和2；多多类类无无序序名名义义变变量量（如如某某病病的的证证1、证证2、证证3）以以及及多多类类有有序序变变量量（等等级级变变量量），编编码码为为1、2、3K；连连续续型型变变量量（计计算算资资料料）填填写写实际观测值。实际观测值。研研究究项项目目应应尽尽量量选选择择客客观观的的测测量量指指标标，即即硬硬数数据据（harddata），尽尽量量减减少少主主观观的的重重复复性性差差的的不不能能准准确确度度量量的的指

37、指标标，即即软软数数据据（softdata）。必必不不可可少少的的软软数数据据（如如症症状状、心心理理学学指指标标等等）收收集集，为为了了统统一一，宜宜做做成成量量表表（scale），要要做做效效度度（真真实实性性、准准确确度度，即即validity）和和信信度度（可可靠靠性性、可可重重复复性性、精精确确度度，即即reliability）检验，可参阅本书有关章节。）检验，可参阅本书有关章节。四、抽样方法及四、抽样方法及SAS统计程序统计程序无论采用何种抽样方法，必须遵守随机抽样的原则。无论采用何种抽样方法，必须遵守随机抽样的原则。因为只有随机样本才能代表总体（自然人群中所有的因为只有随机样本才

38、能代表总体（自然人群中所有的某病某证候者的总体）的各种特征。而且要计算抽样某病某证候者的总体）的各种特征。而且要计算抽样误差（误差（samplingerror），以便对总体参数作出估计，），以便对总体参数作出估计，并对抽样调查结果的真实性（与总体真值的距离，即并对抽样调查结果的真实性（与总体真值的距离，即准确度或效度）和可靠性（可重复性，即精确度或信准确度或效度）和可靠性（可重复性，即精确度或信度）进行评价。度）进行评价。1、单纯随机抽样（、单纯随机抽样（SimplerandomSampling）是是简简单单又又基基本本的的抽抽样样方方法法，从从N个个人人的的总总体体中中随随机机抽抽出出n个个

39、构构成成一一个个含含量量为为n人人的的样样本本。这这里里，基基本本抽抽样单位是样单位是“人人”。笔者建议使用。笔者建议使用SAS程序进行随机抽样程序进行随机抽样单纯随机抽样单纯随机抽样SAS程序举例：程序举例：设设计计例例题题1现现有有有有限限总总体体2000人人，欲欲从从中中随随机机抽抽出出50人人的样本的样本，怎样改写，怎样改写SAS程序？程序？dataebm3.samp;doi=1to50;x=ranuni(20020808);y=int(x*2000);output;end;run;运行结果如下：运行结果如下：OBSIXY110.09759195.答答:欲抽出欲抽出50人，第人，第3行

40、行doi1to50；如欲抽出；如欲抽出80人，人，则把则把50改为改为80；从；从2000人的有限总体中抽样，所以人的有限总体中抽样，所以第第5行行y=int(x*2000)；即随机数先乘；即随机数先乘2,000再取整。如再取整。如有限总体为有限总体为20000人，则把人，则把2000改为改为20000即可。第即可。第4行括号内是种子数，相当于抽样起始点，本例取日期，行括号内是种子数，相当于抽样起始点，本例取日期，读者也可取其它任意正整数。当然首先给有限总体的读者也可取其它任意正整数。当然首先给有限总体的每个抽样单位编号。每个抽样单位编号。上面运行结果的第上面运行结果的第4列（列（Y列）就是在

41、列）就是在2000人的有限人的有限总体内随机抽中的总体内随机抽中的50人随机样本号，它是通过人随机样本号，它是通过x*2000得到的，得到的，x也是随机数。也是随机数。50500.53951107949490.161013222、抽样误差的计算公式：、抽样误差的计算公式：（1）无限总体抽样误差：）无限总体抽样误差：当（当（n/N）5%时，时，N可看作无限总体。例如本例，可看作无限总体。例如本例，n/N=50/200=2.5%，看作无限总体。均数的标准误：，看作无限总体。均数的标准误：率的标准误：率的标准误：式中式中S=样本标准差，样本均数的标准误，样本标准差，样本均数的标准误，p=样本样本率，

42、率，n=样本含量。样本含量。如果如果n/N5%，N看作有限总体。上述看作有限总体。上述2个公式分别乘个公式分别乘上校正数（上校正数（1n/N）的平方根。）的平方根。(2)有限总体的抽样误差：有限总体的抽样误差：均数的标准误：均数的标准误：率的标准误：率的标准误：3、整群抽样（、整群抽样（ClusterSampling）基本抽样单位不是以基本抽样单位不是以“人人”为单位，而是以为单位，而是以“群组群组”为单位，例如城镇的居委会，农村的行政村、学校的为单位，例如城镇的居委会，农村的行政村、学校的班、部队的班、工人的车间班、部队的班、工人的车间等。等。设设计计例例题题2欲欲调调查查北北京京市市200

43、0年年8月月城城区区常常驻驻户户口口居居民民II型糖尿病患病率及型糖尿病患病率及II型糖尿病气阴两虚证的现患率。型糖尿病气阴两虚证的现患率。问：怎样进行整群抽样？问：怎样进行整群抽样？答答：基基本本抽抽样样单单位位是是居居委委会会。按按照照市市公公安安局局户户口口处处的的资资料料，先先把把各各城城区区各各街街道道办办事事处处各各居居委委会会的的名名称称、人人口口列列表表，并并且且按按某某种种人人为为的的顺顺序序对对居居委委会会编编号号，假假定定1到到2000号号。再再计计算算样样本本量量，得得出出平平均均需需要要多多少少个个居居委委会会作作为为样样本本，（样样本本量量计计算算方方法法见见本本章

44、章第第五五节节），假假定定需需要要抽抽出出50个个居居委委会会。最最后后仍仍用用SAS程程序序（数数据据集集名名：ebm3.samp）随随机机抽抽出出50个个居居委委会会。被被抽抽中中的的居居委委会会全全体居民都是调查对象。体居民都是调查对象。整群抽样的抽样误差大于单纯随机抽样，需要增加样整群抽样的抽样误差大于单纯随机抽样，需要增加样本量，一般增加本量，一般增加50左右即可。左右即可。4、分层抽样（、分层抽样（Stratifiedsampling）先按总体人口学特征（年龄组、性别、种族、职业、经先按总体人口学特征（年龄组、性别、种族、职业、经济收水平等）或疾病的病情（轻、中、重、并发症等）济收

45、水平等）或疾病的病情（轻、中、重、并发症等）分成若干层次，并在每层内作单纯随机抽样。一般说，分成若干层次，并在每层内作单纯随机抽样。一般说，当样本会计师足够大时，调查患病率和相关因素，只要当样本会计师足够大时，调查患病率和相关因素，只要单纯随机抽样即可，不一定需要事先分层。但是调查结单纯随机抽样即可，不一定需要事先分层。但是调查结束后处理资料时，分层统计分析是必要的，这是为了控束后处理资料时，分层统计分析是必要的，这是为了控制混杂偏倚（制混杂偏倚（confoundingbias）。）。第四节第四节横断面调查需要的样本含量估计方法横断面调查需要的样本含量估计方法一、调查总体均数时样本含量计算公式

46、（当总体标准一、调查总体均数时样本含量计算公式（当总体标准差和均数未知时，都用样本统计量代替）：差和均数未知时，都用样本统计量代替）：式中，式中，n=样本含量；样本含量；=0.05时，时，s=样本标准差；样本标准差；d=允允许误差，即样本均数与总体均数的容许差值，既可用所许误差，即样本均数与总体均数的容许差值，既可用所求总体均数求总体均数95可信区间间距之半（可信区间间距之半（1.96倍标准误）；倍标准误）；又可根据专业知识提出一个容许差又可根据专业知识提出一个容许差d，通常用后者。，通常用后者。设计例题设计例题3某单位拟抽样调查了解该地区某单位拟抽样调查了解该地区20岁男性身高。岁男性身高。

47、据预调查计算的样本均数为据预调查计算的样本均数为170cm，样本标准差为，样本标准差为30cm，设计者要求容许差值，设计者要求容许差值d=总体均数（用样本均数总体均数（用样本均数代替）的代替）的1/100=170cm/100=1.7cm.问：需要多大样本含量？问：需要多大样本含量？答答：设设定定=0.05，s=30cm，d=1.7cm，t0.05=1.96。将将数数据据代代入入公公式式4-5得得:n=(1.96*30/1.7)2=1196.3。即即决决定至少需要抽样调查定至少需要抽样调查1197名名20岁男性。岁男性。如果想减少样本量必然要增大容许误差如果想减少样本量必然要增大容许误差d，例如

48、令，例如令d=总体均数总体均数(用样本均数代替用样本均数代替)的的2/100=170cm*0.02=3.4cm.将数据代入公式将数据代入公式4-5，得，得：n=(1.96*30/3.4)2=299.1即至少需要抽样调查即至少需要抽样调查300名名20男性。男性。设计例题设计例题4某中医科拟抽样调查了解空腹糖某中医科拟抽样调查了解空腹糖200mg/dL的的II型糖尿病气阴两虚证的证候加权积分均值（指总型糖尿病气阴两虚证的证候加权积分均值（指总体均值，这里的总体是研究结果推论的范围，即全国体均值，这里的总体是研究结果推论的范围，即全国的该证型患者）。的该证型患者）。/15=0.2。据预调查样本计算

49、，样本证候加权分均数为据预调查样本计算，样本证候加权分均数为3分，标准分，标准差约为差约为1分，要求容许误差分，要求容许误差d=总体均数（用样本均数代总体均数（用样本均数代替）的替）的1/15=3分分/15=0.2。问：所需要调查多少人数？问：所需要调查多少人数？答答:设设定定0.05,s=1分分,d=0.2分分，将将数数据据代代入公式入公式45，得：，得：n=(1.96*1/0.2)2=96.04即即至至少少需需调调查查97位位空空腹腹血血糖糖200mg/dL的的II型型糖糖尿尿病病患患者者，才才能能了了解解气气阴阴两两虚虚证证的的加加权权均均数数，其其误差在误差在0.2分范围内。分范围内。

50、二、调查总体阳性率时样本含量计算二、调查总体阳性率时样本含量计算当阳性率当阳性率p在在0.20.8时：时：式中式中n为所需样本含量；当时为所需样本含量；当时0.05时，时，p为总体阳为总体阳性率的估计值；性率的估计值；d为容许误差，一般可取总体率的为容许误差，一般可取总体率的1/10。设计例题设计例题5江南某地区欲用抽样调查了解该地区居民江南某地区欲用抽样调查了解该地区居民血吸虫感染率。总体感染率不知。要求样本感染率与血吸虫感染率。总体感染率不知。要求样本感染率与总体感染率之差不超过总体感染率的总体感染率之差不超过总体感染率的1/10，问需调查，问需调查多少人？多少人？答：因为总体感染率一无所