Simultaneous equivariant estimation of the parameters of linear models

Consider a family of distributions which is invariant under a group of transformations. In this paper, we define an optimality criterion with respect to an arbitrary convex loss function and we prove a characterization theorem for an equivariant estimator to be optimal. Then we consider a linear model Y=Xβ+ε, in which ε has a multivariate distribution with mean vector zero and has a density belonging to a scale family with scale parameter σ. Also we assume that the underlying family of distributions is invariant with respect to a certain group of transformations. First, we find the class of all equivariant estimators of regression parameters and the powers of σ. By using the characterization theorem we discuss the simultaneous equivariant estimation of the parameters of the linear model.     

Tolerance Limits Based on the Multivariate MaxMin Chart

In a previous paper, we have showed how to obtain sequences of number proved random. With this aim, we used sequences of noises y_n such that the conditional probabilities have Lipschitz coefficients not too large. We transformed them using Fibonacci congruences. Then, we obtained sequences x_n which admit the IID model for correct model. This method consisted to value the work of Marsaglia in order to build his CD-ROM. But we did not use Rap Music (as Marsaglia), but texts files. This method also uses an extractor and at the same time the notion of correct models. In this paper, we apply this method to numbers provided by machines or chips. Unfortunately, it is less sure than they have Lipschtiz coefficient not too large. But we can solve this problem: it suffices to use the Central Limit Theorem. We do it modulo 1. In this case, we use a new limit theorem, the XOR Limit theorem : asymptotic distribution of sum of random vectors modulo 1 are asymptotically independent. Then Lipschtiz coefficient of associated sequences are not too large and we can obtain IID sequences by using Fibonacci congruences.     

On tests detecting difference in means and variances simultaneously under normality

In this paper, we propose several tests for detecting difference in means and variances simultaneously between two populations under normality. First of all, we propose a likelihood ratio test. Then we obtain an expression of the likelihood ratio statistic by a product of two functions of random quantities, which can be used to test the two individual partial hypotheses for differences in means and variances. With those individual partial tests, we propose a union-intersection test. Also we consider two optimal tests by combining the p-values of the two individual partial tests. For obtaining null distributions, we apply the permutation principle with the Monte Carlo approach. Then we compare efficiency among the proposed tests with well-known ones through a simulation study. Finally, we discuss some interesting features related to the simultaneous tests and resampling methods as concluding remarks.     

Non-homogeneous Pólya-Aeppli process

Abstract

In this paper we suppose that the intensity parameter of the Pólya-Aeppli process is a function of time t and call the resulting process a non-homogeneous Pólya-Aeppli process (NHPAP). The NHPAP can be represented as a compound non-homogeneous Poisson process with geometric compounding distribution as well as a pure birth process. For this process we give two definitions and show their equivalence. Also, we derive some interesting properties of NHPAP and use simulation the illustrate the process for particular intensity functions. In addition, we introduce the standard risk model based on NHPAP, analyze the ruin probability for this model and include an example of the process under exponentially distributed claims.     

Asymptotic distribution of least square estimators for linear models with dependent errors

In this paper, we consider the usual linear regression model in the case where the error process is assumed strictly stationary. We use a result from Hannan (Central limit theorems for time series regression. Probab Theory Relat Fields. 1973;26(2):157–170), who proved a Central Limit Theorem for the usual least squares estimator under general conditions on the design and on the error process. Whatever the design satisfying Hannan's conditions, we define an estimator of the covariance matrix and we prove its consistency under very mild conditions. As an application, we show how to modify the usual tests on the linear model in this dependent context, in such a way that the type-I error rate remains asymptotically correct, and we illustrate the performance of this procedure through different sets of simulations.     

Using p values to design statistical process control charts

Conventional Phase II statistical process control (SPC) charts are designed using control limits; a chart gives a signal of process distributional shift when its charting statistic exceeds a properly chosen control limit. To do so, we only know whether a chart is out-of-control at a given time. It is therefore not informative enough about the likelihood of a potential distributional shift. In this paper, we suggest designing the SPC charts using p values. By this approach, at each time point of Phase II process monitoring, the p value of the observed charting statistic is computed, under the assumption that the process is in-control. If the p value is less than a pre-specified significance level, then a signal of distributional shift is delivered. This p value approach has several benefits, compared to the conventional design using control limits. First, after a signal of distributional shift is delivered, we could know how strong the signal is. Second, even when the p value at a given time point is larger than the significance level, it still provides us useful information about how stable the process performs at that time point. The second benefit is especially useful when we adopt a variable sampling scheme, by which the sampling time can be longer when we have more evidence that the process runs stably, supported by a larger p value. To demonstrate the p value approach, we consider univariate process monitoring by cumulative sum control charts in various cases.     

Positive-rule stein-type almost unbiased ridge estimator in linear regression model

Abstract

In this article, when it is suspected that regression coefficients may be restricted to a subspace, we discuss the parameter estimation of regression coefficients in a multiple regression model. Then, in order to improve the preliminary test almost ridge estimator, we study the positive-rule Stein-type almost unbiased ridge estimator based on the positive-rule stein-type shrinkage estimator and almost unbiased ridge estimator. After that, quadratic bias and quadratic risk values of the new estimator are derived and compared with some relative estimators. And we also discuss the option of parameter k. Finally, we perform a real data example and a Monte Carlo study to illustrate theoretical results.     

Copula-Based Models for the Power of Independence Tests

ABSTRACT

We consider independence tests and the methods to evaluate their efficiency. First, we observe that many of the most used independence tests are functions of the empirical copula, which is a sufficient statistic. Hence, the power of these tests, such as the tests based on Spearman's ρ, on Kendall's τ, and on Gini's γ, depend solely on the theoretical copula, and not on the marginal distributions. Then, we consider monotone dependence tests and we propose a parametric model to define the power function. Such a model is based on a path of copulas, from the copula of discordance to the copula of concordance, and can be characterized by the copula of the underlying joint distribution. Moreover, we introduce a consistent estimator of the path of copulas. Finally, we provide some examples of applications, and in particular, a bootstrap-plug-in estimator of the power curve, all useful for power comparison.     

An extension of the Anderson–Darling k-sample test to arbitrary sample space partition sizes

In this paper we first show that the k-sample Anderson–Darling test is basically an average of Pearson statistics in 2?×?k contingency tables that are induced by observation-based partitions of the sample space. As an extension, we construct a family of rank test statistics, indexed by c?∈??, which is based on similarly constructed c?×?k partitions. An extensive simulation study, in which we compare the new test with others, suggests that generally very high powers are obtained with the new tests. Finally we propose a decomposition of the test statistic in interpretable components.     

On the property of multivariate generalized hyperbolic distribution and the Stein-type inequality

This note consists of two parts . In the first part, we provide a pedagogic review on the multivariate generalized hyperbolic (MGH) distribution. We show that this probability family is close under margining, conditioning, and linear transforms; however, such property does not hold for its subclasses. In the second part, we obtain the Stein-type inequality in the context of MGH distribution. Moreover, we apply the Stein-type inequality to prove a lower bound for Var[h(X)]. Particularly, we present examples when X belongs to some well-known subclasses in MGH family.     

Central limit theorems of range-based estimators for diffusion models

Abstract

In this article, we consider non parametric range-based estimation procedure for diffusion processes and propose a instantaneous volatility estimator. Under some weak conditions, we certify that the proposed estimator has convergence in probability. Adding some necessary conditions, we prove a central limit theorem. By inference, we reach a conclusion that, with high frequency data in hand, the proposed estimator is more precise than those pure realized instantaneous volatility ones. Numerical simulation illustrates the finite sample properties of the proposed estimator.     

Group sequential x2 statistic for testing the difference of two mean vectors in clinical trials

In this study we discuss the group sequential procedures for comparing two treatments based on multivariate observations in clinical trials. Also we suppose that a response vector on each of two treatments has a multivariate normal distribution with unknown covariance matrix. Then we propose a group sequential x² statistic in order to carry out repeated significance test for hypothesis of no difference between two population mean vectors. In order to realize the group sequential test where average sample number is reduced, we propose another modified group sequential x² statistic by extension of Jennison and Turnbull ( 1991 ). After construction of repeated confidence boundaries for making the repeated significance test, we compare two group sequential procedures based on two statistics regarding the average sample number and the power of the test in the simulations.     

Exchange option pricing in jump-diffusion models based on esscher transform

In the real world, we introduce a dynamic model about the risky asset which is governed by Brownian motion, stationary compound Poisson process and its compensation process. By choosing Esscher transform parameters, we obtain a risk-neural measure Q under which the discounted value of the risky underlying asset is a martingale. Then, we give the pricing formulas of Exchange option by change of numeraire. At last, we analyze the option pricing formula and provide numerical illustrations by introducing BBY stock and SBUX stock.     

Regular automorphism groups on partial geometries

Recently, Ghinelli (Geom Dedicata (1992) 165–174) had studied the generalized quadrangles which admits automorphism groups acting regularly on the points. In this paper, we generalize her idea to partial geometries, pg(s + 1, t + 1, α). Some examples and basic properties are given. In particular, we prove that under certain conditions on the automorphism group and the lines, such a geometry is a translation net. Applying the results to the case when s = t and the automorphism group G is abelian, we find that either the geometry is a translation net or all the lines of the geometry are generated by a subset of G. Also, for this case, we conjecture that the parameter α is either s or s + 1, except (s, α) = (5, 2), and we have checked that it is true for s ⩽ 500.     

Weak convergence for the conditional distribution function in a Koziol–Green model under dependent censoring

In this paper we consider the conditional Koziol–Green model of Braekers and Veraverbeke [2008. A conditional Koziol–Green model under dependent censoring. Statist. Probab. Lett., accepted for publication] in which they generalized the Koziol–Green model of Veraverbeke and Cadarso Suárez [2000. Estimation of the conditional distribution in a conditional Koziol–Green model. Test 9, 97–122] by assuming that the association between a censoring time and a time until an event is described by a known Archimedean copula function. They got in this way, an informative censoring model with two different types of informative censoring. Braekers and Veraverbeke [2008. A conditional Koziol–Green model under dependent censoring. Statist. Probab. Lett., accepted for publication] derived in this model a non-parametric Koziol–Green estimator for the conditional distribution function of the time until an event, for which they showed the uniform consistency and the asymptotic normality. In this paper we extend their results and prove the weak convergence of the process associated to this estimator. Furthermore we show that the conditional Koziol–Green estimator is asymptotically more efficient in this model than the general copula-graphic estimator of Braekers and Veraverbeke [2005. A copula-graphic estimator for the conditional survival function under dependent censoring. Canad. J. Statist. 33, 429–447]. As last result, we construct an asymptotic confidence band for the conditional Koziol–Green estimator. Through a simulation study, we investigate the small sample properties of this asymptotic confidence band. Afterwards we apply this estimator and its confidence band on a practical data set.     

Semiparametric Regression Analysis of Longitudinal Skewed Data

In this paper, we develop a semiparametric regression model for longitudinal skewed data. In the new model, we allow the transformation function and the baseline function to be unknown. The proposed model can provide a much broader class of models than the existing additive and multiplicative models. Our estimators for regression parameters, transformation function and baseline function are asymptotically normal. Particularly, the estimator for the transformation function converges to its true value at the rate n ^{? 1 ∕ 2}, the convergence rate that one could expect for a parametric model. In simulation studies, we demonstrate that the proposed semiparametric method is robust with little loss of efficiency. Finally, we apply the new method to a study on longitudinal health care costs.     

Non‐Gaussian geostatistical modeling using (skew) t processes

We propose a new model for regression and dependence analysis when addressing spatial data with possibly heavy tails and an asymmetric marginal distribution. We first propose a stationary process with t marginals obtained through scale mixing of a Gaussian process with an inverse square root process with Gamma marginals. We then generalize this construction by considering a skew‐Gaussian process, thus obtaining a process with skew‐t marginal distributions. For the proposed (skew) t process, we study the second‐order and geometrical properties and in the t case, we provide analytic expressions for the bivariate distribution. In an extensive simulation study, we investigate the use of the weighted pairwise likelihood as a method of estimation for the t process. Moreover we compare the performance of the optimal linear predictor of the t process versus the optimal Gaussian predictor. Finally, the effectiveness of our methodology is illustrated by analyzing a georeferenced dataset on maximum temperatures in Australia.     

Une application des méthodes séquentielles aux tests de permutations

To carry out a permutation test we have to examine the n! permutations of the observations. In order to make the permutation test feasible, Dwass (1957) proposed to examine only a sample of these permutations. With the help of sequential methods, we obtain a test which is never less efficient than that proposed by Dwass or the permutation test itself, in the sense that it is as powerful and never requires more permutations to make a decision. In practice, we can expect to gain much efficiency.     

A restricted Liu estimator for binary regression models and its application to an applied demand system

In this article, we propose a restricted Liu regression estimator (RLRE) for estimating the parameter vector, β, in the presence of multicollinearity, when the dependent variable is binary and it is suspected that β may belong to a linear subspace defined by Rβ?=?r. First, we investigate the mean squared error (MSE) properties of the new estimator and compare them with those of the restricted maximum likelihood estimator (RMLE). Then we suggest some estimators of the shrinkage parameter, and a simulation study is conducted to compare the performance of the different estimators. Finally, we show the benefit of using RLRE instead of RMLE when estimating how changes in price affect consumer demand for a specific product.     

Kullback–Leibler information of a censored variable and its applications

In this paper, we study the Kullback–Leibler (KL) information of a censored variable, which we will simply call it censored KL information. The censored KL information is shown to have the necessary monotonicity property in addition to inherent properties of nonnegativity and characterization. We also present a representation of the censored KL information in terms of the relative risk and study its relation with the Fisher information in censored data. Finally, we evaluate the estimated censored KL information as a goodness-of-fit test statistic.