基于GARCH族模型的中国旅游酒店板块指数收益率波动分析

运用GARCH族模型分析旅游酒店板块指数日收益率的波动特征,研究表明：旅游酒店板块收益率是一个平稳过程,其波动具有“聚集”现象和“非对称效应”。GARCH（2,1）模型比GARCH（1,1）模型更好地消除了收益率序列的异方差性;TARCH（2,1）模型的拟合效果最好;GARCH—M模型和非对称的CARCH（1,1）模型都不适用于描述收益率的波动特征。     

A note on joint mix random vectors

Abstract

This note studies the dependence of joint mix random vectors from the perspective of covariance matrix. We first provide two useful methods in simulations to construct joint mix for Normal distribution. Then, we propose to characterize joint mix by covariance matrix for general marginal distribution. We present some examples showing that our methodology could provide supplementary results to relevant studies in literature.     

Bayesian Inference for Contact Networks Given Epidemic Data

Abstract. In this article, we estimate the parameters of a simple random network and a stochastic epidemic on that network using data consisting of recovery times of infected hosts. The SEIR epidemic model we fit has exponentially distributed transmission times with Gamma distributed exposed and infectious periods on a network where every edge exists with the same probability, independent of other edges. We employ a Bayesian framework and Markov chain Monte Carlo (MCMC) integration to make estimates of the joint posterior distribution of the model parameters. We discuss the accuracy of the parameter estimates under various prior assumptions and show that it is possible in many scientifically interesting cases to accurately recover the parameters. We demonstrate our approach by studying a measles outbreak in Hagelloch, Germany, in 1861 consisting of 188 affected individuals. We provide an R package to carry out these analyses, which is available publicly on the Comprehensive R Archive Network.     

Extreme-value analysis of canadian wind speeds

Estimates of extreme winds are essential for engineering design, but in preparing such estimates major statistical issues are encountered. In this case study, the analysts were provided with hourly readings on wind speed, wind direction, and barometric pressure at five Canadian stations for observation periods ranging over several recent decades. Their assignment was to calculate point and interval estimates of 10-, 20-, 50-, and 100-year return values (i.e., upper fractiles) for the wind speeds at these stations.     

Perturbation analysis of Markov modulated fluid models

ABSTRACT

We consider perturbations of positive recurrent Markov modulated fluid models. In addition to the infinitesimal generator of the phases, we also perturb the rate matrix, and analyze the effect of those perturbations on the matrix of first return probabilities to the initial level. Our main contribution is the construction of a substitute for the matrix of first return probabilities, which enables us to analyze the effect of the perturbation under consideration.     

上海股票市场的基于小波去噪的混沌性检验

由于经济混沌需要大样本、低噪声的时间序列,所以文章首先利用小波变换对上证指数日收盘价序列进行去噪处理,然后由去噪后的日收盘价序列计算出日收益率序列,姑且称其为去噪后的日收益率序列,并把它同未经过去噪处理得到的日收益率序列进行比较,发现该方法较好地保留了序列自身固有的特性,只是剔除了由于日常细微波动产生的噪声,为有效地探测我国上海证券市场的混沌性打下了基础。最后分别计算去噪前后收益率的关联维数和Lyapunov指数,发现小波去噪并未改变上海证券市场的混沌性,但是去噪后的市场的复杂度要小于去噪前的市场的复杂度。所以进行混沌性探测的时候必须对数据进行去噪处理。     

Reliability estimation from field return data

In this article statistical inference for the failure time distribution of a product from “field return data”, that records the time between the product being shipped and returned for repair or replacement, is described. The problem that is addressed is that the data are not failure times because they also include the time that it took to ship and install the product and then to return it to the manufacturer for repair or replacement. The inference attempts to infer the distribution of time to failure (that is, from installation to failure) from the data when in addition there are separate data on the times from shipping to installation, and from failure to return. The method is illustrated with data from units installed in a telecommunications network. Our collaborator on writing this paper, Ed Lisay of Alcatel-Lucent, passed away suddenly in October 2008. As a tribute, we can state that Ed had an energetic and vigorous charisma in the application of his skills. He brought a sense of fun to his many interests, such as his achievement of becoming a master electrician. Ed is sadly missed by his family, friends and colleagues.     

A monte carlo of plotting positions

A Monte Carlo study was made of the effects of using simple linear regression, on the appropriate probability paper, to estimate parameters, quantiles and cumulative probability for several distributions. These distributions were the Normal, Weibull (shape parameters 1, 2, and 4) and the Type I largest extreme-value distributions. The specific objective was to observe differences arising from choice of plotting positions. Plotting positions used were i/(n+l), (i?3)/(n+.04), (i?.5)/n, either (i?.375)/(n+.25) or (i?.4)/(n+.2), and either F[E(Y_i)] or F[E(£n Y)]. For each combination of 4 sample sizes (n=10(10)(40)), distribution, and plotting position, regression lines were found for each of N =9999 samples. Each regression line was used to estimate: (1) quantiles of 9 specific probabilities, (2) probabilities of 9 specific quantiles, and (3) return periods corresponding to 9 specific quantiles. Compa[rgrave]ison of the mean, variances, mean square error and medians of these estimates and of the regression coefficients confirm some results of Harter [Commun. Statist. A13(13), 1984] and provide further insight.     

Gaussian copula distributions for mixed data,with application in discrimination

The construction of a joint model for mixed discrete and continuous random variables that accounts for their associations is an important statistical problem in many practical applications. In this paper, we use copulas to construct a class of joint distributions of mixed discrete and continuous random variables. In particular, we employ the Gaussian copula to generate joint distributions for mixed variables. Examples include the robit-normal and probit-normal-exponential distributions, the first for modelling the distribution of mixed binary-continuous data and the second for a mixture of continuous, binary and trichotomous variables. The new class of joint distributions is general enough to include many mixed-data models currently available. We study properties of the distributions and outline likelihood estimation; a small simulation study is used to investigate the finite-sample properties of estimates obtained by full and pairwise likelihood methods. Finally, we present an application to discriminant analysis of multiple correlated binary and continuous data from a study involving advanced breast cancer patients.     

Nonparametric maximum likelihood estimation for artificially truncated absence data

In manpower planning it is cornmoniy tue case tnat employees withuraw from active service for a period of time before returning to take up post at a later date. Such periods of absence are frequently of major concern to employers who are anxious to ensure that employees return as soon as possible. The distribution of duration of such periods of absence are therefore of considerable interest as is the probability that such employees will ever return to active service. In this paper we derive a nonparametric estimator for such a lifetime distribution based on renewal data which are subject to various forms of incompleteness, namely right censoring, left and right truncation, and forward recurrence. Artificial truncation is used to ensure that the data are time homogeneous. A nonparametric maximum likelihood estimator for the lifetime.     

中、印、美股市联动差异性研究

基于扩散视角和跳跃视角探究了中、印、美股市联动行为。基于扩散视角,美国和印度股市与中国股市有明显的单向收益溢出效应,中美之间有明显的波动溢出效应,但是中印之间却不存在这种关系。从非对称影响的结果来看,只存在印度股市和美国股市与中国股市单向的非对称影响。基于跳跃视角,中印、中美股市的平均跳跃幅度和平均方差贡献率,与其跳跃强度相比联动性更高,中印联合跳跃比率相关系数和中美联合跳跃比率相关系数都处于较高水平,同时稳健性检验的结果表明结论整体具有一致性。     

Order statistics from overlapping samples: bivariate densities and regression properties

In this paper, we are interested in the joint distribution of two order statistics from overlapping samples. We give an explicit formula for the distribution of such a pair of random variables under the assumption that the parent distribution is absolutely continuous. We are also interested in the question to what extent conditional expectation of one of such order statistic given another determines the parent distribution. In particular, we provide a new characterization by linearity of regression of an order statistic from the extended sample given the one from the original sample, special case of which solves a problem explicitly stated in the literature. It appears that to describe the correct parent distribution it is convenient to use quantile density functions. In several other cases of regressions of order statistics we provide new results regarding uniqueness of the distribution in the sample.     

资本资产定价模型CAPM在中国资本市场中的实证检验

采用中国上海资本市场交易数据对资本资产定价模型（CAPM）的适用性进行了三个方面的检验：资产的风险和收益之间是否存在线性关系;系统风险是否是资产风险的唯一度量;资产的风险和收益是否正相关。结果发现：2003年8月1日至2006年7月31日期间,上海资本市场股票组合的平均超额收益率与其系统风险之间存在正相关关系,与非系统风险不存在显著的线性关系,基本符合标准形式的CAPM。这与国内许多学者对2001年以前中国资本市场CAPM的实证检验结果不太一致。     

Asymptotic Location of Largest Values of a Stationary Random Field

The study of the relationship between extreme values of dependent random fields and their locations has important practical applications, for instance, when dealing with censored data.

In this article we study the asymptotic behavior of the joint locations of the largest order statistics generated by a stationary random field with extremal index as well as the joint limiting distribution of the location of a high level exceedance nearest of the origin and the location of the maximum.     

Evaluation of missing data mechanisms in two and three dimensional incomplete tables

The analysis of incomplete contingency tables is a practical and an interesting problem. In this paper, we provide characterizations for the various missing mechanisms of a variable in terms of response and non-response odds for two and three dimensional incomplete tables. Log-linear parametrization and some distinctive properties of the missing data models for the above tables are discussed. All possible cases in which data on one, two or all variables may be missing are considered. We study the missingness of each variable in a model, which is more insightful for analyzing cross-classified data than the missingness of the outcome vector. For sensitivity analysis of the incomplete tables, we propose easily verifiable procedures to evaluate the missing at random (MAR), missing completely at random (MCAR) and not missing at random (NMAR) assumptions of the missing data models. These methods depend only on joint and marginal odds computed from fully and partially observed counts in the tables, respectively. Finally, some real-life datasets are analyzed to illustrate our results, which are confirmed based on simulation studies.     

Why does “last week” reporting give higher estimates than “last month”?

Abstract

Experiments in various countries with “last week” and “last month” reference periods for reporting of households’ food consumption have generally found that “week”-based estimates are higher. In India the National Sample Survey (NSS) has consistently found that “week”-based estimates are higher than month-based estimates for a majority of food item groups. But why are week-based estimates higher than month-based estimates? It has long been believed that the reason must be recall lapse, inherent in a long reporting period such as a month. But is household consumption of a habitually consumed item “recalled” in the same way as that of an item of infrequent consumption? And why doesn’t memory lapse cause over-reporting (over-assessment) as often as under-reporting? In this paper, we provide an alternative hypothesis, involving a “quantity floor effect” in reporting behavior, under which “week” may cause over-reporting for many items. We design a test to detect the effect postulated by this hypothesis and carry it out on NSS 68th round HCES data. The test results strongly suggest that our hypothesis provides a better explanation of the difference between week-based and month-based estimates than the recall lapse theory.     

A modified two-stage approach for joint modelling of longitudinal and time-to-event data

Joint models for longitudinal and time-to-event data have been applied in many different fields of statistics and clinical studies. However, the main difficulty these models have to face with is the computational problem. The requirement for numerical integration becomes severe when the dimension of random effects increases. In this paper, a modified two-stage approach has been proposed to estimate the parameters in joint models. In particular, in the first stage, the linear mixed-effects models and best linear unbiased predictorsare applied to estimate parameters in the longitudinal submodel. In the second stage, an approximation of the fully joint log-likelihood is proposed using the estimated the values of these parameters from the longitudinal submodel. Survival parameters are estimated bymaximizing the approximation of the fully joint log-likelihood. Simulation studies show that the approach performs well, especially when the dimension of random effects increases. Finally, we implement this approach on AIDS data.     

Overnight and Daytime Stock-Return Dynamites on the London Stock Exchange: The Impacts of “Big Bang” and the 1987 Stock-Market Crash

We explore the time series properties of stock returns on the London Stock Exchange around the 1986 market restructuring (Big Bang) and the 1987 stock-market crash using a modified generalized autoregressive conditional heteroscedasticity model. Using this general dynamic model, which allows (a) intradaily returns to have different impacts and persistence on stock-return volatility, (b) return effects on volatility to be asymmetric, and (c) intradaily returns to follow conditional distributions with different fourth moments, we uncover important changes in return dynamics and conditional fourth moments following Big Bang and the 1987 crash not reported before.     

Some results on residual life and inactivity time at random time

In this article, we enhance the study of residual life at random time (RLRT) and inactivity time at random time (ITRT). To this aim, first we provide some stochastic orderings results among ITRT in two-sample problems when they failed at two different random times. Then, we develop some sufficient conditions which lead to the stochastic comparisons of RLRT and ITRT based on variance residual life order. The results are expected to be useful in reliability theory, forensic science, queue theory, and actuarial science.     

Joint modelling of longitudinal biomarker and gap time between recurrent events: copula-based dependence

In this paper, we will extend the joint model of longitudinal biomarker and recurrent event via copula function for accounting the dependence between the two processes. The general idea of joining separate processes by allowing model-specific random effect may come from different families distribution. It is a main advantage of the proposed method that a copula construction does not constrain the choice of marginal distributions of random effects. A maximum likelihood estimation with importance sampling technique as a simple and easy understanding method is employed to model inference. To evaluate and verify the validation of the proposed joint model, a bootstrapping method as a model-based resampling is developed. Our proposed joint model is also applied to pemphigus disease data for assessing the effect of biomarker trajectory on risk of recurrence.