Finite mixture of regression models for a stratified sample

Despite the popularity and importance, there is limited work on modelling data which come from complex survey design using finite mixture models. In this work, we explored the use of finite mixture regression models when the samples were drawn using a complex survey design. In particular, we considered modelling data collected based on stratified sampling design. We developed a new design-based inference where we integrated sampling weights in the complete-data log-likelihood function. The expectation–maximisation algorithm was developed accordingly. A simulation study was conducted to compare the new methodology with the usual finite mixture of a regression model. The comparison was done using bias-variance components of mean square error. Additionally, a simulation study was conducted to assess the ability of the Bayesian information criterion to select the optimal number of components under the proposed modelling approach. The methodology was implemented on real data with good results.     

Flexible Tweedie regression models for continuous data

Tweedie regression models (TRMs) provide a flexible family of distributions to deal with non-negative right-skewed data and can handle continuous data with probability mass at zero. Estimation and inference of TRMs based on the maximum likelihood (ML) method are challenged by the presence of an infinity sum in the probability function and non-trivial restrictions on the power parameter space. In this paper, we propose two approaches for fitting TRMs, namely quasi-likelihood (QML) and pseudo-likelihood (PML). We discuss their asymptotic properties and perform simulation studies to compare our methods with the ML method. We show that the QML method provides asymptotically efficient estimation for regression parameters. Simulation studies showed that the QML and PML approaches present estimates, standard errors and coverage rates similar to the ML method. Furthermore, the second-moment assumptions required by the QML and PML methods enable us to extend the TRMs to the class of quasi-TRMs in Wedderburn's style. It allows to eliminate the non-trivial restriction on the power parameter space, and thus provides a flexible regression model to deal with continuous data. We provide an R implementation and illustrate the application of TRMs using three data sets.     

Comments on specification of household expenditure functions and equivalence scales by nonparametric regression

This paper considers Lagrange Multiplier (LM) and Likelihood Ratio (LR) tests for determining the cointegrating rank of a vector autoregressive system. n order to deal with outliers and possible fat-tailedness of the error process, non-Gaussian likelihoods are used to carry out the estimation. The limiting distributions of the tests based on these non-Gaussian pseudo-)likelihoods are derived. These distributions depend on nuisance parameters. An operational procedure is proposed to perform inference. It appears that the tests based on non-Gaussian pseudo-likelihoods are much more powerful than their Gaussian counterparts if the errors are fat-tailed. Moreover, the operational LM-type test has a better overall performance than the LR-type test. Copyright O 1998 by Marcel Dekker, Inc.     

An EM algorithm for the estimation of parameters of bivariate generalized exponential distribution under random left censoring

In this paper, we consider the four-parameter bivariate generalized exponential distribution proposed by Kundu and Gupta [Bivariate generalized exponential distribution, J. Multivariate Anal. 100 (2009), pp. 581–593] and propose an expectation–maximization algorithm to find the maximum-likelihood estimators of the four parameters under random left censoring. A numerical experiment is carried out to discuss the properties of the estimators obtained iteratively.     

Inference on cointegrating ranks using lr and lm tests based on pseudo-likelihoods

This paper considers Lagrange Multiplier (LM) and Likelihood Ratio (LR) tests for determining the cointegrating rank of a vector autoregressive system. n order to deal with outliers and possible fat-tailedness of the error process, non-Gaussian likelihoods are used to carry out the estimation. The limiting distributions of the tests based on these non-Gaussian pseudo-)likelihoods are derived. These distributions depend on nuisance parameters. An operational procedure is proposed to perform inference. It appears that the tests based on non-Gaussian pseudo-likelihoods are much more powerful than their Gaussian counterparts if the errors are fat-tailed. Moreover, the operational LM-type test has a better overall performance than the LR-type test. Copyright O 1998 by Marcel Dekker, Inc.     

Heteroscedastic Regression Models and Applications to Off-line Quality Control

We discuss in the present paper the analysis of heteroscedastic regression models and their applications to off-line quality control problems. It is well known that the method of pseudo-likelihood is usually preferred to full maximum likelihood since estimators of the parameters in the regression function obtained are more robust to misspecification of the variance function. Despite its popularity, however, existing theoretical results are difficult to apply and are of limited use in many applications. Using more recent results in estimating equations, we obtain an efficient algorithm for computing the pseudo-likelihood estimator with desirable convergence properties and also derive simple, explicit and easy to apply asymptotic results. These results are used to look in detail at variance minimization in off-line quality control, yielding techniques of inferences for the optimized design parameter. In application of some existing approaches to off-line quality control, such as the dual response methodology, rigorous statistical inference techniques are scarce and difficult to obtain. An example of off-line quality control is presented to discuss the practical aspects involved in the application of the results obtained and to address issues such as data transformation, model building and the optimization of design parameters. The analysis shows very encouraging results, and is seen to be able to unveil some important information not found in previous analyses.     

Modelling columnarity of pyramidal cells in the human cerebral cortex

For modelling the location of pyramidal cells in the human cerebral cortex, we suggest a hierarchical point process in that exhibits anisotropy in the form of cylinders extending along the z-axis. The model consists first of a generalised shot noise Cox process for the xy-coordinates, providing cylindrical clusters, and next of a Markov random field model for the z-coordinates conditioned on the xy-coordinates, providing either repulsion, aggregation or both within specified areas of interaction. Several cases of these hierarchical point processes are fitted to two pyramidal cell data sets, and of these a final model allowing for both repulsion and attraction between the points seem adequate. We discuss how the final model relates to the so-called minicolumn hypothesis in neuroscience.     

利率期限结构模型非线性建模

应用门限模型对利率期限结构模型中漂移项的非线性性进行建模,提出门限(threshold)CKLS模型。用基于自助法(bootstrap)的广义拟似然比检验方法对门限CKLS模型进行了检验。检验结论表明:门限CKLS模型能较好的刻画利率期限结构模型中漂移项的非线性,在0.1的显著水平下优于CKLS模型。     

Diagnostics for repeated measurements in generalized linear mixed effects models

ABSTRACT

As there is an extensive body of research on diagnostics in regression models, various outlier detection methods have been developed. These methods have been extended to mixed effects models and generalized linear models, but there exist intrinsic drawbacks and limitations. This paper presents two-dimensional plots to identify discordant subjects and observations in generalized linear mixed effects models, displaying discordance in two directions. The sTudentized Residual Sum of Squares is not an extension of any regression tools but a new approach designed to efficiently reflect the characteristics of repeated measures. And this noteworthy clustering of outliers is identified in the plot. Applications to real-life examples are presented to illustrate the favorable/beneficial performance of the new tool.     

Space-time estimation of a particle system model

Let X be a discrete time contact process (CP) on ?², as defined by Durrett and Levin (1994, Stochastic spatial models: a user's guide to ecological applications. Philosophical Transactions of the Royal Society of London Series B, 343, 329–350). We study the estimation of the model based on space-time evolution of X, that is, T + 1 successive observations of X on a finite subset S of sites. We consider the maximum marginal pseudo-likelihood (MPL) estimator and show that, when T→∞, this estimator is consistent and asymptotically normal for a non-vanishing supercritical CP. Numerical studies confirm these theoretical ones.