On the importance function in splitting simulation

The splitting method is a simulation technique for the estimation of very small probabilities. In this technique, the sample paths are split into multiple copies, at various stages in the simulation. Of vital importance to the efficiency of the method is the Importance Function (IF). This function governs the placement of the thresholds or surfaces at which the paths are split. We derive a characterisation of the optimal IF and show that for multi-dimensional models the natural choice for the IF is usually not optimal. We also show how nearly optimal splitting surfaces can be derived or simulated using reverse time analysis. Our numerical experiments illustrate that by using the optimal IF, one can obtain a significant improvement in simulation efficiency.

Investigation of electrical conductivity as a function of dopant-ion radius in the systems Zr0.75Ce0.08M0.17O1.92 (M=Nd, Sm, Gd, Dy, Ho, Y, Er, Yb, Sc)

Electrical conductivity versus dopant ionic radius studies in zirconia- and ceria-based, solid oxide fuel cell (SOFC) electrolyte systems have shown that oxygen-ion conductivity is highest when the host and dopant ions are similar in size [J. Am. Ceram. Soc. 48 (1965) 286; Solid State Ionics 37 (1989) 67; Solid State Ionics 5 (1981) 547]. Under these conditions, it is thought that the conduction paths within the crystal lattice become less distorted [Solid State Ionics 8 (1983) 201]. In this study, binary ZrO2-M2O3 unit cells were expanded, via the partial substitution of Ce+4 for Zr+4 into the lattice, in an attempt to identify new, ternary, zirconia/ceria-based electrolyte systems with enhanced electrical conductivity. The compositions Zr0.75Ce0.08M0.17O1.92 (M = Nd, Sm, Gd, Dy, Ho, Y, Yb, Sc) were prepared using traditional solid state techniques. Bulk phase characterisation and precise lattice parameter measurements were performed with X-ray diffraction techniques. Four-probe DC conductivity measurements between 400 and 900 degreesC showed that the dopant-ion radius influenced electrical conductivity. The conductivity versus dopant-ion radius trends previously observed in zirconia-based, binary systems are clearly apparent in the ternary systems investigated in this study. The addition of ceria was found to have a negative influence on the electrical conductivity over the temperature range 400-900 degreesC. It is suggested that distortion of the oxygen-ion conduction path by the presence of the larger M+3 and Ce+4 species (relative to Zr+4) is the reason for the decreasing electrical conductivity as a function of increasing dopant size and ceria addition, respectively. (C) 2002 Elsevier Science B.V. All rights reserved.

Impact of malocclusion and orthognathic reconstruction surgery on resonance and articulatory function: an examination of variability in five cases

Articulatory patterns and nasal resonance were assessed before and 6 months after orthognathic reconstruction surgery in five patients with dentofacial deformities. Perceptual and physiological assessments showed disorders of nasality and articulatory function preoperatively, two patients being hyponasal, and one hypernasal. Four patients had mild articulatory deficits, and four had reduced maximal lip or tongue pressures. Operation resulted in different patterns of change. Nasality deteriorated in three patients and articulatory precision and intelligibility improved in only one patient and showed no change in the other four. Operation improved interlabial pressures in three patients, while its impact on tongue pressures varied, being improved in one case, deteriorating in one, and remaining unchanged in the other three. The variability in the results highlights the need for routine assessment of speech and resonance before and after orthognathic reconstruction. (C) 2002 The British Association of Oral and Maxillofacial Surgeons. Published by Elsevier Science Ltd. All rights reserved.

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.

Empirical evidence for ultrametric structure in multi-layer perceptron error surfaces

Combinatorial optimization problems share an interesting property with spin glass systems in that their state spaces can exhibit ultrametric structure. We use sampling methods to analyse the error surfaces of feedforward multi-layer perceptron neural networks learning encoder problems. The third order statistics of these points of attraction are examined and found to be arranged in a highly ultrametric way. This is a unique result for a finite, continuous parameter space. The implications of this result are discussed.

The impact of genotyping error on haplotype reconstruction and frequency estimation

The choice of genotyping families vs unrelated individuals is a critical factor in any large-scale linkage disequilibrium (LD) study. The use of unrelated individuals for such studies is promising, but in contrast to family designs, unrelated samples do not facilitate detection of genotyping errors, which have been shown to be of great importance for LD and linkage studies and may be even more important in genotyping collaborations across laboratories. Here we employ some of the most commonly-used analysis methods to examine the relative accuracy of haplotype estimation using families vs unrelateds in the presence of genotyping error. The results suggest that even slight amounts of genotyping error can significantly decrease haplotype frequency and reconstruction accuracy, that the ability to detect such errors in large families is essential when the number/complexity of haplotypes is high (low LD/common alleles). In contrast, in situations of low haplotype complexity (high LD and/or many rare alleles) unrelated individuals offer such a high degree of accuracy that there is little reason for less efficient family designs. Moreover, parent-child trios, which comprise the most popular family design and the most efficient in terms of the number of founder chromosomes per genotype but which contain little information for error detection, offer little or no gain over unrelated samples in nearly all cases, and thus do not seem a useful sampling compromise between unrelated individuals and large families. The implications of these results are discussed in the context of large-scale LD mapping projects such as the proposed genome-wide haplotype map.

A curious experiment: the paradigm switch from observation and speculation to experimentation, in the understanding of neuromuscular function and disease

The four-link chain of the motor unit represents the contemporary end-point of some two millennia of evolving knowledge in neuroscience. The paradigm shift in neuromuscular epistemology occurred in the mid-17th century. In 1666, the newly graduated Dutch doctor, Jan Swammerdam (1637-1680) published his former investigations of dissected nerve-muscle preparations. These experiments comprised the quantum leap from observation and speculation, to that of experimentation in the field of neuroanatomy and neurophysiology. In what he termed 'A Curious Experiment' he also described the phenomenon of intrinsic muscle excitability - I cannot observe that the muscle in the living animal ever absolutely ceases from all motion. Eighty years later (1752), von Haller demonstrated experimentally that irritability (contractility) was an intrinsic property of all muscular tissue; and distinguished between the sensibility of nerve impulses and the irritability of muscular contraction. This experimental progression from Swammerdam to von Haller culminated in 1850, when Claude Bernard's studies in experimental pharmacology confirmed that muscle was a functional unit, independent of any electrical innervation via its supplying nerve. This account comprises an audit of Swammerdam's work in the perspective of neuromuscular knowledge. (C) 2002 Elsevier Science B.V. All rights reserved.

Exercise and T-lymphocyte function: a comparison of proliferation in PBMC and NK cell-depleted PBMC culture

This study utilized recently developed microbead technology to remove natural killer (NK) cells from peripheral blood mononuclear cell (PBMC) preparations to determine the effect of acute exercise on T-lymphocyte function, independent of changes in lymphocyte subpopulations. Twelve well-trained male runners completed a 60-min exercise trial at 95% ventilatory threshold and a no-exercise control trial. Six blood samples were taken at each session: before exercise, midexercise, immediately after exercise, and 30, 60, and 90 min after exercise. Isolated PBMC and NK cell-depleted PBMC were stimulated with the mitogen phytohemagglutinin. Cellular proliferation was assessed by using 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide dye uptake. In the PBMC cultures, there was a significantly lower mitogen response to phytohemagglutinin in exercise compared with the control condition immediately postexercise. There were no significant differences between the control and exercise conditions in NK cell-depleted PBMC cultures or in the responses adjusted for the percentage of CD3 cells. The present findings do not support the view that T-lymphocyte function is reduced after exercise.

Coordination and movement pathology: models of structure and function

Here we consider the role of abstract models in advancing our understanding of movement pathology. Models of movement coordination and control provide the frameworks necessary for the design and interpretation of studies of acquired and developmental disorders. These models do not however provide the resolution necessary to reveal the nature of the functional impairments that characterise specific movement pathologies. In addition, they do not provide a mapping between the structural bases of various pathologies and the associated disorders of movement. Current and prospective approaches to the study and treatment of movement disorders are discussed. It is argued that the appreciation of structure-function relationships, to which these approaches give rise, represents a challenge to current models of interlimb coordination, and a stimulus for their continued development. (C) 2002 Elsevier Science B.V. All rights reserved.

Structure and function of plant toxins (with emphasis on cystine knot toxins)

Plant toxins are substances produced and secreted by plants to defend themselves against predators. In a broad sense, this includes all substances that have a toxic effect on targeted organisms, whether they are microbes, other plants, insects, or higher animals. Plant toxins have a diverse range of structures, from small organic molecules through to proteins. This review gives an overview of the various classes of plant toxins but focuses on an interesting class of protein-based plant toxins containing a cystine knot motif. This structural motif confers exceptional stability on proteins containing it and is associated with a wide range of biological activities. The biological activities and structural stability offer many potential applications in the pharmaceutical and agricultural fields. One particularly exciting prospect is in the use of protein-based plant toxins as molecular scaffolds for displaying pharmaceutically important bioactivities. Future applications of plant toxins are likely to involve genetic engineering techniques and molecular pharming approaches.

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.