病例队列研究下带约束的Cox模型中参数的一种加权估计方法

许多大型队列研究的主要预算和成本通常来自昂贵的关键协变量的采集与测量.在有限的预算或者时间下,观测大型队列中所有研究对象的昂贵协变量往往是不可行和低效的.因此,研究人员一直致力于寻找和使用能节约成本并能达到预设效率的抽样设计方法.对于生存数据,病例队列设计正是这样一种具有成本效益的有偏抽样机制.进一步,在病例队列研究中,为了利用更多的数据先验信息来提高研究的效率,可以在统计建模过程中对模型参数进行合理的假设和约束.本文研究病例队列设计下带约束的Cox模型中参数的估计方法.我们提出了一种加权约束估计的方法,并建立了所提出估计的渐近理论.发展了一种新的约束MM算法来实现所提出的加权约束估计的数值计算.通过统计模拟研究评估了所提出方法在有限样本量下的表现.分析了一个肾母细胞瘤的实际数据来展示所提出方法的实际应用价值.     

具有随机删失数据的部分线性模型的拟合优度检验的渐近性质

研究随机删失部分线性回归模型的假设检验问题. 提出了一个检验统计量来检验数据是否满足一个部分线性回归模型, 它是基于残差的cusum过程的平方形式. 研究了零假设下和局部对立假设下检验统计量的渐近分布. 数值模拟表明该检验方法有好的检验功效.     

比例风险模型下病例队列设计中两种加权估计方法及其应用

对于大型队列研究或观察型研究,基于生存数据的病例队列设计是一种能有效节约成本和提高效率的抽样机制.这种抽样设计仅对一个随机抽取的子队列以及子队列之外所有经历了感兴趣事件的病例个体进行关键协变量的测量,具有显著的成本效益.本文研究如何应用比例风险模型拟合病例队列研究数据.探讨逆概率加权和与时间相关加权这两种基于加权估计方程的统计推断方法和其渐近性质等理论结果.通过一系列的统计模拟研究展示了病例队列设计的优良性以及相较于传统简单随机抽样设计的高效性.进一步,应用这两种推断方法分析了两个实际数据,展示了其在实际中的应用价值和前景.     

病例-队列设计下长度偏差数据的比例均值剩余寿命模型的统计推断

在大型队列研究中,病例-队列设计是一种可以有效节约成本的试验设计方法.本文研究了在病例-队列设计下,基于长度偏差数据的比例均值剩余寿命模型的统计推断问题,提出了一种带有时间相依权重的加权混合估计方程方法来估计模型中的回归系数,并证明了在适当条件下,所得到的估计量具有相合性与渐近正态性.模拟结果表明本文所提出的方法在有限样本下的表现不错.最后,我们将所提出的方法应用到了一组实际数据中.     

基于因变量抽样设计下线性回归模型的假设检验问题

在许多大型队列研究中,采用节约成本并能提高效率的抽样机制至关重要,基于因变量的抽样设计正是这样一种有偏抽样机制.这种方法最大的优点在于:能够将资源集中在那些包含有更多的协变量与因变量关系信息的研究群体上.本文研究基于因变量抽样设计下的线性模型中回归方程显著性检验以及回归系数显著性检验问题.基于一种半参数经验轮廓似然的方法,我们分别为回归方程检验与回归系数检验提出了相应的检验统计量,获得了所提出检验统计量的渐近性质.通过模拟研究评估了所提出的检验方法在有限样本下的表现,并应用提出的方法分析了一个孕妇分娩的实际数据.     

基于Bootstrap的厚尾相依序列持久性变点检验

基于修正方差比率函数给出一种检验厚尾序列持久性变点的统计量.在无变点的假设下得到了统计量的渐近分布.为避免检验渐近分布中的厚尾指数,构造Bootstrap抽样方法来确定渐近分布的经验临界值.数值模拟研究结果说明修正方差比率统计量及Bootstrap抽样方法的有效性.     

线性模型中M检验原假设分布的随机加权逼近

在线性模型中M-方法可以用于线性假设检验, 其中M检验、Wald检验和Rao的计分型检验是最常用的检验准则. 但是在计算这些检验的临界值时都涉及到未知参数的估计. 在本文中我们利用随机加权的方法来逼近这些检验的原假设分布. 结果表明在原假设和局部对立假设之下随机加权统计量的渐近分布与原检验统计量在原假设之下的渐近分布相同. 因此我们不需要对冗余参数进行估计,利用随机加权的方法就可以得到这些检验的临界值. 而且在局部对立假设之下可以实现对功效的计算. 当取不同的误差分布和不同的随机权时, 我们对本文的方法进行了蒙特卡洛模拟. 结果表明用随机加权方法来逼近原假设分布是非常精确的.     

线性度量误差模型的序列相关检验

回归模型的序列相关检验是经济和金融数据分析中的一项重要的工作.基于最小二乘残差,作者提出了一个检验统计量以检验线性度量误差模型的误差序列是否存在序列相关性.在零假设下,得到了检验统计量的渐近分布.数值模拟结果表明,这里提出的检验统计量具有良好的有限样本性质.     

面板数据交互固定效应模型的方差分量检验

研究了面板数据交互固定效应模型中方差分量的检验问题.首先依据模型中误差项的估计构造辅助回归模型,然后根据该辅助回归构造检验统计量,对模型中的异方差性进行检验.进一步,通过构造不同的辅助回归模型和检验统计量可以判别异方差的来源.在一定正则条件下,得到了检验统计量在原假设和备择假设下的渐近分布,并说明所提出的检验方法不依赖于误差分布.最后,通过模拟研究对本文的检验方法进行评价,说明所提检验方法是有效的.     

一般偏差右删失数据下剩余寿命分位数回归

剩余寿命是刻画个体预期寿命的一个重要度量,对剩余寿命的早期研究主要集中在剩余均值上.然而当总体生存函数偏态或厚尾时剩余均值函数可能不存在,因此统计学者建议用剩余寿命分位数来刻画预期寿命.在完全数据和右删失数据下,剩余寿命分位数的建模和理论已经很完善.但是,在实际的调查研究中经常会遇到偏差抽样数据.例如,临床医学中的左截断数据,流行病学中的病例队列抽样数据,医学大型队列研究中的长度偏差抽样数据等等.忽略抽样偏差会导致参数估计有偏和不合理的推断结果.本文考虑一般偏差右删失数据下剩余寿命分位数回归的统计推断问题.首先,我们提出了一个一般偏差右删失数据下的剩余寿命分位数回归模型,并利用一般估计方程方法对模型中的参数进行了估计.针对已有文献常用的删失变量与协变量独立性假设,本文重点考虑了删失变量依赖于协变量场合.其次,由于估计量的渐近方差中涉及非参密度函数,在估计渐近方差时,本文采用Bootstrap方法.最后,数值模拟显示本文提出的方法有限样本性质表现很好.     

Additive hazards regression under generalized case-cohort sampling

Case-cohort design usually requires the disease rate to be low in large cohort study,although it has been extensively used in practice.However,the disease with high rate is frequently observed in many clinical studies.Under such circumstances,it is desirable to consider a generalized case-cohort design,where only a fraction of cases are sampled.In this article,we propose the inference procedure for the additive hazards regression under the generalized case-cohort sampling.Asymptotic properties of the proposed estimators for the regression coefcients are established.To demonstrate the efectiveness of the generalized case-cohort sampling,we compare it with simple random sampling in terms of asymptotic relative efciency.Furthermore,we derive the optimal allocation of the subsamples for the proposed design.The fnite sample performance of the proposed method is evaluated through simulation studies.     

Additive-multiplicative hazards model for case-cohort studies with multiple disease outcomes

Case-cohort design is an efficient and economical design to study risk factors for diseases with expensive measurements, especially when the disease rate is low. When several diseases are of interest, multiple case-cohort design studies may be conducted using the same subcohort. To study the association between risk factors and each disease occurrence or death, we consider a general additive-multiplicative hazards model for case-cohort designs with multiple disease outcomes. We present an estimation procedure for the regression parameters of the additive-multiplicative hazards model, and show that the proposed estimator is consistent and asymptotically normal. Large sample approximation works well in finite sample studies in simulation. Finally, we apply the proposed method to a real data example for illustration.     

Information and asymptotic efficiency of the case-cohort sampling design in Cox's regression model

Efficiencies of the maximum pseudolikelihood estimator and a number of related estimators for the case-cohort sampling design in the proportional hazards regression model are studied. The asymptotic information and lower bound for estimating the parametric regression parameter are calculated based on the effective score, which is obtained by determining the component of the parametric score orthogonal to the space generated by the infinite-dimensional nuisance parameter. The asymptotic distributions of the maximum pseudolikelihood and related estimators in an i.i.d. setting show that these estimators do not achieve the computed asymptotic lower bound. Simple guidelines are provided to determine in which instances such estimators are close enough to efficient for practical purposes.     

基于加速失效时间模型模拟生存时间应用于检测生存性状位点

对于考察预指定情形下的统计模型的性能、性质及适应性,模拟研究是非常重要的统计工具.作为生存分析中两个最受欢迎的模型之一,由于加速失效时间模型中的因变量是生存时间的对数,且此模型能够以线性形式回归带有易解释的参数的协变量,从而加速失效模型比COX比例风险模型更便于拟合生存数据.首先提出了关于带有广义F-分布的加速失效模型的模拟研究中生成生存时间的方法,然后给出了描述加速失效时间模型的误差分布和相应的生存时间之间的一般的关系式,并给出了广义F-分布是如何生成生存时间的.最后,为证实所建议模拟技术的性能和有效性,将此方法应用于检测生存性状位点的模型中.     

Weighted Kaplan-Meier Tests for Umbrella Alternatives

The problem of testing for umbrella alternatives in a one-way layout with right-censored survival data is considered. Testing procedures based on the two-sample weighted Kaplan-Meier statistics suggested by Pepe and Fleming (1989, Biometrics, 45, 497–507; 1991, J. Roy. Statist. Soc. Ser. B, 53, 341–352) are suggested for both cases when the peak of the umbrella is known or unknown. The asymptotic relative efficiency of the weighted Kaplan-Meier test and the weighted logrank test proposed by Chen and Wolfe (2000, Statist. Sinica, 10, 595–612) is computed for the umbrella peak-known setting where the piecewise exponential survival distributions have the proportional or crossing hazards, or the related hazards differ at early or late times. Moreover, the results of a Monte Carlo study are presented to investigate the level and power performances of the umbrella tests. Finally, application of the proposed procedures to an appropriated data set is illustrated.     

Asymptotics on Semiparametric Analysis of Multivariate Failure Time Data Under the Additive Hazards Model

Many survival studies record the times to two or more distinct failures on each subject. The failures may be events of different natures or may be repetitions of the same kind of event. In this article, we consider the regression analysis of such multivariate failure time data under the additive hazards model. Simple weighted estimating functions for the regression parameters are proposed, and asymptotic distribution theory of the resulting estimators are derived. In addition, a class of generalized Wald and generalized score statistics for hypothesis testing and model selection are presented, and the asymptotic properties of these statistics are examined.     

Empirical likelihood method for linear transformation models

Empirical likelihood inferential procedure is proposed for right censored survival data under linear transformation models, which include the commonly used proportional hazards model as a special case. A log-empirical likelihood ratio test statistic for the regression coefficients is developed. We show that the proposed log-empirical likelihood ratio test statistic converges to a standard chi-squared distribution. The result can be used to make inference about the entire regression coefficients vector as well as any subset of it. The method is illustrated by extensive simulation studies and a real example.     

左截断双删失数据下可加可乘风险率模型的估计

在生存分析中,可加可乘风险率模型常用来研究协变量对初始事件和终止事件之间持续时间的影响效应。在本文中,我们考虑了在初始事件存在部分区间删失,同时终止事件存在左截断右删失的情形下,持续时间的可加可乘风险率模型的估计问题。我们提出了一个两阶段估计过程来估计模型的回归参数。并通过模拟分析验证了估计的大样本性质。最后利用该方法分析了恶性黑色素瘤手术治疗数据。     

����ɾʧ�»������Ӻ����IJ���ģ�͵�ͳ�Ʒ���

In collecting clinical data, data would be censored due to competing risks or patient withdrawal. The statistical inference for censoring data is always based on the assumption that the failure time and censoring time is independent. But in practice the failure time and censoring time are often dependent. Dependent censoring make the job to deal with censoring data more complicated. In this paper, we assume that the joint distribution of the failure time variable and censoring time variable is a function of their marginal distributions. This function is called a copula. Under prespecified copulas, the maximum likelihood estimators for cox proportional hazards models are worked out. Statistical analysis results are carried by simulations. When dependent censoring happens, the proposed method will do better than the traditional method used in independent situations. Simulation results show that the proposed method can get efficient estimations.     

Adaptive Penalized Weighted Least Absolute Deviations Estimation for the Accelerated Failure Time Model

The accelerated failure time model always offers a valuable complement to the traditional Cox proportional hazards model due to its direct and meaningful interpretation. We propose a variable selection method in the context of the accelerated failure time model for survival data, which can simultaneously complete variable selection and parameter estimation. Meanwhile, the proposed method can deal with the potential outliers in survival times as well as heteroscedastic model errors, which are frequently encountered in practice. Specifically, utilizing the general nonconvex penalty, we propose the adaptive penalized weighted least absolute deviation estimator for the accelerated failure time model. Under some regularity conditions, we show that the proposed method yields consistent estimator and possesses the oracle property. In addition, we propose a new algorithm to compute the estimate in the high dimensional settings, and evaluate the practical utility of the proposed method through extensive simulation studies and two real examples.