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.

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.

Model specification tests in nonparametric stochastic regression models

In this paper, we consider testing for additivity in a class of nonparametric stochastic regression models. Two test statistics are constructed and their asymptotic distributions are established. We also conduct a small sample study for one of the test statistics through a simulated example. (C) 2002 Elsevier Science (USA).

Cross-national comparison of internal migration: issues and measures

Our objectives are to identify the issues that researchers encounter when measuring internal migration in different countries and to propose key indicators that analysts can use to compare internal migration at the 'national' level. We establish the benefits to be gained by a rigorous approach to cross-national comparisons of internal migration and discuss issues that affect such comparisons. We then distinguish four dimensions of internal migration on which countries can be compared and, for each dimension, identify a series of summary measures. We illustrate the issues and measures proposed by comparing migration in Australia and Great Britain.

Estimating and evaluating the statistics of gapped local-alignment scores

We present a novel maximum-likelihood-based algorithm for estimating the distribution of alignment scores from the scores of unrelated sequences in a database search. Using a new method for measuring the accuracy of p-values, we show that our maximum-likelihood-based algorithm is more accurate than existing regression-based and lookup table methods. We explore a more sophisticated way of modeling and estimating the score distributions (using a two-component mixture model and expectation maximization), but conclude that this does not improve significantly over simply ignoring scores with small E-values during estimation. Finally, we measure the classification accuracy of p-values estimated in different ways and observe that inaccurate p-values can, somewhat paradoxically, lead to higher classification accuracy. We explain this paradox and argue that statistical accuracy, not classification accuracy, should be the primary criterion in comparisons of similarity search methods that return p-values that adjust for target sequence length.

Determinants of tissue Doppler measures of regional diastolic function during dobutamine stress echocardiography

Background Diastolic dysfunction induced by ischemia may alter transmitral blood flow, but this reflects global ventricular function, and pseudonormalization may occur with increased preload. Tissue Doppler may assess regional diastolic function and is relatively load-independent, but limited data exist regarding its application to stress testing. We sought to examine the stress response of regional diastolic parameters to dobutomine echocardiography (DbE). Methods Sixty-three patients underwent study with DbE: 20 with low probability of coronary artery disease (CAD) and 43 with CAD who underwent angiography. A standard DbE protocol was used, and segments were categorized as ischemic, scar, or normal. Color tissue Doppler was acquired at baseline and peak stress, and waveforms in the basal and mid segments were used to measure early filling (Em), late filling (Am), and E deceleration time. Significant CAD was defined by stenoses >50% vessel diameter. Results Diastolic parameters had limited feasibility because of merging of Em and Am waves at high heart rates and limited reproducibility. Nonetheless, compared with normal segments, segments subtended with significant stenoses showed a lower Em velocity at rest (6.2 +/- 2.6 cm/s vs 4.8 +/- 2.2 cm/s, P < .0001) and peak (7.5 +/- 4.2 cm/s vs 5.1 +/- 3.6 cm/s, P < .0001), Abnormal segments also showed a shorter E deceleration time (51 +/- 27 ms vs 41 +/- 27 ms, P = .0001) at base and peak. No changes were documented in Am. The same pattern was seen with segments identified as ischemic with wall motion score. However, in the absence of ischemia, segments of patients with left ventricular hypertrophy showed a lower Em velocity, with blunted Em responses to stress. Conclusion Regional diastolic function is sensitive to ischemia. However, a number of practical limitations limit the applicability of diastolic parameters for the quantification of stress echocardiography.

Development of a fully quantitative approach to the interpretation of stress echocardiography using radial and longitudinal myocardial velocities

Background: Tissue Doppler may be used to quantify regional left ventricular function but is limited by segmental variation of longitudinal velocity from base to apex and free to septal walls. We sought to overcome this by developing a composite of longitudinal and radial velocities. Methods and Results. We examined 82 unselected patients undergoing a standard dobutamine echocardiogram. Longitudinal velocity was obtained in the basal and mid segments of each wall using tissue Doppler in the apical views. Radial velocities were derived in the same segments using an automated border detection system and centerline method with regional chords grouped according to segment location and temporally averaged. In 25 patients at low probability of coronary disease, the pattern of regional variation in longitudinal velocity (higher in the septum) was the opposite of radial velocity (higher in the free wall) and the combination was homogenous. In 57 patients undergoing angiography, velocity in abnormal segments was less than normal segments using longitudinal (6.0 +/- 3.6 vs 9.0 +/- 2.2 cm/s, P = .01) and radial velocity (6.0 +/- 4.0 vs 8.0 +/- 3.9 cm/s, P = .02). However, the composite velocity permitted better separation of abnormal and normal segments (13.3 +/- 5.6 vs 17.5 +/- 4.2 cm/s, P = .001). There was no significant difference between the accuracy of this quantitative approach and expert visual wall motion analysis (81% vs 84%, P = .56). Conclusion: Regional variation of uni-dimensional myocardial velocities necessitates site-specific normal ranges, probably because of different fiber directions. Combined analysis of longitudinal and radial velocities allows the derivation of a composite velocity, which is homogenous in all segments and may allow better separation of normal and abnormal myocardium.

Quantum dynamics of two coupled qubits

We investigate the difference between classical and quantum dynamics of coupled magnetic dipoles. We prove that in general the dynamics of the classical interaction Hamiltonian differs from the corresponding quantum model, regardless of the initial state. The difference appears as nonpositive-definite diffusion terms in the quantum evolution equation of an appropriate positive phase-space probability density. Thus, it is not possible to express the dynamics in terms of a convolution of a positive transition probability function and the initial condition as can be done in the classical case. It is this feature that enables the quantum system to evolve to an entangled state. We conclude that the dynamics are a quantum element of nuclear magnetic resonance quantum-information processing. There are two limits where our quantum evolution coincides with the classical one: the short-time limit before spin-spin interaction sets in and the long-time limit when phase diffusion is incorporated.

Effects of paternal drinking, conduct disorder and childhood home environment on the development of alcohol use disorders in a Thai population

Aims To identify influences on the development of alcohol use disorders in a Thai population, particularly parental drinking and childhood environment. Design Case-control study. Setting A university hospital, a regional hospital and a community hospital in southern Thailand. Participants Ninety-one alcohol-dependents and 177 hazardous/harmful drinkers were recruited as cases and 144 non-or infrequent drinkers as controls. Measurements Data on parental drinking, family demographic characteristics, family activities, parental disciplinary practice, early religious life and conduct disorder were obtained using a structured interview questionnaire. The main outcome measure was the subject's classification as alcohol-dependent, hazardous/harmful drinker or non-/infrequent drinker. Findings A significant relationship was found between having a drinking father and the occurrence of hazardous/harmful drinking or alcohol dependence in the subjects. Childhood factors (conduct disorder and having been a temple boy, relative probability ratios, RPRs and 95% CI: 6.39, 2.81-14.55 and 2.21, 1.19-4.08, respectively) also significantly predicted alcohol dependence, while perceived poverty and ethnic alienation was reported less frequently by hazardous/harmful drinkers and alcohol-dependents (RPRS and 95% CIs = 0.34, 0.19-0.62 and 0.59, 0.38-0.93, respectively) than the controls. The relative probability ratio for the effect of the father's infrequent drinking on the son's alcohol dependence was 2.92 (95% CI = 1.42-6.02) and for the father's heavy or dependent drinking 2.84 (95% CI=1.31-6.15). Conclusions Being exposed to a light-drinking, father increases the risk of a son's alcohol use disorders exhibited either as hazardous-harmful or dependent drinking. However, exposure to a heavy- or dependent-drinking father is associated more uniquely with an increased risk of his son being alcohol-dependent. The extent to which this is seen in other cultures is worthy of exploration.

An atomistic simulation of solid state sintering using Monte Carlo methods

This paper presents results on the simulation of the solid state sintering of copper wires using Monte Carlo techniques based on elements of lattice theory and cellular automata. The initial structure is superimposed onto a triangular, two-dimensional lattice, where each lattice site corresponds to either an atom or vacancy. The number of vacancies varies with the simulation temperature, while a cluster of vacancies is a pore. To simulate sintering, lattice sites are picked at random and reoriented in terms of an atomistic model governing mass transport. The probability that an atom has sufficient energy to jump to a vacant lattice site is related to the jump frequency, and hence the diffusion coefficient, while the probability that an atomic jump will be accepted is related to the change in energy of the system as a result of the jump, as determined by the change in the number of nearest neighbours. The jump frequency is also used to relate model time, measured in Monte Carlo Steps, to the actual sintering time. The model incorporates bulk, grain boundary and surface diffusion terms and includes vacancy annihilation on the grain boundaries. The predictions of the model were found to be consistent with experimental data, both in terms of the microstructural evolution and in terms of the sintering time. (C) 2002 Elsevier Science B.V. All rights reserved.