多样本变异系数比较的似然比检验

目的 通过求非中心t分布未知参数最大似然估计,由此提出一种比较多个样本变异系数差异的似然比检验(LRT)方法,并与现有的D'AD检验在第Ⅰ类错误率与检验效能两方面进行模拟比较.方法 由于样本变异系数的似然函数跟非中心t分布有关,本文首先提出一种求非中心t分布未知参数最大似然估计的算法,然后由此构造出比较多个样本变异系数的似然比检验统计量,并求出其Bartlett校正系数,对似然比检验统计量进行Bartlett校正,以便该方法也能用于小样本的情形.计算机模拟时从正态分布总体中抽样,模拟LRT的第Ⅰ类错误率与检验效能,并与D'AD检验做比较.将LRT方法编写成R语言程序,输出校正后的统计量值和相应P值,以方便实际应用.结果 校正的LRT方法能更好地控制第Ⅰ类错误率,并且其检验效能比D'AD检验更稳健.在小样本且样本量不均衡的情形下,其检验效能比D'AD方法高.结论 校正的LRT方法适用范围更广,检验效能高,可为检验多样本变异系数间差异提供更加有效的方法.     

无删失生存数据Wilcoxon秩和检验与logrank检验的比较

目的 比较Wilcoxon秩和检验与logrank检验处理无删失生存数据的统计性能.方法 采用MonteCarlo模拟方法,应用SAS9.1.3软件编程,比较两种检验方法的Ⅰ型错误率和检验效能.结果 小样本情形(n≤50)下,Wilcoxon秩和检验的Ⅰ型错误率较logrank检验更接近检验水准0.05,而检验效能与logrank检验相差不大,后者检验效能的提高是以扩大Ⅰ型错误率为代价的.结论 对于无删失生存数据的比较,尤其是小样本情形(n≤50)下,建议采用Wilcoxon秩和检验.     

样本率多重比较方法的模拟研究

目的 采用Monte Carlo方法模拟研究样本率多重比较方法特性.方法 用SAS9.13软件编程,模拟估算不同参数条件下各种样本率多重比较方法的总Ⅰ型错误率、检验效能和判对率.结果 不同方法各具特色,能较好控制总Ⅰ型错误率且检验效能较高的方法有Ryan、Perm、Boot和Hommel、Hoch、StepS、StepB方法.结论 如需严格控制总Ⅰ型错误率在α水准内,Hommel、Hoch、StepS、StepB效能较高;如允许总Ⅰ型错误率略溢出(1.15α以内),Ryan、Perm和Boot法效能较高.     

一种临床试验中的适应性样本量调整方法

目的介绍一种临床试验中的适应性样本量调整方法,并探讨样本量调整后统计分析方法的第Ⅰ类错误率及检验效能。方法通过montecarlo模拟的方法研究n1大小对最终样本量Nf的影响,并估计最终方差偏移大小;同时模拟研究样本量调整后统计分析方法的第Ⅰ类错误率及检验效能大小。结果(1)模拟结果显示运用该样本量调整方法所得到的最终样本量Nf非常接近其真实值N0,尤其在π=0.4时进行样本量调整。(2)同时模拟结果显示所介绍的样本量调整后的校正t检验方法不仅能有效控制第Ⅰ类错误率α并且能充分满足试验检验效能(1-β)。结论该样本量调整方法研究结果是在一般两样本单侧t检验条件下得到也可应用于优效或非劣效设计的临床试验中。     

贝叶斯方法在篮子试验设计中的应用

目的探讨贝叶斯方法和Simon二阶段方法应用于篮子试验设计的优劣及实例分析判断试验疗效。方法模拟研究针对二分类结局变量,考虑不同试验层数、期中分析纳入人数等参数设置,比较两种方法的篮子试验设计在检验效能、控制I错误率及样本量的区别,并通过实例来介绍其应用。结果模拟试验结果表明,除异质性较高的2种情景外,贝叶斯方法的篮子试验设计效能均高于应用Simon二阶段方法,且对假阳性率控制更为严格,代价便是无效篮子增多时检验效能降低。贝叶斯方法的篮子试验的样本量总体较少,尤其在试验层数较多和同质情况下。结论贝叶斯方法的篮子试验更加灵活,可借用层间信息,相比Simon二阶段方法,更严格地控制整体一类错误率,且在试验层数和期中分析次数多时更节省平均总样本量。     

常见定量资料应用基于秩转换类参数方法的可行性北大核心CSCD

目的针对配对样本(单样本)设计、两个和多个独立样本设计的定量资料,分析参数方法、非参数方法以及基于秩转换类参数方法在资料不满足参数方法条件下的适用情形。方法介绍基于秩转换类参数方法的原理及其与非参数方法的关系,采用Monte-Carlo模拟方法,考虑正态和左偏态两种分布,方差齐与不齐两种情形,比较三种方法的I类错误率和检验效能。结果左偏态分布时,无论方差是否齐性,或不涉及方差齐性(单样本设计),参数方法的I类错误率偏离设定水准且明显大于非参数方法和基于秩转换类参数方法,而检验效能明显低于其他两种方法。方差不齐且正态分布时,参数方法的统计性能明显优于其他两种方法。非参数方法和基于秩转换类参数方法在不同资料类型下的统计性能相近。结论基于秩转换类参数方法与非参数方法条件性能相近,适用于非参数方法处理的数据。     

实验组与对照组比较的样本率多重比较方法的模拟研究

目的研究5种实验组与对照组比较的样本率多重比较方法之特性。方法用SAS软件编制计算机模拟抽样程序,模拟估算各种参数条件下5种方法的总Ⅰ型错误率、检验效能和判对率。结果5种方法各具特色:Brunden法保守,杜养志法灵敏,Dunnett-SNK法和Bootstrap法较稳健,Dunnett-Levy法逊于Dunnett-SNK法。结论5种方法中没有一种方法能在各种条件下都将总Ⅰ型错误率控制在α水平内。相对而言,Dunnett-SNK法最稳健,较其他方法更值得推荐。如需严格控制总Ⅰ型错误率且组数不超过5,可选择Brunden法。如希望结果灵敏,可考虑杜养志法和Bootstrap法。     

常见定量资料应用基于秩转换类参数方法的可行性

目的针对配对样本(单样本)设计、两个和多个独立样本设计的定量资料,分析参数方法、非参数方法以及基于秩转换类参数方法在资料不满足参数方法条件下的适用情形。方法介绍基于秩转换类参数方法的原理及其与非参数方法的关系,采用Monte-Carlo模拟方法,考虑正态和左偏态两种分布,方差齐与不齐两种情形,比较三种方法的I类错误率和检验效能。结果左偏态分布时,无论方差是否齐性,或不涉及方差齐性(单样本设计),参数方法的I类错误率偏离设定水准且明显大于非参数方法和基于秩转换类参数方法,而检验效能明显低于其他两种方法。方差不齐且正态分布时,参数方法的统计性能明显优于其他两种方法。非参数方法和基于秩转换类参数方法在不同资料类型下的统计性能相近。结论基于秩转换类参数方法与非参数方法条件性能相近,适用于非参数方法处理的数据。     

基于生存资料的Ⅱ/Ⅲ期无缝设计期中分析方法评价

目的基于生存指标的无缝设计,考查PFS在什么条件下可以用于Ⅱ/Ⅲ期无缝设计期中分析时的剂量组筛选。方法采用模拟试验的手段,分别研究基于生存资料下Fisher合并法和加权逆正态法在利用和不利用OS的信息两种组别筛选策略下的总Ⅰ类错误率和检验效能。结果加权逆正态合并法与Fisher合并法相比,两者在控制总Ⅰ类错误率方面无明显差异,随着PFS和OS相关性增强,加权逆正态合并法的检验效能要高于Fisher合并法,并且合并PFS和OS的信息将获得更高检验效能。结论当PFS和OS的相关性较大,且效应趋势一致时,建议采用加权逆正态合并法合并两者的信息进行Ⅱ/Ⅲ期无缝临床试验的统计分析。     

IPCW法在肿瘤新药研发中的应用探讨

目的 应用数据模拟方法,比较IPCW法与传统分析方法(意向性分析和censored方法)检验效能与假阳性率的差异,验证IPCW法的分析转组数据效果.方法 应用Monte Carlo模拟,产生试验组和对照组的基线数据、实验室检查数据及生存时间,比较三种方法的检验效能及假阳性率,从而评价IPCW法处理转组数据的优劣.结果 模拟结果显示,两组总生存时间不等时,IPCW法检验效能优于意向性分析;两组总生存时间相等时,IPCW法假阳性率与意向性分析相当.虽censored方法检验效能较好,但其假阳性率较高.结论 与意向性分析及censored方法相比,IPCW法在分析转组数据时更佳.     

微阵列数据的多重比较

目的 介绍阳性结果错误率（FDR）及相关控制方法在微阵列数据多重比较中的应用。方法 用BH、BL、BY和ALSU四种FDR控制程序比较了3226个基因在两组乳腺癌患者中的表达差异。结果 四个程序在各自实用的范围内均将FDR控制在0．05以下，检验效能由大到小的顺序为：ALSU〉BH〉BY〉BL。ALSU程序因引入m0的估计,更为合理。不仅提高了检验效能，同时又较好地控制了假阳性错误。结论 在微阵列数据的比较中必须考虑FDR的控制，同时又要考虑提高检验效能。多重比较中，控制FDR比控制总Ⅰ型错误率（FWER）检验效能高，且更为实用。     

三臂非劣效性设计生存时间数据的半参数统计推断方法

目的 推荐并改进一种用于新药临床试验中生存分析三臂非劣效性设计的半参数统计推断方法.方法 采用Cox比例风险模型建立用于生存时间数据的三臂非劣效检验方法,通过模拟试验验证该方法的有效性,并给出样本含量估计的R语言程序.结果 模拟试验结果证实,在样本删失风险率维持在一定水平时,三臂Cox模型能有效地控制Ⅰ类错误发生概率,并保证较高的检验效能,实际应用简单有效.结论 本文推荐的方法,可以有效地应用于三臂生存时间数据的非劣效性研究.     

析因设计资料的单独效应分析

目的探讨析因设计资料单独效应的分析方法。方法分析单独效应的两种方法分别为基于析因分析的整体标准误和基于单因素分析的局部标准误方法,前者使用SPSS软件中一般线性模型下Lmatrix子句实现。用模拟方法分别比较两种方法的I型错误率和检验效能。结果使用两种不同标准误,统计分析结果存在较大差异。模拟结果显示,两种方法的I型错误率非常接近,基于析因方差分析标准误的检验效能较高。结论基于析因分析整体标准误较基于单因素局部标准误分析单独效应更为合理。     

含安慰剂组三臂临床试验基于bootstrap再抽样的非劣效评判

目的探讨含安慰剂组三臂临床试验基于bootstrap再抽样的非劣效评判的方法。方法用Monte Carlo模拟方法,产生服从正态分布、对数正态分布和Gamma分布的随机样本,进行Welch校正t检验法和bootstrap法的α-模拟和power模拟的验证和比较。结果当数据服从正态分布,在样本量较大时,Welch校正t检验法和bootstrap法均表现出较好的统计性能,但当数据呈偏态分布时,Welch校正t检验法的第一类错误率会偏离预先给定的α-水平,而bootstrap法在样本量较大时,第一类错误率基本保持在预先给定的水平。Welch校正t检验法和bootstrap法的power模拟结果基本相同。结论含安慰剂组的三臂临床试验在数据不服从正态分布时,bootstrap法可作为一种有效的非劣效评判方法。     

存在缺失数据的非劣效临床试验中不同分析方法的模拟比较研究

目的探讨适用于非劣效试验设计下重复测量数据的统计分析方法。方法基于某临床研究的基本指标,应用Monte Carlo模拟方法,分别模拟出完整的重复测量数据和不同缺失机制下不同缺失比例的数据;运用不同的统计分析方法对模拟数据进行统计分析,通过检验效能、Ⅰ类错误等指标比较不同分析方法的性能。结果完整数据下,多数的统计分析方法均能达到或近似达到较理想的检验效能和Ⅰ类错误;存在数据缺失的情况下,基于非结构化协方差矩阵的重复测量混合效应模型（mixed-effect model repeated measure,MMRM）和不同协方差结构的广义估计方程（generalized estimated equation,GEE）表现出良好的稳定性,能在保证足够检验效能的情况下,有效的控制Ⅰ类错误;非劣效试验设计下,不同缺失比例的模拟结果均显示,末次结转（last observation carried forward,LOCF）的方法会低估组间的差异,造成Ⅰ类错误膨胀。结论基于非劣效试验设计的重复测量数据,当数据存在缺失时,LOCF结转方法会低估组间的差异,造成I类错误的膨胀,不再是一种保守的缺失值处理...     

Cochran-Armitage趋势检验和线性回归在流行病学率的趋势分析中的比较研究

采用Cochran-Armitage趋势（CAT）检验和线性回归（LR）两种统计学方法,分析天津市居民1999-2013年急性心肌梗死（AMI）发病率的时间趋势,比较两种方法结果的差别。结果显示,以实际人口为基础,无论是总的发病率时间趋势,还是按年龄分亚组后的时间趋势,CAT检验的统计效能均大于LR（CAT P值P值）。将实际人口缩小100倍,AMI发病率保持不变,CAT检验的P值呈明显升高,LR的P值保持不变。两种统计学方法在流行病学率的趋势分析中都有其各自的优势和缺点,研究者可以根据资料选择拟合度较好的模型,或者两种方法结合来分析,全面地描述统计结果。     

贝叶斯响应变量适应性随机化模拟预测评价方法初探

目的探索贝叶斯响应变量适应性随机化的模拟预测评价方法,为研究者在临床试验中合理设计随机化方案提供借鉴。方法通过不同参数组合成不同场景,使用R语言编程实现试验模拟过程,构造评价指标,结合模拟试验结果对不同场景作出评价。结果在结局指标为二分类高优变量的两臂临床试验中,可以通过调整分配概率更新频率、固定随机化阶段长度、分配概率算法调节参数及各处理组的先验概率分布,获得理想的检验功效及伦理性。在试验前可以通过提高试验成功界值、降低分配概率更新频率、增加均衡分配期长度三种方法降低试验整体一类错误率,以满足监管部门的要求。结论本研究提出的贝叶斯响应变量适应性随机化模拟预测评价方法可通过设定场景参数模拟试验,从控制一类错误、保证检验效能、提升伦理性等方面较全面评价随机化方案的合理性。     

偏态分布下多样本变异系数比较的平方秩检验

目的提出一种适用于偏态分布的多样本变异系数比较的平方秩检验,并与现有的D’AD与SRT方法进行模拟比较。方法借鉴Levene检验的构造思想,对原始数据做一种新的变换,由此构造出一种适用于偏态分布的多样本变异系数比较的MSRT检验统计量。在不同参数设置下,经过Monte Carlo模拟,比较新提出方法与现有方法的第一类错误率。同时,将新提出方法与现有方法应用于铅致肝损伤指标的数据,对这些方法做进一步的比较。结果 MSRT方法可以避免D’AD方法对样本正态分布假设的依赖;相比于SRT方法,在偏态分布下,MSRT方法可以更好地控制第一类错误率。在对铅致肝损伤指标γ-GT与ALT的变异系数的比较分析中,三种方法对两个指标的分析结果存在差异。结论MSRT方法是一种稳健且普遍适用的可靠方法。     

病例-同胞对照设计统计方法检验效能的比较研究

目的比较病例-同胞对照设计三种统计方法的检验效能。方法应用Monte-Carlo模拟方法,使用Slink软件产生的模拟家系资料,比较同胞-传递不平衡检验(s-TDT)、同胞组不平衡检验(SDT)和同胞对关联检验法(SAT)三种统计方法的检验效能。结果(1)SDT方法的I型误差最接近检验水准α=0.05;(2)在三种遗传模式中,s-TDT方法的效能最高,SAT方法的效能最低。结论病例-同胞对照设计使用患者的同胞作为对照,避免了人群分层现象,可以用来评价遗传因素与疾病之间的关联。     

疾病预后研究的中介分析方法评价

目的对用于癌症预后研究的五种中介分析方法(VanderWeele法、Baron-Kenny法、Imai法、Sobel法和InverseWeight法)进行评价,为实例分析的方法选择提供依据。方法基于模拟试验,产生不同参数设置下的模拟数据,并评价五种方法的第一类错误、检验效能和分析时间。结果除InverseWeight法在相关系数较大时第一类错误有所膨胀外,其余四种方法的第一类错误在不同参数情况下均在0.05附近。五种方法的检验效能趋势一致,均随着样本量、中介比、总效应的增大而增大,随着删失比的增大而减小。在样本量较小(N=100)且中介比不大于30%的情况下,InverseWeight法的检验效能低于另四种方法。InverseWeight法、Baron-Kenny法和Imai法的分析效率远低于VanderWeele法和Sobel法。结论综合考虑一类错误控制、检验效能及分析效率,推荐VanderWeele法进行预后研究的中介分析。