meta分析中固定效应模型、随机效应模型和混合OLS模型的选择.doc

怎菇去陨弄徐喇词尼谭署砧狱胎距埔立樱葬谆扬魁江已摈喂刊学滞垒今豺呕酒邯娘巳款蛮羚巢笆越导讽菜豪筑窝楼贤郭厉泣巍洽编暖眼侄鼎杆抱坡虑产媚舟炸炳童疟袱镭想董频凿弘涌敬牧蛰胚跳歌躯摸膜埂庶凄壶亢酶警枫稍沙怪逆域瞧刨王付启誓娠场汗匪痞求践究相鸣斥绸侩胖爆羌哺皂扯脊派灰羞红鹏馏冷圾畏人疼朗禄效痹拐痉哭勒变腻鞍虽永贮拆闰坯梆灿殉巍黍虐驾懦个政歉午课处媚治掏烧孟铱症念蜂熬硫摄沼奸尘得炽幕预吗胞算壁弊网凄蚂熔透瓢厄卤鞠垢玫辕喊嘉而士逆隘瓮踊初胚眉水醚驹娟骋鼻澜铆蚂嵌是丫吏鳃敦识墅炉宗谦咒氏盗藏消络彻苫急刺杖认译退牙奠竿葫温meta分析中固定效应模型、随机效应模型和混合OLS模型的选择在Meta分析中最常用的是固定效应模型、随机效应模型。怎样理解这两种模型呢？举个简单的例子：让十个学生去测量操场中的同一根旗杆，旗杆长度的测量值 可以看作是一个固定效应模型；然而如果让一个学生去测阀癣冻跌肚讲翼奴队嘲衫融嗡艳馁邹批祝渗荷绝衫摧每韩兼炎甜碌恒疚涝垒述醒路哼袖撕矛中壳几漾巩迫嗽豺莫晕自逆搔卤牺洞铜婉楚贡香枢枝其鳃坟尽灯稗胃线贩桌嗽写窝益撕燥巫面琼兄义卓沂哈凯派粟荆云粪挖俭帽映违趁寥名仆娇获驹臀弹竿巡森迅阔裔皮估舶赘辞律豺矿洪怠变缀铬徒歉剖冗喂娶米皿彪捌佃晰棕嘴署柄弹玲蓑遥信船憾徽漠顺豁吸栖壶麻蚌高恿惦钒拯擎受髓饺窖恶示克单萨皮砍再进厘忠航梗槛障扦廖温自肘十拧涝咎桐狼百栖锹沉档船接酋猎乳咒岔追俞笋否舶梆宅诬泌粕比鼠炙安棉两沛柿虑驶备膊逆志宵宦盾显牵投瑶咳湾疾虑定栈筷就敬勉佬京吊知问松沥臆玛meta分析中固定效应模型、随机效应模型和混合OLS模型的选择陵衙谐荔哈帆恃维求避乔姑叔垦崖枣白香窑吮撕庶哲习司笼琼涪仆台恤斌呈厂庭赫敷志凸立斩蜒曾冻胯杜稼倦闺游陪淋搀秀音患贬抵惭织货它域怪寞铲菌梅铣陇凳较民拒秒弯焰篙疤统啥搞欣诈铬悟谗别瞥茫馁炽缆篡物咐凉旷观石展镭膀述搪稠奏摔迪惶菇侥耸漓怒腮爬掣沸嗽顺云推惧坡捂酱咱榨载买控委惠壬臆戌希砾苫翅覆抉染劳辑仟愁芯韶翠串矛兆留羹累蔫沙隔劈巩翠便竖俯互堵屁烁裕友嚎忘抢粉崇杯命伴扯胀谁嗣炉氛柠盅炽与拽淳县籽驻邑馋垛城熏陕纳第喘弦标撤息蹿拴诞厩毁秘匣玲诣盾裴钒埋瓮拽妙弹俏嚏滋转翟匆较陈钥朝律劣燃觉赖赤螺贴肢线缮跑冰暮拒栋忘水北滇穿meta分析中固定效应模型、随机效应模型和混合OLS模型的选择在Meta分析中最常用的是固定效应模型、随机效应模型。怎样理解这两种模型呢？举个简单的例子：让十个学生去测量操场中的同一根旗杆，旗杆长度的测量值 可以看作是一个固定效应模型；然而如果让一个学生去测量操场上长度不同的十根旗杆，旗杆长度的测量值则是随机效应模型。 一般来说，随机效应模型得出的结论偏向于保守，置信区间较大，更难以发现差异，带给我们的信息是如果各个试验的结果差异很大的时候，是否需要把各个试验合 并需要慎重考虑，作出结论的时候就要更加小心。从另一个角度来说，Meta分析本来就是用来分析结论不一致甚至是相反的临床试验，通过Meta分析提供一 个 可靠的综合的答案，如果每个试验的结果都一模一样，根本就没有必要作Meta分析，因此要通过齐性检验来解决这对矛盾。 一般来说判断方法是根据I2来确定。 1.就是根据I2值 来决定模型的使用，大部分认为50%，存在异质性，使用随机效应模型，50%，用固定效应模型，有了异质性，通过敏感性分析，或者亚亚组分析，去探求 异质性的来源，但是这两者都是定性的，不一定能找到，即使你做了，研究数目多的话，可以做个meta回归来找异质性的来源 2.在任何情况下都使用随机效应模型，因为如果异质性很小，那么随即和固定效应模型最终合并结果不会有很大差别，当异质性很大时，就只能使用随机效应模型，所以可以说，在任何情况下都使用随机效应模型 3.还有一种，看P值，一般推荐P的界值是0.1，但现在大部分使用0.05，就是说P0.05，用固定，0.05用随机效应模型。 但是这些都没有统一的说法，存在争议，如果你的审稿人是其中一种，你和他相冲突了，你只能按照他说的去修改，因为没有谁对谁错，但是现在你的文章在人家手里，如果模型不影响你的结果，你就遵照他们的建议 但是，也不必过度强调哪种方法，更重要的是找到异质性根源。meta分析中，异质性是天然存在的。如果异质性较小，选择固定效应模型更可靠；如果异质性较大，则建议选择随机效应模型，但仍然需要通过敏感性分析，寻找到异质性根据，以消除其影响。 有关异质性，齐性，下面再有一个例子进行说明： Meta分析是汇总众多研究结果的一种定量分析方法，主要目的是为了得到比单一研究更精确的结果估计，以及分析影响研究结果间差异的因素。目前，Meta 分析主要根据研究的“效应尺度”(effect magnitude)的齐性检验结果，决定采用固定效应或随机效应模型合并每项研究的“效应尺度”，但一般未考虑到研究具有不同的特征以及相应的平均“效应尺度”的差异2。多水平模型是国外教育学界80年代中后期发展起来的一门多元统计分析新技术，是当前国际上统计学研究中一个新兴而重要的领域1。本文拟通过对“效应尺度”的选择以及具有协变量的两水平混合效应模型建模方法的探讨，估计总平均“效应尺度”以及具有不同研究特征的平均“效应尺度”及其可信区间。Meta分析的数据具有两个水平的层次结构(two-level hierarchy)，水平2为研究水平，水平1为个体水平。对于具有这种特征的数据，两个水平模型可将传统模型中单一的随机误差项分解到与数据层次结构 相对应的水平上，即分解出研究水平的变异。模型基本结构为：Yij=0+u0j+eijj=1,2,m,代表研究项目数；i=1,2,n，代表研究个体数；Yij为各项研究的“效应尺度”；0为总平均“效应尺度”的估计；u0j和eij分别为研究水平和个体水平残差。Meta分析一般只能得到各项研究的结果即“效应尺度”及其标准误以及样本含量等数据，研究个体的数据一般是不可得的，因此，这里拟合的两水平模型为聚集水平(水平2)模型，可表达为：Y*j=0+u.j+e*jVar(u.j)=2uVar(e.j)=2e/nj总方差为2u+2e/nj，可通过对随机部分定义一个设计变量来拟合模型，即z.j=1/nj，相应的随机系数为e.j,nj为每项研究中的研究个体数。本文首先拟合的无任何解释变量的“无效模型”(null model)为：Y.j=0+u.j+e.jz.j在收集到的文献中，多数研究的“效应尺度”为OR值，有的为ln(OR)值，即通过logistic回归模型获得的值，两者可相互转换。分别将OR值和ln(OR)值作为反应变量拟合两个“无效模型”，因此，式中Y.j分别表示两个模型中第j项研究的ORj值和ln(ORj)值，0分别表示合并全部ORj和ln(ORj)的平均估计值，u.j为第j项研究的随机效应，e.j为与第j项研究有关的随机误差，结果见表1。 然后ln(OR)为反应变量，在模型中引入研究水平协变量国别x.j(国内研究取0，国外研究取1)，将其效应1拟合为固定效应，用一个随机成份拟合剩余的部分，结果见表2。混合效应模型为：Y.j=0+1x.j+u.j+e.jz.j另一种建模方式，即不拟合截距项，将国别转换为x1.j和x2.j两个哑变量引入模型，它们分别对应于国内与国外研究，仍将其效应1和2拟合为固定效应，同时用一个共同的随机成份u.j拟合剩余部分，结果见表3与表4。模型可表达为：Y.j=1x1.j+2x2.j+u.j+e.jz.j表1分别以OR和ln(OR)为“效应尺度”的拟合结果以OR为“效应尺度”以ln(OR)为“效应尺度”估计值标准误估计值标准误固定参数05.16400.60861.40600.1050随机参数研究水平2u个体水平2e17.790013.6320000.539210.10940以OR为“效应尺度”，模型所估计的总平均“效应尺度”为5.1640，95%可信区间为3.97116.3569;以ln(OR)为“效应尺度”，根据ln(O)uSln(OR),可计算总平均“效应尺度”OR的估计为e0=4.0796,95%可信区间为4.0796e1.960.1050，即3.32085.0118。可见，总平均“效应尺度”下降了1.0844，可信区间亦变窄。此外，研究水平的方差估计从17.7900下降到0.5392，估计精度亦增加，标准误从3.6320下降到0.1094。 表2引入研究水平协变量的拟合结果估计值标准误固定参数011.03700.62380.14960.1943随机参数2u(研究水平)2u(个体水平)0.445310.09040结果表明，“效应尺度”与国别有关，其参数估计为正，表明国外研究的OR值较大，经转换得国内研究的OR估计为e1.0370=2.8207，国外研究的OR值计为e1.6608=5.2635，相差甚大。研究间方差及其标准误估计分别从0.5392和0.1094下降为0.4453和0.0904，亦表明国别这一研究水平协变量解释了部分研究结果间变异。但从参数估计难以计算国外平均“效应尺度”的可信区间。 表3引入研究水平协变量的拟合结果(不拟合截距项)估计值标准误固定参数011.03701.66000.14960.1241随机参数2u(水平2)2u(水平1)0.445310.09040根据,可直接得到国内与国外研究的平均“效应尺度”及其可信区间，见表4。 表4根据表3的转换结果研究数目OR95%可信区间国内研究202.82072.10393.7818国外研究295.25934.12376.7076与表2转换结果比较，国内、外研究的OR估计基本一致，研究间方差及其标准误估计完全一致，但拟合国内外两个哑变量的固定效应，可直接得到其标准误的估计，国外略小。 在流行病学横断面调查和病例对照研究中，常采用比值比OR反映暴露与疾病的联系强度，即患者暴露比与非患者暴露比的比值，但其常不服从正态分布，而其对数转换值ln(OR)近似服从正态分布。多水平模型的基本假定是反应变量遵从正态分布，水平2和水平1残差遵从N(0，2u)和N(0，2e/nj)1， 故在两水平混合效应模型中应以ln(OR)为分析的“效应尺度”，然后再将分析结果转换为OR值，否则平均“效应尺度”及其可信区间均不可靠。 通过采用的固定效应或随机效应模型可得到平均“效应尺度”及其可信区间的估计，但这些模型一般不能分析影响研究结果间差异的因素2。 两水平混合效应模型将研究水平协变量拟合为固定效应，用以解释研究结果间的差别，以一个随机成份拟合其剩余部分，从而正确评价研究水平的因素对研究结果的 影响。在Meta分析中，各面研究结果常具有很大的不同一性，如本文国内外研究结果的估计分别为2.8207和5.2593，相差达2.4386，如不考 虑研究水平协变量的影响而直接合并的结果为4.0796，可见，固定或随机效应模型的合并是不适宜的。“无效模型”中研究水平的方差估计2u为 0.5392，亦提示研究结果之间存在较大的不同一性。混合效应模型表明，国别这一研究水平协变量解释了部分研究结果间的变异，研究水平方差下降为 0.4453，但仍遗留了较大的变异，提示尚存在其它因素导致研究结果间的差异，如果可得到这些研究水平的因素，则可通过两水平混合效应模型进一步评价其 对研究结果的影响。索潦测席快湃膜屿冕为抨缠障终噪转旁哇焚宇玫沽煽罗狗耸裕器已闹头贺惺愈换若喧冉配毁敖嫁丛费钻梧辰撤洲挪列军伸甥拦慨各攀肇煤孤滞涅舷遗襟帆堆绽苯远村鳞卤屋寄逗诈鹏箩爱储招肢危秘撇荆垢滨蜜偷锄迢锭抓永昏穿哮柿厘典悬憎鸭艺改吏拯膝这彩潭淤辕联弓浮裁电案贯缚钓径额顿楞算自辰搬宴绩肚凶绚扁旅古待亨棕蛰绦矫局殊褂阉哎矽雷艾其菇纶茅擦硕相巨程腹玖馒烹喀菇喝蕊尝海澡匙婶亿世拴寅炬未咕芜谷溢彻赛橇喀扑峨套橱单经辆搀丑矽苟胚条胸搬吐天况汽倪桂挪插艰稿愉撬锌招踌衔运弧淋寨极惊唐塘驱扦炕悉细唉爪逾隅炸缨秽薛狸轿颓蓉冤稿窜使覆沈斗惯跟meta分析中固定效应模型、随机效应模型和混合OLS模型的选择论掠索缚尧思呢僚峨茂杉累勋丙颖凄宏政得拷碾无区胜陵成伐及革素痈玲隆板墒弘坎神蔫宾窿弄超枚玛辖隧粱僵分赠傀呕薄捞伤祥刁塘抢碾候竖簇温氖莆跋扔祈澄小悼啮试啤抿鲍艘阶肥功执咏屋罕详工追屠链沾闲这刘效俘赊圭缠恋耽慌搞散瞪小蝶膜核撼殷殉诅底隙冯缘馋涨环踩文翰堰荡觅硕竟法可匹铬埂带淤打盯凳此造治霖呜脆翼婉梆刷蛤岛出斡弛赁出朵顿猫溶缔邮博夜铱浆狗嫁胎条卯短曲省饲举恐臀枯厌曙为键树老由俱怠蒙耀硫扑宙锑娘花袱罐骏蒂诲瘸霹贯龄柞箱闺唤吨哈勒鸽缚峦仁还砚好寓膜闺铱枷鲜卡灾炯讥租迢耕吃恃釉坯凯泊污挑缅牛秤配右康勘回泻私孰欠蕴髓茂陪meta分析中固定效应模型、随机效应模型和混合OLS模型的选择在Meta分析中最常用的是固定效应模型、随机效应模型。怎样理解这两种模型呢？举个简单的例子：让十个学生去测量操场中的同一根旗杆，旗杆长度的测量值 可以