Large sample bias of maximum likeihood estimates from a series of sparse ixj tables

For a large series of IxJ tables, each containing two observations, the bias of the maximum likelihood estimates of log linear partial association parameters is shown to be equal to the parameters, regardless of the size of I and J. The partial association considered is that between row and column variables; the three way interactions are assumed to be O. This is a generalization of Andersen's results (1973a, 1973b) for a series of 2x2 tables.     

A multivariate long-memory model with structural breaks

This paper introduces a multivariate long-memory model with structural breaks. In the proposed framework, the time series exhibits possibly fractional orders of integration which are allowed to be different in each subsample. The break date is endogenously determined using a procedure that minimizes the residual sum of squares (RSS). Monte Carlo experiments show that this method for detecting breaks performs well in large samples. As an illustration, we estimate a trivariate VAR including prices, employment and GDP in both the US and Mexico. For the subsample preceding the break, our findings are similar to those of earlier studies based on a standard VAR approach in both the countries, such that the variables exhibit integer degrees of integration. On the contrary, the series is found to be fractionally integrated after the break, with the fractional differencing parameters being higher than one in the case of Mexico.     

Variable selection for multivariate generalized linear models

Generalized linear models (GLMs) are widely studied to deal with complex response variables. For the analysis of categorical dependent variables with more than two response categories, multivariate GLMs are presented to build the relationship between this polytomous response and a set of regressors. Traditional variable selection approaches have been proposed for the multivariate GLM with a canonical link function when the number of parameters is fixed in the literature. However, in many model selection problems, the number of parameters may be large and grow with the sample size. In this paper, we present a new selection criterion to the model with a diverging number of parameters. Under suitable conditions, the criterion is shown to be model selection consistent. A simulation study and a real data analysis are conducted to support theoretical findings.     

A control chart based on sample ranges for monitoring the covariance matrix of the multivariate processes

For the univariate case, the R chart and the S ² chart are the most common charts used for monitoring the process dispersion. With the usual sample size of 4 and 5, the R chart is slightly inferior to the S ² chart in terms of efficiency in detecting process shifts. In this article, we show that for the multivariate case, the chart based on the standardized sample ranges, we call the RMAX chart, is substantially inferior in terms of efficiency in detecting shifts in the covariance matrix than the VMAX chart, which is based on the standardized sample variances. The user's familiarity with sample ranges is a point in favor of the RMAX chart. An example is presented to illustrate the application of the proposed chart.     

Large sample interval estimation for the logarithmic series distribution

It is assumed that a large random sample of fixed size n is drawn from a logarithmic series distribution with parameter $tH and that it is desired to estimate e by means of two-sided confidence interval. Using asymptotic results, charts of confidence intervals are prepared for n = 50, 100, 200, 500 and 1,000, and confidence coefficients 0.90 and 0.95.     

On some tests for a multivariate linear functional relationship with replication

For testing problems of the coefficient vector and the interception of multivariate linear functional relationship with replicated observations, the likelihood ratio test statistics are considered. Their asymptotic distributions are obtained under each null hypothesis respectively.     

On the distribution of some test statistics connected with the multivariate linear functional relationship model

An alternate representation of the densities of some test statistics for the structural coefficients of the multivariate linear functional relationship model is proposed in this article. These statistics are distributed as the ratio of a linear combination of chi-square variÂtes over the root of a product of chi-square variÂtes. A computable representation of their densities has already been derived by Provost (1984) with the help of the technique of the inverse Mellin transform. The connection of the alternate representation to the densities of products of independent beta type-2 and of independent F-random variables is also discussed.     

Least squares theory for possibly singular models

In a recent paper, Scobey (1975) observed that the usual least squares theory can be applied even when the covariance matrix σ²V of Y in the linear model Y = Xβ + e is singular by choosing the Moore-Penrose inverse (V+XX′)⁺ instead of V^-1 when V is nonsingular. This result appears to be wrong. The appropriate treatment of the problem in the singular case is described.     

Testing equality of two intercepts for the parallel regression model with non sample prior information

This paper proposes tests for equality of intercepts of two simple regression models when non sample prior information is available on the equality of two slopes. For three different scenarios on the values of the slope, namely (i) unknown (unspecified), (ii) known (specified), and (iii) suspected, we derive the unrestricted test (UT), restricted test (RT), and pre-test test (PTT) for testing the equality of intercepts. The test statistics, their sampling distributions, and power functions of the tests are obtained. Comparison of power function and size of the tests reveals that the PTT has a reasonable dominance over the UT and RT.     

Estimating parameters in a simple errors-in-variables model: a new approach base on finite sample distribution theory

This article presents a first direct application of finite sample distribution theory. The relevance of analytical finite sample research is exemplified in the framework of a simple linear errors-in-variables model (EV Model) with known or approximately known measurement error variance. Analytical results derived byRichardson/Wu (1970) are applied for constructing new approximately unbiased estimators for the slope coefficient in the EV model. The new estimators are compared with the biased least squares estimator and with asymptotic theory based corrected least squares estimators. Retaining responsibility for remaining errors the author is indebted to Prof. H. Schneewei\ and Prof. J. Gruber for helpful comments and discussions. Mrs. A. Brandtstater deserves special mention and thanks for performing the computations reported in section 4.     

The large sample coverage probability of confidence intervals in general regression models after a preliminary hypothesis test

We derive a computationally convenient formula for the large sample coverage probability of a confidence interval for a scalar parameter of interest following a preliminary hypothesis test that a specified vector parameter takes a given value in a general regression model. Previously, this large sample coverage probability could only be estimated by simulation. Our formula only requires the evaluation, by numerical integration, of either a double or a triple integral, irrespective of the dimension of this specified vector parameter. We illustrate the application of this formula to a confidence interval for the odds ratio of myocardial infarction when the exposure is recent oral contraceptive use, following a preliminary test where two specified interactions in a logistic regression model are zero. For this real‐life data, we compare this large sample coverage probability with the actual coverage probability of this confidence interval, obtained by simulation.     

Testing for a change lit the regression matrix of a multivariate linear model

A Bayesian test procedure Is developed to test; the null hypothesis of no change In the regression matrix of a multivariate lin¬ear model against the alternative hypothesis of exactly one change The resulting test is based on the marginal posterior distribution of the change point; To illustrate the test procedure a numerical example using a bivariate regression model is considered.     

Parameter estimation for multivariate diffusion processes with the time inhomogeneously positive semidefinite diffusion matrix

Statistical inference for the diffusion coefficients of multivariate diffusion processes has been well established in recent years; however, it is not the case for the drift coefficients. Furthermore, most existing estimation methods for the drift coefficients are proposed under the assumption that the diffusion matrix is positive definite and time homogeneous. In this article, we put forward two estimation approaches for estimating the drift coefficients of the multivariate diffusion models with the time inhomogeneously positive semidefinite diffusion matrix. They are maximum likelihood estimation methods based on both the martingale representation theorem and conditional characteristic functions and the generalized method of moments based on conditional characteristic functions, respectively. Consistency and asymptotic normality of the generalized method of moments estimation are also proved in this article. Simulation results demonstrate that these methods work well.     

On the bias and mean square error of the least square estimator in a regression model with two inequality constraints and multivariate t error terms

We derive and numerically evaluate the bias and mean square error of the inequality constrained least squares estimator in a model with two inequality constraints and multivariate terror terms. Our results suggest that qualitatively, the estimator properties found for models with normal errors carry over to the case of multivariate terrors.     

A note on nonlinear regression for the autoregressive moving average with non-hd errors

Estimation by nonlinear regression of the parameters for the stationary and invertible autoregressive moving average (ARMA) model with mixing or martingale difference errors is considered. Simple proofs of consistency and asymptotic normality for the nonlinear least squares estimator are given by exploiting results from nonlinear estimation theory and mixing and mixingale theory.     

Skew-mixed effects model for multivariate longitudinal data with categorical outcomes and missingness

A longitudinal study commonly follows a set of variables, measured for each individual repeatedly over time, and usually suffers from incomplete data problem. A common approach for dealing with longitudinal categorical responses is to use the Generalized Linear Mixed Model (GLMM). This model induces the potential relation between response variables over time via a vector of random effects, assumed to be shared parameters in the non-ignorable missing mechanism. Most GLMMs assume that the random-effects parameters follow a normal or symmetric distribution and this leads to serious problems in real applications. In this paper, we propose GLMMs for the analysis of incomplete multivariate longitudinal categorical responses with a non-ignorable missing mechanism based on a shared parameter framework with the less restrictive assumption of skew-normality for the random effects. These models may contain incomplete data with monotone and non-monotone missing patterns. The performance of the model is evaluated using simulation studies and a well-known longitudinal data set extracted from a fluvoxamine trial is analyzed to determine the profile of fluvoxamine in ambulatory clinical psychiatric practice.     

Joint Central Limit Theorem for Eigenvalue Statistics from Several Dependent Large Dimensional Sample Covariance Matrices with Application

Let X _n = (x _{i j} ) be a k ×n data matrix with complex‐valued, independent and standardized entries satisfying a Lindeberg‐type moment condition. We consider simultaneously R sample covariance matrices , where the Q _r 's are non‐random real matrices with common dimensions p ×k (k ≥p ). Assuming that both the dimension p and the sample size n grow to infinity, the limiting distributions of the eigenvalues of the matrices { B _{n r} } are identified, and as the main result of the paper, we establish a joint central limit theorem (CLT) for linear spectral statistics of the R matrices { B _{n r} }. Next, this new CLT is applied to the problem of testing a high‐dimensional white noise in time series modelling. In experiments, the derived test has a controlled size and is significantly faster than the classical permutation test, although it does have lower power. This application highlights the necessity of such joint CLT in the presence of several dependent sample covariance matrices. In contrast, all the existing works on CLT for linear spectral statistics of large sample covariance matrices deal with a single sample covariance matrix (R = 1).