张楠利用MonteCarlo方法研究匹配条件变化对.ppt

《张楠利用MonteCarlo方法研究匹配条件变化对.ppt》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、1利用MonteCarlo方法研究匹配条件变化对于条件逻辑回归结果的影响河北医科大学第四医院信息中心日本新潟大学医学部医疗信息研究室张楠背景条件逻辑回归分析模型（Conditionallogisticregressionmodel）是利用对照数据（Matchedcase-controldataset）基于逻辑回归模型（Logisticregressionmodel）进行统计分析的一种方法。条件逻辑回归模型在医学相关研究中得到广泛应用，经常被用于回顾研究和前瞻研究。数据匹配标准（1）1对照数据中的作为数据匹配的变量分布或者比例应当与病人组被用于匹配的变量分布应当近似。2被选为数

2、据匹配的变量应当为独立变量，不会受到其他变量的影响。例如:年龄，性别。数据匹配标准（2）数据对照1年龄性别癌症患者60男正常人60男数据对照2年龄性别癌症患者60男正常人58男在数据匹配过程中，正常人的数据通过随机方法取得。常见的数据匹配方法1在每一对数据中，病人组的匹配变量数值与对照组匹配变量数值网全相同。2在每一对数据中，病人组的匹配变量数值与对照组匹配变量数值存在一定程度的浮动（±1~±5）。3将匹配变量分组化，然后按照分组进行数据匹配。研究目的及条件逻辑回归模型我们设计了模拟程序用来发现当对照数据匹配标准发生变化时，条件逻辑回归模型的分析结果是如何变化的。分析的模型基

3、于如下模型：采用MonteCarlo方法的必要性Time=1因素OddsRatio95%ConfidencelimitsPvalue吸烟0.7780.489-1.2390.291饮酒4.5352.684-7.622<0.001饮茶0.7301.095-2.7430.019Time=2因素OddsRatio95%ConfidencelimitsPvalue吸烟0.8660.527-1.4210.568饮酒4.5662.714-7.684<0.001饮茶0.2571.464-3.820<0.001Time=3因素OddsRatio95%ConfidencelimitsPvalue

4、吸烟1.1670.743-1.8300.503饮酒3.5102.169-5.680<0.001饮茶0.4130.954-2.4470.078方法1利用模拟程序生成符合事先指定分布及比例的随机数，并根据匹配标准自动生成匹配数据。为了观察不同数据量下匹配标准变化的情况，我们模拟了4组不同数量的模拟数据。通过比较条件逻辑回归模型的系数以及P值等统计值观察匹配标准变化时对分析结果的影响。2我们将程序应用于真实的临床数据分析上，以观察对于匹配标准对于实际应用上的影响。本次研究的数据模拟以及条件逻辑回归分析在R下完成。模拟数据参数变量模拟病理组模拟对照组年龄（标准偏差）60.00±5.0

5、050.00±5.00性别（男/女）70%/30%50%/50%男性吸烟与否（吸/不吸）70%/30%40%/60%女性吸烟与否（吸/不吸）3%/97%1%/99%男性饮酒与否（饮/不饮）50%/50%30%/70%女性饮酒与否（饮/不饮）1%/99%1%/99%合计n*10n*n*=(150,250,350,450)模拟程序结果-回归参数变化趋势吸烟变量回归参数变化趋势饮酒变量回归参数变化趋势模拟程序结果-P值变化趋势吸烟变量P值变化趋势饮酒变量P值变化趋势随机数据匹配率模拟数据量匹配对数匹配率1509865.3%25017268.8%35024469.7%45031670

6、.2%数据匹配过程均由程序根据指定的匹配标准随机选择数据自动完成。真实数据的分布特点变量病例组对照组年龄（标准偏差）60.67±9.1851.03±9.74性别（男/女）420/1481064/639吸烟与否（吸/不吸）292/276654/1049饮酒与否（饮/不饮）211/357425/1278饮茶与否（饮/不饮）130/438476/1227受教育水平(0/1/2/3)43/108/175/124299/626/754/24合计5681703临床数据严格对应分析结果系数标准差95%CIORP-value吸烟0.7500.2071.25-3.752.1620.028饮酒0

7、.8910.2301.37-4.602.5060.012饮茶-1.6330.2291.13-1.460.200<0.001受教肓水平2.2040.1635.84-14.499.188<0.001匹配标准为±0时临床数据浮动匹配分析结果系数标准差95%CIORP-value吸烟0.3770.1520.93-2.331.4750.178饮酒0.7350.1751.28-3.492.1170.016饮茶-1.2770.2031.15-1.440.285<0.001受教肓水平1.3210.0882.97-4.7