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.

Control of chaotic instability in a dual-spin spacecraft with dissipation using energy methods

Control of chaotic instability in a rotating multibody system in the form of a dual-spin spacecraft with an axial nutational damper is achieved using an algorithm derived using energy methods. The control method is implemented on two realistic spacecraft parameter configurations which have been found to exhibit chaotic instability when a sinusoidally varying torque is applied to the spacecraft for a range of forcing amplitudes and frequencies. Such a torque, in practice, may arise under malfunction of the control system or from an unbalanced rotor. Chaotic instabilities arising from these torques could introduce uncertainties and irregularities into a spacecraft's attitude and consequently impair pointing accuracy. The control method is formulated from nutational stability results derived using an energy sink approximation for a dual-spin spacecraft with an asymmetric platform and axisymmetric rotor. The effectiveness of the control method is shown numerically and the results are studied by means of time history, phase space, Poincare map, Lyapunov characteristic exponents and Bifurcation diagrams.

The measurement of actual acidity in acid sulfate soils and the determination of sulfidic acidity in suspension after peroxide oxidation

Improvements to the routine methods for the determination of actual acidity in suspension for acid sulfate soils (ASS) are introduced. The titratable sulfidic acidity (TSA) results using an improved peroxide-based method were compared with the theoretical acidity predicted by the chromium reducible sulfur method for 9 acid sulfate soils. The regression between these 2 measures of sulfidic acidity was highly significant, the slope of the regression line not significantly different from unity (P = 0.05) and the intercept not significantly different from zero. This contrasts with results of other workers using earlier peroxide oxidation methods, where TSA substantially underestimated the theoretical acidity predicted by reduced inorganic sulfur analysis. Comparison was made between the 2 principal measurements from the improved peroxide method (TSA and S-POS), with S-POS converted to theoretical sulfidic acidity to allow comparison. The relationship between these 2 measurements was highly significant. The effects of titration in suspension, as well as raising titration end points to pH 6.5, were investigated, principally with respect to the titratable actual acidity (TAA) result. TAA results obtained by KCl extraction were compared with those obtained using BaCl2, MgCl2, and water extraction. TAA in 1 M KCl suspensions titrated to pH 6.5 agreed well with titratable actual acidity measured using the 25-h extraction approach of the Lin et al. (2000a) BaCl2 method. Both BaCl2 and KCl solutions were ineffective at fully recovering acidity from synthetic jarosite without repeated extraction and titration. The application of correction factors for the estimation of total actual acidity in ASS is not supported by the results of this investigation. Acid sulfate soils that contain substantial quantities of jarosite or other acid-producing but relatively insoluble sulfate minerals continue to prove problematic to chemically analyse; however, an approach for estimating this component is discussed.

Improvements to peroxide oxidation methods for analysing sulfur in acid sulfate soils

Improvements to peroxide oxidation methods for analysing acid sulfate soils (ASS) are introduced. The soil solution ratio has been increased to 1 : 40, titrations are performed in suspension, and the duration of the peroxide digest stage is substantially shortened. For 9 acid sulfate soils, the peroxide oxidisable sulfur value obtained using the improved method was compared with the reduced inorganic sulfur result obtained using the chromium reducible sulfur method. Their regression was highly significant, the slope of the regression line was not significantly different (P = 0.05) from unity, and the intercept not significantly different from zero. A complete sulfur budget for the improved method showed there was no loss of sulfur as has been reported for earlier peroxide oxidation techniques. When soils were very finely ground, efficient oxidation of sulfides was achieved, despite the milder digestion conditions. Highly sulfidic and organic soils were shown to be the most difficult to analyse using either the improved method or the chromium method. No single analytical method can be universally applied to all ASS, rather a suite of methods is necessary for a thorough understanding of many ASS. The improved peroxide method, in combination with the chromium method and the 4 M HCl extraction, form a sound platform for informed decision making on the management of acid sulfate soils.

Predictor-corrector methods of Runge-Kutta type for stochastic differential equations

In this paper we construct predictor-corrector (PC) methods based on the trivial predictor and stochastic implicit Runge-Kutta (RK) correctors for solving stochastic differential equations. Using the colored rooted tree theory and stochastic B-series, the order condition theorem is derived for constructing stochastic RK methods based on PC implementations. We also present detailed order conditions of the PC methods using stochastic implicit RK correctors with strong global order 1.0 and 1.5. A two-stage implicit RK method with strong global order 1.0 and a four-stage implicit RK method with strong global order 1.5 used as the correctors are constructed in this paper. The mean-square stability properties and numerical results of the PC methods based on these two implicit RK correctors are reported.

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.

Comparison of two population pharmacokinetic programs, NONMEM and P-PHARM, for tacrolimus

Objectives: To compare the population modelling programs NONMEM and P-PHARM during investigation of the pharmacokinetics of tacrolimus in paediatric liver-transplant recipients. Methods: Population pharmacokinetic analysis was performed using NONMEM and P-PHARM on retrospective data from 35 paediatric liver-transplant patients receiving tacrolimus therapy. The same data were presented to both programs. Maximum likelihood estimates were sought for apparent clearance (CL/F) and apparent volume of distribution (V/F). Covariates screened for influence on these parameters were weight, age, gender, post-operative day, days of tacrolimus therapy, transplant type, biliary reconstructive procedure, liver function tests, creatinine clearance, haematocrit, corticosteroid dose, and potential interacting drugs. Results: A satisfactory model was developed in both programs with a single categorical covariate - transplant type - providing stable parameter estimates and small, normally distributed (weighted) residuals. In NONMEM, the continuous covariates - age and liver function tests - improved modelling further. Mean parameter estimates were CL/F (whole liver) = 16.3 1/h, CL/F (cut-down liver) = 8.5 1/h and V/F = 565 1 in NONMEM, and CL/F = 8.3 1/h and V/F = 155 1 in P-PHARM. Individual Bayesian parameter estimates were CL/F (whole liver) = 17.9 +/- 8.8 1/h, CL/F (cutdown liver) = 11.6 +/- 18.8 1/h and V/F = 712 792 1 in NONMEM, and CL/F (whole liver) = 12.8 +/- 3.5 1/h, CL/F (cut-down liver) = 8.2 +/- 3.4 1/h and V/F = 221 1641 in P-PHARM. Marked interindividual kinetic variability (38-108%) and residual random error (approximately 3 ng/ml) were observed. P-PHARM was more user friendly and readily provided informative graphical presentation of results. NONMEM allowed a wider choice of errors for statistical modelling and coped better with complex covariate data sets. Conclusion: Results from parametric modelling programs can vary due to different algorithms employed to estimate parameters, alternative methods of covariate analysis and variations and limitations in the software itself.

Population pharmacokinetics of tacrolimus in adult kidney transplant recipients

Objectives: The aims of this study were to investigate the population pharmacokinetics of tacrolimus in adult kidney transplant recipients and to identify factors that explain variability. Methods: Population analysis was performed on retrospective data from 70 patients who received oral tacrolimus twice daily. Morning blood trough concentrations were measured by liquid chromatography-tandem mass spectrometry. Maximum likelihood estimates were sought for apparent clearance (CL/F) and apparent volume of distribution (V/F), with the use of NONMEM (GloboMax LLC, Hanover, Md). Factors screened for influence on these parameters were weight, age, gender, postoperative day, days of tacrolimus therapy, liver function tests, creatinine clearance, hematocrit fraction, corticosteroid dose, and potential interacting drugs. Results. CL/F was greater in patients with abnormally low hematocrit fraction (data from 21 patients only), and it decreased with increasing days of therapy and AST concentrations (P