Time of treatment influences the appearance of drug-resistant parasites in Plasmodium falciparum infections

A deterministic mathematical model which predicts the probability of developing a new drug-resistant parasite population within the human host is reported, The model incorporates the host's specific antibody response to PfEMP1, and also investigates the influence of chemotherapy on the probability of developing a viable drug-resistant parasite population within the host. Results indicate that early, treatment, and a high antibody threshold coupled with a long lag time between antibody stimulation and activity, are risk factors which increase the likelihood of developing a viable drug-resistant parasite population. High parasite mutation rates and fast PfEMP1 var gene switching are also identified as risk factors. The model output allows the relative importance of the various risk factors as well as the relationships between them to be established, thereby increasing the understanding of the conditions which favour the development of a new drug-resistant parasite population.

Primary fallopian tube carcinoma: the Queensland experience

The purpose of this study was to review the experience with fallopian tube carcinoma in Queensland and to compare it with previously published data. Thirty-six patients with primary fallopian tube carcinoma treated at the Queensland Gynaecological Cancer Center from 1988 to 1999 were reviewed in a retrospective clinicopathologic study. All patients had primary surgery and 31/36 received chemotherapy postoperatively. Abnormal vaginal bleeding (15/36) and abdominal pain (14/36) were the most common presenting symptoms at the time of diagnosis. Median follow-up was 70.3 months and the median overall survival was 68.1 months. Surgical stage I disease (P = 0.02) and the absence of residual tumor after operation (P = 0.03) were the only factors associated with improved survival. Twenty of the 36 patients (55%) presented with stage I disease and survival was 62.7% at 5 years. No patient with postoperative residual tumor survived. The majority of the patients with fallopian tube carcinoma present with stage I disease at diagnosis, but their survival probability is low compared with that of other early stage gynecological malignancies. If primary surgical debulking cannot achieve macroscopic tumor clearence, the chance of survival is extremely low.

A risk-based system to penalize and reward line management for occupational safety and health performance

Penalizing line management for the occurrence of lost time injuries has in some cases had unintended negative consequences. These are discussed. An alternative system is suggested that penalizes line management for accidents where the combination of the probability of recurrence and the maximum reasonable consequences such a recurrence may have exceeds an agreed limit. A reward is given for prompt effective control of the risk to below the agreed risk limit. The reward is smaller than the penalty. High-risk accidents require independent investigation by a safety officer using analytical techniques. Two case examples are given to illustrate the system. Continuous safety improvement is driven by a planned reduction in the agreed risk limit over time and reward for proactive risk assessment and control.

The design of hazard risk assessment matrices for ranking occupational health risks and their application in mining and minerals processing

Two hazard risk assessment matrices for the ranking of occupational health risks are described. The qualitative matrix uses qualitative measures of probability and consequence to determine risk assessment codes for hazard-disease combinations. A walk-through survey of an underground metalliferous mine and concentrator is used to demonstrate how the qualitative matrix can be applied to determine priorities for the control of occupational health hazards. The semi-quantitative matrix uses attributable risk as a quantitative measure of probability and uses qualitative measures of consequence. A practical application of this matrix is the determination of occupational health priorities using existing epidemiological studies. Calculated attributable risks from epidemiological studies of hazard-disease combinations in mining and minerals processing are used as examples. These historic response data do not reflect the risks associated with current exposures. A method using current exposure data, known exposure-response relationships and the semi-quantitative matrix is proposed for more accurate and current risk rankings.

Brief note on the variation of constants formula for fuzzy differential equations

This note gives a theory of state transition matrices for linear systems of fuzzy differential equations. This is used to give a fuzzy version of the classical variation of constants formula. A simple example of a time-independent control system is used to illustrate the methods. While similar problems to the crisp case arise for time-dependent systems, in time-independent cases the calculations are elementary solutions of eigenvalue-eigenvector problems. In particular, for nonnegative or nonpositive matrices, the problems at each level set, can easily be solved in MATLAB to give the level sets of the fuzzy solution. (C) 2002 Elsevier Science B.V. All rights reserved.

Patchy populations in stochastic environments: Critical number of patches for persistence

We introduce a model for the dynamics of a patchy population in a stochastic environment and derive a criterion for its persistence. This criterion is based on the geometric mean (GM) through time of the spatial-arithmetic mean of growth rates. For the population to persist, the GM has to be greater than or equal to1. The GM increases with the number of patches (because the sampling error is reduced) and decreases with both the variance and the spatial covariance of growth rates. We derive analytical expressions for the minimum number of patches (and the maximum harvesting rate) required for the persistence of the population. As the magnitude of environmental fluctuations increases, the number of patches required for persistence increases, and the fraction of individuals that can be harvested decreases. The novelty of our approach is that we focus on Malthusian local population dynamics with high dispersal and strong environmental variability from year to year. Unlike previous models of patchy populations that assume an infinite number of patches, we focus specifically on the effect that the number of patches has on population persistence. Our work is therefore directly relevant to patchily distributed organisms that are restricted to a small number of habitat patches.

A zero-augmented gamma mixed model for longitudinal data with many zeros

In many occupational safety interventions, the objective is to reduce the injury incidence as well as the mean claims cost once injury has occurred. The claims cost data within a period typically contain a large proportion of zero observations (no claim). The distribution thus comprises a point mass at 0 mixed with a non-degenerate parametric component. Essentially, the likelihood function can be factorized into two orthogonal components. These two components relate respectively to the effect of covariates on the incidence of claims and the magnitude of claims, given that claims are made. Furthermore, the longitudinal nature of the intervention inherently imposes some correlation among the observations. This paper introduces a zero-augmented gamma random effects model for analysing longitudinal data with many zeros. Adopting the generalized linear mixed model (GLMM) approach reduces the original problem to the fitting of two independent GLMMs. The method is applied to evaluate the effectiveness of a workplace risk assessment teams program, trialled within the cleaning services of a Western Australian public hospital.

A mixture model-based approach to the clustering of microarray expression data

Motivation: This paper introduces the software EMMIX-GENE that has been developed for the specific purpose of a model-based approach to the clustering of microarray expression data, in particular, of tissue samples on a very large number of genes. The latter is a nonstandard problem in parametric cluster analysis because the dimension of the feature space (the number of genes) is typically much greater than the number of tissues. A feasible approach is provided by first selecting a subset of the genes relevant for the clustering of the tissue samples by fitting mixtures of t distributions to rank the genes in order of increasing size of the likelihood ratio statistic for the test of one versus two components in the mixture model. The imposition of a threshold on the likelihood ratio statistic used in conjunction with a threshold on the size of a cluster allows the selection of a relevant set of genes. However, even this reduced set of genes will usually be too large for a normal mixture model to be fitted directly to the tissues, and so the use of mixtures of factor analyzers is exploited to reduce effectively the dimension of the feature space of genes. Results: The usefulness of the EMMIX-GENE approach for the clustering of tissue samples is demonstrated on two well-known data sets on colon and leukaemia tissues. For both data sets, relevant subsets of the genes are able to be selected that reveal interesting clusterings of the tissues that are either consistent with the external classification of the tissues or with background and biological knowledge of these sets.

A simulated annealing algorithm for finding consensus sequences

Motivation: A consensus sequence for a family of related sequences is, as the name suggests, a sequence that captures the features common to most members of the family. Consensus sequences are important in various DNA sequencing applications and are a convenient way to characterize a family of molecules. Results: This paper describes a new algorithm for finding a consensus sequence, using the popular optimization method known as simulated annealing. Unlike the conventional approach of finding a consensus sequence by first forming a multiple sequence alignment, this algorithm searches for a sequence that minimises the sum of pairwise distances to each of the input sequences. The resulting consensus sequence can then be used to induce a multiple sequence alignment. The time required by the algorithm scales linearly with the number of input sequences and quadratically with the length of the consensus sequence. We present results demonstrating the high quality of the consensus sequences and alignments produced by the new algorithm. For comparison, we also present similar results obtained using ClustalW. The new algorithm outperforms ClustalW in many cases.

Efficient simulation of a Tandem Jackson Network

The two-node tandem Jackson network serves as a convenient reference model for the analysis and testing of different methodologies and techniques in rare event simulation. In this paper we consider a new approach to efficiently estimate the probability that the content of the second buffer exceeds some high level L before it becomes empty, starting from a given state. The approach is based on a Markov additive process representation of the buffer processes, leading to an exponential change of measure to be used in an importance sampling procedure. Unlike changes of measures proposed and studied in recent literature, the one derived here is a function of the content of the first buffer. We prove that when the first buffer is finite, this method yields asymptotically efficient simulation for any set of arrival and service rates. In fact, the relative error is bounded independent of the level L; a new result which is not established for any other known method. When the first buffer is infinite, we propose a natural extension of the exponential change of measure for the finite buffer case. In this case, the relative error is shown to be bounded (independent of L) only when the second server is the bottleneck; a result which is known to hold for some other methods derived through large deviations analysis. When the first server is the bottleneck, experimental results using our method seem to suggest that the relative error is bounded linearly in L.

Reconciling paleodistribution models and comparative phylogeography in the Wet Tropics rainforest land snail Gnarosophia bellendenkerensis (Brazier 1875)

Comparative phylogeography has proved useful for investigating biological responses to past climate change and is strongest when combined with extrinsic hypotheses derived from the fossil record or geology. However, the rarity of species with sufficient, spatially explicit fossil evidence restricts the application of this method. Here, we develop an alternative approach in which spatial models of predicted species distributions under serial paleoclimates are compared with a molecular phylogeography, in this case for a snail endemic to the rainforests of North Queensland, Australia. We also compare the phylogeography of the snail to those from several endemic vertebrates and use consilience across all of these approaches to enhance biogeographical inference for this rainforest fauna. The snail mtDNA phylogeography is consistent with predictions from paleoclimate modeling in relation to the location and size of climatic refugia through the late Pleistocene-Holocene and broad patterns of extinction and recolonization. There is general agreement between quantitative estimates of population expansion from sequence data (using likelihood and coalescent methods) vs. distributional modeling. The snail phylogeography represents a composite of both common and idiosyncratic patterns seen among vertebrates, reflecting the geographically finer scale of persistence and subdivision in the snail. In general, this multifaceted approach, combining spatially explicit paleoclimatological models and comparative phylogeography, provides a powerful approach to locating historical refugia and understanding species' responses to them.

Laying date and laying sequence influence the sex ratio of crimson rosella broods

We examine the patterns of sex allocation in crimson rosellas Platycercus elegans, a socially monogamous Australian parrot. Overall, 41.8% of nestlings were male, a significant female bias. However underlying this population-level bias were non-random patterns of sex allocation within broods. Broods produced early in the season were female-biased, but the proportion of males in a brood increased as the breeding season progressed. Female rosellas may obtain greater fitness benefits from early-fledging daughters than sons because daughters can breed as 1-year-olds whereas sons do not breed until they are at least 2 years old. Laying date and laying sequence also interacted to influence the sex ratio of eggs. The sex of early-laid eggs strongly followed the brood level pattern, whereas the sex of middle- and late-laid eggs did not change significantly as the season progressed. Nevertheless, late-laid eggs were very unlikely to be male at the end of the season. We argue these differing seasonal patterns reflect the relative costs and benefits to producing early-hatched males and females at different times of the season. Female rosellas appear to maximise the probability that daughters are able to breed early but to minimise competitive asymmetries within the brood. In particular, late-hatched male chicks are disadvantaged if their oldest sibling is male, explaining the dearth of broods containing late-hatched males at the end of the breeding season.

Acute symptoms, not rectally administered sucralfate, predict for late radiation proctitis: Longer term follow-up of a phase III trial - trans-Tasman radiation oncology group

Image : To assess the potential for sucralfate administered rectally to reduce the risk of late rectal morbidity in patients undergoing nonconformal radiotherapy (RT) for carcinoma of the prostate and to study the variables potentially contributing to late rectal morbidity and particularly to explore the relationship between acute and late toxicity. Image : Eighty-six patients with localized prostate carcinoma were randomized in a double-blind, placebo-controlled study to a daily enema of 3 g of sucralfate in a 15-mL suspension or the same suspension without sucralfate. The enema began the first day of RT and was continued for 2 weeks after treatment completion. The primary end point of the study was acute Radiation Therapy Oncology Group (RTOG)/European Organization for Research and Treatment of Cancer (EORTC) toxicity; however, the patients were followed for an additional 5 years on a 6-month basis. The evaluation included late RTOG/EORTC toxicity and a patient self-assessment questionnaire. Image : With a median follow-up of 5 years, the Kaplan-Meier probability of late Grade 2 RTOG/EORTC toxicity was 12% (95% confidence interval [CI] 2–22%) for placebo and 5% (95% CI 0–12%) for sucralfate (p = 0.26). The probability of late rectal bleeding was 59% (95% CI 45–73%) for placebo and 54% (95% CI 40–68%) for sucralfate. No statistically significant difference was found between the treatment arms for the peak incidence of any of the other patient self-assessment variables. Cox proportional hazards modeling indicated acute RTOG/EORTC toxicity of Grade 2 or greater was associated with a hazard ratio of 2.74 (95% CI 1.31–5.73) for the development of late toxicity of Grade 1 or greater. Substituting the patient self-assessment variables for acute RTOG/EORTC toxicity revealed that rectal pain of a moderate or severe grade during RT was the best predictor of the subsequent development of late toxicity, with a hazard ratio of 3.44 (95% CI 1.68–7). Image : The results of this study do not support the use of sucralfate administered rectally as a method for reducing the late toxicity of nonconformal RT for prostate cancer. There appears to be an association between the development of acute and subsequent late toxicity, although the nature of this association remains to be determined

Two-dimensional balanced sampling plans excluding contiguous units

A balanced sampling plan excluding contiguous units (or BSEC for short) was first introduced by Hedayat, Rao and Stufken in 1988. These designs can be used for survey sampling when the units are arranged in one-dimensional ordering and the contiguous units in this ordering provide similar information. In this paper, we generalize the concept of a BSEC to the two-dimensional situation and give constructions of two-dimensional BSECs with block size 3. The existence problem is completely solved in the case where lambda = 1.

Do female rainbowfish (Melanotaenia spp.) prefer to shoal with familiar individuals under predation pressure?

Shoaling with familiar individuals may have many benefits including enhanced escape responses or increased foraging efficiency. This study describes the results of two complimentary experiments. The first utilised a simple binary choice experiment to determine if rainbowfish (Melanotaenia spp.) preferred to shoal with familiar individuals or with strangers. The second experiment used a free range situation where familiar and unfamiliar individuals were free to intermingle and were then exposed to a predator threat. Like many other small species of fish, rainbowfish were capable of identifying and distinguishing between individuals and choose to preferentially associate with familiar individuals as opposed to strangers. Contrary to expectations. however. rainbowrish did not significantly increase their preference for familiar individuals in the presence of a stationary predator model. Griffiths [J Fish Biol (1997) 51:489-4951 conducted similar studies under semi-natural conditions examining, the shoaling preferences of European minnows and showed similar results. Both the current study and that of Griffiths were conducted using predator wary populations of fish. It is suggested that, in predator sympatric populations, the benefits of shoaling with familiar individuals are such that it always pays to stay close to familiar individuals even when the probability If predator attack is remote.