Maximum Likelihood Estimation of Mixture Densities for Binned and Truncated Multivariate Data

Binning and truncation of data are common in data analysis and machine learning. This paper addresses the problem of fitting mixture densities to multivariate binned and truncated data. The EM approach proposed by McLachlan and Jones (Biometrics, 44: 2, 571-578, 1988) for the univariate case is generalized to multivariate measurements. The multivariate solution requires the evaluation of multidimensional integrals over each bin at each iteration of the EM procedure. Naive implementation of the procedure can lead to computationally inefficient results. To reduce the computational cost a number of straightforward numerical techniques are proposed. Results on simulated data indicate that the proposed methods can achieve significant computational gains with no loss in the accuracy of the final parameter estimates. Furthermore, experimental results suggest that with a sufficient number of bins and data points it is possible to estimate the true underlying density almost as well as if the data were not binned. The paper concludes with a brief description of an application of this approach to diagnosis of iron deficiency anemia, in the context of binned and truncated bivariate measurements of volume and hemoglobin concentration from an individual's red blood cells.

Integrable coupling in a model for Josephson tunneling between non-identical BCS systems

We extend a recent construction for an integrable model describing Josephson tunneling between identical BCS systems to the case where the BCS systems have different single particle energy levels. The exact solution of this generalized model is obtained through the Bethe ansatz.

Field quantization, photons and non-Hermitean modes

Field quantization in unstable optical systems is treated by expanding the vector potential in terms of non-Hermitean (Fox-Li) modes. We define non-Hermitean modes and their adjoints in both the cavity and external regions and make use of the important bi-orthogonality relationships that exist within each mode set. We employ a standard canonical quantization procedure involving the introduction of generalized coordinates and momenta for the electromagnetic (EM) field. Three-dimensional systems are treated, making use of the paraxial and monochromaticity approximations for the cavity non-Hermitean modes. We show that the quantum EM field is equivalent to a set of quantum harmonic oscillators (QHOs), associated with either the cavity or the external region non-Hermitean modes, and thus confirming the validity of the photon model in unstable optical systems. Unlike in the conventional (Hermitean mode) case, the annihilation and creation operators we define for each QHO are not Hermitean adjoints. It is shown that the quantum Hamiltonian for the EM field is the sum of non-commuting cavity and external region contributions, each of which can be expressed as a sum of independent QHO Hamiltonians for each non-Hermitean mode, except that the external field Hamiltonian also includes a coupling term responsible for external non-Hermitean mode photon exchange processes. The non-commutativity of certain cavity and external region annihilation and creation operators is associated with cavity energy gain and loss processes, and may be described in terms of surface integrals involving cavity and external region non-Hermitean mode functions on the cavity-external region boundary. Using the essential states approach and the rotating wave approximation, our results are applied to the spontaneous decay of a two-level atom inside an unstable cavity. We find that atomic transitions leading to cavity non-Hermitean mode photon absorption are associated with a different coupling constant to that for transitions leading to photon emission, a feature consequent on the use of non-Hermitean mode functions. We show that under certain conditions the spontaneous decay rate is enhanced by the Petermann factor.

Existence of a semicontinuous or continuous utility function: a unified approach and an elementary proof

In this paper, we present a new unified approach and an elementary proof of a very general theorem on the existence of a semicontinuous or continuous utility function representing a preference relation. A simple and interesting new proof of the famous Debreu Gap Lemma is given. In addition, we prove a new Gap Lemma for the rational numbers and derive some consequences. We also prove a theorem which characterizes the existence of upper semicontinuous utility functions on a preordered topological space which need not be second countable. This is a generalization of the classical theorem of Rader which only gives sufficient conditions for the existence of an upper semicontinuous utility function for second countable topological spaces. (C) 2002 Elsevier Science B.V. All rights reserved.

Lower Bounds on the State Complexity of Geometric Goppa Codes

We reinterpret the state space dimension equations for geometric Goppa codes. An easy consequence is that if deg G less than or equal to n-2/2 or deg G greater than or equal to n-2/2 + 2g then the state complexity of C-L(D, G) is equal to the Wolf bound. For deg G is an element of [n-1/2, n-3/2 + 2g], we use Clifford's theorem to give a simple lower bound on the state complexity of C-L(D, G). We then derive two further lower bounds on the state space dimensions of C-L(D, G) in terms of the gonality sequence of F/F-q. (The gonality sequence is known for many of the function fields of interest for defining geometric Goppa codes.) One of the gonality bounds uses previous results on the generalised weight hierarchy of C-L(D, G) and one follows in a straightforward way from first principles; often they are equal. For Hermitian codes both gonality bounds are equal to the DLP lower bound on state space dimensions. We conclude by using these results to calculate the DLP lower bound on state complexity for Hermitian codes.

Factorizations of and by powers of complete graphs

Let K-k(d) denote the Cartesian product of d copies of the complete graph K-k. We prove necessary and sufficient conditions for the existence of a K-k(r)-factorization of K-pn(s), where p is prime and k > 1, n, r and s are positive integers. (C) 2002 Elsevier Science B.V. All rights reserved.

Propagation of nonstationary curved and stretched premixed flames with multistep reaction mechanisms

The propagation speed of a thin premixed flame disturbed by an unsteady fluid flow of a larger scale is considered. The flame may also have a general shape but the reaction zone is assumed to be thin compared to the flame thickness. Unlike in preceding publications, the presented asymptotic analysis is performed for a general multistep reaction mechanism and, at the same time, the flame front is curved by the fluid flow. The resulting equations define the propagation speed of disturbed flames in terms of the properties of undisturbed planar flames and the flame stretch. Special attention is paid to the near-equidiffusion limit. In this case, the flame propagation speed is shown to depend on the effective Zeldovich number Z(f) , and the flame stretch. Unlike the conventional Zeldovich number, the effective Zeldovich number is not necessarily linked directly to the activation energies of the reactions. Several examples of determining the effective Zeldovich number for reduced combustion mechanisms are given while, for realistic reactions, the effective Zeldovich number is determined from experiments. Another feature of the present approach is represented by the relatively simple asymptotic technique based on the adaptive generalized curvilinear system of coordinates attached to the flame (i.e., intrinsic disturbed flame equations [IDFE]).

Lead stimulus modality change and the attentional modulation of the acoustic and electrical blink reflex

Two experiments investigated the effects of the sensory modality of the lead and of the blink-eliciting stimulus during lead stimulus modality change on blink modulation at lead intervals of 2500 and 3500 ins. Participants were presented with acoustic, visual, or tactile change stimuli after habituation training with lead stimuli from the same or a different sensory modality. In Experiment 1, latency and magnitude of the acoustic blink were facilitated during a change to acoustic or visual lead stimuli, but not during a change to tactile lead stimuli. After habituation to acoustic lead stimuli, blink magnitude was smaller during tactile change stimuli than during habituation stimuli. The latter finding was replicated in Experiment 2 in which blink was elicited by electrical stimulation of the trigeminal nerve. The consistency of the findings across different combinations of lead stimulus and blink-eliciting stimulus modalities does not support a modality-specific account of attentional blink modulation. Rather, blink modulation during generalized orienting reflects modality non-specific processes, although modulation may not always be found during tactile lead stimuli. (C) 2002 Elsevier Science B.V. All rights reserved.

The epileptology of Theodore Herpin (1799-1865)

The Frenchman, Theodore Herpin (1799-1865), in Des Acces Incomplets d'Epilepsie, published posthumously in 1867, provided a very detailed account of a wide range of the possible manifestations of nonconvulsive epileptic seizures. However, he did not note the presence of absence seizures in any of his 300 patients who had experienced, at least in some of their attacks, what he considered were incomplete manifestations of epilepsy, the word epilepsy being taken to refer to full generalized tonic-clonic seizures. In the one patient, Herpin recognized that all epileptic seizures, whether complete or incomplete, began in the same way, and deduced that they must originate in the same place in that patient's brain. He did not develop the latter idea further. His observations, and his interpretation of them, seem to have preceded John Hughlings Jackson's independent development of similar concepts, but Jackson's more extensive intellectual exploration of the implications of his observations made him a more important figure than Herpin in the history of epileptology.

Doramectin as a treatment for canine and feline demodicosis

Twenty-three dogs and three cats with generalised demodicosis were treated with weekly subcutaneous injections of doramectin at a dose of 600mcg/kg body weight. All dogs and cats responded to treatment by going into remission. The median time until skin scrape results were negative was eight weeks (range five to 20 weeks). Ten dogs remained in remission after the first treatment, five were lost to follow up and seven needed a second course of doramectin or were maintained in remission by monthly injections. Time to when skin scrape results were negative for the cats was two to three weeks with one cat still in remission four years later. The other two cats were euthanased due to their underlying disease whilst only four and six months into remission. Weekly subcutaneous doramectin injections are useful in the treatment of generalised demodicosis in dogs and cats. No systemic side effects of the drug were seen, the injection was painless and no reactions were seen at the sites.

Overlarge sets of 2-(11, 5, 2) designs and related configurations

We consider the construction of several configurations, including: • overlarge sets of 2-(11,5,2) designs, that is, partitions of the set of all 5-subsets of a 12-set into 72 2-(11,5,2) designs; • an indecomposable doubly overlarge set of 2-(11,5,2) designs, that is, a partition of two copies of the set of all 5-subsets of a 12-set into 144 2-(11,5,2) designs, such that the 144 designs can be arranged into a 12 × 12 square with interesting row and column properties; • a partition of the Steiner system S(5,6,12) into 12 disjoint 2-(11,6,3) designs arising from the diagonal of the square; • bidistant permutation arrays and generalized Room squares arising from the doubly overlarge set, and their relation to some new strongly regular graphs.

Conformationally constrained macrocycles that mimic tripeptide beta-strands in water and aprotic solvents

The beta-strand conformation is unknown for short peptides in aqueous solution, yet it is a fundamental building block in proteins and the crucial recognition motif for proteolytic enzymes that enable formation and turnover of all proteins. To create a generalized scaffold as a peptidomimetic that is preorganized in a beta-strand, we individually synthesized a series of 15-22-membered macrocyclic analogues of tripeptides and analyzed their structures. Each cycle is highly constrained by two trans amide bonds and a planar aromatic ring with a short nonpeptidic linker between them. A measure of this ring strain is the restricted rotation of the component tyrosinyl aromatic ring (DeltaG(rot) 76.7 kJ mol(-1) (16-membered ring), 46.1 kJ mol(-1) (17-membered ring)) evidenced by variable temperature proton NMR spectra (DMF-d(7), 200-400 K). Unusually large amide coupling constants ((3)J(NH-CHalpha) 9-10 Hz) corresponding to large dihedral angles were detected in both protic and aprotic solvents for these macrocycles, consistent with a high degree of structure in solution. The temperature dependence of all amide NH chemical shifts (Deltadelta/T7-12 ppb/deg) precluded the presence of transannular hydrogen bonds that define alternative turn structures. Whereas similar sized conventional cyclic peptides usually exist in solution as an equilibrium mixture of multiple conformers, these macrocycles adopt a well-defined beta-strand structure even in water as revealed by 2-D NMR spectral data and by a structure calculation for the smallest (15-membered) and most constrained macrocycle. Macrocycles that are sufficiently constrained to exclusively adopt a beta-strand-mimicking structure in water may be useful pre-organized and generic templates for the design of compounds that interfere with beta-strand recognition in biology.

Fitting genetic models to twin data with binary and ordered categorical responses: A comparison of structural equation modelling and Bayesian hierarchical models

We compare Bayesian methodology utilizing free-ware BUGS (Bayesian Inference Using Gibbs Sampling) with the traditional structural equation modelling approach based on another free-ware package, Mx. Dichotomous and ordinal (three category) twin data were simulated according to different additive genetic and common environment models for phenotypic variation. Practical issues are discussed in using Gibbs sampling as implemented by BUGS to fit subject-specific Bayesian generalized linear models, where the components of variation may be estimated directly. The simulation study (based on 2000 twin pairs) indicated that there is a consistent advantage in using the Bayesian method to detect a correct model under certain specifications of additive genetics and common environmental effects. For binary data, both methods had difficulty in detecting the correct model when the additive genetic effect was low (between 10 and 20%) or of moderate range (between 20 and 40%). Furthermore, neither method could adequately detect a correct model that included a modest common environmental effect (20%) even when the additive genetic effect was large (50%). Power was significantly improved with ordinal data for most scenarios, except for the case of low heritability under a true ACE model. We illustrate and compare both methods using data from 1239 twin pairs over the age of 50 years, who were registered with the Australian National Health and Medical Research Council Twin Registry (ATR) and presented symptoms associated with osteoarthritis occurring in joints of the hand.

Bridging the gap between different statistical approaches: An integrated framework for modelling

This paper proposes a template for modelling complex datasets that integrates traditional statistical modelling approaches with more recent advances in statistics and modelling through an exploratory framework. Our approach builds on the well-known and long standing traditional idea of 'good practice in statistics' by establishing a comprehensive framework for modelling that focuses on exploration, prediction, interpretation and reliability assessment, a relatively new idea that allows individual assessment of predictions. The integrated framework we present comprises two stages. The first involves the use of exploratory methods to help visually understand the data and identify a parsimonious set of explanatory variables. The second encompasses a two step modelling process, where the use of non-parametric methods such as decision trees and generalized additive models are promoted to identify important variables and their modelling relationship with the response before a final predictive model is considered. We focus on fitting the predictive model using parametric, non-parametric and Bayesian approaches. This paper is motivated by a medical problem where interest focuses on developing a risk stratification system for morbidity of 1,710 cardiac patients given a suite of demographic, clinical and preoperative variables. Although the methods we use are applied specifically to this case study, these methods can be applied across any field, irrespective of the type of response.

Indirect evidence of density-dependent population regulation in Aponomma hydrosauri (Acari : Ixodidae), an ectoparasite of reptiles

The extent to which density-dependent processes regulate natural populations is the subject of an ongoing debate. We contribute evidence to this debate showing that density-dependent processes influence the population dynamics of the ectoparasite Aponomma hydrosauri (Acari: Ixodidae), a tick species that infests reptiles in Australia. The first piece of evidence comes from an unusually long-term dataset on the distribution of ticks among individual hosts. If density-dependent processes are influencing either host mortality or vital rates of the parasite population, and those distributions can be approximated with negative binomial distributions, then general host-parasite models predict that the aggregation coefficient of the parasite distribution will increase with the average intensity of infections. We fit negative binomial distributions to the frequency distributions of ticks on hosts, and find that the estimated aggregation coefficient k increases with increasing average tick density. This pattern indirectly implies that one or more vital rates of the tick population must be changing with increasing tick density, because mortality rates of the tick's main host, the sleepy lizard, Tiliqua rugosa, are unaffected by changes in tick burdens. Our second piece of evidence is a re-analysis of experimental data on the attachment success of individual ticks to lizard hosts using generalized linear modelling. The probability of successful engorgement decreases with increasing numbers of ticks attached to a host. This is direct evidence of a density-dependent process that could lead to an increase in the aggregation coefficient of tick distributions described earlier. The population-scale increase in the aggregation coefficient is indirect evidence of a density-dependent process or processes sufficiently strong to produce a population-wide pattern, and thus also likely to influence population regulation. The direct observation of a density-dependent process is evidence of at least part of the responsible mechanism.