一般线性模型ppt课件ppt课件.ppt

一般线性模型 GeneralLinearModel菜单 1 方差分析中的术语 1 因素与处理 Factor 因素是影响因变量变化的客观条件 处理是影响因变量变化的人为条件 也可通称为因素 用分类变量表示 取有限的离散值 2 水平 Level 因素的不同等级称作水平 水平值取有限的离散值 如 性别因素有两个水平0 1 男 女 等3 单元 Cell 指各因素的水平之间的每个组合 如性别 0 1 和年龄 10 11 12 的六种组合 2 方差分析中的术语 续 4 因素的主效应和因素间的交互效应 如药物A B的主效应及AB的交互效应 5 均值比较 均值的相对比较是比较各因素对因变量的效应大小的相对比较 如研究A B的单独效应之和是否等于它们的交互效应 或A B的效应是否相等 均值的多重比较是研究因素单元对因变量的影响之间是否存在显著性差异 如A B的疗效是否存在显著性差异 6 单元均值 边际均值 在多因素方差分析中 每种因素水平组合的因变量均值称为单元均值 一个因素水平的因变量均值称为边际均值 MarginalMeans 3 方差分析中的术语 续 7 协方差分析 在一般进行方差分析时 要求除研究的因素外应该保证其他条件的一致 作动物实验往往采用同一胎动物分组给予不同的处理 研究不同处理对研究对象的影响就是这个道理 如研究身高与体重的关系时要求按性别分别进行分析 以消除性别因素的影响 要消除其他因素的影响 应采用协方差分析 8 重复测量 组内变异的主要的原因是实验对象之间的个体差异 由于个体差异存在 即使实验对象受到相同的处理 他们的因变量值也可能相当不同 重复测量设计的方差分析也是像协方差分析一样 是在研究中减少个体差异带来的误差方差的一种有效方法 而且由于对相同个体进行重复测量 在一定程度上降低了人力 物力 财力的消耗 如果重复测量是在一段时间内或一个温度间隔内进行的 还可以研究因变量对时间 温度等自变量的变化趋势 这种重复测量研究称为趋势研究 4 方差分析过程 1 One Way过程 单因素简单方差分析过程 在CompareMeans菜单项中 可以进行单因素方差分析 均值多重比较和相对比较 2 GeneralLinearModel 简称GLM 过程 GLM过程由Analyze菜单直接调用 这些过程可以完成简单的多因素方差分析和协方差分析 不但可以分析各因素的主效应 还可以分析各因素间的交互效应 5 一般线性模型 GeneralLinearModel菜单 6 一般线性模型 单应变量的检验 多应变量的检验 重复测量数据的检验 方差分量分析 7 Univariate子菜单 提供回归分析和一个因变量和一个或几个因素变量的方差分析 Multivariate子菜单 可进行多因变量的多因素分析RepetedMeasures子菜单 可进行重复测量方差分析VarianceComponents子菜单 可进行方差成分分析 通过计算方差估计值 可以帮助我们分析如何减小方差 8 成组设计的方差分析 即单因素方差分析 配伍设计的方差分析 即两因素方差分析 交叉设计的方差分析析因设计的方差分析重复测量的方差分析协方差分析等 9 某湖水不同季节氯化物含量测定值如 分析不同季节氯化物含量有无差别 10 方法1 comparemeans onewayANOVA 11 一般线性模型 12 13 14 例对小白鼠喂以A B C三种不同的营养素 目的是了解不同营养素增重的效果 采用随机区组设计方法 以窝别作为划分区组的特征 以消除遗传因素对体重增长的影响 现将同品系同体重的24只小白鼠分为8个区组 每个区组3只小白鼠 三周后体重增量结果 克 列于下表 问小白鼠经三种不同营养素喂养后所增体重有无差别 问题1 影响增重的因素有几个 问题2 要设几个变量 15 16 应变量 固定因素 随机因素 协变量 限选1个变量 且为数值变量 分组变量 可选1个或多个 加权变量 17 用于设置在模型中包含哪些主效应和交互因子 默认情况为Fullfactorial 即分析所有的主效应和交互作用 定义模型 全因素模型 自定义模型 18 进入模型的因素交互作用级别 即是分析主效应 两阶交互 三阶交互 还是全部分析 方差分析模型类别 是否在模型中包括截距 因素及其类型 F 固定因素 R 随机因素 C 协变量 模型方式 19 20 Contrast钮 弹出Contrast对话框 用于对精细趋势检验和精确两两比较的选项进行定义 使用频率少 Plots钮 用于指定用模型的某些参数作图 比如用food和group来作图 用的也比较少 指国内 因为它主要是用来做模型诊断用的 PostHoc钮 该按钮弹出的两两比较对话框 本题对food作两两比较 方法为SNK法 Save钮 将模型拟合时产生的中间结果或参数保存为新变量供继续分析时用 如预测值 残差 诊断用指标等 Options钮 定义选项 可以定义输出哪些指标的估计均数 并做所选择的两两比较 还有其他一些输出 如常用描述指标 方差齐性检验等 21 首先是所用方差分析模型的检验 F值为11 517 P小于0 05 因此所用的模型有统计学意义 可以用它来判断模型中系数有无统计学意义 第二行是截距 它在我们的分析中没有实际意义 忽略即可 第三行是变量GROUP 可见它也有统计学意义 不过我们关心的也不是他 第四行是真正要分析的FOOD P值为0 084 没有统计学意义 结论是 尚不能认为三种营养素喂养的小白鼠体重增量有差别 22 随机区组设计的方差分析 例9 2为研究不同卡坏对牙齿的固定效果 以10颗取自新鲜尸体的牙齿为实验对象 每颗牙齿同时随机在不同部位固定三种卡环 测试抗拉强度 分析3中卡环的固定效果有无差异 23 Group变量Teeth区组变量Pull因变量 两因素 24 25 26 27 28 不能用fullfactorial 两因素间不存在交互作用 29 30 31 32 对于初学者来说可能有用的几个问题 需要分析的影响因素可以都选入fixedfactor框 如果不是复杂的模型 一般分析结果不会有误 方差分析模型多数情况下要选modelIII 但这在数据存在缺失值 设计不平衡等情况下要慎重考虑 因为此时往往会要求模型进行详细的设置 model的设置对分析是非常重要的 如果设置不正确 可能什么都做不出来 比如无重复数据的方差分析纳入了交互作用 析因设计的方差分析纳入了设计中不存在的因素 就会做不出结果 一般线性模型的复杂性是超出大家想象的 一但有存在疑问的内容 一定要查阅有关统计书籍 33 重复测量资料数据形式 表12 3 g 1m 4n 8 单组重复测量资料 34 重复测量设计 方法 重复测量的方差分析目的 推断处理 时间 处理 时间对试验对象的试验指标的作用资料 处理因素分g个水平 每组随机分配n个试验对象 共ng个 g 1时间因素分m个水平 m个时点 每个对象有m个时点上的测量值 共gnm个 m 2特例 g 1 单组重复测量资料m 2 前后重复测量资料 35 实验操作方法 重复测量数据的两因素多水平设计 两因素包括一个干预因素 A因素 和测量时间因素 B因素 多水平指干预 A因素 有g 2 个水平 测量时间 B因素 有m 2 个水平 测量时间点 随机化分组采用完全随机设计的分组方式 将gn个观察对象随机分配到g个处理组中 数据收集在m个时间点上进行 每一个观察对象在完全相同的时间点上重复进行m次测量 36 表12 2高血压患者治疗前后的舒张压 mmHg g 2m 2n 10 37 表12 2数据的统计学分析问题 计算前后测量数据的差值 上述数据即可转化为完全随机设计 两组 的资料形式 一般情况下 针对前后测量数据差值的成组t检验方法是可取的 但应注意其应用条件 即方差齐性的问题 38 对象内Mk 对象间Bj g 3m 5n 5 Ai Tij 39 重复测量设计资料的统计分析方法 对于重复测量数据 临床上常称纵向监测数据 实质上每个受试对象的观察结果是多次重复测量结果的连线 统计分析的目的是比较这些连线变化趋势的特征 重复测量试验数据的方差分析需要考虑两个因素 一是处理分组 二是测量时间 可采用的统计分析方法 1 多元方差分析方法2 重复测量数据的方差分析 40 变异分解思路 重复测量数据的变异由两大部分组成 一是观察对象间差异 二是重复测量间差异 观察对象间差异包括处理组间差异和观察对象个体间变异两部分 重复测量间差异包括测量时间之间差异 处理与测量时间的交互作用和组内误差三个部分 因此 重复测量数据的总变异可分解为处理组 测量时间 处理组与测量时间的交互作用 观察对象间随机误差以及重复测量误差等五个部分 41 42 43 44 定义重复测量的变量名 重复测量的模型 应变量被重复测量了几次 分别存放在几个变量中 所以我们这里要自行定义应变量 定义重复测量次数 添加 修改和删除重复测量的变量 定义嵌套的重复测量变量 45 46 对象内变量 对象间变量 47 五次重复测量的变量名 分组变量情况 48 对组内变量以及它和分组变量交互作用的多元方差分析 四种多元检验方法 49 球形检验结果 满足球