Population balance model of in vivo neutrophil formation following bone marrow rescue therapy

In this paper, we develop a simple four parameter population balance model of in vivo neutrophil formation following bone marrow rescue therapy. The model is used to predict the number and type of neutrophil progenitors required to abrogate the period of severe neutropenia that normally follows a bone marrow transplant. The estimated total number of 5 billion neutrophil progenitors is consistent with the value extrapolated from a human trial. The model provides a basis for designing ex vivo expansion protocols.

A wavelet visible difference predictor

In this paper, we describe a model of the human visual system (HVS) based on the wavelet transform. This model is largely based on a previously proposed model, but has a number of modifications that make it more amenable to potential integration into a wavelet based image compression scheme. These modifications include the use of a separable wavelet transform instead of the cortex transform, the application of a wavelet contrast sensitivity function (CSP), and a simplified definition of subband contrast that allows us to predict noise visibility directly from wavelet coefficients. Initially, we outline the luminance, frequency, and masking sensitivities of the HVS and discuss how these can be incorporated into the wavelet transform. We then outline a number of limitations of the wavelet transform as a model of the HVS, namely the lack of translational invariance and poor orientation sensitivity. In order to investigate the efficacy of this wavelet based model, a wavelet visible difference predictor (WVDP) is described. The WVDP is then used to predict visible differences between an original and compressed (or noisy) image. Results are presented to emphasize the limitations of commonly used measures of image quality and to demonstrate the performance of the WVDP, The paper concludes with suggestions on bow the WVDP can be used to determine a visually optimal quantization strategy for wavelet coefficients and produce a quantitative measure of image quality.

Phase diagram of a Heisenberg spin-Peierls model with quantum phonons

Using a new version of the density-matrix renormalization group we determine the phase diagram of a model of an antiferromagnetic Heisenberg spin chain where the spins interact with quantum phonons. A quantum phase transition from a gapless spin-fluid state to a gapped dimerized phase occurs at a nonzero value of the spin-phonon coupling. The transition is in the same universality class as that of a frustrated spin chain, to which the model maps in the diabatic limit. We argue that realistic modeling of known spin-Peierls materials should include the effects of quantum phonons.

Extended integrability regime for the supersymmetric U model

An extension of the supersymmetric U model for correlated electrons is given and integrability is established by demonstrating that the model can he constructed through the quantum inverse scattering method using an R-matrix without the difference property. Some general symmetry properties of the model are discussed and from the Bethe ansatz solution an expression for the energies is presented.

Integrable Kondo impurities in the one-dimensional supersymmetric extended Hubbard model

An integrable Kondo problem in the one-dimensional supersymmetric extended Hubbard model is studied by means of the boundary graded quantum inverse scattering method. The boundary K-matrices depending on the local moments of the impurities are presented as a non-trivial realization of the graded reflection equation algebras in a two-dimensional impurity Hilbert space. Further, the model is solved by using the algebraic Bethe ansatz method and the Bethe ansatz equations are obtained.

On modifications to the long-term survival mixture model in the presence of competing risks

A mixture model for long-term survivors has been adopted in various fields such as biostatistics and criminology where some individuals may never experience the type of failure under study. It is directly applicable in situations where the only information available from follow-up on individuals who will never experience this type of failure is in the form of censored observations. In this paper, we consider a modification to the model so that it still applies in the case where during the follow-up period it becomes known that an individual will never experience failure from the cause of interest. Unless a model allows for this additional information, a consistent survival analysis will not be obtained. A partial maximum likelihood (ML) approach is proposed that preserves the simplicity of the long-term survival mixture model and provides consistent estimators of the quantities of interest. Some simulation experiments are performed to assess the efficiency of the partial ML approach relative to the full ML approach for survival in the presence of competing risks.

Using Lindenmayer systems to model morphogenesis in a tropical pasture legume Stylosanthes scabra

This paper summarizes the processes involved in designing a mathematical model of a growing pasture plant, Stylosanthes scabra Vog. cv. Fitzroy. The model is based on the mathematical formalism of Lindenmayer systems and yields realistic computer-generated images of progressive plant geometry through time. The processes involved in attaining growth data, retrieving useful growth rules, and constructing a virtual plant model are outlined. Progressive output morphological data proved useful for predicting total leaf area and allowed for easier quantification of plant canopy size in terms of biomass and total leaf area.

New integrable model of correlated electrons with off-diagonal long-range order from so(5) symmetry

We present a new integrable model for correlated electrons which is based on so(5) symmetry. By using an eta-pairing realization we construct eigenstates of the Hamiltonian with off-diagonal long-range order. It is also shown that these states lie in the ground state sector. We exactly solve the model on a one-dimensional lattice by the Bethe ansatz.

Protein sequence threading, the alignment problem, and a two-step strategy

Conventionally, protein structure prediction via threading relies on some nonoptimal method to align a protein sequence to each member of a library of known structures. We show how a score function (force field) can be modified so as to allow the direct application of a dynamic programming algorithm to the problem. This involves an approximation whose damage can be minimized by an optimization process during score function parameter determination. The method is compared to sequence to structure alignments using a more conventional pair-wise score function and the frozen approximation. The new method produces results comparable to the frozen approximation, but is faster and has fewer adjustable parameters. It is also free of memory of the template's original amino acid sequence, and does not suffer from a problem of nonconvergence, which can be shown to occur with the frozen approximation. Alignments generated by the simplified score function can then be ranked using a second score function with the approximations removed. (C) 1999 John Wiley & Sons, Inc.

Ameloblast apoptosis and IGF-1 receptor expression in the continuously erupting rat incisor model

Enamel-producing cells (ameloblasts) pass through several phenotypic and functional stages during enamel formation. In the transition between secretory and maturation stages, about one quarter of the ameloblasts suddenly undergo apoptosis. We have studied this phenomenon using the continuously erupting rat incisor model. A special feature of this model is that all stages of ameloblast differentiation are presented within a single longitudinal section of the developing tooth. This permits investigation of the temporal sequence of gene and growth factor receptor expression during ameloblast differentiation and apoptosis. We describe the light and electron microscopic morphology of ameloblast apoptosis and the pattern of insulin-like growth factor-1 receptor expression by ameloblasts in the continuously erupting rat incisor model. In the developing rat incisor, ameloblast apoptosis is associated with downregulated expression of the insulin-like growth factor-1 receptor. These data are consistent with the hypothesis that ameloblasts are hard wired for apoptosis and that insulin-like growth factor-1 receptor expression is required to block the default apoptotic pathway. Possible mechanisms of insulin-like growth factor-1 inhibition of ameloblast apoptosis are presented. The rat incisor model may be useful in studies of physiological apoptosis as it presents apoptosis in a predictable pattern in adult tissues.

RWMODEL II: Computer simulation of the Rescorla-Wagner model of Pavlovian conditioning

RWMODEL II simulates the Rescorla-Wagner model of Pavlovian conditioning. It is written in Delphi and runs under Windows 3.1 and Windows 95. The program was designed for novice and expert users and can be employed in teaching, as well as in research. It is user friendly and requires a minimal level of computer literacy but is sufficiently flexible to permit a wide range of simulations. It allows the display of empirical data, against which predictions from the model can be validated.

Class attainment among British men - A multivariate extension of the CASMIN model of intergenerational class mobility

The CASMIN Project is arguably the most influential contemporary study of class mobility in the world. However, CASMIN results with respect to weak vertical status effects on class mobility have been extensively criticized. Drawing on arguments about how to model vertical mobility, Hout and Hauser (1992) show that class mobility is strongly determined by vertical socioeconomic differences. This paper extends these arguments by estimating the CASMIN model while explicitly controlling for individual determinants of socioeconomic attainment. Using the 1972 Oxford Mobility Data and the 1979 and 1983 British Election Studies, the paper employs mixed legit models to show how individual socioeconomic factors and categorical differences between classes shape intergenerational mobility. The findings highlight the multidimensionality of class mobility and its irreducibility to vertical movement up and down a stratification hierarchy.

Determining forest structural attributes using an inverted geometric-optical model in mixed eucalypt forests, Southeast Queensland, Australia

The Montreal Process indicators are intended to provide a common framework for assessing and reviewing progress toward sustainable forest management. The potential of a combined geometrical-optical/spectral mixture analysis model was assessed for mapping the Montreal Process age class and successional age indicators at a regional scale using Landsat Thematic data. The project location is an area of eucalyptus forest in Emu Creek State Forest, Southeast Queensland, Australia. A quantitative model relating the spectral reflectance of a forest to the illumination geometry, slope, and aspect of the terrain surface and the size, shape, and density, and canopy size. Inversion of this model necessitated the use of spectral mixture analysis to recover subpixel information on the fractional extent of ground scene elements (such as sunlit canopy, shaded canopy, sunlit background, and shaded background). Results obtained fron a sensitivity analysis allowed improved allocation of resources to maximize the predictive accuracy of the model. It was found that modeled estimates of crown cover projection, canopy size, and tree densities had significant agreement with field and air photo-interpreted estimates. However, the accuracy of the successional stage classification was limited. The results obtained highlight the potential for future integration of high and moderate spatial resolution-imaging sensors for monitoring forest structure and condition. (C) Elsevier Science Inc., 2000.

Marine protected areas for spatially structured exploited stocks

Many harvested marine and terrestrial populations have segments of their range protected in areas free from exploitation. Reasons for areas being protected from harvesting include conservation, tourism, research, protection of breeding grounds, stock recovery, harvest regulation, or habitat that is uneconomical to exploit. In this paper we consider the problem of optimally exploiting a single species local population that is connected by dispersing larvae to an unharvested local population. We define a spatially-explicit population dynamics model and apply dynamic optimization techniques to determine policies for harvesting the exploited patch. We then consider how reservation affects yield and spawning stock abundance when compared to policies that have not recognised the spatial structure of the metapopulation. Comparisons of harvest strategies between an exploited metapopulation with and without a harvest refuge are also made. Results show that in a 2 local population metapopulation with unidirectional larval transfer, the optimal exploitation of the harvested population should be conducted as if it were independent of the reserved population. Numerical examples suggest that relative source populations should be exploited if the objective is to maximise spawning stock abundance within a harvested metapopulation that includes a protected local population. However, this strategy can markedly reduce yield over a sink harvested reserve system and may require strict regulation for conservation goals to be realised. If exchange rates are high, results indicate that spawning stock abundance can be less in a reserve system than in a fully exploited metapopulation. In order to maximise economic gain in the reserve system, results indicate that relative sink populations should be harvested. Depending on transfer levels, loss in harvest through reservation can be minimal, and is likely to be compensated by the potential environmental and economic benefits of the reserve.