Spatial design matrices and associated quadratic forms: structure and properties

The paper provides significant simplifications and extensions of results obtained by Gorsich, Genton, and Strang (J. Multivariate Anal. 80 (2002) 138) on the structure of spatial design matrices. These are the matrices implicitly defined by quadratic forms that arise naturally in modelling intrinsically stationary and isotropic spatial processes. We give concise structural formulae for these matrices, and simple generating functions for them. The generating functions provide formulae for the cumulants of the quadratic forms of interest when the process is Gaussian, second-order stationary and isotropic. We use these to study the statistical properties of the associated quadratic forms, in particular those of the classical variogram estimator, under several assumptions about the actual variogram.     

Analysis of two-sample truncated data using generalized logistic models

Parallel to Cox's [JRSS B34 (1972) 187-230] proportional hazards model, generalized logistic models have been discussed by Anderson [Bull. Int. Statist. Inst. 48 (1979) 35-53] and others. The essential assumption is that the two densities ratio has a known parametric form. A nice property of this model is that it naturally relates to the logistic regression model for categorical data. In astronomic, demographic, epidemiological, and other studies the variable of interest is often truncated by an associated variable. This paper studies generalized logistic models for the two-sample truncated data problem, where the two lifetime densities ratio is assumed to have the form exp{α+φ(x;β)}. Here φ is a known function of x and β, and the baseline density is unspecified. We develop a semiparametric maximum likelihood method for the case where the two samples have a common truncation distribution. It is shown that inferences for β do not depend the nonparametric components. We also derive an iterative algorithm to maximize the semiparametric likelihood for the general case where different truncation distributions are allowed. We further discuss how to check goodness of fit of the generalized logistic model. The developed methods are illustrated and evaluated using both simulated and real data.     

A triangular map with homoclinic orbits and no infinite ω-limit set containing periodic points

Recently, Forti, Paganoni and Smítal constructed an example of a triangular map of the unite square, F(x,y)=(f(x),g(x,y)), possessing periodic orbits of all periods and such that no infinite ω-limit set of F contains a periodic point. In this note we show that the above quoted map F has a homoclinic orbit. As a consequence, we answer in the negative the problem presented by A.N. Sharkovsky in the eighties whether, for a triangular map of the square, existence of a homoclinic orbit implies the existence of an infinite ω-limit set containing a periodic point. It is well known that, for a continuous map of the interval, the answer is positive.     

Unavoidable Traces Of Set Systems

Sauer, Shelah, Vapnik and Chervonenkis proved that if a set system on n vertices contains many sets, then the set system has full trace on a large set. Although the restriction on the size of the groundset cannot be lifted, Frankl and Pach found a trace structure that is guaranteed to occur in uniform set systems even if we do not bound the size of the groundset. In this note we shall give three sequences of structures such that every set system consisting of sufficiently many sets contains at least one of these structures with many sets.     

Higher Order Peak Algebras

Using the theory of noncommutative symmetric functions, we introduce the higher order peak algebras (Sym(N))_N≥1, a sequence of graded Hopf algebras which contain the descent algebra and the usual peak algebra as initial cases (N=1 and N=2). We compute their Hilbert series, introduce and study several combinatorial bases, and establish various algebraic identities related to the multisection of formal power series with noncommutative coefficients. Received November 19, 2004     

On the Computational Complexity of the Rooted Subtree Prune and Regraft Distance

The graph-theoretic operation of rooted subtree prune and regraft is increasingly being used as a tool for understanding and modelling reticulation events in evolutionary biology. In this paper, we show that computing the rooted subtree prune and regraft distance between two rooted binary phylogenetic trees on the same label set is NP-hard. This resolves a longstanding open problem. Furthermore, we show that this distance is fixed parameter tractable when parameterised by the distance between the two trees.Received March 16, 2004     

Cochran's statistical theorem for outer inverses of matrices and matrix quadratic forms

We extend the matrix version of Cochran's statistical theorem to outer inverses of a matrix. As applications, we investigate the Wishartness and independence of matrix quadratic forms for Kronecker product covariance structures.     

Admissible slopes for monotone and convex interpolation

Summary In many applications, interpolation of experimental data exhibiting some geometric property such as nonnegativity, monotonicity or convexity is unacceptable unless the interpolant reflects these characteristics. This paper identifies admissible slopes at data points of variousC ¹ interpolants which ensure a desirable shape. We discuss this question, in turn for the following function classes commonly used for shape preserving interpolations: monotone polynomials,C ¹ monotone piecewise polynomials, convex polynomials, parametric cubic curves and rational functions.     

EPR and XPS measurements of polymeric catalysts doped with hereopolyacids in oxygen adsorption studies

The goal of our studies was to determine the nature of interactions between catalytic active samples, such as polypyrrole (PPy) doped with: chloride anions, heteropolyacids such as H₅PMo₁₀V₂O₄₀, or H₄SiW₁₂O₄₀, and oxygen particles. In order to reveal the mechanism of the linkage between the catalysts and oxygen we provided the electron paramagnetic resonance (EPR) measurements of synthesized samples in O₂ and N₂ flow. Moreover, the X-ray photoelectron spectroscopy (XPS) measurements were performed to show the properties of doped heteropolyacids and the state of nitrogen in polypyrrole matrix. The results of EPR studies of PPy(Cl), PPy(H₄SiW₁₂O₄₀), PPy(H₅PMo₁₀V₂O₄₀) samples, show that the polypyrrole doped with chloride anions interacts much easier with oxygen than polypyrrole doped with heteropolyacids (in the above given order, of preference with the most active polypyrrole system first). It correlates with XPS results, which show that a molybdenium-vanadium anion is more reduced than a silicon-tungsten one. The redox reactions of heteropolyacids involve the oxidation-reduction of the conjugated polymer chain, leaving no or little place for interactions with oxygen. Polypyrrole doped with H₄SiW₁₂O₄₀ shows some oxygen sensitivity as observed in the EPR studies, as opposed to polypyrrole doped with H₅PMo₁₀V₂O₄₀, which is in line with the XPS results.