A Score Test for Zero-Inflation in Correlated Count Data

To account for the preponderance of zero counts and simultaneous correlation of observations, a class of zero-inflated Poisson mixed regression models is applicable for accommodating the within-cluster dependence. In this paper, a score test for zero-inflation is developed for assessing correlated count data with excess zeros. The sampling distribution and the power of the test statistic are evaluated by simulation studies. The results show that the test statistic performs satisfactorily under a wide range of conditions. The test procedure is further illustrated using a data set on recurrent urinary tract infections. Copyright (c) 2005 John Wiley & Sons, Ltd.

Comparison of blast furnace raceway size with theory

A previously developed one-dimensional mathematical model, to explain raceway hysteresis, is used to predict the raceway diameter in operating blast furnaces and hot models. Raceway size obtained from the open literature under various conditions for various blast furnaces are compared with computed predictions. In addition the predictions are also compared with published outcomes from other hot models. Simulated results on raceway diameter are in very good agreement with published operating blast furnace and hot model data. The effect of various parameters such as tuyere and hearth diameter, coke size and density, void fraction and bed height on raceway diameter has been studied.

Improved algorithms for rare event simulation with heavy tails

The estimation of P(S-n > u) by simulation, where S, is the sum of independent. identically distributed random varibles Y-1,..., Y-n, is of importance in many applications. We propose two simulation estimators based upon the identity P(S-n > u) = nP(S, > u, M-n = Y-n), where M-n = max(Y-1,..., Y-n). One estimator uses importance sampling (for Y-n only), and the other uses conditional Monte Carlo conditioning upon Y1,..., Yn-1. Properties of the relative error of the estimators are derived and a numerical study given in terms of the M/G/1 queue in which n is replaced by an independent geometric random variable N. The conclusion is that the new estimators compare extremely favorably with previous ones. In particular, the conditional Monte Carlo estimator is the first heavy-tailed example of an estimator with bounded relative error. Further improvements are obtained in the random-N case, by incorporating control variates and stratification techniques into the new estimation procedures.

Isomeric distribution and catalyzed isomerization of cobalt(III) complexes with pentadentate macrocyclic ligands. Importance of hydrogen bonding

We have investigated the isomeric distribution and rearrangement of complexes of the type [CoXLn](2+,3+) (where X = Cl-, OH-, H2O, and L-n represents a pentadentate 13-, 14-, and 15-membered tetraaza or diaza-dithia (N-4 or N2S2) macrocycle bearing a pendant primary amine). The preparative procedures for chloro complexes produced almost exclusively kinetically preferred cis isomers (where the pendant primary amine is cis to the chloro ligand) that can be separated by careful cation-exchange chromatography. For L-13 and L-14 the so-called cis-V isomer is isolated as the kinetic product, and for L-15 the cis-VI form (an N-based diastereomer) is the preferred, while for the L-14(S) complex both cis-V and trans-I forms are obtained. All these complexes rearrange to form stable trans isomers in which the pendent primary amine is trans to the monodentate aqua or hydroxo ligand, depending on pH and the workup procedure. In total 11 different complexes have been studied. From these, two different trans isomers of [CoCIL14S](2+) have been characterized crystallographically for the first time in addition to a new structure of cis-V-[CoCIL14S](2+); all were isolated as their chloride perchlorate salts. Two additional isomers have been identified and characterized by NMR as reaction intermediates. The remaining seven forms correspond to the complexes already known, produced in preparative procedures. The kinetic, thermal, and baric activation parameters for all the isomerization reactions have been determined and involve large activation enthalpies and positive volumes of activation. Activation entropies indicate a very important degree of hydrogen bonding in the reactivity of the complexes, confirmed by density functional theory studies on the stability of the different isomeric forms. The isomerization processes are not simple and even some unstable intermediates have been detected and characterized as part of the above-mentioned 11 forms of the complexes. A common reaction mechanism for the isomerization reactions has been proposed for all the complexes derived from the observed kinetic and solution behavior.

Multi-level zero-inflated Poisson regression modelling of correlated count data with excess zeros

Count data with excess zeros relative to a Poisson distribution are common in many biomedical applications. A popular approach to the analysis of such data is to use a zero-inflated Poisson (ZIP) regression model. Often, because of the hierarchical Study design or the data collection procedure, zero-inflation and lack of independence may occur simultaneously, which tender the standard ZIP model inadequate. To account for the preponderance of zero counts and the inherent correlation of observations, a class of multi-level ZIP regression model with random effects is presented. Model fitting is facilitated using an expectation-maximization algorithm, whereas variance components are estimated via residual maximum likelihood estimating equations. A score test for zero-inflation is also presented. The multi-level ZIP model is then generalized to cope with a more complex correlation structure. Application to the analysis of correlated count data from a longitudinal infant feeding study illustrates the usefulness of the approach.

The role of habitat disturbance and recovery in metapopulation persistence

Classical metapopulation theory assumes a static landscape. However, empirical evidence indicates many metapopulations are driven by habitat succession and disturbance. We develop a stochastic metapopulation model, incorporating habitat disturbance and recovery, coupled with patch colonization and extinction, to investigate the effect of habitat dynamics on persistence. We discover that habitat dynamics play a fundamental role in metapopulation dynamics. The mean number of suitable habitat patches is not adequate for characterizing the dynamics of the metapopulation. For a fixed mean number of suitable patches, we discover that the details of how disturbance affects patches and how patches recover influences metapopulation dynamics in a fundamental way. Moreover, metapopulation persistence is dependent not only oil the average lifetime of a patch, but also on the variance in patch lifetime and the synchrony in patch dynamics that results from disturbance. Finally, there is an interaction between the habitat and metapopulation dynamics, for instance declining metapopulations react differently to habitat dynamics than expanding metapopulations. We close, emphasizing the importance of using performance measures appropriate to stochastic systems when evaluating their behavior, such as the probability distribution of the state of the. metapopulation, conditional on it being extant (i.e., the quasistationary distribution).

Toy model for a relational formulation of quantum theory

In the absence of an external frame of reference-i.e., in background independent theories such as general relativity-physical degrees of freedom must describe relations between systems. Using a simple model, we investigate how such a relational quantum theory naturally arises by promoting reference systems to the status of dynamical entities. Our goal is twofold. First, we demonstrate using elementary quantum theory how any quantum mechanical experiment admits a purely relational description at a fundamental. Second, we describe how the original non-relational theory approximately emerges from the fully relational theory when reference systems become semi-classical. Our technique is motivated by a Bayesian approach to quantum mechanics, and relies on the noiseless subsystem method of quantum information science used to protect quantum states against undesired noise. The relational theory naturally predicts a fundamental decoherence mechanism, so an arrow of time emerges from a time-symmetric theory. Moreover, our model circumvents the problem of the collapse of the wave packet as the probability interpretation is only ever applied to diagonal density operators. Finally, the physical states of the relational theory can be described in terms of spin networks introduced by Penrose as a combinatorial description of geometry, and widely studied in the loop formulation of quantum gravity. Thus, our simple bottom-up approach (starting from the semiclassical limit to derive the fully relational quantum theory) may offer interesting insights on the low energy limit of quantum gravity.

Air Entrainment Processes in a Full-Scale Rectangular Dropshaft at Large Flows

Rectangular dropshafts, commonly used in sewers and storm water systems, are characterised by significant flow aeration. New detailed air-water flow measurements were conducted in a near-full-scale dropshaft at large discharges. In the shaft pool and outflow channel, the results demonstrated the complexity of different competitive air entrainment mechanisms. Bubble size measurements showed a broad range of entrained bubble sizes. Analysis of streamwise distributions of bubbles suggested further some clustering process in the bubbly flow although, in the outflow channel, bubble chords were in average smaller than in the shaft pool. A robust hydrophone was tested to measure bubble acoustic spectra and to assess its field application potential. The acoustic results characterised accurately the order of magnitude of entrained bubble sizes, but the transformation from acoustic frequencies to bubble radii did not predict correctly the probability distribution functions of bubble sizes.

Liquid distribution as a means to describing the granule growth mechanism

Copper concentrate (chalcopyrite) was granulated in a rotating drum with a diameter of 0.3 m and a length of 0.2 m. Water was used as the binder and it was sprayed onto the powder bed with a nozzle. This material exhibited induction type behaviour, which was defined by Iveson and Litster [AIChE J. 44 (1998) 1510]. Induction type behaviour is characterized by the occurrence of an induction stage, during which the granules are gradually being compacted and little or no growth occurs. At the end of this induction stage, binder liquid is squeezed from the interior of the granules onto the granule surface and the granules are then surface-wet. This results in a rapid growth rate of the granules. Different types of experiments were conducted. The influence of the nozzle pressure and the distance from the nozzle to the powder bed on the growth behaviour of the granules as well as on the binder distribution was examined. The results of these experiments led to the postulation of a modified mechanism for induction type behaviour: it was found that after the binder was delivered, there were large granules containing a high amount of binder and small granules containing less binder. During the induction stage, the granules are compacted and binder liquid continuously appears at the surface of the large granules. These wet spots that are continuously being formed pick up the dry and small granules. When all the small granules have been picked up, further expulsion of binder liquid onto the granules' surface results in granules that remain surface-wet. This phenomenon marks the end of the induction stage and it coincides with the disappearance of the small granules. The hypothesis was tested by selectively removing the smaller granules during an experiment. As expected, this resulted in a shorter induction time.

Expertise and the perception of kinematic and situational probability information

Two experiments were conducted on the nature of expert perception in the sport of squash. In the first experiment, ten expert and fifteen novice players attempted to predict the direction and force of squash strokes from either a film display (occluded at variable time periods before and after the opposing player had struck the ball) or a matched point-light display (containing only the basic kinematic features of the opponent's movement pattern). Experts outperformed the novices under both display conditions, and the same basic time windows that characterised expert and novice pick-up of information in the film task also persisted in the point-light task. This suggests that the experts' perceptual advantage is directly related to their superior pick-up of essential kinematic information. In the second experiment, the vision of six expert and six less skilled players was occluded by remotely triggered liquid-crystal spectacles at quasi-random intervals during simulated match play. Players were required to complete their current stroke even when the display was occluded and their prediction performance was assessed with respect to whether they moved to the correct half of the court to match the direction and depth of the opponent's stroke. Consistent with experiment 1, experts were found to be superior in their advance pick-up of both directional and depth information when the display was occluded during the opponent's hitting action. However, experts also remained better than chance, and clearly superior to less skilled players, in their prediction performance under conditions where occlusion occurred before any significant pre-contact preparatory movement by the opposing player was visible. This additional source of expert superiority is attributable to their superior attunement to the information contained in the situational probabilities and sequential dependences within their opponent's pattern of play.

It's all very well, in theory: Theoretical perspectives and their applications in contemporary pedagogical research

Current debates about educational theory are concerned with the relationship between knowledge and power and thereby issues such as who possesses a truth and how have they arrived at it, what questions are important to ask, and how should they best be answered. As such, these debates revolve around questions of preferred, appropriate, and useful theoretical perspectives. This paper overviews the key theoretical perspectives that are currently used in physical education pedagogy research and considers how these inform the questions we ask and shapes the conduct of research. It also addresses what is contested with respect to these perspectives. The paper concludes with some cautions about allegiances to and use of theories in line with concerns for the applicability of educational research to pressing social issues.

Quantum noise in optical fibers. I. Stochastic equations

We analyze the quantum dynamics of radiation propagating in a single-mode optical fiber with dispersion, nonlinearity, and Raman coupling to thermal phonons. We start from a fundamental Hamiltonian that includes the principal known nonlinear effects and quantum-noise sources, including linear gain and loss. Both Markovian and frequency-dependent, non-Markovian reservoirs are treated. This treatment allows quantum Langevin equations, which have a classical form except for additional quantum-noise terms, to be calculated. In practical calculations, it is more useful to transform to Wigner or 1P quasi-probability operator representations. These transformations result in stochastic equations that can be analyzed by use of perturbation theory or exact numerical techniques. The results have applications to fiber-optics communications, networking, and sensor technology.