ON CHARACTERIZATIONS THROUGH MIXED SUMS

Let X, X‘ and X“ be independent and have the same law, and let U be uniform on [0,1] and independent of X’ and XI”. A problem which has received considerable attention is to identify the laws of X for which X =c U(X‘+X“), the equality being in law. If X is positive it has an exponential law, and mixture solutions exist without the sign restriction. The problem has been extended to allow more summands on either side and to let U have a special beta law. In the positive case X then has a gamma law. This problem is tackled in full generality by considering an integral equation for the characteristic function of X. Uniqueness results and necessary conditions are given. The general case is reduced to that where X is positive, and the solution of the corresponding integral equation is constructed under the necessary conditions. This is shown to characterize all possible laws of X when there is a single summand on the left and U has an arbitrary law on [0,1]. Examples are given and known particular results are sharpened.     

A COMPARISON OF MIXED MODEL SPLINES FOR CURVE FITTING

Three types of polynomial mixed model splines have been proposed: smoothing splines, P‐splines and penalized splines using a truncated power function basis. The close connections between these models are demonstrated, showing that the default cubic form of the splines differs only in the penalty used. A general definition of the mixed model spline is given that includes general constraints and can be used to produce natural or periodic splines. The impact of different penalties is demonstrated by evaluation across a set of functions with specific features, and shows that the best penalty in terms of mean squared error of prediction depends on both the form of the underlying function and the signal:noise ratio.     

LIMIT THEOREMS FOR ASYMPTOTICALLY MINIMAX ESTIMATION OF A DISTRIBUTION WITH INCREASING FAILURE RATE UNDER A RANDOM MIXED CENSORSHIP/TRUNCATION MODEL

ABSTRACT

The search for optimal non-parametric estimates of the cumulative distribution and hazard functions under order constraints inspired at least two earlier classic papers in mathematical statistics: those of Kiefer and Wolfowitz[1] Kiefer, J. and Wolfowitz, J. 1976. Asymptotically Minimax Estimation of Concave and Convex Distribution Functions. Z. Wahrsch. Verw. Gebiete, 34: 73–85. [Crossref], [Web of Science ®] , [Google Scholar] and Grenander[2] Grenander, U. 1956. On the Theory of Mortality Measurement. Part II. Scand. Aktuarietidskrift J., 39: 125–153.  [Google Scholar] respectively. In both cases, either the greatest convex minorant or the least concave majorant played a fundamental role. Based on Kiefer and Wolfowitz's work, Wang3-4 Wang, J.L. 1986. Asymptotically Minimax Estimators for Distributions with Increasing Failure Rate. Ann. Statist., 14: 1113–1131. Wang, J.L. 1987. Estimators of a Distribution Function with Increasing Failure Rate Average. J. Statist. Plann. Inference, 16: 415–427.   found asymptotically minimax estimates of the distribution function F and its cumulative hazard function Λ in the class of all increasing failure rate (IFR) and all increasing failure rate average (IFRA) distributions. In this paper, we will prove limit theorems which extend Wang's asymptotic results to the mixed censorship/truncation model as well as provide some other relevant results. The methods are illustrated on the Channing House data, originally analysed by Hyde.5-6 Hyde, J. 1977. Testing Survival Under Right Censoring and Left Truncation. Biometrika, 64: 225–230. Hyde, J. 1980. “Survival Analysis with Incomplete Observations”. In Biostatistics Casebook 31–46. New York: Wiley Series in Probability and Mathematical Statistics: Applied Probability and Statistics.       

LOG-LINEAR MODELS FOR MEAN AND DISPERSION IN MIXED POISSON REGRESSION MODELS

This paper is concerned with the analysis of repeated measures count data overdispersed relative to a Poisson distribution, with the overdispersion possibly heterogeneous. To accommodate the overdispersion, the Poisson random variable is compounded with a gamma random variable, and both the mean of the Poisson and the variance of the gamma are modelled using log linear models. Maximum likelihood estimates (MLE) are then obtained. The paper also gives extended quasi-likelihood estimates for a more general class of compounding distributions which are shown to be approximations to the MLEs obtained for the gamma case. The theory is illustrated by modelling the determination of asbestos fibre intensity on membrane filters mounted on microscope slides.     

BAYESIAN ANALYSIS OF THE ADDITIVE MIXED MODEL FOR RANDOMIZED BLOCK DESIGNS

This paper deals with the Bayesian analysis of the additive mixed model experiments. Consider b randomly chosen subjects who respond once to each of t treatments. The subjects are treated as random effects and the treatment effects are fixed. Suppose that some prior information is available, thus motivating a Bayesian analysis. The Bayesian computation, however, can be difficult in this situation, especially when a large number of treatments is involved. Three computational methods are suggested to perform the analysis. The exact posterior density of any parameter of interest can be simulated based on random realizations taken from a restricted multivariate t distribution. The density can also be simulated using Markov chain Monte Carlo methods. The simulated density is accurate when a large number of random realizations is taken. However, it may take substantial amount of computer time when many treatments are involved. An alternative Laplacian approximation is discussed. The Laplacian method produces smooth and very accurate approximates to posterior densities, and takes only seconds of computer time. An example of a pipeline cracks experiment is used to illustrate the Bayesian approaches and the computational methods.     

FAMILY BUILDING DIFFERENCES BETWEEN SAME SEX AND MIXED SEX FAMILIES IN AUSTRALIA

From this investigation it seems that an aversion among some couples to a same sex family has some relevance in family building. However, although there are significant differences in the proportions of same sex and mixed sex families having an additional child, the overall effect on average family size of the (apparent) attempt to achieve a child of each sex appears to be slight—an increase of less than 3% in the case of the Melbourne data. Despite the evidence that a significantly greater proportion of women with same sex families have another child, relatively few women admit that the sex structure of the children was the main factor in such a decision. Another insight into the relationship between sex structure and family building was that the desire for a child of each sex may have a negative effect on fertility in that couples who have already achieved such a configuration may decide not to have an additional child that they had originally planned to have. Obviously, it is difficult to separate attitudes to the importance of having both sons and daughters from rationalisations associated with the sex structure of the respondent's own family. Nevertheless the desire for children of each sex seems to be related to traditional attitudes to woman's role, Southern European background, an earlier year of birth, and a lower level of education. However, at the same time women in these categories, particularly the last, seem to be relatively less likely to exhibit behaviour indicative of a controlled response to the sex structure of the family; thus among such groups the probability of same sex families having another child is found to be relatively similar to the proportion of mixed sex families having another child. In other words, they are less likely than other women, who care less about having both sons and daughters, to stop at two (or three) children when at least one son and one daughter have been attained. What of the future relationship between proportions of same sex and mixed sex families having another child? Although it would seem that the combined effect of better educational opportunities and less traditional attitudes of each new generation would further reduce the importance attached to having both sons and daughters, at the same time one would expect an increase in a couple's ability to stop family building when two or three children of the desired sex structure had been achieved. Another factor is that possibly an increasing preference for smaller families will outweigh the wish to have a third (or fourth) child to achieve one child of each sex or for any other reason.     

INFERENCE IN DIRICHLET PROCESS MIXED GENERALIZED LINEAR MODELS BY USING MONTE CARLO EM

Generalized linear mixed models are widely used for describing overdispersed and correlated data. Such data arise frequently in studies involving clustered and hierarchical designs. A more flexible class of models has been developed here through the Dirichlet process mixture. An additional advantage of using such mixture models is that the observations can be grouped together on the basis of the overdispersion present in the data. This paper proposes a partial empirical Bayes method for estimating all the model parameters by adopting a version of the EM algorithm. An augmented model that helps to implement an efficient Gibbs sampling scheme, under the non‐conjugate Dirichlet process generalized linear model, generates observations from the conditional predictive distribution of unobserved random effects and provides an estimate of the average number of mixing components in the Dirichlet process mixture. A simulation study has been carried out to demonstrate the consistency of the proposed method. The approach is also applied to a study on outdoor bacteria concentration in the air and to data from 14 retrospective lung‐cancer studies.     

BAYESIAN PREDICTION FOR SPATIAL GENERALISED LINEAR MIXED MODELS WITH CLOSED SKEW NORMAL LATENT VARIABLES

Spatial generalised linear mixed models are used commonly for modelling non‐Gaussian discrete spatial responses. In these models, the spatial correlation structure of data is modelled by spatial latent variables. Most users are satisfied with using a normal distribution for these variables, but in many applications it is unclear whether or not the normal assumption holds. This assumption is relaxed in the present work, using a closed skew normal distribution for the spatial latent variables, which is more flexible and includes normal and skew normal distributions. The parameter estimates and spatial predictions are calculated using the Markov Chain Monte Carlo method. Finally, the performance of the proposed model is analysed via two simulation studies, followed by a case study in which practical aspects are dealt with. The proposed model appears to give a smaller cross‐validation mean square error of the spatial prediction than the normal prior in modelling the temperature data set.     

ASYMPTOTIC PROPERTIES OF RESTRICTED MAXIMUM LIKELIHOOD (REML) ESTIMATES FOR HIERARCHICAL MIXED LINEAR MODELS

This paper explores the asymptotic distribution of the restricted maximum likelihood estimator of the variance components in a general mixed model. Restricting attention to hierarchical models, central limit theorems are obtained using elementary arguments with only mild conditions on the covariates in the fixed part of the model and without having to assume that the data are either normally or spherically symmetrically distributed. Further, the REML and maximum likelihood estimators are shown to be asymptotically equivalent in this general framework, and the asymptotic distribution of the weighted least squares estimator (based on the REML estimator) of the fixed effect parameters is derived.     

SIMULATED MAXIMUM LIKELIHOOD APPLIED TO NON-GAUSSIAN AND NONLINEAR MIXED EFFECTS AND STATE–SPACE MODELS

The paper presents an overview of maximum likelihood estimation using simulated likelihood, including the use of antithetic variables and evaluation of the simulation error of the resulting estimates. It gives a general purpose implementation of simulated maximum likelihood and uses it to re‐visit four models that have previously appeared in the published literature: a state–space model for count data; a nested random effects model for binomial data; a nonlinear growth model with crossed random effects; and a crossed random effects model for binary salamander‐mating data. In the case of the last three examples, this appears to be the first time that maximum likelihood fits of these models have been presented.     

Some new partial orderings on Sn and Sn×Sn

A new class of partial orderings on S_n, the set of permutations of {1,…,n}, is given. Each of these partial orderings is shown to be a subordering of a recently described partial ordering on S_n (Block, Chhetry, Fang and Sampson (1990)) which is related to Schriever's (1987) more associated ordering on bivariate distributions. Also given is an extension of three known partial orders on S_n to partial orders on S_n×S_n. These extensions facilitate the study of functions from S_n×S_n into , which preser these partial orderings, thereby, providing a methodology for extending the notion of arrangement increasing functions.     

A formula for P(bi≤Ri≤ai, 1≤i≤m) under a general class of alternatives

Let X_i≤?≤X_m and Y_i≤?≤Y_n be two sets of independent order statistics from continous distributions with distribution functions F and G respectively. Let R_i denote the rank of X_i in the combined order sample. Steck (1980) has found an expression for P(b_i≤R_i≤a_i, all i) when F = h(G), h being the incomplete beta function with parameters (α,β?α+1). An alternative expression for the same probability is obtained which is computationally a substantial improvement on Steck's result.     

Lp-norm spherical distribution

In this paper, the spherical distribution is generalized to the L_p-norm spherical distribution, S(n, p), through the stochastic representation. Then the form of its p.d.f. characterizing the S(n, p) distribution is derived. Marginal, conditional distributions and the moments of S(n, p) distribution are obtained. The order statistics and the invariance property of S(n,p) distribution are also investigated.     

Robustness of weighted L p–depth and L p–median

Summary: L ^p –norm weighted depth functions are introduced and the local and global robustness of these weighted L ^p –depth functions and their induced multivariate medians are investigated via influence function and finite sample breakdown point. To study the global robustness of depth functions, a notion of finite sample breakdown point is introduced. The weighted L ^p –depth functions turn out to have the same low breakdown point as some other popular depth functions. Their influence functions are also unbounded. On the other hand, the weighted L ^p –depth induced medians are globally robust with the highest possible breakdown point for any reasonable estimator. The weighted L ^p –medians are also locally robust with bounded influence functions for suitable weight functions. Unlike other existing depth functions and multivariate medians, the weighted L ^p depth and medians are easy to calculate in high dimensions. The price for this advantage is the lack of affine invariance and equivariance of the weighted L ^p depth and medians, respectively.*The author thanks the referees for their very insightful and constructive comments and suggestions which led to corrections and substantial improvements. Supported in part by NSF Grants DMS-0071976 and DMS-0134628.     

Bayesian analysis of Er/M/1 and Er/M/c queues

Bayesian inference and prediction tasks for Er/M/1 and Er/M/c queues are undertaken. Equilibrium probabilities of the queue size and waiting time distributions are estimated using conditional Monte-Carlo simulation methods. We illustrate that some standard queueing measures do not exist when independent priors are used for the arrival and service rates of a G/M/1 queue.     

Lattice paths combinatorics applied to transient queue length distribution of C2/M/1 queues and busy period analysis of bulk queues C2/M/1

This paper aims at deriving explicit transient queue length distribution for GI/M/1 system and busy period analysis of bulk queue GI^b/M/1 through lattice paths (LPs) combinatorics. The general interarrival time distribution is approximated by two-phase Cox distribution, C₂, that has Markovian property, enabling us to represent the processes by two-dimensional LPs. As distributions C₂ cover a wide class of distributions that have rational Laplace–Stieltjes transforms (LSTs) with square coefficient of variation lying in , the results obtained are applicable to a large class of real life situations. Some numerical results for the C₂^b/M/1 model are also given.     

An improved Cp criterion for spline smoothing

Spline smoothing is a popular technique for curve fitting, in which selection of the smoothing parameter is crucial. Many methods such as Mallows’ C_p, generalized maximum likelihood (GML), and the extended exponential (EE) criterion have been proposed to select this parameter. Although C_p is shown to be asymptotically optimal, it is usually outperformed by other selection criteria for small to moderate sample sizes due to its high variability. On the other hand, GML and EE are more stable than C_p, but they do not possess the same asymptotic optimality as C_p. Instead of selecting this smoothing parameter directly using C_p, we propose to select among a small class of selection criteria based on Stein's unbiased risk estimate (SURE). Due to the selection effect, the spline estimate obtained from a criterion in this class is nonlinear. Thus, the effective degrees of freedom in SURE contains an adjustment term in addition to the trace of the smoothing matrix, which cannot be ignored in small to moderate sample sizes. The resulting criterion, which we call adaptive C_p, is shown to have an analytic expression, and hence can be efficiently computed. Moreover, adaptive C_p is not only demonstrated to be superior and more stable than commonly used selection criteria in a simulation study, but also shown to possess the same asymptotic optimality as C_p.     

A modified Cp statistic in a system-of-equations model

A new statistic, SΓ(p), is developed for variable selection in a system-of-equations model. The standardized total mean square error in the SΓ(p)statistic is weighted by the covariance matrix of dependent variables instead of the error covariance matrix of the true model as in the original definition. The new statistic can be also used for model selection in the non-nested models. The estimate of SΓ(p), SC(p), is derived and shown to become SC_ε(p) in the similar form of C_p in a single-equation model when the covariance matrix of sampled dependent variables is replaced by the error covariance matrix under the full model.