Modelling the Distribution of Ischaemic Stroke-Specific Survival Time Using an EM-based Mixture Approach with Random Effects Adjustment

A two-component survival mixture model is proposed to analyse a set of ischaemic stroke-specific mortality data. The survival experience of stroke patients after index stroke may be described by a subpopulation of patients in the acute condition and another subpopulation of patients in the chronic phase. To adjust for the inherent correlation of observations due to random hospital effects, a mixture model of two survival functions with random effects is formulated. Assuming a Weibull hazard in both components, an EM algorithm is developed for the estimation of fixed effect parameters and variance components. A simulation study is conducted to assess the performance of the two-component survival mixture model estimators. Simulation results confirm the applicability of the proposed model in a small sample setting. Copyright (C) 2004 John Wiley Sons, Ltd.

MULTIDIMENSIONAL SCALING OF MUSICAL TIME ESTIMATIONS

The aim of this study was to identify the psycho-musical factors that govern time evaluation in Western music from baroque, classic, romantic, and modern repertoires. The excerpts were previously found to represent variability in musical properties and to induce four main categories of emotions. 48 participants (musicians and nonmusicians) freely listened to 16 musical excerpts (lasting 20 sec. each) and grouped those that seemed to have the same duration. Then, participants associated each group of excerpts to one of a set of sine wave tones varying in duration from 16 to 24 sec. Multidimensional scaling analysis generated a two-dimensional solution for these time judgments. Musical excerpts with high arousal produced an overestimation of time, and affective valence had little influence on time perception. The duration was also overestimated when tempo and loudness were higher, and to a lesser extent, timbre density. In contrast, musical tension had little influence.

Real-time detection of pneumothorax using electrical impedance tomogyaphy

Objectives: Pneumothorax is a frequent complication during mechanical ventilation. Electrical impedance tomography (EIT) is a noninvasive tool that allows real-time imaging of regional ventilation. The purpose of this study was to 1) identify characteristic changes in the EIT signals associated with pneumothoraces; 2) develop and fine-tune an algorithm for their automatic detection; and 3) prospectively evaluate this algorithm for its sensitivity and specificity in detecting pneumothoraces in real time. Design: Prospective controlled laboratory animal investigation. Setting: Experimental Pulmonology Laboratory of the University of Sao Paulo. Subjects: Thirty-nine anesthetized mechanically ventilated supine pigs (31.0 +/- 3.2 kg, mean +/- SD). Interventions. In a first group of 18 animals monitored by EIT, we either injected progressive amounts of air (from 20 to 500 mL) through chest tubes or applied large positive end-expiratory pressure (PEEP) increments to simulate extreme lung overdistension. This first data set was used to calibrate an EIT-based pneumothorax detection algorithm. Subsequently, we evaluated the real-time performance of the detection algorithm in 21 additional animals (with normal or preinjured lungs), submitted to multiple ventilatory interventions or traumatic punctures of the lung. Measurements and Main Results: Primary EIT relative images were acquired online (50 images/sec) and processed according to a few imaging-analysis routines running automatically and in parallel. Pneumothoraces as small as 20 mL could be detected with a sensitivity of 100% and specificity 95% and could be easily distinguished from parenchymal overdistension induced by PEEP or recruiting maneuvers, Their location was correctly identified in all cases, with a total delay of only three respiratory cycles. Conclusions. We created an EIT-based algorithm capable of detecting early signs of pneumothoraces in high-risk situations, which also identifies its location. It requires that the pneumothorax occurs or enlarges at least minimally during the monitoring period. Such detection was operator-free and in quasi real-time, opening opportunities for improving patient safety during mechanical ventilation.

Modeling the optical constants of solids using acceptance-probability-controlled simulated annealing with an adaptive move generation procedure

The acceptance-probability-controlled simulated annealing with an adaptive move generation procedure, an optimization technique derived from the simulated annealing algorithm, is presented. The adaptive move generation procedure was compared against the random move generation procedure on seven multiminima test functions, as well as on the synthetic data, resembling the optical constants of a metal. In all cases the algorithm proved to have faster convergence and superior escaping from local minima. This algorithm was then applied to fit the model dielectric function to data for platinum and aluminum.

Analysis and implementation of a new constant acceleration subcycling algorithm

An algorithm for explicit integration of structural dynamics problems with multiple time steps is proposed that averages accelerations to obtain subcycle states at a nodal interface between regions integrated with different time steps. With integer time step ratios, the resulting subcycle updates at the interface sum to give the same effect as a central difference update over a major cycle. The algorithm is shown to have good accuracy, and stability properties in linear elastic analysis similar to those of constant velocity subcycling algorithms. The implementation of a generalised form of the algorithm with non-integer time step ratios is presented. (C) 1997 by John Wiley & Sons, Ltd.

The application of adaptive partitioned random search in feature selection problem

Feature selection is one of important and frequently used techniques in data preprocessing. It can improve the efficiency and the effectiveness of data mining by reducing the dimensions of feature space and removing the irrelevant and redundant information. Feature selection can be viewed as a global optimization problem of finding a minimum set of M relevant features that describes the dataset as well as the original N attributes. In this paper, we apply the adaptive partitioned random search strategy into our feature selection algorithm. Under this search strategy, the partition structure and evaluation function is proposed for feature selection problem. This algorithm ensures the global optimal solution in theory and avoids complete randomness in search direction. The good property of our algorithm is shown through the theoretical analysis.

Rhipicephalus (Boophilus) microplus: Clotting time in tick-infested skin varies according to local inflammation and gene expression patterns in tick salivary glands

Ticks deposit saliva at the site of their attachment to a host in order to inhibit haemostasis, inflammation and innate and adaptive immune responses. The anti-haemostatic properties of tick saliva have been described by many studies, but few show that tick infestations or its anti-haemostatic components exert systemic effects in vivo. In the present study, we extended these observations and show that, compared with normal skin, bovine hosts that are genetically susceptible to tick infestations present an increase in the clotting time of blood collected from the immediate vicinity of haemorrhagic feeding pools in skin infested with different developmental stages of Rhipicepahlus microplus; conversely, we determined that clotting time of tick-infested skin from genetically resistant bovines was shorter than that of normal skin. Coagulation and inflammation have many components in common and we determined that in resistant bovines, eosinophils and basophils, which are known to contain tissue factor, are recruited in greater numbers to the inflammatory site of tick bites than in susceptible hosts. Finally, we correlated the observed differences in clotting times with the expression profiles of transcripts for putative anti-haemostatic proteins in different developmental stages of R. microplus fed on genetically susceptible and resistant hosts: we determined that transcripts coding for proteins similar to these molecules are overrepresented in salivary glands from nymphs and males fed on susceptible bovines. Our data indicate that ticks are able to modulate their host`s local haemostatic reactions. In the resistant phenotype, larger amounts of inflammatory cells are recruited and expression of anti-coagulant molecules is decreased tick salivary glands, features that can hamper the tick`s blood meal. (C) 2010 Elsevier Inc. All rights reserved.

Cytokine Signatures of Innate and Adaptive Immunity in 17DD Yellow Fever Vaccinated Children and Its Association With the Level of Neutralizing Antibody

Background. The live attenuated yellow fever (YF) vaccines have been available for decades and are considered highly effective and one of the safest vaccines worldwide. Methods. The impact of YF-17DD-antigens recall on cytokine profiles of YF-17DD-vaccinated children were characterized using short-term cultures of whole blood samples and single-cell flow cytometry. This study enrolled seroconverters and nonseroconverters after primovaccination (PV-PRNT(+) and PV-PRNT(-)), seroconverters after revaccination (RV-PRNT(+)), and unvaccinated volunteers (UV-PRNT(-)). Results. The analysis demonstrated in the PV-PRNT(+) group a balanced involvement of pro-inflammatory/regulatory adaptive immunity with a prominent participation of innate immunity pro-inflammatory events (IL-12(+) and TNF-alpha(+) NEU and MON). Using the PV-PRNT(+) cytokine signature as a reference profile, PV-PRNT(+) presented a striking lack of innate immunity proinflammatory response along with an increased adaptive regulatory profile (IL-4(+) CD4(+) T cells and IL-10(+) and IL-5(+) CD8(+) T cells). Conversely, the RV-PRNT(+) shifted the overall cytokine signatures toward an innate immunity pro-inflammatory profile and restored the adaptive regulatory response. Conclusions. The data demonstrated that the overall cytokine signature was associated with the levels of PRNT antibodies with a balanced innate/adaptive immunity with proinflammatory/regulatory profile as the hallmark of PV-PRNT(MEDIUM+), whereas a polarized regulatory response was observed in PV-PRNT(-) and a prominent proinflammatory signature was the characteristic of PV-PRNT(HIGH+).

Anisotropy-based performance analysis of linear discrete time invariant control systems

The anisotropic norm of a linear discrete-time-invariant system measures system output sensitivity to stationary Gaussian input disturbances of bounded mean anisotropy. Mean anisotropy characterizes the degree of predictability (or colouredness) and spatial non-roundness of the noise. The anisotropic norm falls between the H-2 and H-infinity norms and accommodates their loss of performance when the probability structure of input disturbances is not exactly known. This paper develops a method for numerical computation of the anisotropic norm which involves linked Riccati and Lyapunov equations and an associated special type equation.

Tolerance or immunity to a tumor antigen expressed in somatic cells can be determined by systemic proinflammatory signals at the time of first antigen exposure

Mice transgenic for the E7 tumor Ag of human papillomavirus type 16, driven from a keratin 14 promoter, express E7 in keratinocytes but not dendritic cells. Grafted E7-transgenic skin is not rejected by E7-immunized mice that reject E7-transduced transplantable tumors. Rejection of recently transplanted E7-transgenic skin grafts, but not of control nontransgenic grafts or of established E7-transgenic grafts, is induced by systemic administration of live or killed Listeria monocytogenes or of endotoxin. Graft recipients that reject an E7 graft reject a subsequent E7 graft more rapidly and without further L. monocytogenes exposure, whereas recipients of an E7 graft given without L. monocytogenes do not reject a second graft, even if given with L. monocytogenes. Thus, cross-presentation of E7 from keratinocytes to the adaptive immune system occurs with or without a proinflammatory stimulus, but proinflammatory stimuli at the time of first cross-presentation of Ag can determine the nature of the immune response to the Ag. Furthermore, immune effector mechanisms responsible for rejection of epithelium expressing a tumor Ag in keratinocytes are different from those that reject an E7-expressing transplantable tumor. These observations have implications for immunotherapy for epithelial cancers.

Time-dependent master equation simulation of complex elementary reactions in combustion: Application to the reaction of (CH2)-C-1 with C2H2 from 300-2000 K

Computational simulations of the title reaction are presented, covering a temperature range from 300 to 2000 K. At lower temperatures we find that initial formation of the cyclopropene complex by addition of methylene to acetylene is irreversible, as is the stabilisation process via collisional energy transfer. Product branching between propargyl and the stable isomers is predicted at 300 K as a function of pressure for the first time. At intermediate temperatures (1200 K), complex temporal evolution involving multiple steady states begins to emerge. At high temperatures (2000 K) the timescale for subsequent unimolecular decay of thermalized intermediates begins to impinge on the timescale for reaction of methylene, such that the rate of formation of propargyl product does not admit a simple analysis in terms of a single time-independent rate constant until the methylene supply becomes depleted. Likewise, at the elevated temperatures the thermalized intermediates cannot be regarded as irreversible product channels. Our solution algorithm involves spectral propagation of a symmetrised version of the discretized master equation matrix, and is implemented in a high precision environment which makes hitherto unachievable low-temperature modelling a reality.

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.

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.

Control of chaotic instabilities in a spinning spacecraft with dissipation using Lyapunov's method

Control of chaotic instability in a simplified model of a spinning spacecraft with dissipation is achieved using an algorithm derived using Lyapunov's second method. The control method is implemented on a realistic spacecraft parameter configuration which has been found to exhibit chaotic instability for a range of forcing amplitudes and frequencies when a sinusoidally varying torque is applied to the spacecraft. Such a torque, may arise in practice from an unbalanced rotor or from vibrations in appendages. Numerical simulations are performed and the results are studied by means of time history, phase space, Poincare map, Lyapunov characteristic exponents and bifurcation diagrams. (C) 2002 Elsevier Science Ltd. All rights reserved.

Efficient time-independent wave packet scattering calculations within a Lanczos subspace: H+O-2 (J=0) state-to-state reaction probabilities

An efficient Lanczos subspace method has been devised for calculating state-to-state reaction probabilities. The method recasts the time-independent wave packet Lippmann-Schwinger equation [Kouri , Chem. Phys. Lett. 203, 166 (1993)] inside a tridiagonal (Lanczos) representation in which action of the causal Green's operator is affected easily with a QR algorithm. The method is designed to yield all state-to-state reaction probabilities from a given reactant-channel wave packet using a single Lanczos subspace; the spectral properties of the tridiagonal Hamiltonian allow calculations to be undertaken at arbitrary energies within the spectral range of the initial wave packet. The method is applied to a H+O-2 system (J=0), and the results indicate the approach is accurate and stable. (C) 2002 American Institute of Physics.