Modelling High-Dimensional Data by Mixtures of Factor Analyzers

We focus on mixtures of factor analyzers from the perspective of a method for model-based density estimation from high-dimensional data, and hence for the clustering of such data. This approach enables a normal mixture model to be fitted to a sample of n data points of dimension p, where p is large relative to n. The number of free parameters is controlled through the dimension of the latent factor space. By working in this reduced space, it allows a model for each component-covariance matrix with complexity lying between that of the isotropic and full covariance structure models. We shall illustrate the use of mixtures of factor analyzers in a practical example that considers the clustering of cell lines on the basis of gene expressions from microarray experiments. (C) 2002 Elsevier Science B.V. All rights reserved.

The embryonic development of the bodywall and nervous system of the cestode flatworm Hymenolepis diminuta

Cestodes (tapeworms) are a derived, parasitic clade of the phylum Platyhelminthes (flatworms). The cestode body wall represents an adaptation to its endoparasitic lifestyle. The epidermis forms a nonciliated syncytium, and both muscular and nervous system are reduced. Morphological differences between cestodes and free-living flatworms become apparent already during early embryogenesis. Cestodes have a complex life cycle that begins with an infectious larva, called the oncosphere. In regard to cell number, cestode oncospheres are among the simplest multicellular organisms, containing in the order of 50-100 cells. As part of our continuing effort to analyze embryonic development in flatworms, we describe here the staining pattern obtained with acTub in embryos and larvae of the cestode Hymenolepis diminuta and, briefly, the monogenean Neoheterocotyle rhinobatidis. In addition, we labeled the embryonic musculature of Hymenolepis with phalloidin. In Hymenolepis embryos, two different cell types that we interpret as neurons and epidermal gland cells express acTub. There exist only two neurons that develop close to the midline at the anterior pole of the embryo. The axons of these two neurons project posteriorly into the center of the oncosphere, where they innervate the complex of muscles that is attached to the booklets. In addition to neurons, acTub labels a small and invariant set of epidermal gland cells that develop at superficial positions, anteriorly adjacent to the neurons, in the dorsal midline, and around the posteriorly located hooklets. During late stages of embryogenesis they spread and form a complete covering of the embryo. We discuss these data in the broader context of platyhelminth embryology.

Initial effects of elbow taping on pain-free grip strength and pressure pain threshold

Study design: Single-blind, placebo control, randomized, crossover, experimental Study with repeated measures, Objective: To determine the initial effects of a taping technique on grip strength and pain in individuals with lateral epicondylalgia. Background: Taping techniques are advocated for chronic musculoskeletal conditions such as lateral epicondylalgia, a prevalent disorder with significant impact on the individual and community. Little evidence exists supporting the effects of taping techniques on musculoskeletal pain. Methods and Measures: Sixteen participants (mean age +/- SD, 45.8 +/- 10.2 years) with chronic lateral epicondylalgia (rnean duration +/- SD, 13.1 +/- 9.9 months) participated in a placebo control study of an elbow taping technique. Outcome measures were pain-free grip strength and pressure pain threshold taken before, immediately after, and 30 minutes after application of tape. Results: The taping technique significantly improved pain-free grip strength by 24% from baseline (P = .028). The treatment effect was greater than that for placebo and control conditions. Changes in pressure pain threshold (19%), although positive, were not statistically significant. Conclusion: This preliminary study demonstrated an initial ameliorative effect of a taping technique for lateral epicondylalgia and suggests that it should be considered as an adjunct in the management of this condition.

Hypoalgesia induced by elbow manipulation in lateral epicondylalgia does not exhibit tolerance

Previous studies have demonstrated that the initial hypoalgesic effect of spinal manipulative therapy was not antagonized by naloxone and did not exhibit tolerance with repeated applications. The implication is that endogenous opioid mechanisms of pain relief are probably not at play in spinal manipulative therapy. The role of endogenous opioid peptides in manipulation of the peripheral joints has not been investigated. The aim of this study was to evaluate whether the initial hypoalgesic effect of a peripheral manipulative technique (mobilization-with-movement treatment for the elbow) demonstrated a tolerance to repeated applications (ie, reduction in magnitude of effect over repeated applications). Twenty-four participants with unilateral chronic lateral epicondylalgia participated in the study. A repeated measures study was conducted to examine the effect of repeated applications of the mobilization-with-movement treatment for the elbow on 6 separate treatment occasions at least 2 days apart. Pain-free grip strength and pressure pain threshold were chosen as the pain-related outcome measures. Changes in the percent maximum possible effect scores of measures of hypoalgesia were evaluated across the 6 treatment sessions by using linear trend analysis. The results showed no significant difference for the hypoalgesic effect of the treatment technique between sessions (P >.05). This peripheral manipulative therapy treatment technique appeared to have a similar effect profile to previously studied spinal manipulative therapy techniques, thereby contributing to the body of knowledge that indicates that manipulative therapy most likely induces a predominant non-opioid form of analgesia.

Product form approximations for highly linear loss networks with trunk reservation

Admission controls, such as trunk reservation, are often used in loss networks to optimise their performance. Since the numerical evaluation of performance measures is complex, much attention has been given to finding approximation methods. The Erlang Fixed-Point (EFP) approximation, which is based on an independent blocking assumption, has been used for networks both with and without controls. Several more elaborate approximation methods which account for dependencies in blocking behaviour have been developed for the uncontrolled setting. This paper is an exploratory investigation of extensions and synthesis of these methods to systems with controls, in particular, trunk reservation. In order to isolate the dependency factor, we restrict our attention to a highly linear network. We will compare the performance of the resulting approximations against the benchmark of the EFP approximation extended to the trunk reservation setting. By doing this, we seek to gain insight into the critical factors in constructing an effective approximation. (C) 2003 Elsevier Ltd. All rights reserved.

Model-based clustering in gene expression microarrays: An application to breast cancer data

In microarray studies, the application of clustering techniques is often used to derive meaningful insights into the data. In the past, hierarchical methods have been the primary clustering tool employed to perform this task. The hierarchical algorithms have been mainly applied heuristically to these cluster analysis problems. Further, a major limitation of these methods is their inability to determine the number of clusters. Thus there is a need for a model-based approach to these. clustering problems. To this end, McLachlan et al. [7] developed a mixture model-based algorithm (EMMIX-GENE) for the clustering of tissue samples. To further investigate the EMMIX-GENE procedure as a model-based -approach, we present a case study involving the application of EMMIX-GENE to the breast cancer data as studied recently in van 't Veer et al. [10]. Our analysis considers the problem of clustering the tissue samples on the basis of the genes which is a non-standard problem because the number of genes greatly exceed the number of tissue samples. We demonstrate how EMMIX-GENE can be useful in reducing the initial set of genes down to a more computationally manageable size. The results from this analysis also emphasise the difficulty associated with the task of separating two tissue groups on the basis of a particular subset of genes. These results also shed light on why supervised methods have such a high misallocation error rate for the breast cancer data.

Explicit/implicit partitioning and a new explicit form of the generalized alpha method

Most finite element packages use the Newmark algorithm for time integration of structural dynamics. Various algorithms have been proposed to better optimize the high frequency dissipation of this algorithm. Hulbert and Chung proposed both implicit and explicit forms of the generalized alpha method. The algorithms optimize high frequency dissipation effectively, and despite recent work on algorithms that possess momentum conserving/energy dissipative properties in a non-linear context, the generalized alpha method remains an efficient way to solve many problems, especially with adaptive timestep control. However, the implicit and explicit algorithms use incompatible parameter sets and cannot be used together in a spatial partition, whereas this can be done for the Newmark algorithm, as Hughes and Liu demonstrated, and for the HHT-alpha algorithm developed from it. The present paper shows that the explicit generalized alpha method can be rewritten so that it becomes compatible with the implicit form. All four algorithmic parameters can be matched between the explicit and implicit forms. An element interface between implicit and explicit partitions can then be used, analogous to that devised by Hughes and Liu to extend the Newmark method. The stability of the explicit/implicit algorithm is examined in a linear context and found to exceed that of the explicit partition. The element partition is significantly less dissipative of intermediate frequencies than one using the HHT-alpha method. The explicit algorithm can also be rewritten so that the discrete equation of motion evaluates forces from displacements and velocities found at the predicted mid-point of a cycle. Copyright (C) 2003 John Wiley Sons, Ltd.

The first strand transfer reaction of HIV-1 reverse transcription is more efficient in infected cells than in cell-free natural endogenous reverse transcription reactions

Background: In the presence of dNTPs, intact HIV-1 virions are capable of reverse transcribing at least part of their genome, a process known as natural endogenous reverse transcription (NERT). PCR analysis of virion DNA produced by NERT revealed that the first strand transfer reaction (1stST) was inefficient in intact virions, with minus strand (-) strong stop DNA (ssDNA) copy numbers up to 200 times higher than post-1stST products measured using primers in U3 and U5. This was in marked contrast to the efficiency of 1stST observed in single-round cell infection assays, in which (-) ssDNA and U3-U5 copy numbers were indistinguishable. Objectives: To investigate the reasons for the discrepancy in first strand transfer efficiency between intact cell-free virus and the infection process. Study design: Alterations of both NERT reactions and the conditions of cell infection were used to test whether uncoating and/or entry play a role in the discrepancy in first strand transfer efficiency. Results and Conclusions: The difference in 1stST efficiency could not be attributed simply to viral uncoating, since addition of very low concentrations of detergent to NERT reactions removed the viral envelope without disrupting the reverse transcription complex, and these conditions resulted in no improvement in 1stST efficiency. Virus pseudotyped with surface glycoproteins from either vesicular stomatitis virus or amphotrophic murine leukaemia virus also showed low levels of 1stST in low detergent NERT assays and equivalent levels of (-) ssDNA and 1stST in single-round infections of cells, demonstrating that the gp120-mediated infection process did not select for virions capable of carrying out 1stST. These data indicate that a post-entry event or factor may be involved in efficient HIV-1 reverse transcription in vivo. (C) 2002 Elsevier Science B.V. All rights reserved.

Modeling diffraction in free-space optical interconnects by the mode expansion method

Free-space optical interconnects (FSOIs), made up of dense arrays of vertical-cavity surface-emitting lasers, photodetectors and microlenses can be used for implementing high-speed and high-density communication links, and hence replace the inferior electrical interconnects. A major concern in the design of FSOIs is minimization of the optical channel cross talk arising from laser beam diffraction. In this article we introduce modifications to the mode expansion method of Tanaka et al. [IEEE Trans. Microwave Theory Tech. MTT-20, 749 (1972)] to make it an efficient tool for modelling and design of FSOIs in the presence of diffraction. We demonstrate that our modified mode expansion method has accuracy similar to the exact solution of the Huygens-Kirchhoff diffraction integral in cases of both weak and strong beam clipping, and that it is much more accurate than the existing approximations. The strength of the method is twofold: first, it is applicable in the region of pronounced diffraction (strong beam clipping) where all other approximations fail and, second, unlike the exact-solution method, it can be efficiently used for modelling diffraction on multiple apertures. These features make the mode expansion method useful for design and optimization of free-space architectures containing multiple optical elements inclusive of optical interconnects and optical clock distribution systems. (C) 2003 Optical Society of America.

Transdermal penetration of vasoconstrictors - Present understanding and assessment of the human epidermal flux and retention of free bases and ion-pairs

Purpose. As reductions in dermal clearance increase the residence time of solutes in the skin and underlying tissues we compared the topical penetration of potentially useful vasoconstrictors (VCs) through human epidermis as both free bases and ion-pairs with salicylic acid (SA). Methods. We determined the in vitro epidermal flux of ephedrine, naphazoline, oxymetazoline, phenylephrine, and xylometazoline applied as saturated solutions in propylene glycol: water (1: 1) and of ephedrine, naphazoline and tetrahydrozoline as 10% solutions of 1: 1 molar ratio ion-pairs with SA in liquid paraffin. Results. As free bases, ephedrine had the highest maximal flux, Jmax = 77.4 +/- 11.7 mug/cm(2)/h, being 4-fold higher than tetrahydrozoline and xylometazoline, 6-fold higher than phenylephrine, 10-fold higher than naphazoline and 100-fold higher than oxymetazoline. Stepwise regression of solute physicochemical properties identified melting point as the most significant predictor of flux. As ion-pairs with SA, ephedrine and naphazoline had similar fluxes (11.5 +/- 2.3 and 12.0 +/- 1.6 mug/cm(2)/h respectively), whereas tetrahydrozoline was approximately 3-fold slower. Corresponding fluxes of SA from the ion-pairs were 18.6 +/- 0.6, 7.8 +/- 0.8 and 1.1 +/- 0.1 respectively. Transdermal transport of VC's is discussed. Conclusions. Epidermal retention of VCs and SA did not correspond to their molar ratio on application and confirmed that following partitioning into the stratum corneum, ion-pairs separate and further penetration is governed by individual solute characteristics.

A review of drainage and spontaneous rupture in free standing thin films with tangentially immobile interfaces

A review of spontaneous rupture in thin films with tangentially immobile interfaces is presented that emphasizes the theoretical developments of film drainage and corrugation growth through the linearization of lubrication theory in a cylindrical geometry. Spontaneous rupture occurs when corrugations from adjacent interfaces become unstable and grow to a critical thickness. A corrugated interface is composed of a number of waveforms and each waveform becomes unstable at a unique transition thickness. The onset of instability occurs at the maximum transition thickness, and it is shown that only upper and lower bounds of this thickness can be predicted from linear stability analysis. The upper bound is equivalent to the Freakel criterion and is obtained from the zeroth order approximation of the H-3 term in the evolution equation. This criterion is determined solely by the film radius, interfacial tension and Hamaker constant. The lower bound is obtained from the first order approximation of the H-3 term in the evolution equation and is dependent on the film thinning velocity A semi-empirical equation, referred to as the MTR equation, is obtained by combining the drainage theory of Manev et al. [J. Dispersion Sci. Technol., 18 (1997) 769] and the experimental measurements of Radoev et al. [J. Colloid Interface Sci. 95 (1983) 254] and is shown to provide accurate predictions of film thinning velocity near the critical thickness of rupture. The MTR equation permits the prediction of the lower bound of the maximum transition thickness based entirely on film radius, Plateau border radius, interfacial tension, temperature and Hamaker constant. The MTR equation extrapolates to Reynolds equation under conditions when the Plateau border pressure is small, which provides a lower bound for the maximum transition thickness that is equivalent to the criterion of Gumerman and Homsy [Chem. Eng. Commun. 2 (1975) 27]. The relative accuracy of either bound is thought to be dependent on the amplitude of the hydrodynamic corrugations, and a semiempirical correlation is also obtained that permits the amplitude to be calculated as a function of the upper and lower bound of the maximum transition thickness. The relationship between the evolving theoretical developments is demonstrated by three film thickness master curves, which reduce to simple analytical expressions under limiting conditions when the drainage pressure drop is controlled by either the Plateau border capillary pressure or the van der Waals disjoining pressure. The master curves simplify solution of the various theoretical predictions enormously over the entire range of the linear approximation. Finally, it is shown that when the Frenkel criterion is used to assess film stability, recent studies reach conclusions that are contrary to the relevance of spontaneous rupture as a cell-opening mechanism in foams. (C) 2003 Elsevier Science B.V. All rights reserved.

A kinetic study of the thermally initiated free radical reactions of styrene with octanedithiol as a model for the reactions of multifunctional thiols with divinylbenzene

The kinetics of chain reactions of octanedithiol with styrene, thermally initiated with TX29B50 (a 50:50 wt% solution of TX29 diperoxy initiator in a phthalate plasticizer), have been studied over a range of initiator concentrations, a range of mixture formulations and a range of temperatures. This system has been investigated as a model system for the reactions of polyfunctional thiols with divinyl benzene. The reactions have been shown to follow first-order kinetics for both the thiol and the ene species and to be characterized by a dependence on the initiator concentration to the power of one half. The kinetic rate parameters have been shown to adhere to Arrhenius behaviour. A kinetic model for the chain reactions for this system has been proposed. (C) 2003 Society of Chemical Industry.