Singular matrix variate beta distribution

In this paper, we determine the symmetrised density of doubly noncentral singular matrix variate beta type I and II distributions under different definitions. As particular cases we obtain the noncentral singular matrix variate beta type I and II distributions and the corresponding joint density of the nonnull eigenvalues. In addition, we propose an alternative approach to find the corresponding nonsymmetrised densities. From the latter, we solve the integral proposed by Constantine [Noncentral distribution problems in multivariate analysis, Ann. Math. Statist. 34 (1963) 1270-1285] and Khatri [A note on Mitra's paper “A density free approach to the matrix variate beta distribution”, Sankhyā A 32 (1970) 311-318] and reconsidered in Farrell [Multivariate Calculation: Use of the Continuous Groups, Springer Series in Statistics, Springer, New York, 1985, p. 191], see also Díaz-García and Gutiérrez-Jáimez [Noncentral matrix variate beta distribution, Comunicación Técnica, No. I-06-06 (PE/CIMAT), Guanajuato, México, 2006, 〈http://www.cimat.mx/biblioteca/RepTec/index.html?m=2〉], for the singular and nonsingular cases.     

Wishart ratios with dependent structure: New members of the bimatrix beta type IV

In multivariate statistics under normality, the problems of interest are random covariance matrices (known as Wishart matrices) and “ratios” of Wishart matrices that arise in multivariate analysis of variance (MANOVA) (see 24). The bimatrix variate beta type IV distribution (also known in the literature as bimatrix variate generalised beta; matrix variate generalization of a bivariate beta type I) arises from “ratios” of Wishart matrices. In this paper, we add a further independent Wishart random variate to the “denominator” of one of the ratios; this results in deriving the exact expression for the density function of the bimatrix variate extended beta type IV distribution. The latter leads to the proposal of the bimatrix variate extended F distribution. Some interesting characteristics of these newly introduced bimatrix distributions are explored. Lastly, we focus on the bivariate extended beta type IV distribution (that is an extension of bivariate Jones’ beta) with emphasis on P(X₁<X₂) where X₁ is the random stress variate and X₂ is the random strength variate.     

Complex bimatrix variate generalised beta distributions

In this paper, the study of bivariate generalised beta types I and II distributions is extended to the complex matrix variate case, for which the corresponding density functions are found. In addition, for complex bimatrix variate beta type I distributions, several basic properties, including the joint eigenvalue density and the maximum eigenvalue distribution, are studied.     

Moments on Catalan numbers

By combining inverse series relations with binomial convolutions and telescoping method, moments of Catalan numbers are evaluated, which resolves a problem recently proposed by Gutiérrez et al. [J.M. Gutiérrez, M.A. Hernández, P.J. Miana, N. Romero, New identities in the Catalan triangle, J. Math. Anal. Appl. 341 (1) (2008) 52-61].     

Exact nonnull distribution of Wilks’ statistic: The ratio and product of independent components

The study of the noncentral matrix variate beta type distributions has been sidelined because the final expressions for the densities depend on an integral that has not been resolved in an explicit way. We derive an exact expression for the nonnull distribution of Wilks’ statistic and precise expressions for the densities of the ratio and product of two independent components of matrix variates where one matrix variate has the noncentral matrix variate beta type I distribution and the other has the matrix variate beta type I distribution. We provide the expressions for the densities of the determinant of the ratio and the product of these two components. These distributions play a fundamental role in various areas of statistics, for example in the criteria proposed by Wilks.     

On Jörgens, Calabi, and Pogorelov type theorem and isolated singularities of parabolic Monge-Ampère equations

In the paper, we extend Jörgens, Calabi, and Pogorelov's theorem on entire solutions of elliptic Monge-Ampère equations to parabolic equations associated with Gauss curvature flows. Our results include Gutiérrez and Huang's previous work as a special case. Besides, we also treat the isolated singularities for parabolic Monge-Ampère equations that was firstly studied by Jörgens for elliptic case in two dimensions.     

A certain critical density property for invariant Harnack inequalities in H-type groups

We consider second order linear degenerate elliptic operators which are elliptic with respect to horizontal directions generating a stratified algebra of H-type. Extending a result by Gutiérrez and Tournier (2011) for the Heisenberg group, we prove a critical density estimate by assuming a condition of Cordes–Landis type. We then deduce an invariant Harnack inequality for the non-negative solutions from a result by Di Fazio, Gutiérrez, and Lanconelli (2008).     

Recurrence relations for a Newton-like method in Banach spaces

The convergence of iterative methods for solving nonlinear operator equations in Banach spaces is established from the convergence of majorizing sequences. An alternative approach is developed to establish this convergence by using recurrence relations. For example, the recurrence relations are used in establishing the convergence of Newton's method [L.B. Rall, Computational Solution of Nonlinear Operator Equations, Robert E. Krieger, New York, 1979] and the third order methods such as Halley's, Chebyshev's and super Halley's [V. Candela, A. Marquina, Recurrence relations for rational cubic methods I: the Halley method, Computing 44 (1990) 169–184; V. Candela, A. Marquina, Recurrence relations for rational cubic methods II: the Halley method, Computing 45 (1990) 355–367; J.A. Ezquerro, M.A. Hernández, Recurrence relations for Chebyshev-type methods, Appl. Math. Optim. 41 (2000) 227–236; J.M. Gutiérrez, M.A. Hernández, Third-order iterative methods for operators with bounded second derivative, J. Comput. Appl. Math. 82 (1997) 171–183; J.M. Gutiérrez, M.A. Hernández, Recurrence relations for the Super–Halley method, Comput. Math. Appl. 7(36) (1998) 1–8; M.A. Hernández, Chebyshev's approximation algorithms and applications, Comput. Math. Appl. 41 (2001) 433–445 [10]].     

Inexact Newton-type methods

We introduce the new idea of recurrent functions to provide a semilocal convergence analysis for an inexact Newton-type method, using outer inverses. It turns out that our sufficient convergence conditions are weaker than in earlier studies in many interesting cases (Argyros, 2004 [5] and [6], Argyros, 2007 [7], Dennis, 1971 [14], Deuflhard and Heindl, 1979 [15], Gutiérrez, 1997 [16], Gutiérrez et al., 1995 [17], Häubler, 1986 [18], Huang, 1993 [19], Kantorovich and Akilov, 1982 [20], Nashed and Chen, 1993 [21], Potra, 1982 [22], Potra, 1985 [23]).     

Singular random matrix decompositions: distributions

Assuming that Y has a singular matrix variate elliptically contoured distribution with respect to the Hausdorff measure, the distributions of several matrices associated to QR, modified QR, SV and polar decompositions of matrix Y are determined, for central and non-central, non-singular and singular cases, as well as their relationship to the Wishart and pseudo-Wishart generalized singular and non-singular distributions. Some of these results are also applied to two particular subfamilies of elliptical distributions, the singular matrix variate normal distribution and the singular matrix variate symmetric Pearson type VII distribution.     

Distributions of the compound and scale mixture of vector and spherical matrix variate elliptical distributions

Several matrix variate hypergeometric type distributions are derived. The compound distributions of left-spherical matrix variate elliptical distributions and inverted hypergeometric type distributions with matrix arguments are then proposed. The scale mixture of left-spherical matrix variate elliptical distributions and univariate inverted hypergeometric type distributions is also derived as a particular case of the compound distribution approach.     

Third-order iterative methods with applications to Hammerstein equations: A unified approach

The geometrical interpretation of a family of higher order iterative methods for solving nonlinear scalar equations was presented in [S. Amat, S. Busquier, J.M. Gutiérrez, Geometric constructions of iterative functions to solve nonlinear equations. J. Comput. Appl. Math. 157(1) (2003) 197-205]. This family includes, as particular cases, some of the most famous third-order iterative methods: Chebyshev methods, Halley methods, super-Halley methods, C-methods and Newton-type two-step methods. The aim of the present paper is to analyze the convergence of this family for equations defined between two Banach spaces by using a technique developed in [J.A. Ezquerro, M.A. Hernández, Halley’s method for operators with unbounded second derivative. Appl. Numer. Math. 57(3) (2007) 354-360]. This technique allows us to obtain a general semilocal convergence result for these methods, where the usual conditions on the second derivative are relaxed. On the other hand, the main practical difficulty related to the classical third-order iterative methods is the evaluation of bilinear operators, typically second-order Fréchet derivatives. However, in some cases, the second derivative is easy to evaluate. A clear example is provided by the approximation of Hammerstein equations, where it is diagonal by blocks. We finish the paper by applying our methods to some nonlinear integral equations of this type.     

Factorization of entire mappings of nuclear bounded type

In this paper we study the relationships between the spaces of entire mappings of bounded type, entire mappings of nuclear bounded type, entire mappings of Pietsch integral bounded type, and entire mappings of Grothendieck integral bounded type. Several results due to Alencar (Proc. Roy. Irish Acad.85A (1985) 131-138) and Cilia and Gutiérrez (J. Aust. Math. Soc.76 (2004) 269-280) for homogeneous polynomials are extended to entire mappings. In the main result we prove that an entire mapping is of nuclear bounded type if and only if it factors through an entire mapping of Pietsch integral bounded type.     

From Archimedean to Liouville copulas

We use a recent characterization of the d-dimensional Archimedean copulas as the survival copulas of d-dimensional simplex distributions (McNeil and Nešlehová (2009) [1]) to construct new Archimedean copula families, and to examine the relationship between their dependence properties and the radial parts of the corresponding simplex distributions. In particular, a new formula for Kendall’s tau is derived and a new dependence ordering for non-negative random variables is introduced which generalises the Laplace transform order. We then generalise the Archimedean copulas to obtain Liouville copulas, which are the survival copulas of Liouville distributions and which are non-exchangeable in general. We derive a formula for Kendall’s tau of Liouville copulas in terms of the radial parts of the corresponding Liouville distributions.     

Maximum entropy characterizations of the multivariate Liouville distributions

A random vector X=(X₁,X₂,…,X_n) with positive components has a Liouville distribution with parameter θ=(θ₁,θ₂,…,θ_n) if its joint probability density function is proportional to , θ_i>0 [R.D. Gupta, D.S.P. Richards, Multivariate Liouville distributions, J. Multivariate Anal. 23 (1987) 233-256]. Examples include correlated gamma variables, Dirichlet and inverted Dirichlet distributions. We derive appropriate constraints which establish the maximum entropy characterization of the Liouville distributions among all multivariate distributions. Matrix analogs of the Liouville distributions are considered. Some interesting results related to I-projection from a Liouville distribution are presented.     

On Newton’s method for solving equations containing Fréchet-differentiable operators of order at least two

We provide sufficient conditions for the semilocal convergence of Newton’s method to a locally unique solution of a nonlinear operator equation containing operators that are Fréchet-differentiable of order at least two, in a Banach space setting. Numerical examples are also provided to show that our results apply to solve nonlinear equations in cases earlier ones cannot [J.M. Gutiérrez, A new semilocal convergence theorem for Newton’s method, J. Comput. Appl. Math. 79(1997) 131-145; Z. Huang, A note of Kantorovich theorem for Newton iteration, J. Comput. Appl. Math. 47 (1993) 211-217; F.A. Potra, Sharp error bounds for a class of Newton-like methods, Libertas Mathematica 5 (1985) 71-84].     

Nonparametric estimation of competing risks models with covariates

In competing risks model, several failure times arise potentially. The smallest failure time and its index only are observed. Without specific assumptions, the joint or even the marginal distribution functions of the underlying failure times are not identifiable (A. Tsiatis, Proc. Natl. Acad. Sci. USA 72 (1975) 20). Nonetheless, if each individual is characterized by a “sufficiently informative” set of covariates, these distributions are identifiable under some conditions of regularity (J.J. Heckman and B. Honoré, Biometrika 76 (1989) 325). In this paper, nonparametric kernel estimators of the joint distribution function of failure times conditional on the covariates are proposed. Their weak and strong consistency are discussed.     

On Mardia's and Song's measures of kurtosis in elliptical distributions

The main objective of this paper is the calculation and the comparative study of two general measures of multivariate kurtosis, namely Mardia's measure β_2,p and Song's measure S(f). In this context, general formulas for the said measures are derived for the broad family of the elliptically contoured symmetric distributions and also for specific members of this family, like the multivariate t-distribution, the multivariate Pearson type II, the multivariate Pearson type VII, the multivariate symmetric Kotz type distribution and the uniform distribution in the unit sphere. Analytic expressions for computing Shannon and Rényi entropies are obtained under the elliptic family. The behaviour of Mardia's and Song's measures, their similarities and differences, possible interpretations and uses in practice are investigated by comparing them in specific members of the elliptic family of multivariate distributions. An empirical estimator of Song's measure is moreover proposed and its asymptotic distribution is investigated under the elliptic family of multivariate distributions.     

Nonparametric estimation of the dependence function for a multivariate extreme value distribution

Understanding and modeling dependence structures for multivariate extreme values are of interest in a number of application areas. One of the well-known approaches is to investigate the Pickands dependence function. In the bivariate setting, there exist several estimators for estimating the Pickands dependence function which assume known marginal distributions [J. Pickands, Multivariate extreme value distributions, Bull. Internat. Statist. Inst., 49 (1981) 859-878; P. Deheuvels, On the limiting behavior of the Pickands estimator for bivariate extreme-value distributions, Statist. Probab. Lett. 12 (1991) 429-439; P. Hall, N. Tajvidi, Distribution and dependence-function estimation for bivariate extreme-value distributions, Bernoulli 6 (2000) 835-844; P. Capéraà, A.-L. Fougères, C. Genest, A nonparametric estimation procedure for bivariate extreme value copulas, Biometrika 84 (1997) 567-577]. In this paper, we generalize the bivariate results to p-variate multivariate extreme value distributions with p?2. We demonstrate that the proposed estimators are consistent and asymptotically normal as well as have excellent small sample behavior.