多水平模型简介

MLwiN worksheet for practice “satisfaction.ws”
Variables Patid: Patient id Sat: raw score of satisfaction Phys: physician id Pt_age-60: patients’ age – 60 Md_age-44: physicians’ age – 44 Yrs_prac-20: years of physician’s working experience – 20 Const: constant New_sat: transformed sat score Ptage_ave: average age of patients who attended the same clinic
示例二
1482名Ⅱ型糖尿病患者对69名社区医生
服务的满意度评价，反应变量为满意度分数
(取对数)。病人年龄31－91岁，社区病人平 均年龄54－74岁。
来自百度文库
考查社区年龄对个人满意度的影响
表1.2 Ⅱ型糖尿病患者对社会医疗服务的满意程度模型比较
单水平线性回归 单水平集合线性回 归 两水平线性回归
(SE）
2. 模型可准确地估计高水平变量对个人水平 上地反应变量的影响的程度
MLwiN多水平模型分析软件
1. 使用迭代广义最小二乘估计算法(IGLS)和 MCMC算法 2. 能处理任意水平的层次结构数据 3. 窗口工具和宏程序功能 4. 拟合模型包括以下类型的数据： 正态、Binary/Binomial、Poisson，多项 式、负二项分布、生存时间、集合数据
1289.02
0.1354(0.0051)
处理非独立数据传统方法的弊端
1.S.E.的稳健估计、回归系数的校正法
－治标不治本
2.GEE估计
－将数据聚集原因处理为干扰因子而弃之
多水平模型对传统方法的补充
1. 在模型中进一步分解由不同原因导致数据
非独立所致的方差分量，进而获得准确的
回归系数估计值及其标准误
反应变量相似性(非独立性)图示
非独立数据的一般结构
社区
个人 个人
动物窝别
幼仔 幼仔
测量1
个体
测量2
二水平层次结构：水平一单位在水平二内聚集
忽略多水平层次结构的后果
1.回归模型中的参数估计值尤其是其标准误有
偏差
2.残差方差偏大，即模型拟合优度差
3.损失环境因素(水平二变量)对反应结果的影 响信息
1.正态性不满足 －数据变换，增加样本含量 2.方差非齐性 －增加协变量 －数据变换 －广义线性模型或非线性模型
假设不满足时的处理
3.独立性不满足 －S.E.的稳健估计 －基于组内相关系数校正回归系数的校正法 －GEE估计方法 －拟合非独立性来源的模型
非独立性来源
1.生活、工作、社区环境对反应变量的影响
－卫生服务区域的资源、社会经济条件和政策
会影响对病人的服务质量
－高血压发病率可能有地区聚集性，取决于经
济文化背景和居民饮食习惯
非独立性来源
2.同一观察对象、某一观察结果的重复测量 无疑具有强相关 －分子生物学研究中重复测量数据处理中的问
题
3.区组设计和多中心试验 －卫生毒理实验研究中同窝动物的相似性 －同中心内病人病情、病种相似性
MLwiN主要窗口工具
1. 数据处理：
－Excel文件和TXT文件的输入，TXT文件输出
－缺失值处理，结果显示精度
－数据修改，筛选和显示 －构造数据的层次结构
MLwiN主要窗口工具
2. 模型定义窗口 － Equations 3. 算法选择窗口 － Estimates 4. 显著性检验窗口 － Intervals and tests 5. 残差计算窗口 － Residuals 6. 模型预测窗口 － Predictions 7. 结果图示解释 － Graphs
个人年龄 社区平均年龄 0.0032（0.0009) -0.0298(0.0029）
(SE）
--------0.0231（0.0039）
（SE）
0.0032（0.0008） -0.0288（0.0040）
-2对数似然值
残差标准误平方
1311.8
0.1419
-99.81
0.01378（0.0024）
多水平模型简介
杨珉 英国伦敦大学玛丽王后医学院 二ΟΟ六年七月二日
单水平模型复习
1,2,...,i,...n个观察对象
y i 0 1 x i e i , e i ~ N (0, e )
2
模型假设： 正态性、独立性、残差方差齐同性 协变量的影响在所有对象不变
假设不满足时的处理
示例一
全国第三次卫生服务调查在四川农村抽
样的部分数据2239个调查对象来自647个住
户，反应变量两周内患病与否有家庭聚集性。
单水平模型与二水平模型结果对比
表1.1 两周患病及影响因素的单水平正态模型和两水平正态线性模型的比较
变量
单水平线性模型 （SE）
两水平线性模型 (SE）
0.075（0.026）
1. 系统数据文件 “ satisfaction .ws” 熟悉主要窗口的界面和功能
2. 整理自己的数据，使用数据处理将其转为
MLwiN系统文件.
年龄≥65
0.085（0.026）
女 经常饮酒
患有慢性病 -2对数似然值 残差标准误平方
0.036（0.018） 0.097（0.025）
0.474（0.026） 2215.2 0.157（0.005）
0.035（0.017） 0.088（0.025）
0.463（0.026） 2180.87 0.139（0.005）