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.

Expression of galectin-1 in the olfactory nerve pathway of rat

The olfactory neuroepithelium is characterised by the mosaic distribution of primary olfactory neurons that express different odorant receptors and cell surface glycoconjugates. Carbohydrates are believed to form a glycocode that mediates sorting out and fasciculation of primary olfactory axons through interactions with carbohydrate-binding proteins such as galectin-1. In the present study, we describe in detail the expression pattern of galectin-1 in the developing and adult rat olfactory system. We demonstrate that galectin-1 is expressed by olfactory ensheathing cells both in olfactory nerve and within the nerve fibre layer of the olfactory bulb of the embryonic and adult rat. In the adult rat, galectin-1 was preferentially expressed by olfactory ensheathing cells in the nerve fibre layer of the ventromedial and lateral surfaces of the olfactory bulb. Galectin-1 was also expressed by subsets of periglomerular cells and granule cells, particularly in the ventromedial region of the olfactory bulb. In adult rat, the galectin-1 ligand, N-acetyl-lactosamine, was expressed by primary olfactory axons that terminated in glomeruli present in the ventromedial and lateral olfactory bulb. These results suggest that expression of galectin-1 may provide a mechanism for the sorting of subpopulations of axons in the nerve fibre layer of the olfactory bulb during development as well as play a role in the postnatal maintenance of specific glomerular connections. (C) 1999 Elsevier Science B.V. All rights reserved.

The mechanical properties of Saccharomyces cerevisiae

Cell-wall mechanical properties play an integral part in the growth and form of Saccharomyces cerevisiae, In contrast to the tremendous knowledge on the genetics of S. cerevisiae, almost nothing is known about its mechanical properties. We have developed a micromanipulation technique to measure the force required to burst single cells and have recently established a mathematical model to extract the mechanical properties of the cell wall from such data, Here we determine the average surface modulus of the S, cerevisiae cell wall to be 11.1 +/- 0.6 N/m and 12.9 +/- 0.7 N/m in exponential and stationary phases, respectively, giving corresponding Young's moduli of 112 +/- 6 MPa and 107 +/- 6 MPa, This result demonstrates that yeast cell populations strengthen as they enter stationary phase by increasing wall thickness and hence the surface modulus, without altering the average elastic properties of the cell-wall material. We also determined the average breaking strain of the cell wall to be 82% +/- 3% in exponential phase and 80% +/- 3% in stationary phase, This finding provides a failure criterion that can be used to predict when applied stresses (e,g,, because of fluid flow) will lead to wall rupture, This work analyzes yeast compression experiments in different growth phases by using engineering methodology.

Determination of activation energy distributions for chemisorption of oxygen on carbon: an improved approach

The present paper proposes an approach to obtaining the activation energy distribution for chemisorption of oxygen onto carbon surfaces, while simultaneously allowing for the activation energy dependence of the pre-exponential factor of the rate constant. Prior studies in this area have considered this factor to be uniform, thereby biasing estimated distributions. The results show that the derived activation energy distribution is not sensitive to the chemisorption mechanism because of the step function like property of the coverage. The activation energy distribution is essentially uniform for some carbons, and has two or possibly more discrete stages, suggestive of at least two types of sites, each with its own uniform distribution. The pre-exponential factors of the reactions are determined directly from the experimental data, and are found not to be constant as assumed in earlier work, but correlated with the activation energy. The latter results empirically follow an exponential function, supporting some earlier statistical and experimental work. The activation energy distribution obtained in the present paper permits improved correlation of chemisorption data in comparison to earlier studies. (C) 2000 Elsevier Science Ltd. All rights reserved.

Computational immunology: The coming of age

The explosive growth in biotechnology combined with major advancesin information technology has the potential to radically transformimmunology in the postgenomics era. Not only do we now have readyaccess to vast quantities of existing data, but new data with relevanceto immunology are being accumulated at an exponential rate. Resourcesfor computational immunology include biological databases and methodsfor data extraction, comparison, analysis and interpretation. Publiclyaccessible biological databases of relevance to immunologists numberin the hundreds and are growing daily. The ability to efficientlyextract and analyse information from these databases is vital forefficient immunology research. Most importantly, a new generationof computational immunology tools enables modelling of peptide transportby the transporter associated with antigen processing (TAP), modellingof antibody binding sites, identification of allergenic motifs andmodelling of T-cell receptor serial triggering.

Molecular evidence for definition of genera in the Oxylobium group (Fabaceae : Mirbelieae)

The circumscription of Oxylobium and related genera has been problematic for nearly 200 years. Traditional definitions of genera in the group have relied on morphological features of the leaves, flower, and fruit that overlap extensively between genera. Therefore sequences of cpDNA (trnL-F intron and spacer) and nrDNA (ITS) were used to estimate the phylogeny of the group in an attempt to redefine the genera as monophyletic groups. Oxylobium sens. str. was found to be a well supported clade in both data sets, with the inclusion of Mirbelia oxylobioides. No other genus in the group was supported by these data, except Gastrolobium sens. lat. Some species groups within Chorizema, Mirbelia, and Podolobium were supported but relationships among these, Oxylobium and Gastrolobium differed significantly between the chloroplast and nuclear data sets. No group supported by the molecular data had a morphological synapomorphy, not even Oxylobium or Gastrolobium. Therefore it may be necessary to adopt a much broader generic concept in this group than has been done previously. Incongruence between the two molecular data sets, and very short internal basal branches in both, suggest a rapid early radiation in this group, possibly combined with hybridization and lineage sorting.

Infiltration from surface and buried point sources: The average wetting water content

The assumption in analytical solutions for flow from surface and buried point sources of an average water content, (θ) over bar, behind the wetting front is examined. Some recent work has shown that this assumption fitted some field data well. Here we calculated (θ) over bar using a steady state solution based on the work by Raats [1971] and an exponential dependence of the diffusivity upon the water content. This is compared with a constant value of (θ) over bar calculated from an assumption of a hydraulic conductivity at the wetting front of 1 mm day(-1) and the water content at saturation. This comparison was made for a wide range of soils. The constant (θ) over bar generally underestimated (θ) over bar at small wetted radii and overestimated (θ) over bar at large radii. The crossover point between under and overestimation changed with both soil properties and flow rate. The largest variance occurred for coarser texture soils at low-flow rates. At high-flow rates in finer-textured soils the use of a constant (θ) over bar results in underestimation of the time for the wetting front to reach a particular radius. The value of (θ) over bar is related to the time at which the wetting front reaches a given radius. In coarse-textured soils the use of a constant value of (θ) over bar can result in an error of the time when the wetting front reaches a particular radius, as large as 80% at low-flow rates and large radii.

Use of bulk density profiles from X-radiography to examine structural crust models

Exponential and sigmoidal functions have been suggested to describe the bulk density profiles of crusts. The present work aims to evaluate these conceptual models using high resolution X-radiography. Repacked seedbeds from two soil materials, air-dried or prewetted by capillary rise, were subjected to simulated rain, which resulted in three types of structural crusts, namely, slaking, infilling, and coalescing. Bulk density distributions with depth were generated using high-resolution (70 mum), calibrated X-ray images of slices from the resin-impregnated crusted seedbeds. The bulk density decreased progressively with depth, which supports the suggestion that a crust should be considered as a nonuniform layer. For the slaking and the coalescing crusts, the exponential function underestimated the strong change in bulk density across the morphologically defined transition between the crust and the underlying material; the sigmoidal function provided a better description. Neither of these crust models effectively described the shape of the bulk density profiles through the whole seedbed. Below the infilling and slaking crusts, bulk density increased linearly with depth as a result of slumping. In the coalescing crusted seedbed, the whole seedbed uniformly collapsed and most of the bulk density change within the crust could be ascribed to slumping (0.33 g cm(-3)) rather than to crusting (0.12 g cm(-3)). Finally, (i) X-radiography appears as a unique tool to generate high resolution bulk density profiles and (ii) in structural crusts, bulk density profiles could be modeled using the existing exponential and sigmoidal crusting models, provided a slumping model would be coupled.

Cumulative sedimentation analysis of Escherichia coli debris size

A new method to measure Escherichia coil cell debris size after homogenization is presented. It is based on cumulative sedimentation analysis under centrifugal force, coupled with Sodium Dodecyl Sulfate-Polyacrylamide Gel Electrophoresis (SDS-PAGE) analysis of sedimented proteins. The effects that fermentation and homogenization conditions have on the resulting debris distributions were investigated using this method. Median debris size decreased significantly from approximately 0.5 mu m to 0.3 mu m as the number of homogenization passes increased from 2 to 10. Under identical homogenization conditions, uninduced host cells in stationary phase had a larger debris size than exponential cells after 5 homogenizer passes. This difference was not evident after 2 or in passes, possibly because of confounding intact cells and the existence of a minimum debris size for the conditions investigated. Recombinant cells containing protein inclusion bodies had the smallest debris size following homogenization. The method was also used to measure the size distribution of inclusion bodies. This result compared extremely well with an independent determination using centrifugal disc photosedimentation (CDS), thus validating the method. This is the first method that provides accurate size distributions of E. coli debris without the need for sample pretreatment, theoretical approximations (e.g. extinction coefficients), or the separation of debris and inclusion bodies prior to analysis. (C) 1997 John Wiley & Sons, Inc.

The subcycled Newmark algorithm

The popular Newmark algorithm, used for implicit direct integration of structural dynamics, is extended by means of a nodal partition to permit use of different timesteps in different regions of a structural model. The algorithm developed has as a special case an explicit-explicit subcycling algorithm previously reported by Belytschko, Yen and Mullen. That algorithm has been shown, in the absence of damping or other energy dissipation, to exhibit instability over narrow timestep ranges that become narrower as the number of degrees of freedom increases, making them unlikely to be encountered in practice. The present algorithm avoids such instabilities in the case of a one to two timestep ratio (two subcycles), achieving unconditional stability in an exponential sense for a linear problem. However, with three or more subcycles, the trapezoidal rule exhibits stability that becomes conditional, falling towards that of the central difference method as the number of subcycles increases. Instabilities over narrow timestep ranges, that become narrower as the model size increases, also appear with three or more subcycles. However by moving the partition between timesteps one row of elements into the region suitable for integration with the larger timestep these the unstable timestep ranges become extremely narrow, even in simple systems with a few degrees of freedom. As well, accuracy is improved. Use of a version of the Newmark algorithm that dissipates high frequencies minimises or eliminates these narrow bands of instability. Viscous damping is also shown to remove these instabilities, at the expense of having more effect on the low frequency response.

Performance on tests sensitive to impaired executive ability in schizophrenia, mania and well controls: acute and subacute phases

Aims: To compare the performance of schizophrenia, mania and well control groups on tests sensitive to impaired executive ability, and to assess the within-group stability of these measures across the acute and subacute phases of psychoses. Method: Recently admitted patients with schizophrenia (n=36), mania (n=18) and a well control group (n=20) were assessed on two occasions separated by 4 weeks. Tests included: the Controlled Oral Word Association Test, the Stroop Test, the Wisconsin Card Sort Test, and the Trail Making Test. Results: The two patient groups were significantly impaired on the Stroop Test at both time points compared to the control group. Significant group differences were also found for the Trail Making Test at Time 1 and for the Wisconsin Card Sort Test at Time 2. When controlled for practice effect, significant improvements over time were found on the Stroop and Trail Making tests in the schizophrenia group and on WCST Categories Achieved in the mania group. Discussion: Compared to controls, the patient groups were impaired on measures related to executive ability. The pattern of improvement on test scores between the acute and subacute phases differed between patients with schizophrenia versus patients with mania. (C) 1997 Elsevier Science B.V.

Metastable chaos in the ammonia ring laser

We report experimental studies of metastable chaos in the far-infrared ammonia ring: laser. When the laser pump power is switched from above chaos threshold to slightly below, chaotic intensity pulsations continue for a varying time afterward before decaying to either periodic or cw emission. The behavior is in good qualitative agreement with that predicted by the Lorenz equations, previously used to describe this laser. The statistical distribution of the duration of the chaotic transient is measured and shown to be in excellent agreement with the Lorenz equations in showing a modified exponential distribution. We also give a brief numerical analysis and graphical visualization of the Lorenz equations in phase space illustrating the boundary between the metastable chaotic and the stable fixed point basins of attraction. This provides an intuitive understanding of the metastable dynamics of the Lorenz equations and the experimental system.

Substrate-bound carbohydrates stimulate signal transduction and neurite outgrowth in an olfactory neuron cell line

SUBPOPULATIONS of olfactory receptor neurons, which are dispersed throughout the olfactory neuroepithelium, express specific cell surface carbohydrates and project to discrete regions of the olfactory bulb. Cell surface carbohydrates such as N-acetyl-lactosamine have been postulated to mediate sorting and selective fasciculation of discrete axon subpopulations during development of the olfactory pathway. Substrate-bound N-acetyl-lactosamine promotes neurite outgrowth by both clonal olfactory receptor neuron cell lines and olfactory receptor neurons in vitro, indicating that cell surface carbohydrates may be ligands for receptor-mediated stimulation of axon growth in vivo. In the present study, the role of transmembrane signaling in N-acetyl-lactosamine-stimulated neurite outgrowth was examined in the clonal olfactory neuron cell line 4.4.2. Substrate-bound N-acetyl-lactosamine stimulated neurite outgrowth which was specifically inhibited by antagonists to N- and L-type calcium channels and to tyrosine kinase phosphorylation. These results indicate that N-acetyl-lactosamine can evoke transmembrane receptor-mediated responses capable of influencing neurite outgrowth.

Modeling of subcutaneous absorption kinetics of infusion solutions in the elderly using technetium

Absorption kinetics of solutes given with the subcutaneous administration of fluids is ill-defined. The gamma emitter, technitium pertechnetate, enabled estimates of absorption rate to be estimated independently using two approaches. In the first approach, the counts remaining at the site were estimated by imaging above the subcutaneous administration site, whereas in the second approach, the plasma technetium concentration-time profiles were monitored up to 8 hr after technetium administration. Boluses of technetium pertechnetate were given both intravenously and subcutaneously on separate occasions with a multiple dosing regimen using three doses on each occasion. The disposition of technetium after iv administration was best described by biexponential kinetics with a V-ss of 0.30 +/- 0.11 L/kg and a clearance of 30.0 +/- 13.1 ml/min. The subcutaneous absorption kinetics was best described as a single exponential process with a half-life of 18.16 +/- 3.97 min by image analysis and a half-life of 11.58 +/- 2.48 min using plasma technetium time data. The bioavailability of technetium by the subcutaneous route was estimated to be 0.96 +/- 0.12. The absorption half-life showed no consistent change with the duration of the subcutaneous infusion. The amount remaining at the absorption site with time was similar when analyzed using image analysis, and plasma concentrations assuming multiexponential disposition kinetics and a first-order absorption process. Profiles of fraction remaining at the absorption sire generated by deconvolution analysis, image analysis, and assumption of a constant first-order absorption process were similar. Slowing of absorption from the subcutaneous administration site is apparent after the last bolus dose in three of the subjects and can De associated with the stopping of the infusion. In a fourth subject, the retention of technetium at the subcutaneous site is more consistent with accumulation of technetium near the absorption site as a result of systemic recirculation.

On the decay rates of buffers in continuous flow lines

Consider a tandem system of machines separated by infinitely large buffers. The machines process a continuous flow of products, possibly at different speeds. The life and repair times of the machines are assumed to be exponential. We claim that the overflow probability of each buffer has an exponential decay, and provide an algorithm to determine the exact decay rates in terms of the speeds and the failure and repair rates of the machines. These decay rates provide useful qualitative insight into the behavior of the flow line. In the derivation of the algorithm we use the theory of Large Deviations.