Random weighting method for Cox’s proportional hazards model

Variance of parameter estimate in Cox’s proportional hazards model is based on asymptotic variance. When sample size is small, variance can be estimated by bootstrap method. However, if censoring rate in a survival data set is high, bootstrap method may fail to work properly. This is because bootstrap samples may be even more heavily censored due to repeated sampling of the censored observations. This paper proposes a random weighting method for variance estimation and confidence interval estimation for proportional hazards model. This method, unlike the bootstrap method, does not lead to more severe censoring than the original sample does. Its large sample properties are studied and the consistency and asymptotic normality are proved under mild conditions. Simulation studies show that the random weighting method is not as sensitive to heavy censoring as bootstrap method is and can produce good variance estimates or confidence intervals.     

沪深股市大盘指数收益率分布尾部的实证研究

通过H ill估计的改进方法对上证综合指数和深圳成分指数的收益率分布的尾部指数进行了参数估计,用χ2检验验证了指数的稳定性及其置信区间.在此基础上提出用尾部指数估计尾概率,达到风险控制的目的.实证研究表明,沪深大盘指数收益率分布具有肥尾的特征,但并不服从无限方差分布.     

Generalized Levinson-Durbin sequences, binomial coefficients and autoregressive estimation

For a discrete time second-order stationary process, the Levinson-Durbin recursion is used to determine the coefficients of the best linear predictor of the observation at time k+1, given k previous observations, best in the sense of minimizing the mean square error. The coefficients determined by the recursion define a Levinson-Durbin sequence. We also define a generalized Levinson-Durbin sequence and note that binomial coefficients form a special case of a generalized Levinson-Durbin sequence. All generalized Levinson-Durbin sequences are shown to obey summation formulas which generalize formulas satisfied by binomial coefficients. Levinson-Durbin sequences arise in the construction of several autoregressive model coefficient estimators. The least squares autoregressive estimator does not give rise to a Levinson-Durbin sequence, but least squares fixed point processes, which yield least squares estimates of the coefficients unbiased to order 1/T, where T is the sample length, can be combined to construct a Levinson-Durbin sequence. By contrast, analogous fixed point processes arising from the Yule-Walker estimator do not combine to construct a Levinson-Durbin sequence, although the Yule-Walker estimator itself does determine a Levinson-Durbin sequence. The least squares and Yule-Walker fixed point processes are further studied when the mean of the process is a polynomial time trend that is estimated by least squares.     

Poisson分布下基于鞍点逼近的慢性病风险差的置信区间构造

风险差是流行病学中重要的指标之一,常用来比较两种治疗或两种诊断的有效性.因此,风险差区间的精确估计对流行病病情的诊断以及治疗方案的选择有很重要的意义.结合Poisson抽样的优点以及慢性病发病周期长和发病率低的特点,利用鞍点逼近方法来构造了Poisson分布下风险差的置信区间.同时,通过实例和Monte Carlo模拟对传统的四种区间构造方法进行评价.模拟结果表明:在小样本情况下,鞍点逼近方法得到的置信区间大多数能保证覆盖率近似于期望的置信水平并且使得区间长度最短,是一种很好的置信区间构造方法.     

一类重尾风险因子的模拟及其投资高风险值和置信区间的估计

由于金融市场中的日周期或短周期对数回报率的样本数据多数呈现胖尾分布，于是现有的正态或对数正态分布模型都在不同程度上失效，为了准确模拟这种胖尾分布和提高投资风险估计及金融管理，本文引进了一种可根据实际金融市场数据作出调正的蒙特卡洛模拟方法．这个方法可以有效地复制金融产品价格的日周期对数回报率数据的胖尾分布．结合非参数估计方法，利用该模拟方法还得到投资高风险值以及高风险置信区间的准确估计。     

Self-crossing Points of a Line Segment Process

This paper is devoted to planar stationary line segment processes. The segments are assumed to be independent, identically distributed, and independent of the locations (reference points). We consider a point process formed by self-crossing points between the line segments. Its asymptotic variance is explicitly expressed for Poisson segment processes. The main result of the paper is the central limit theorem for the number of intersection points in expanding rectangular sampling window. It holds not only for Poisson processes of reference points but also for stationary point processes satisfying certain conditions on absolute regularity (β-mixing) coefficients. The proof is based on the central limit theorem for β-mixing random fields. Approximate confidence intervals for the intensity of intersections can be constructed.     

区间数据情形下线性模型的经验似然推断

§1Introduction Instatisticalapplications,weoftenencounterintervalcensoreddatawhenafailure timeYcannotbeobserved,butcanonlybedeterminedtolieinanintervalobtainedfroma sequenceofexaminationtimes.Forinstance,themaximumdosagewhichpatientscan endureisconcerned.LetYibethemaximumdosagewhichtheithpatientcanendure,Ui,j(j=1,2,...,k)bethedosagewhichthepatienthasbeentested.ItisobviousthatYiis unobservable.SupposetheithpatientisnormalwhenthedosageisUi,j,andhe(orshe)is abnormalwhenthedosageisUi,j+1.Then…     

Panel模型中的精确检验和置信区间

利用广义p-值和广义置信区间的概念,研究了Panel模型中未知参数的检验和置信区间问题.对于回归系数,分别考虑了单个情形和多个线性无关情形下的检验和置信区间问题,得到了精确检验和置信区间.对于方差分量,研究了其任意线性组合的检验和置信区间问题,建立了精确检验和置信区间.基于广义p-值和广义置信区间,获取精确检验和置信区间的方法具有计算方便、易应用于小样本问题的特点.最后,分别从理论和数值上研究了这些精确检验和置信区间的统计性质.     

Complementary moments and the summation of series

In previous papers the author has treated the summation of power and other series whose coefficients re moments of a function f(u), usually over the interval 0 </ u </ 1. The present paper defines “complementary moments” of f(u) for the same interval, and demonstrates their use in convergence acceleration and in the analytic summation of series having somewhat complicated coefficients of a rather special form. Simple algorithms are established, and these enable one, for example, to obtain exact Laplace-transform inverses of unusual functions of the Laplace-transform operator p.     

Empirical likelihood analysis of the Buckley-James estimator

The censored linear regression model, also referred to as the accelerated failure time (AFT) model when the logarithm of the survival time is used as the response variable, is widely seen as an alternative to the popular Cox model when the assumption of proportional hazards is questionable. Buckley and James [Linear regression with censored data, Biometrika 66 (1979) 429-436] extended the least squares estimator to the semiparametric censored linear regression model in which the error distribution is completely unspecified. The Buckley-James estimator performs well in many simulation studies and examples. The direct interpretation of the AFT model is also more attractive than the Cox model, as Cox has pointed out, in practical situations. However, the application of the Buckley-James estimation was limited in practice mainly due to its illusive variance. In this paper, we use the empirical likelihood method to derive a new test and confidence interval based on the Buckley-James estimator of the regression coefficient. A standard chi-square distribution is used to calculate the P-value and the confidence interval. The proposed empirical likelihood method does not involve variance estimation. It also shows much better small sample performance than some existing methods in our simulation studies.     

Exact confidence coefficients of simultaneous confidence intervals for multinomial proportions

Simultaneous confidence intervals for multinomial proportions are useful in many areas of science. Since 1964, approximate simultaneous 1-α confidence intervals have been proposed for multinomial proportions. Although at each point in the parameter space, these confidence sets have asymptotic 1-α coverage probability, the exact confidence coefficients of these simultaneous confidence intervals for a fixed sample size are unknown before.In this paper, we propose a procedure for calculating exact confidence coefficients for simultaneous confidence intervals of multinomial proportions for any fixed sample size. With this methodology, exact confidence coefficients can be clearly derived, and the point at which the infimum of the coverage probability occurs can be clearly identified.     

基于区间数的证券组合投资模型研究

提出了证券组合投资的区间数线性规划模型.通过引入区间数线性规划问题中的目标函数优化水平α和约束水平β将目标函数和约束条件均为区间数的线性规划问题转化为确定型的线性规划问题.投资者可以根据自己的风险喜好程度和客观情况,对这两个参数做出不同的估计,从而得到相应情况下的有效投资方案,使证券组合投资决策更具柔性.最后通过实例分析说明了该模型的可行性.     

Technology scoring model for reflecting evaluator’s perception within confidence limits

The main purpose of this study is to propose a new technology scoring model for reflecting the total perception scoring phenomenon which happens often in many evaluation settings. A base model used is a logistic regression for non-default prediction of a firm. The point estimator used to predict the probability for non-default based on this model does not consider the risk involved in the estimation error. We propose to update the point estimator within its confidence interval using the evaluator’s perception. The proposed approach takes into account not only the risk involved in the estimation error of the point estimator but also the total perception scoring phenomenon. Empirical evidence of a better prediction ability of the proposed model is displayed in terms of the area under the ROC curves. Additionally, we showed that the proposed model can take advantage when it is applied to smaller data size. It is expected that the proposed approach can be applied to various technology related decision-makings such as R&D investment, alliance, transfer, and loan.     

The structure of unrevealed Bads in Good/Bad risk scores

In calculating risk scores for making predictions and decisions about loan defaults, it is common practice to base assessments on a population of individuals whose loans have not yet attained a final status or trapped state of Good (G: paid in full) or Bad (B: default, bankrupt, written off, no response, etc). When active accounts are examined prior to end of loan term, we describe them as Contaminated Goods (CG) because they contain some Bads that default at a later time. In such cases, one can easily misestimate or misinterpret the eventual population odds and scores because the CG to B odds at any point in time is larger than G to B at the end of the loan. It is shown that if the risk score is a sufficient statistic and if the Information Odds score for Goods at the end-of-term is normal with variance σ² in a population of terminated loan accounts, then so also is the conditional score distribution for Bads; surprisingly, the theoretical means are ±0.5σ². When active accounts are contaminated by unrevealed Bads not yet classified as such, the conditional score distribution is a mixture of normal distributions with a variance larger than σ²; thus, variances of Active (CG) and Bad (B) accounts are unequal and the log of fitted odds versus score is convex, departing from the traditional assumption of a linear fit.     

Correlation and the time interval in multiple regression models

In this paper we investigate the time interval effect of multiple regression models in which some of the variables are additive and some are multiplicative. The effect on the partial regression and correlation coefficients is influenced by the selected time interval. We find that the partial regression and correlation coefficients between two additive variables approach one-period values as n increases. When one of the variables is multiplicative, they will approach zero in the limit. We also show that the decreasing speed of the n-period correlation coefficients between both multiplicative variables is faster than others, except that a one-period correlation has a higher positive value. The results of this paper can be widely applied in various fields where regression or correlation analyses are employed.     

Sequential Estimation for a Functional of the Spectral Density of a Gaussian Stationary Process

Integral functional of the spectral density of stationary process is an important index in time series analysis. In this paper we consider the problem of sequential point and fixed-width confidence interval estimation of an integral functional of the spectral density for Gaussian stationary process. The proposed sequential point estimator is based on the integral functional replaced by the periodogram in place of the spectral density. Then it is shown to be asymptotically risk efficient as the cost per observation tends to zero. Next we provide a sequential interval estimator, which is asymptotically efficient as the width of the interval tends to zero. Finally some numerical studies will be given.     

随机设计非线性混合模型的统计分析

本文研究了个体观察次数为随机的非线性 混合效应模型中参数的点估计以及区间估计. 在仅给出适当的矩条件下, 给出了固定效应、随机效应的方差阵以及误差方差的矩估计, 并证明了估计量的相合性及渐近正态性. 为给出误差方差以及随机效应方差分量的置信区间, 本文也给出了误差及随机效应的四阶矩估计. 随机模拟说明了方法的有效性.     

证券组合投资的多目标区间数线性规划模型

本文提出了证券组合投资的多目标区间数线性规划模型，引入了收益——风险偏好参数和优化水平参数。投资者可以根据对风险的喜好程度和金融市场的客观情况，适当估计这两个参数，从而得到相应情况下的有效投资方案，使投资过程更具柔性，而且更接近于实际情况。     

A note on interval estimation of the standard deviation of a gamma population with applications to statistical quality control

Traditional process control charts for a measurement standard deviation are based on the assumption of normality, which may not always be valid. Assuming that measurements follow a gamma distribution, we have obtained an approximate distribution of the sample variance, scaled appropriately. This approximate distribution, which happens to be another gamma model, is used to derive an interval estimate of the population standard deviation. Further, the above approximate gamma model for the sample variance can be used to develop a process control chart as demonstrated by a simulated data set.     

Fourier expansion along geodesics on Riemann surfaces

For an eigenfunction of the Laplacian on a hyperbolic Riemann surface, the coefficients of the Fourier expansion are described as intertwining functionals. All intertwiners are classified. A refined growth estimate for the coefficients is given and a summation formula is proved.