医学科研设计与统计分析错误辨析与释疑.ppt

题目 医学科研设计与统计分析错误辨析与释疑单位 军事医学科学院生物医学统计咨询中心 1 问候 向您们 重视科研工作科学性与严谨性的人们表示衷心地感谢 希望我的讲课能给您们带来一些欢乐和帮助 2 开场白 别人喜欢以美酒会友本人习惯以问题酬宾首先请大家思考并回答下面的问题 3 问题1 错在哪 为了将某种治疗失眠的药物与安慰剂的疗效进行比较 根据失眠程度打分 1 轻度 2 中度 3 重度 4 极重度 很多人用t检验比较试验组与对照组失眠的平均值 若有5个访视点 仍对每个时间点进行t检验 对吗 4 问题2 错在哪 某单位试图比较6种降压药 含安慰剂 对飞行员的降压效果 在一段时间内让他们每人随机地依次服6种药 每服一种药要重复观测7个数据 即从每位飞行员身上要观测42个数据 然后分析各药物的疗效 5 问题3 用t检验处理下面的资料错在哪 6 心功能不全 与溶栓的关系 组别有 例数无 例数溶栓 n 216 75141未溶栓 n 356 162194计算结果 2 6 422 P 0 01 原作者的结论 溶栓与心功能不全的相关分析 发现两者的联系有统计学差异 问题4 这样下结论错在哪 7 问题的总结 在实际工作者运用统计学的过程中 从基本概念 科研设计 统计分析 到结论的陈述 错误无处不在 80 的误用率是十分保守的估计 我们生物医学统计咨询中心试图通过我们的努力 使80 10 8 介绍 军事医学科学院生物医学统计咨询中心成果3项 著作7部 论文百篇举办全国培训班数十期 业务范围详见下页 9 服务范围 1 大 中型科研课题协作 2 新药临床试验数据分析 3 生物医学科研人员和临床大夫科研素质的提高 科研设计与统计分析能力的培养 10 敬告 科研要创新 创新不容易 理念很重要 设计定乾坤 资料要准确 分析应合理 推理有依据 结论方可信 劝君别着急 我能帮助您 11 培训班消息 科研设计和统计分析及统计软件 SAS SPSS 应用培训班时间 2004年08月12 17日地点 北京市太平路27号主讲 胡良平 12 与胡良平联系方式 电话及传真 010 66931130E mail hulp lphu812 网站 100850北京市太平路27号军事医学科学院生物医学统计咨询中心 13 讲授计划 一 统计分析方法合理选用的预备知识与要领二 统计分析常见错误辨析与释疑 14 专题一 统计分析方法合理选用的预备知识与要领 15 一 预备知识 1 统计学包括哪些内容 2 什么叫因素和水平 3 资料类型如何划分 4 什么叫交互作用 16 一 预备知识 5 什么叫多因素统计分析 6 什么叫多元统计分析 7 什么叫实验设计类型 17 1 统计学包括哪些内容 统计研究设计 统计资料的表达与描述 定量与定性资料统计分析 相关与回归分析 多元统计分析 生存资料统计分析 时间序列分析 18 2 1 什么叫因素和水平 因素 研究者关心的实验条件 水平 因素的具体表现形式 例 温度 30 60 90 药物种类 A B C 药 药物剂量 1 5 10 mg 19 2 1 实验因素与区组因素 实验因素 研究者特别关心的实验条件 如 温度 药物种类 药物剂量 等 区组因素 即重要的非实验因素 如 性别 年龄 体重 20 3 资料类型如何划分 定量 计量 计数 定性 名义 有序 计量资料 如成人血液中钾含量 4 6mmol L 计数资料 如脉搏 72次 min 名义资料 如血型 A B O AB 有序资料 如疗效 痊愈 好转 无效 死亡 21 4 什么叫交互作用 假定某药对某病患者的疗效如下 药物剂量 小中大A优中差B中优差C差中优三种药对某病的疗效将随着剂量的改变而改变 22 5 什么叫多因素统计分析 一般指 同时考察影响结果的原因变量有多个 而结果变量只有一个 如 多因素设计资料的方差分析 多元线性回归分析 多元Logistic回归分析 等 23 6 什么叫多元统计分析 一般指 同时考察影响结果的原因变量有一个或多个 而结果变量也有多个 如 多因素设计资料的多元方差分析 多重线性回归分析 主成分分析 聚类分析 等 24 7 什么叫实验设计类型 在设计实验和分析定量资料时 都不可避免地要涉及到 实验设计类型 的概念 简单地说 实验设计类型就是合理安排实验因素和区组因素的具体方法 例如 25 7 什么叫实验设计类型 例1 成组设计 分别测量学术厅内男性与女性的身高所得的实验资料 就叫来自 成组设计的定量资料 例2 配对设计 将在场的所有人按性别 年龄 身高配对 然后用两台血压计分别测各对人的血压所得的资料叫来自 条件相近者配对设计资料 26 最常见的设计类型 单因素设计 包括 单组设计 配对设计 成组设计 单因素多水平设计 析因设计 具有重复测量设计 27 最常见的设计类型 三类常见实验设计类型举例 28 例1 已知正常人乙酰胆碱脂酶的平均值为1 44单位 现测得13例慢性气管炎患者的乙酰胆碱脂酶分别为 1 50 2 19 2 32 2 41 2 11 2 54 2 20 2 36 1 42 2 17 1 84 1 96 2 39问 慢性气管炎患者与正常人的平均乙酰胆碱脂酶之间的差别有无显著性意义 单组设计举例 29 例2 从八窝大白鼠中分别选出同性别 体重相近的两只 喂以水解蛋白和酪蛋白的饲料 四周后测定其体重增加情况 结果如下 窝编号12345678含酪蛋白组8266747882767390含水解蛋白组1528292824382137问 两种饲料对大白鼠体重增加量之间的差别的影响有无显著性意义 配对设计举例 30 例3 有两组小儿肝炎患者 其一为甲 乙肝炎病毒混合感染 其二为重叠感染 治愈者的治愈天数如下 混合感染 n1 9 552514264223402756重叠感染 n2 9 2034173278816577725试比较两组疗效之间的差别有无显著性意义 成组设计举例 31 例4 研究单味中药对小白鼠细胞免疫机能的影响 把40只小白鼠随机分为四组 每组10只 雌雄各半 用药15d后 进行E 玫瑰结形成率测定 结果如下 对照组1410121613141210139党参组21241817221918232018黄芪组24202218172118221923淫羊藿组35273329314035302836 单因素四水平设计举例 32 例5 某医科大学病理生理学教研室研究三种因素 小鼠种别A 体重B和性别C 对皮下移植SRS瘤细胞生长特性影响的结果 A B C三因素各有两个水平 A分为A1 昆明种 A2 沪白1号 B分为B1 24 25克 B2 13 15克 C分为C1 雄性 C2 雌性 共选了24只小鼠 在接种后第8天测得肿瘤体积见表1 请问 这是一种什么设计类型 应当对此资料作怎样的分析 33 表1三因素影响下小鼠第8天肿瘤体积 因素肿瘤体积 cm3 C因素A与B A1 B1B2 A2 B1B2 C10 70691 08380 06280 47120 78540 94250 09420 08800 35810 33350 04710 1759C20 07850 50270 01260 22460 18850 95500 01260 25130 34030 92150 00940 36760 25030 85140 01250 1327 三因素析因设计举例 34 表2接受不同处理的家兔血浆中K 含量的测定结果 处理组家兔血浆中K 含量 ppm 因素A 编号缺氧时间T h T1 0 T2 0 5 T3 1 T4 2 A1 适应组 1154 5129 8122 7171 72173 0124 0170 4168 63186 0131 0137 0138 04161 0154 0178 0128 05187 0158 0162 0152 0A2 平原组 6144 4135 5149 9129 17147 1134 2138 3146 18183 6189 3190 5227 39181 7160 3163 3163 310166 7136 8134 6142 5A3 急性缺11173 3231 4293 7401 4氧组 12155 1199 6191 8203 413177 9153 6158 7240 314158 2146 7135 2228 715180 4170 8226 2267 7 具有一个重复测量的两因素设计 35 例6 某研究者为了探讨两种不同处理 因素A 和不同剂量 因素B 作用下不同时间点 因素T 上NFS 60细胞凋亡情况 将全部样品随机地均分为6个组 每组含3个样品 对各组中的每个样品在处理后5个不同时间点上观测200个细胞 数出凋亡细胞个数 因素A分为 即A1 单纯照射 A2 照射 过钒酸钠 Perv 各处理组的照射剂量 因素B Gy 又分为0 3 10 设计格式和资料见表3 试判断此资料所取自的实验设计类型 并对资料作相应的统计分析 36 表3不同处理和不同剂量作用下不同时间点上NFS 60细胞凋亡情况 处理组照射剂量样品200个细胞中凋亡细胞数 因素A Gy 编号时间T h T1 12 T2 24 T3 36 T4 48 T5 60 A1 照射 01377101725357183568141431471020332341422373359193110161838831012416429112035184592117A2 加Perv 014641214224614173373111931577141723591219376615161015101633492412123744339143652 具有一个重复测量的三因素设计 37 例7 用贲门癌患者的标本制成液体 在三种不同处理条件下观测鸡胚背根神经节与鸡胚交感神经节中长出突起的神经节的比例 现有贲门癌患者10例 将每人的标本均分成三份 分别给予三种不同的处理 因素A 即A1 加入100ng ml神经生长因子 A2 加入200ng ml神经生长因子 和A3 单用贲门癌培养液 并对每种处理后的标本中的两种类型的神经节 因素B 即B1 背根神经节 与B2 交感神经节 观测长出突起的神经节的比例 Y 设计格式和资料见表4 试判断此资料所取自的实验设计类型 并对资料作相应的统计分析 38 表4贲门癌患者的标本经三种处理后两种神经节中长出突起的神经节的比例 病例号Y 长出突起的神经节的比例 A1 B1B2 A2 B1B2 A3 B1B2 10 500 430 500 430 800 5020 630 380 550 500 710 6330 500 500 540 500 830 6740 430 430 500 380 700 5750 500 400 500 570 700 7160 440 380 630 500 770 6370 430 500 500 380 690 5080 630 440 400 380 700 5090 440 500 500 440 670 50100 440 430 450 400 600 75 具有两个重复测量的两因素设计 39 二 合理选用统计分析方法的要领 1 定量资料统计分析的要领 2 定性资料统计分析的要领 3 相关与回归分析的要领 4 多元资料统计分析的要领 40 1 定量资料 统计分析方法合理选用的要领 1 正确判定设计类型2 检查资料前提条件 41 定量资料的参数与非参数检验 参数检验 u检验 t检验 方差分析 T2检验 Wilks s 检验 非参数检验 符号秩检验 秩和检验 等 42 参数检验前提条件 独立性正态性方差齐性 43 可用t检验处理的定量资料所对应的设计类型及前提条件 独立性1 单组设计 正态性 2 配对设计 差量正态性 3 成组设计 正态性 方差齐性 44 可用方差分析处理的定量资料所对应的设计类型及前提条件 设计类型 单因素多水平设计其他多因素设计前提条件 独立性正态性方差齐性 45 定量资料实例 要想合理选用统计分析方法处理下面的定量资料 应当如何着手考虑 46 表1甲 乙两种药物治疗高血脂病前后血清胆固醇 mmol L 变化的测定结果 编血清胆固醇值 mmol L 号甲药 治疗前治疗后乙药 治疗前治疗后 16 786 276 805 4427 386 486 505 4036 625 777 345 0346 155 557 155 0456 005 306 105 9767 107 106 354 6876 455 807 506 0886 156 276 155 719 6 006 9510 6 455 0111 6 706 1012 6 176 45 47 本例的设计名称 具有一个重复测量的两因素设计 一个 指重复的方向上有一个因素 即 时间 两个 指整个实验涉及到两个实验因素 即 时间与药物 48 2 定性资料 统计分析方法合理选用的要领 1 认准列联表类型2 检查资料的前提条件3 弄清分析的目的 49 列联表类型 2 2表有四型 R C表有四型 高维表有三型 每一类仅举一例 50 第一类表 介绍2 2表中的第一型 即横断面研究设计的2 2表 51 横断面研究设计2 2表举例 手术时感染未感染间 h 例数例数 513229 5756 52 横断面研究设计2 2表计算方法 1 N 40且理论频数T 5 可用一般卡方检验 2 N 40但至少有一个1 T 5 可用校正的卡方检验 3 N 40或有一个T 1 应选用Fisher的精确检验 4 若用软件计算 均可选用Fisher的精确检验结果 53 另外三类四格表 1 来自队列研究设计的四格表 2 来自病例 对照研究设计的四格表 3 来自配对设计的四格表 实例和分析方法从略 54 第二类表 介绍R C表中的第二型 即结果变量为有序变量的R C表 55 单向有序R C表举例 56 结果变量为有序变量的单向有序R C表计算方法 1 秩和检验 2 Ridit分析 3 有序变量的Logistic回归分析 57 另外三型R C表 1 双向无序的R C表 2 双向有序且属性相同的R C表 3 双向有序且属性不同的R C表 实例与方法选用从略 58 第三类表 介绍高维表中的第一型 即结果变量为二值变量的高维表 59 结果变量为二值变量的三维列联表举例 60 高维表的计算方法 结果变量为二值变量的高维列联表的计算方法 1 加权卡方检验 2 多元Logistic回归分析 3 对数线性模型 61 另外两型高维表 1 结果变量为多值有序变量的高维列联表 2 结果变量为多值名义变量的高维列联表 实例与方法从略 62 3 相关回归分析要领 统计分析方法合理选用的要领 1 以专业知识为依据2 绘制变量间散布图 63 3 相关回归分析要领 例如 测量教室中10名男性的身高与体重的数据 可否进行直线相关与回归分析 若测量教室中10名男性身高 测量10名女性体重的数据 可否进行直线相关与回归分析 64 3 相关回归分析要领 测量教室中10名男性的身高与血压的数据 可否进行直线相关与回归分析 若测量10只动物用某药后的血糖随时间变化的数据 可否进行直线相关与回归分析 65 4 多元资料 统计分析方法合理选用的要领 1 以专业知识为依据2 了解变量基本关系3 弄清统计分析目的 66 专题二 统计分析常见错误辨析与释疑 67 专题二提纲 第一部分实验设计第二部分定量资料统计分析第三部分定性资料统计分析第四部分直线相关回归分析 68 专题二提纲 第五部分多元统计分析第六部分资料的综合分析第七部分结果的解释第八部分统计描述与表达 69 第一部分 与实验设计有关的问题 请辨析问题中与实验设计有关的错误 并分析产生错误的根本原因 危害性和消除错误的对策与建议 70 例1 1 买得起马 买不起鞍 2 缺胳膊少腿我不怕 因为我身穿盔甲 3 迫着脑袋定方案 随心所欲高科研 71 例1 的具体内容 为了考察阿仑膦酸防治牙槽骨吸收的效果 有人设计了如下的实验 取3 4个月龄的SD大鼠45只 将45只大鼠随机分入六个组 其中第一组至第五组均为8只大鼠 第六组为5只大鼠 各组含义如下 72 组别形成各组实验的条件 1给大鼠去势 结扎 用药2给大鼠结扎 用药3给大鼠去势 结扎 不用药4给大鼠结扎 不用药5给大鼠去势 不用药6空白对照 即假手术 组 错在哪 73 结构变形分析法 结构变形现原型 结扎与否去势与否用药与否 用不用结扎去势A组B组结扎不去势C组D组不结扎去势X1E组不结扎不去势X2F组由于原作者少设计了两组 导致一个非平衡的多因素设计的出现 组间可比性差 结论可靠性差 74 例2 独生子女成了 小皇帝 为了观察甲紫注入某种动物正常腮腺后组织病理变化情况 有人选择6个月龄 体重20 25kg的中国实验用小型猪15只 雄性9只 雌性6只 每只动物任选一侧腮腺为实验侧 另一侧作为正常对照 以消除个体差异及增龄对实验结果的影响 75 按注入药物后1周 2周 1个月 3个月及6个月将15只动物随机分为5组 每组3只 每个组的3只动物分别随机注入0 6ml 1 0ml及4 0ml1 某药物溶液 然后观察组织病理变化情况 76 例3 动物不会提抗议 任我将它们来摆布 某研究者为研究bFGF对肌成纤维细胞的作用及其对创面愈合的影响 利用大鼠30 深 烫伤模型 观察的指标有 平滑肌肌动蛋白等定量指标 77 原作者设计的三个组如下 1 正常对照组 将大鼠背部皮肤置于37 水中8s 2 单纯损伤组 造成30 TBSA深 烫伤 在伤后3h切取创面组织检测多项定量观测指标的数值 3 在 2 的基础上 用bFGF治疗 在伤前 伤后3h 6h 1d 3d 7d 14d切取创面组织检测多项定量观测指标的数值 将第 3 组产生的7个小组数据分别与第 1 2 组比较 用t检验 问 这样的实验设计科学吗 78 例4 单枪匹马斗群雄 分4批实验分别确定4个实验因素的所谓 最佳水平 1 摸索 培养基含水量 对产物的影响 2 摸索 起始pH值 对产物的影响 3 摸索 温度 对产物的影响 2 摸索 光线 对产物的影响 将四批实验的最好结果组合在一起 就一定是最优的实验条件吗 79 第二部分 与定量分析有关的问题 所谓 与定量分析有关的问题 就是在实验研究中 影响因素为定性的 有时也包含定量的 观测结果是定量 希望比较各实验条件下定量观测结果的平均值之间的差别是否具有统计学意义 80 例1 单腿 踢 向全世界 用多次t检验处理此资料合适吗 81 例2 组别 一词真厉害它使 博导 成 驳倒 某研究者在研究 N 甲基 D 天冬氨酸及其非竞争性受体阻滞剂和钙离子对培养的鼠视网膜神经细胞的作用 时 做了一个较复杂的实验 其设计格式和资料如表1所示 82 主要实验结果见下表 用t检验处理下面的资料正确吗 83 对差错的辨析与释疑 该实验研究中涉及到3个实验因素 NMDA N 甲基 D 天冬氨酸 用否 MK801 非竞争性受体阻滞剂 用否 Ca2 用否 84 对差错的辨析与释疑 何设计 怎样分析 85 例3 夏天穿棉袄 反常 某人在研究 准分子激光治疗性角膜切削术治疗角膜浅层病变的实验和临床研究 时 收集到表1资料 86 例3 续 87 该资料的正确表达格式 88 例4 问题虽然不复杂 却能难倒科学家 您若不信请尝试 吃一堑来长一智 89 表1不同药物对小鼠迟发超敏反应的影响结果 药物剂量鼠数耳肿重量P值 g kg 只 mg 与对照比与Cy比 对照 1021 2 2 7补肾药51022 3 3 5 0 05补肾药101018 8 3 1 0 05补肾药201016 5 2 4 0 01Cy0 0251011 2 1 5Cy 补肾药0 025 51014 3 2 9 0 01Cy 补肾药0 025 101018 6 3 6 0 01Cy 补肾药0 025 201019 2 3 4 0 01 注 补肾药全称为补肾益寿胶囊 请问 如何正确分析此资料 90 表2表1资料的第一种变形结果 药物耳肿重量 mg 种类药物剂量 g kg 00 02551020 补肾药21 2 2 7 22 3 3 518 8 3 116 5 2 4Cy药21 2 2 711 2 1 5 注 各组均有10只小鼠 表示补肾药未用的剂量 表示Cy药未用的剂量 91 表3表1资料的第二种变形结果 Cy药剂量耳肿重量 mg g kg 补肾药剂量 g kg 051020 021 2 2 722 3 3 518 8 3 116 5 2 40 02511 2 1 514 3 2 918 6 3 619 2 3 4 注 各组均有10只小鼠 92 例5 我只问一个 简单问题 某临床医生前来进行统计咨询 她说 我有一个简单的问题 其大意如下 给30名某病患者服用一种药物后测量其体内的放射性物质的含量 具体做法如下 每天在离每位患者不同距离处 5种距离 测量其体内的放射性物质的含量 连续共测了8天 观测指标是定量的 请