Continuity of f-projections and applications to the iterative proportional fitting procedure*

This paper proves continuity of f-projections and the continuous dependence of the limit matrix of the iterative proportional fitting procedure (IPF procedure) on the given matrix as well as the given marginals under certain regularity constraints. For finite spaces, the concept of f-projections of finite measures on a compact and convex set is introduced and continuity of f-projections is proven. This result is applied to the IPF procedure. Given a nonnegative matrix as well as row and column marginals the IPF procedure generates a sequence of matrices, called the IPF sequence, by alternately fitting rows and columns to match their respective marginals. The procedure is equivalent to cyclic f-projections. If the IPF sequence converges, the application of the continuity of f-projections yields the continuous dependence of the limit matrix on the given matrix. By generalized convex programming and under some constraints, it is shown that the limit matrix of the IPF sequence continuously depends not only on the given matrix but also on the marginals.     

Regression for non-Euclidean data using distance matrices

Regression methods for common data types such as measured, count and categorical variables are well understood but increasingly statisticians need ways to model relationships between variable types such as shapes, curves, trees, correlation matrices and images that do not fit into the standard framework. Data types that lie in metric spaces but not in vector spaces are difficult to use within the usual regression setting, either as the response and/or a predictor. We represent the information in these variables using distance matrices which requires only the specification of a distance function. A low-dimensional representation of such distance matrices can be obtained using methods such as multidimensional scaling. Once these variables have been represented as scores, an internal model linking the predictors and the responses can be developed using standard methods. We call scoring as the transformation from a new observation to a score, whereas backscoring is a method to represent a score as an observation in the data space. Both methods are essential for prediction and explanation. We illustrate the methodology for shape data, unregistered curve data and correlation matrices using motion capture data from an experiment to study the motion of children with cleft lip.     

A perspective of composite sampling

Composite samples are formed by physically mixing samples. Usually, composite samples are used to reduce the overall cost associated with analytical procedures that must be performed on each sample, but they can also be used to protect the privacy of individuals.

Composite sampling can reduce the cost of identifying individual cases that have a certain trait, such as those with a rare disease or those exceeding pollution-level standards. Not much is lost by applying this method as long as the trait is relatively rare.

Composite sampling can reduce the cost of estimating the mean of some process. When samples are composited, the ability to estimate the variance is lost. In spite of this, the potential savings are so great that composite samples have been used.

Much of this paper deasl with the variance of estimators based on composite sampling when the porportions of hte original samples comprising the composite sample are actually random. Taking repeated samples and measurements on several composite samples complicates the prodcedure, but allows the estimation of between and within variation as well as measurement error.     

Weighted Stochastic Restricted Estimation in Linear Measurement Error Models

This article discusses the consistent estimation of the parameters in a linear measurement error model when stochastic linear restrictions on regression coefficients are available. We propose some methodologies to obtain the consistent estimation when either the covariance matrix of the measurement errors or the reliability matrix of independent variables is known. Their finite- and large-sample properties are derived with not necessarily normal errors. A Monte Carlo simulation is carried out to study the the finite properties of the estimators.     

Finding environmentally important industry clusters: Multiway cut approach using nonnegative matrix factorization

This paper proposes an optimal combinatorial method for finding groups of industries with relatively large CO₂ emissions through industrial relations. Using an economic input–output table, we estimated a non-symmetric matrix describing how much CO₂ is emitted in producing the commodity of industry i, which was purchased to produce commodity of industry j, to meet the final demand for a specific commodity. A symmetric strength of relations matrix describing the CO₂ emissions associated with the industrial relations was further estimated using the non-symmetric matrix. The strength of relations matrix can be viewed as a representation of the supply-chain network of the final commodity. In this study, we estimated the strength of relations matrix associated with the final demand for automobiles and applied the multiway cut approach using nonnegative matrix factorization to the matrix in order to find environmentally important industry clusters in the Japanese automobile supply chain. According to our empirical results, the optimal number of industry clusters is 19, and 4 industry clusters are playing a key role in CO₂ emission reduction.     

On the simple step-stress model for two-parameter exponential distribution

In this paper, we consider the simple step-stress model for a two-parameter exponential distribution, when both the parameters are unknown and the data are Type-II censored. It is assumed that under two different stress levels, the scale parameter only changes but the location parameter remains unchanged. It is observed that the maximum likelihood estimators do not always exist. We obtain the maximum likelihood estimates of the unknown parameters whenever they exist. We provide the exact conditional distributions of the maximum likelihood estimators of the scale parameters. Since the construction of the exact confidence intervals is very difficult from the conditional distributions, we propose to use the observed Fisher Information matrix for this purpose. We have suggested to use the bootstrap method for constructing confidence intervals. Bayes estimates and associated credible intervals are obtained using the importance sampling technique. Extensive simulations are performed to compare the performances of the different confidence and credible intervals in terms of their coverage percentages and average lengths. The performances of the bootstrap confidence intervals are quite satisfactory even for small sample sizes.     

A simplified proof of a theorem on the difference of the moore–penrose inverses of two positive semi–difinte matrices

Simplified proofs are given of a standard result that establishes positive semi–definiteness of the difference of the inverses of two non–singular matrices, and of the extension of this result by Milliken and Akdeniz (1977) to the difference of the Moore–Penrose inverse of two singular matrices.     

On square ordinal contingency tables: a comparison of social class and income mobility for the same individuals

Summary.  Square contingency tables with matching ordinal rows and columns arise in particular as empirical transition matrices and the paper considers these in the context of social class and income mobility tables. Such tables relate the socio-economic position of parents to the socio-economic position of their child in adulthood. The level of association between parental and child socio-economic position is taken as a measure of mobility. Several approaches to analysis are described and illustrated by UK data in which interest focuses on comparisons of social class and income mobility tables that are derived from the same individuals. Account is taken of the use of the same individuals in the two tables. Additionally comparisons over time are considered.     

Optimum chemical balance weighing designs for estimating the total weight

This deals with chemical balance weighing designs which attain a lower bound for the variance of the estimated total weight, The results extend those of Chacko Dey (1978).     

A generalized scheffé method of multiple comparisons

Standard resulrs on the extrema of quotients of quadratic forms are extended to the non-negative definite case. The maximum and the set over which it is achieved are characterized explicitly both in terms of generalized inverse matrices and generalized eigenvalues. These results become the basis of Scheffe type multiple comparisons in the usual way. To demonstrate their application to statistics with singular covariance matrices, the method is detailed for Mantel-Haenszel, Breslow, and Cox statistics. An example is presented illustrating a situation where the proposed Scheffe type comparisons may be better than the pairwise method.