On the type-B integral equation and the distribution of wilks statist for testing independence of several groups of variables

In this paper, WILKS'type-B integral equation is solved in the general form of a series of beta functions and a series of weighted gamma functions as proposed by WALD and BROOKNER 1941. The coefficients in both representations can be obtained by explicit recurrence relartions, therefore the results solve many distributional problems and have the fewest computational difficulties of any representation that has surfaced to date. The radius of convergence of the second series representation is given, whereas the convergence property of the first series representation is given, whereas the convergence property of the first series representation was studied by WALD and Brookner. The exact null distributions of WILKS' statistic A for testing the independence of several groups of variables and of V = -log A are given. The coefficients in all the series representration can be computed recursilvely and hence can be obtained easily with the help of modern computatinal facilities     

Joint distribution of rises,falls, and number of runs in random sequences

By using the matrix formulation of the two-step approach to the distributions of runs, a recursive relation and an explicit expression are derived for the generating function of the joint distribution of rises and falls for multivariate random sequences in terms of generating functions of individual letters, from which the generating functions of the joint distribution of rises, falls, and number of runs are obtained. An explicit formula for the joint distribution of rises and falls with arbitrary specification is also obtained.     

On the Expressions of Estimability in Testing General Linear Hypotheses

In testing a general linear hypothesis of the form K ′β ? ( W ′) under a general linear model, an equivalent hypothesis involving only estimable parametric functions is provided, and then an explicit test statistic in terms of the model matrices is given. The corresponding results are expanded to the case of a general linear model with a restriction and are illustrated by an example.     

Estimable functions in the nonlinear model

When the method of least squares is used to estimate the parameters in a general model and the generated system of normal equations is linearly dependent, the estimate of the vector of parameters which satisfies the criterion is not unique. However, there exist certain functions of the estimated vector of parameters which are invariant to the least squares solution obtained from the normal equations. We define those invariant functions to be estimable, and present a technique to determine the functions of the parameters which are estimable for the general model. The method results in solving either a linear first order partial differential equation or a system of linear first order partial differential equations corresponding, respectively, to a single or multiple dependency between columns of the Jacobian matrix of the mean of the model. The usual results concerning estimability for linear models are a special case of the general results developed.     

On distributions Of generalized order statistics

In a wide subclass of generalized order statistics, representations of marginal density and distribution functions are developed. The results are applied to obtain several relations, such as recurrence relations, and explicit expressions for the moments of generalized order statistics from Pareto, power function and Weibull distributions Moreover, characterizations of exponential distributions are shown by means of a distributional identity as well as by* an identity of expectations involving a subrange and a corresponding generalized order statistic.     

Estimating the common hazard rate of several exponential distributions

Abstract

In the present communication, we consider the estimation of the common hazard rate of several exponential distributions with unknown and unequal location parameters with a common scale parameter under a general class of bowl-shaped scale invariant loss functions. We have shown that the best affine equivariant estimator (BAEE) is inadmissible by deriving a non smooth improved estimator. Further, we have obtained a smooth estimator which improves upon the BAEE. As an application, we have obtained explicit expressions of improved estimators for special loss functions. Finally, a simulation study is carried out for numerically comparing the risk performance of various estimators.     

Bivariate generalized gamma distributions of Kibble's type

In this paper, a new type of bivariate generalized gamma (BGG) distribution derived from the bivariate gamma distribution of Kibble [Two-variate gamma-type distribution. Sankh?a 1941;5:137–150] by means of a power transformation is presented. The explicit expressions of statistical properties of the BGG distribution are presented. The estimation of marginal and dependence parameters using the method of moments and the method of inference functions for margins are discussed, and their performance through a Monte Carlo simulation study is assessed. Finally, an example is given to illustrate the applicability of the distributions introduced here.     

The star-shaped Λ-coalescent and Fleming–Viot process

The star-shaped Λ-coalescent and corresponding Λ-Fleming–Viot process, where the Λ measure has a single atom at unity, are studied in this article. The transition functions and stationary distribution of the Λ-Fleming–Viot process are derived in a two-type model with mutation. The distribution of the number of non-mutant lines back in time in the star-shaped Λ-coalescent is found. Extensions are made to a model with d types, either with parent-independent mutation or general Markov mutation, and an infinitely-many-types model, when d → ∞. An eigenfunction expansion for the transition functions is found, which has polynomial right eigenfunctions and left eigenfunctions described by hyperfunctions. A further star-shaped model with general frequency-dependent change is considered and the stationary distribution in the Fleming–Viot process derived. This model includes a star-shaped Λ-Fleming–Viot process with mutation and selection. In a general Λ-coalescent explicit formulae for the transition functions and stationary distribution, when there is mutation, are unknown. However, in this article, explicit formulae are derived in the star-shaped coalescent.     

On the Product of Inverse Wishart and Normal Distributions with Applications to Discriminant Analysis and Portfolio Theory

Abstract. In this article we analyse the product of the inverse Wishart matrix and a normal vector. We derive the explicit joint distribution of the components of the product. Furthermore, we suggest several exact tests of general linear hypothesis about the elements of the product. We illustrate the developed techniques on examples from discriminant analysis and from portfolio theory.     

Quotient spaces and statistical models

The purpose of this paper is to draw attention to the widespread occurrence of quotient spaces in statistical work. Quotient spaces are intrinsic to probability distributions, residuals and interaction in linear models, covariance functions and variograms of stochastic processes, etc. The theme is that explicit recognition of the quotient space can offer surprising conceptual simplification. The advantages of working directly with the quotient space are hard to describe in general. As the examples demonstrate, the answer lies partly in directness of approach.     

Robust inference for bivariate point processes

In the analysis of recurrent events where the primary interest lies in studying covariate effects on the expected number of events occurring over a period of time, it is appealing to base models on the cumulative mean function (CMF) of the processes (Lawless & Nadeau 1995). In many chronic diseases, however, more than one type of event is manifested. Here we develop a robust inference procedure for joint regression models for the CMFs arising from a bivariate point process. Consistent parameter estimates with robust variance estimates are obtained via unbiased estimating functions for the CMFs. In most situations, the covariance structure of the bivariate point processes is difficult to specify correctly, but when it is known, an optimal estimating function for the CMFs can be obtained. As a convenient model for more general settings, we suggest the use of the estimating functions arising from bivariate mixed Poisson processes. Simulation studies demonstrate that the estimators based on this working model are practically unbiased with robust variance estimates. Furthermore, hypothesis tests may be based on the generalized Wald or generalized score tests. Data from a trial of patients with bronchial asthma are analyzed to illustrate the estimation and inference procedures.     

Posterior analysis of the compound Rayleigh distribution under balanced loss functions for censored data

This paper develops Bayesian analysis in the context of progressively Type II censored data from the compound Rayleigh distribution. The maximum likelihood and Bayes estimates along with associated posterior risks are derived for reliability performances under balanced loss functions by assuming continuous priors for parameters of the distribution. A practical example is used to illustrate the estimation methods. A simulation study has been carried out to compare the performance of estimates. The study indicates that Bayesian estimation should be preferred over maximum likelihood estimation. In Bayesian estimation, the balance general entropy loss function can be effectively employed for optimal decision-making.     

Multi-step forecasts from threshold ARMA models using asymmetric loss functions

The forecasts generation from nonlinear time series models is investigated under general loss functions. After presenting the main results and some relevant features of these functions, the Linex loss has been used to generate multi-step forecasts from threshold autoregressive moving average models showing their main properties and some results connected to a proper transformation of the forecast errors. A simulation exercise highlights interesting properties of the proposed predictors, both in terms of their bias and their distribution, further clarifying how the Linex predictor can be helpful in empirical applications.     

ESTIMATION FOR THE GENERAL SAMPLE SELECTION MODELS

Consider a general regression model with an arbitrary and unknown link function and a stochastic selection variable that determines whether the outcome variable is observable or missing. The paper proposes U-statistics that are based on kernel functions as estimators for the directions of the parameter vectors in the link function and the selection equation, and shows that these estimators are consistent and asymptotically normal.     

Analysis of Linear Models with One Factor Having Both Fixed and Random Levels

ABSTRACT

A general theory for a case where a factor has both fixed and random effect levels is developed under one-way treatment structure model. Estimation procedures for the fixed effects and variance components are consider for the model. The testing of fixed effects is considered when the variance–covariance matrix is known and unknown. Confidence intervals for estimable functions and prediction intervals for predictable functions are constructed. The computational procedures are illustrated using data from an on-farm trial.     

Cointegrated systems II

The three invited papers in this special issue of Econometric Reviews on "Cointegrated Systems II" complement the previous special issue of the journal. The paper by Eric Zivot and Peter Phillips provides a comprehensive Bayesian analysis of trend determination in economic time series. Two interesting comments on some aspects of current research involving cointegration and the modelling of dynamic economic relationships are provided by Clive Granger and Denzil Fiebig.     

Estimation in Epidemics with Incomplete Observations

The construction of estimating equations by martingale methods is generalized to yield estimators with explicit expressions for the parameters of the birth-and-death and the general epidemic processes when only partial observations are available. (For the birth-and-death process the death process is observed but the number of births is observed only at the end and for the general epidemic process only the removal process is observed.) For large populations, the use of the martingale central limit theorem yields asymptotic confidence regions for the parameters. Explicit expressions are derived for estimators of the variances of the large sample distributions. The range of validity and usefulness of the new estimators is determined by simulation.     

Simple and Explicit Estimating Functions for a Discretely Observed Diffusion Process

We consider a one-dimensional diffusion process X , with ergodic property, with drift b ( x , θ) and diffusion coefficient a ( x , θ) depending on an unknown parameter θ that may be multidimensional. We are interested in the estimation of θ and dispose, for that purpose, of a discretized trajectory, observed at n equidistant times t_i = iΔ , i = 0, ..., n . We study a particular class of estimating functions of the form ∑ f (θ, X _{t i −1}) which, under the assumption that the integral of f with respect to the invariant measure is null, provide us with a consistent and asymptotically normal estimator. We determine the choice of f that yields the estimator with minimum asymptotic variance within the class and indicate how to construct explicit estimating functions based on the generator of the diffusion. Finally the theoretical study is completed with simulations.     

A Note on Continuous Spatial-Temporal Dynamics of Stochastic Processes

We propose a general form to analyze the space-time interdependency of continuous space-time stochastic processes. We present a new space-time approach based on the intensity function of the underlying point process. These formulations can be, to some extent, analytically solved to obtain explicit formulae of interest. We define a general function that controls the space-time interaction and allows for closed forms depending on the particular choice of several mathematical tools playing a role in this interaction function. In particular, we make use of copulas and Laplace transforms to provide interesting examples of the dynamics of the random intensity function and, in turn, of the number of points contained in a given region.     

Robust estimation for functional coefficient regression models with spatial data

A global smoothing procedure is developed using B-spline function approximation for estimating the unknown functions of a functional coefficient regression model with spatial data. A general formulation is used to treat mean regression, median regression, quantile regression and robust mean regression in one setting. The global convergence rates of the estimators of unknown coefficient functions are established. Various applications of the main results, including estimating conditional quantile coefficient functions and robustifying the mean regression coefficient functions are given. Finite sample properties of our procedures are studied through Monte Carlo simulations. A housing data example is used to illustrate the proposed methodology.