Analysis of a continuous-time proportional hazard model using discrete duration data

Many economic duration variables are often available only up to intervals, and not up to exact points. However, continuous time duration models are conceptually superior to discrete ones. Hence, in duration analyses, one faces a situation with discrete data and a continuous model. This paper discusses (i) the asymptotic bias of a conventional approximation procedure in which a discrete duration is treated as an exact observation; and (ii) the efficiency of a correct maximum likelihood estimator which appropriately accounts for the discrete nature of the data.     

A comparison of two approaches for power and sample size calculations in logistic regression models

Whittemore (1981) proposed an approach for calculating the sample size needed to test hypotheses with specified significance and power against a given alternative for logistic regression with small response probability. Based on the distribution of covariate, which could be either discrete or continuous, this approach first provides a simple closed-form approximation to the asymptotic covariance matrix of the maximum likelihood estimates, and then uses it to calculate the sample size needed to test a hypothesis about the parameter. Self et al. (1992) described a general approach for power and sample size calculations within the framework of generalized linear models, which include logistic regression as a special case. Their approach is based on an approximation to the distribution of the likelihood ratio statistic. Unlike the Whittemore approach, their approach is not limited to situations of small response probability. However, it is restricted to models with a finite number of covariate configurations. This study compares these two approaches to see how accurate they would be for the calculations of power and sample size in logistic regression models with various response probabilities and covariate distributions. The results indicate that the Whittemore approach has a slight advantage in achieving the nominal power only for one case with small response probability. It is outperformed for all other cases with larger response probabilities. In general, the approach proposed in Self et al. (1992) is recommended for all values of the response probability. However, its extension for logistic regression models with an infinite number of covariate configurations involves an arbitrary decision for categorization and leads to a discrete approximation. As shown in this paper, the examined discrete approximations appear to be sufficiently accurate for practical purpose.     

N-stage splitting for maximum penalized likelihood estimation

An n-stage splitting algorithm for the solution of maximum penalized likelihood estimation (MPLE) problems is compared to the one-step-late (OSL) algorithm. General conditions under which the asymptotic rate of convergence of this splitting algorithm. exceeds that of the OSL algorithm are given. A one-dimensional positive data example, illustrates the comparison of the rates of convergence of these two algorithms.     

A pattern-mixture odds ratio model for incomplete categorical data

Most models for incomplete data are formulated within the selection model framework. Pattern-mixture models are increasingly seen as a viable alternative, both from an interpretational as well as from a computational point of view (Little 1993, Hogan and Laird 1997, Ekholm and Skinner 1998). Whereas most applications are either for continuous normally distributed data or for simplified categorical settings such as contingency tables, we show how a multivariate odds ratio model (Molenberghs and Lesaffre 1994, 1998) can be used to fit pattern-mixture models to repeated binary outcomes with continuous covariates. Apart from point estimation, useful methods for interval estimation are presented and data from a clinical study are analyzed to illustrate the methods.     

A New Discrete Distribution Related to Generalized Gamma Distribution and Its Properties

In this article, a new discrete distribution related to the generalized gamma distribution (Stacy, 1962) is derived from a statistical mechanical setup. This new distribution can be seen as generalization of two-parameter discrete gamma distribution (Chakraborty and Chakravarty, 2012) and encompasses discrete version of many important continuous distributions. Some basic distributional and reliability properties, parameter estimation by different methods, and their comparative performances using simulation are investigated. Two-real life data sets are considered for data modeling and likelihood ratio test for illustrating the advantages of the proposed distribution over two-parameter discrete gamma distribution.     

A note on properties of spatial yule-walker estimators

For two-dimensional spatial autoregressive (AR) models, asymptotic properties of the spatial Yule-Walker (YW) estimators (Tjøstheim, 1978) are studied. These estimators although consistent, are shown to be asymptotically biased. Estimators from the first-order spatial bilateral AR model are looked at in more detail and the spatial YW estimators for this model are compared with the exact maximum likelihood estimators. Small sample properties of both estimators are also discussed briefly and some simulation results are presented.     

Beta Regression for Modelling Rates and Proportions

This paper proposes a regression model where the response is beta distributed using a parameterization of the beta law that is indexed by mean and dispersion parameters. The proposed model is useful for situations where the variable of interest is continuous and restricted to the interval (0, 1) and is related to other variables through a regression structure. The regression parameters of the beta regression model are interpretable in terms of the mean of the response and, when the logit link is used, of an odds ratio, unlike the parameters of a linear regression that employs a transformed response. Estimation is performed by maximum likelihood. We provide closed-form expressions for the score function, for Fisher's information matrix and its inverse. Hypothesis testing is performed using approximations obtained from the asymptotic normality of the maximum likelihood estimator. Some diagnostic measures are introduced. Finally, practical applications that employ real data are presented and discussed.     

The beta distribution as a latent response model for ordinal data (I): Estimation of location and dispersion parameters

Ordinal data are often modeled using a continuous latent response distribution, which is partially observed through windows of adjacent intervals defined by cutpoints. In this paper we propose the beta distribution as a model for the latent response. The beta distribution has several advantages over the other common distributions used, e.g. , normal and logistic. In particular, it enables separate modeling of location and dispersion effects which is essential in the Taguchi method of robust design. First, we study the problem of estimating the location and dispersion parameters of a single beta distribution (representing a single treatment) from ordinal data assuming known equispaced cutpoints. Two methods of estimation are compared: the maximum likelihood method and the method of moments. Two methods of treating the data are considered: in raw discrete form and in smoothed continuousized form. A large scale simulation study is carried out to compare the different methods. The mean square errors of the estimates are obtained under a variety of parameter configurations. Comparisons are made based on the ratios of the mean square errors (called the relative efficiencies). No method is universally the best, but the maximum likelihood method using continuousized data is found to perform generally well, especially for estimating the dispersion parameter. This method is also computationally much faster than the other methods and does not experience convergence difficulties in case of sparse or empty cells. Next, the problem of estimating unknown cutpoints is addressed. Here the multiple treatments setup is considered since in an actual application, cutpoints are common to all treatments, and must be estimated from all the data. A two-step iterative algorithm is proposed for estimating the location and dispersion parameters of the treatments, and the cutpoints. The proposed beta model and McCullagh's (1980) proportional odds model are compared by fitting them to two real data sets.     

On Discrete Gamma Distribution

A new discrete counterpart of gamma distribution for modelling discrete life data is defined based on similar mathematical form and properties of the continuous version. The main statistical and reliability properties of this distribution are derived and it is shown that this model can deal with both over and under-dispersed data. Geometric variables and finite sum of geometric variables, i.e., negative binomial are shown to be special cases of the proposed discrete gamma. Also, the size-biased discrete gamma distribution is derived and discussed. Moreover, different estimation methods of the underlying parameters of this distribution are utilized and comparisons of their performance have been made. Finally, an application in real-life data is used to elucidate the earlier results of this article.     

时变系数空间自回归面板数据模型的极大似然估计

本文对扰动项存在跨时期的异方差、但不存在序列相关的时变系数空间自回归模型提出了极大似然的估计方法,并证明了该估计量的一致性,同时,证明了该估计量渐进服从正态分布,由此说明该估计量具有优良的大样本性质。同时,我们还对本文所提出估计量的小样本性质进行了数值模拟。本文研究表明,估计量虽然在N较小时偏差较大,但是随着N的不断增加,估计量偏差减小,体现了比较优良的渐进性质。同时,估计量的偏差会随着时期数的增加而变大,这说明本文所提出的估计方法适用于个体数较多、时期数较少的短面板数据。     

A Likelihood Ratio Test for Idiosyncratic Unit Roots in the Exact Factor Model with Integrated Factors

We consider an exact factor model with unobservable common stochastic trends imposed by non-stationary factors, and study, by simulation, the power of the likelihood ratio test for unit roots in the idiosyncratic components. The power of the test is compared with the analogous Lagrange multiplier test and the Fisher-type test proposed by Bai and Ng. The results suggest that the benefit of the likelihood ratio test is in panels with a small cross-section.     

Maximum likelihood estimation in the proportional hazards model of random censorship

The maximum likelihood estimator (MLE) for the survival function S_Tunder the proportional hazards model of censorship is derived and shown to differ from the Abdushukurov-Cheng-Lin estimator when the class of allowable distributions includes all continuous and discrete distributions. The estimators are compared via an example. The MLE is calculated using a Newton-Raphson iterative procedure and implemented via a FORTRAN algorithm.     

Latent variable models with mixed continuous and polytomous data

Owing to the nature of the problems and the design of questionnaires, discrete polytomous data are very common in behavioural, medical and social research. Analysing the relationships between the manifest and the latent variables based on mixed polytomous and continuous data has proven to be difficult. A general structural equation model is investigated for these mixed outcomes. Maximum likelihood (ML) estimates of the unknown thresholds and the structural parameters in the covariance structure are obtained. A Monte Carlo–EM algorithm is implemented to produce the ML estimates. It is shown that closed form solutions can be obtained for the M-step, and estimates of the latent variables are produced as a by-product of the analysis. The method is illustrated with a real example.     

Finite sample properties of maximum likelihood and quasi-maximum likelihood estimators of egarch models

In this paper I examine finite sample properties of the maximum likelihood and quasi-maximum likelihood estimators of EGARCH(1,1) processes using Monte Carlo methods. I use response surface methodology to summarize the results of a wide array of experiments which suggest that the maximum likelihood estimator has reasonable finite sample properties. The Gaussian quasi-maximum likelihood estimator has poor finite sample properties when the data generating process has conditional excess kurtosis. Some of these poor properties appear to be asymptotic in nature.     

Finite sample properties of maximum likelihood and quasi-maximum likelihood estimators of egarch models

In this paper I examine finite sample properties of the maximum likelihood and quasi-maximum likelihood estimators of EGARCH(1,1) processes using Monte Carlo methods. I use response surface methodology to summarize the results of a wide array of experiments which suggest that the maximum likelihood estimator has reasonable finite sample properties. The Gaussian quasi-maximum likelihood estimator has poor finite sample properties when the data generating process has conditional excess kurtosis. Some of these poor properties appear to be asymptotic in nature.     

An exact density-based empirical likelihood ratio test for paired data

The Wilcoxon rank-sum test and its variants are historically well-known to be very powerful nonparametric decision rules for testing no location difference between two groups given paired data versus a shift alternative. In this title, we propose a new alternative empirical likelihood (EL) ratio approach for testing the equality of marginal distributions given that sampling is from a continuous bivariate population. We show that in various shift alternative scenarios the proposed exact test is superior to the classic nonparametric procedures, which may break down completely or are frequently inferior to the density-based EL ratio test. This is particularly true in the cases where there is a nonconstant shift under the alternative or the data distributions are skewed. An extensive Monte Carlo study shows that the proposed test has excellent operating characteristics. We apply the density-based EL ratio test to analyze real data from two medical studies.     

Estimation of the mean of the exponential distribution using moving extremes ranked set sampling

Moving Extremes Ranked Set Sampling (MERSS) is a useful modification of Ranked Set Sampling (RSS). Unlike RSS, MERSS allows for an increase of set size without introducing too much ranking error. The method is considered parametrically under exponential distribution. Maximum likelihood estimator (MLE), and a modified MLE are considered and their properties are studied. The method is studied under both perfect and imperfect ranking (with error in ranking). It appears that these estimators can be real competitors to the MLE using the usual simple random sampling (SRS).     

Testing the Existence of Change-Point in NHPP Software Reliability Models

In the software testing process, the nature of the failure data is affected by many factors, such as the testing environment, testing strategy, and resource allocation. These factors are unlikely to all be kept stable during the entire process of software testing. As a result, the statistical structure of the failure data is likely to experience major changes. Recently, some useful non homogeneous Poisson process (NHPP) models with change-point are proposed. However, in many realistic situations, whether a change-point exists is unknown. Furthermore, some real data seem to have two or more change-points. In this article we propose test statistics to test the existence of change-point(s). The experimental results of real data show that our tests perform well.     

Saddlepoint Approximations for Stopped-Sum Distributions

One of the common used classes of distributions is the stopped-sum class. This class includes Hermite distribution, Polya–Aeppli distribution, Poisson-Gamma distribution, and Neyman type A. This article introduces the saddlepoint approximations to the stopped-sum class in continuous and discrete settings. We discuss approximations for mass/density and cumulative distribution functions of stopped-sum distributions. Examples of continuous and discrete distributions from the Poisson stopped-sum class are presented. Comparisons between saddlepoint approximations and the exact calculations show the great accuracy of the saddlepoint methods.     

Negative moments of a modified power series distribution and bias of the maximum likelihood estimator

The class of Modified Power Series distributions (MPSD) containing Lagrangian Poisson (LPD) (Consul and Jain, 1973) and Lagrangian binomial distributions (LBD) (Jain and Consul, 1971) was studied by Gupta (1974). We investigate the problem of finding the negative momentsE[X^-r ], of displaced and decapitated Modified Power Series Distributions. We derive the relationship between rand (r-1) negative moments. The negative moments of the decapitated and displaced LPD are obtained. These results are, then, used to find the exact amount of bias in the ML estimators of the parameters in the LPD and the LBD. We have also given the variances of the ML estimator and the minimum variance unbiased estimator of the parameter in the LPD.