Some tests for the allometric extension model in regression

The allometric extension model is a multivariate regression model recently proposed by Tarpey and Ivey (2006 Tarpey, T., Ivey, C.T. (2006). Allometric extension for multivariate regression. J. Data Sci. 4:479–495. [Google Scholar]). This model holds when the matrix of covariances between the variables in the response vector y and the variables in the vector of regressors x has a particular structure. In this paper, we consider tests of hypotheses for this structure when (y′, x′)′ has a multivariate normal distribution. In particular, we investigate the likelihood ratio test and a Wald test.     

Principal points of univariate continuous distributions

The K principal points of a p-variate random variable X are defined as those points ₁,..., _K which minimize the expected squared distance of X from the nearest of the _k . This paper reviews some of the theory of principal points and presents a method of determining principal points of univariate continuous distributions. The method is applied to the uniform distribution, to the normal distribution and to the exponential distribution.     

Inference for a skew extension of the Grubbs model

In this paper, we discuss inferential aspects for the Grubbs model when the unknown quantity x (latent response) follows a skew-normal distribution, extending early results given in Arellano-Valle et al. (J Multivar Anal 96:265–281, 2005b). Maximum likelihood parameter estimates are computed via the EM-algorithm. Wald and likelihood ratio type statistics are used for hypothesis testing and we explain the apparent failure of the Wald statistics in detecting skewness via the profile likelihood function. The results and methods developed in this paper are illustrated with a numerical example.     

Probability computations for an extension of the geometric distribution

An extension of the stochastic process associated with the geometric distribution is presented. Combinatorial arguments are used to derive probabilities for various events of interest. Probabilities are approximated by evaluating truncated series. Bounds on the errors of approximation are developed. An example is presented and some additional applications are noted.     

Optimal allocation of time points for the random effects model

In this article the problem of the optimal selection and allocation of time points in repeated measures experiments is considered. D‐ optimal designs for linear regression models with a random intercept and first order auto‐regressive serial correlations are computed numerically and compared with designs having equally spaced time points. When the order of the polynomial is known and the serial correlations are not too small, the comparison shows that for any fixed number of repeated measures, a design with equally spaced time points is almost as efficient as the D‐ optimal design. When, however, there is no prior knowledge about the order of the underlying polynomial, the best choice in terms of efficiency is a D‐ optimal design for the highest possible relevant order of the polynomial. A design with equally‐spaced time points is the second best choice     

Semiparametric method for detecting multiple change points model in financial time series

The problem of multiple change points has been discussed in these years on the background of financial shocks. In order to decrease the damage, it is worthy to find a more available model for the problem as precise as possible by the information from data set. This paper proposes the problem of detecting the change points by semiparametric test. The change points estimations are obtained by empirical likelihood method. Then some asymptotic results for multiple change points are obtained by loglikelihood ratio test and law of large numbers. Furthermore, the consistency of change points estimations is presented. Indeed, the method and steps to find the change points are derived. The simulation experiments prove that the semiparametric test is more efficient than nonparametric test. The diagnosis with simulation and the applications for multiple change points also illustrates the proposed model well.     

A 2-dimensional extension of the Bradley–Terry model for paired comparisons

The Bradley–Terry model is widely and often beneficially used to rank objects from paired comparisons. The underlying assumption that makes ranking possible is the existence of a latent linear scale of merit or equivalently of a kind of transitiveness of the preference. However, in some situations such as sensory comparisons of products, this assumption can be unrealistic. In these contexts, although the Bradley–Terry model appears to be significantly interesting, the linear ranking does not make sense. Our aim is to propose a 2-dimensional extension of the Bradley–Terry model that accounts for interactions between the compared objects. From a methodological point of view, this proposition can be seen as a multidimensional scaling approach in the context of a logistic model for binomial data. Maximum likelihood is investigated and asymptotic properties are derived in order to construct confidence ellipses on the diagram of the 2-dimensional scores. It is shown by an illustrative example based on real sensory data on how to use the 2-dimensional model to inspect the lack-of-fit of the Bradley–Terry model.     

Estimation of principal points for a multivariate binary distribution using a log-linear model

This article proposes a method for estimating principal points for a multivariate binary distribution, assuming a log-linear model for the distribution. Through numerical simulation studies, the proposed parametric estimation method using a log-linear model is compared with a nonparametric estimation method.     

Bayesian analysis for a skew extension of the multivariate null intercept measurement error model

Skew-normal distribution is a class of distributions that includes the normal distributions as a special case. In this paper, we explore the use of Markov Chain Monte Carlo (MCMC) methods to develop a Bayesian analysis in a multivariate, null intercept, measurement error model [R. Aoki, H. Bolfarine, J.A. Achcar, and D. Leão Pinto Jr, Bayesian analysis of a multivariate null intercept error-in-variables regression model, J. Biopharm. Stat. 13(4) (2003b), pp. 763–771] where the unobserved value of the covariate (latent variable) follows a skew-normal distribution. The results and methods are applied to a real dental clinical trial presented in [A. Hadgu and G. Koch, Application of generalized estimating equations to a dental randomized clinical trial, J. Biopharm. Stat. 9 (1999), pp. 161–178].     

On the effect of the number of quadrature points in a logistic random effects model: an example

Although generalized linear mixed models are recognized to be of major practical importance, it is also known that they can be computationally demanding. The problem is the evaluation of the integral in calculating the marginalized likelihood. The straightforward method is based on the Gauss–Hermite technique, based on Gaussian quadrature points. Another approach is provided by the class of penalized quasi-likelihood methods. It is commonly believed that the Gauss–Hermite method works relatively well in simple situations but fails in more complicated structures. However, we present here a strikingly simple example of a logistic random-intercepts model in the context of a longitudinal clinical trial where the method gives valid results only for a high number of quadrature points ( Q ). As a consequence, this result warns the practitioner to examine routinely the dependence of the results on Q . The adaptive Gaussian quadrature, as implemented in the new SAS procedure NLMIXED, offered the solution to our problem. However, even the adaptive version of Gaussian quadrature needs careful handling to ensure convergence.     

An extension of geisser's discrimination model to proportional covariance matrices

The predictive discrimination model of Geisser is derived for the case in which the covariance matrices of the different populations are assumed proportional, with unknown constants of proportionality. It is shown that this model provides a better fit to a set of data than does the usual homoscedastic model.     

On uniformly distributed random points in an n-Ball

Consider a set of r+1 independently and identically and uniformly distributed random points X₀, X₁,…,X_r in RⁿThese points determine almost surely via their convex hull a unique r-simplex in R^e This article deals with the exact density of the r-content of this random r-simplex when the points are such that p of them are in the interior and r+l?p of them are on the surface of a unit n-ball. This problem is transformed into a distribution problem connected with multivariate test statistics. Various possible representations of the exact density in the general case, are also pointed out.     

Estimation of critical points in the mixture inverse Gaussian model

The maximum likelihood estimation for the critical points of the failure rate and the mean residual life function are presented in the case of mixture inverse Gaussian model. Several important data sets are analyzed from this point of view. For each of the data sets, Bootstrapping is used to construct confidence intervals of the critical points.     

Weibull extension of bivariate exponential regression model with different frailty distributions

We propose bivariate Weibull regression model with frailty in which dependence is generated by a gamma or positive stable or power variance function distribution. We assume that the bivariate survival data follows bivariate Weibull of Hanagal (Econ Qual Control 19:83–90, 2004; Econ Qual Control 20:143–150, 2005a; Stat Pap 47:137–148, 2006a; Stat Methods, 2006b). There are some interesting situations like survival times in genetic epidemiology, dental implants of patients and twin births (both monozygotic and dizygotic) where genetic behavior (which is unknown and random) of patients follows known frailty distribution. These are the situations which motivate to study this particular model. David D. Hanagal is on leave from Department of Statistics, University of Pune, Pune 411007, India.     

An extension of the Freund's bivariate distribution to model cascading failures

The notion of cascading failures is a common phenomenon we observe around us. Here the initial failure alters the structure function of the system, which leads to subsequent failures within a short period of time referred to as threshold time. The concept of cascading failures within the framework of reliability theory and the Freund bivariate exponential distribution to model cascading failures has been studied by few authors. The Freund bivariate exponential distribution allows modelling a parallel redundant system with two components. In this system, the lifetimes of the two components behave as if they are independent, until one of the components fail, after which the remaining component suffers an increased/decreased stress. In this article, we further generalize this model to accommodate cascading failures. Various properties of the model are investigated and statistical inference procedures are developed using L-moments and method of moments. A practical application of this model is illustrated using data from www.espncricinfo.com. Also well analysed Diabetic Retinopathy Study (DRS) data is further analysed using this model and our findings are presented.     

Representative points for location-biased datasets

Representative points (RPs) are a set of points that optimally represents a distribution in terms of mean square error. When the prior data is location biased, the direct methods such as the k-means algorithm may be inefficient to obtain the RPs. In this article, a new indirect algorithm is proposed to search the RPs based on location-biased datasets. Such an algorithm does not constrain the parameter model of the true distribution. The empirical study shows that such algorithm can obtain better RPs than the k-means algorithm.     

Principal varying coefficient estimator for high-dimensional models

ABSTRACT

We consider principal varying coefficient models in the high-dimensional setting, combined with variable selection, to reduce the effective number of parameters in semiparametric modelling. The estimation is based on B-splines approach. For the unpenalized estimator, we establish non-asymptotic bounds of the estimator and then establish the (asymptotic) local oracle property of the penalized estimator, as well as non-asymptotic error bounds. Monte Carlo studies reveal the favourable performance of the estimator and an application on a real dataset is presented.     

A Diagnostic Method for Principal Components Tests

The local influence method is adapted to testing hypotheses about principal components for investigating the influence of observations on the test statistic. Simultaneous perturbations on all observations are considered. The main diagnostic is the direction vector of the maximum slope of the surface formed by the perturbed test statistic. A perturbation is constructed whose result is the same as that of the influence function method. An example is given for illustration.