Models involving two inverse Weibull distributions

In this paper, we look at three models (mixture, competing risk and multiplicative) involving two inverse Weibull distributions. We study the shapes of the density and failure-rate functions and discuss graphical methods to determine if a given data set can be modelled by one of these models. (C) 2001 Elsevier Science Ltd. All rights reserved.

Fitting mixtures of Kent distributions to aid in joint set identification

When examining a rock mass, joint sets and their orientations can play a significant role with regard to how the rock mass will behave. To identify joint sets present in the rock mass, the orientation of individual fracture planer can be measured on exposed rock faces and the resulting data can be examined for heterogeneity. In this article, the expectation-maximization algorithm is used to lit mixtures of Kent component distributions to the fracture data to aid in the identification of joint sets. An additional uniform component is also included in the model to accommodate the noise present in the data.

Development and implementation of the population Fisher information matrix for the evaluation of population pharmacokinetic designs.

In population pharmacokinetic studies, the precision of parameter estimates is dependent on the population design. Methods based on the Fisher information matrix have been developed and extended to population studies to evaluate and optimize designs. In this paper we propose simple programming tools to evaluate population pharmacokinetic designs. This involved the development of an expression for the Fisher information matrix for nonlinear mixed-effects models, including estimation of the variance of the residual error. We implemented this expression as a generic function for two software applications: S-PLUS and MATLAB. The evaluation of population designs based on two pharmacokinetic examples from the literature is shown to illustrate the efficiency and the simplicity of this theoretic approach. Although no optimization method of the design is provided, these functions can be used to select and compare population designs among a large set of possible designs, avoiding a lot of simulations.

Relationship of glycosaminoglycan and matrix changes to vascular smooth cell phenotype modulation in rabbit arteries after acute injury

Purpose: The phenotype of vascular smooth muscle cells (SMCs) is altered in several arterial pathologies, including the neointima formed after acute arterial injury. This study examined the time course of this phenotypic change in relation to changes in the amount and distribution of matrix glycosaminoglycans. Methods: The immunochemical staining of heparan sulphates (HS) and chondroitin sulphates (CS) in the extracellular matrix of the arterial wall was examined at early points after balloon catheter injury of the rabbit carotid artery. SMC phenotype was assessed by means of ultrastructural morphometry of the cytoplasmic volume fraction of myofilaments. The proportions of cell and matrix components in the media were analyzed with similar morphometric techniques. Results: HS and CS were shown in close association with SMCs of the uninjured arterial media as well as being more widespread within the matrix. Within 6 hours after arterial injury, there was loss of the regular pericellular distribution of both HS and CS, which was associated with a significant expansion in the extracellular space. This preceded the change in ultrastructural phenotype of the SMCs. The glycosaminoglycan loss was most exaggerated at 4 days, after which time the HS and CS reappeared around the medial SMCs. SMCs of the recovering media were able to rapidly replace their glycosaminoglycans, whereas SMCs of the developing neointima failed to produce HS as readily as they produced CS. Conclusions: These studies indicate that changes in glycosaminoglycans of the extracellular matrix precede changes in SMC phenotype after acute arterial injury. In the recovering arterial media, SMCs replace their matrix glycosaminoglycans rapidly, whereas the newly established neointima fails to produce similar amounts of heparan sulphates.

Calculation of product state distributions from resonance decay via Lanczos subspace filter diagonalization: Application to HO2

Resonance phenomena associated with the unimolecular dissociation of HO2 have been investigated quantum-mechanically by the Lanczos homogeneous filter diagonalization (LHFD) method. The calculated resonance energies, rates (widths), and product state distributions are compared to results from an autocorrelation function-based filter diagonalization (ACFFD) method. For calculating resonance wave functions via ACFFD, an analytical expression for the expansion coefficients of the modified Chebyshev polynomials is introduced. Both dissociation rates and product state distributions of O-2 show strong fluctuations, indicating the dissociation of HO2 is essentially irregular. (C) 2001 American Institute of Physics.

Lanczos subspace filter diagonalization: Homogeneous recursive filtering and a low-storage method for the calculation of matrix elements

We develop a new iterative filter diagonalization (FD) scheme based on Lanczos subspaces and demonstrate its application to the calculation of bound-state and resonance eigenvalues. The new scheme combines the Lanczos three-term vector recursion for the generation of a tridiagonal representation of the Hamiltonian with a three-term scalar recursion to generate filtered states within the Lanczos representation. Eigenstates in the energy windows of interest can then be obtained by solving a small generalized eigenvalue problem in the subspace spanned by the filtered states. The scalar filtering recursion is based on the homogeneous eigenvalue equation of the tridiagonal representation of the Hamiltonian, and is simpler and more efficient than our previous quasi-minimum-residual filter diagonalization (QMRFD) scheme (H. G. Yu and S. C. Smith, Chem. Phys. Lett., 1998, 283, 69), which was based on solving for the action of the Green operator via an inhomogeneous equation. A low-storage method for the construction of Hamiltonian and overlap matrix elements in the filtered-basis representation is devised, in which contributions to the matrix elements are computed simultaneously as the recursion proceeds, allowing coefficients of the filtered states to be discarded once their contribution has been evaluated. Application to the HO2 system shows that the new scheme is highly efficient and can generate eigenvalues with the same numerical accuracy as the basic Lanczos algorithm.

Complex natural resonances of conducting planar objects buried in a dielectric half-space

A scheme is presented to incorporate a mixed potential integral equation (MPIE) using Michalski's formulation C with the method of moments (MoM) for analyzing the scattering of a plane wave from conducting planar objects buried in a dielectric half-space. The robust complex image method with a two-level approximation is used for the calculation of the Green's functions for the half-space. To further speed up the computation, an interpolation technique for filling the matrix is employed. While the induced current distributions on the object's surface are obtained in the frequency domain, the corresponding time domain responses are calculated via the inverse fast Fourier transform (FFT), The complex natural resonances of targets are then extracted from the late time response using the generalized pencil-of-function (GPOF) method. We investigate the pole trajectories as we vary the distance between strips and the depth and orientation of single, buried strips, The variation from the pole position of a single strip in a homogeneous dielectric medium was only a few percent for most of these parameter variations.

Freeform fabrication of aluminum metal-matrix composites

A series of metal-matrix composites were formed by extrusion freeform, fabrication of a sinterable aluminum alloy in combination with silicon carbide particles and whiskers, carbon fibers, alumina particles, and hollow flyash cenospheres. Silicon carbide particles were most successful in that the composites retained high density with up to 20 vol% of reinforcement and the strength approximately doubles over the strength of the metal matrix alone. Comparison with simple models suggests that this unexpectedly high degree of reinforcement can be attributed to the concentration of small silicon carbide particles around the larger metal powder. This fabrication method also allows composites to be formed with hollow spheres that cannot be formed by other powder or melt methods.

Effects of moist chilling and solid matrix priming on germination of loblolly pine (Pinus taeda L.) seeds

Loblolly pine ( Pinus taeda L.) seeds from sources with a mild climate under maritime influence (North Carolina) required shorter moist chilling to achieve maximum germination vigor than seeds from sources with a harsher continental climate (Oklahoma). Solid matrix priming (SMP) for 6 d achieved as much as 60 d of moist chilling to improve rapidity, synchrony and completeness of germination for three of the four families studied. SMP after moist chilling increased the rapidity, synchrony and completeness of germination. The benefit of SMP was greatest for non-stratified seeds and the benefit decreased with length of moist chilling. In general, delaying planting for one week after SMP had minor effects on germination when seeds were kept in the SMP matrix at 4 degreesC. Delayed planting after SMP can increase germination rapidity and synchrony of seeds that have received long moist chilling and reduce the benefit of SMP in non-moist-chilled seeds.

Calculation of resonances and product state distributions for the unimolecular dissociation of H2S

Resonance phenomena associated with the unimolecular dissociation of H2S --> SH + H have been investigated quantum mechanically by the Lanczos homogeneous filter diagonalization method using a newly developed potential energy surface (J. Chem. Phys. 2001, 114, 320). Resonance energies, widths (rates), and product state distributions have been obtained. Both dissociation rates and product state distributions of SH show, strong fluctuations, indicating that the dissociation of H2S is essentially irregular. Statistical analysis of neighboring level spacing and width distributions also confirms this behavior. The dissociation rates and product state distributions are compared to the predictions of quantum phase space theory.

Analysis of genetic diversity in Cassia brewsteri with randomly amplified DNA fingerprints (RAFS)

Genetic diversity in Cassia brewsteri (F. Muell.) F. Muell. ex Benth. was assessed with Randomly Amplified DNA Fingerprints (RAFs). Thirty accessions of C. brewsteri collected from throughout its natural distribution were analysed with three random decamer primers, along with three accessions of C. tomentella (Benth.) Domin and a single accession of each of C. queenslandica C. T. White and C. marksiana (F. M. Bailey) Domin. The three primers yielded a reproducible amplification profile of 265 scorable polymorphic fragments for the 35 accessions. These molecular markers were used to calculate Nei and Li similarity coefficients between each pair of individuals. A matrix of dissimilarity of each pair of individuals was examined by multidimensional scaling (MDS). The analysis supports the division of C. brewsteri into two subspecies and the suggestion that intergradation of C. brewsteri and C. tomentella can occur where the distributions of these species meet.

Inferring genome trees by using a filter to eliminate phylogenetically discordant sequences and a distance matrix based on mean normalized BLASTP scores

Darwin's paradigm holds that the diversity of present-day organisms has arisen via a process of genetic descent with modification, as on a bifurcating tree. Evidence is accumulating that genes are sometimes transferred not along lineages but rather across lineages. To the extent that this is so, Darwin's paradigm can apply only imperfectly to genomes, potentially complicating or perhaps undermining attempts to reconstruct historical relationships among genomes (i.e., a genome tree). Whether most genes in a genome have arisen via treelike (vertical) descent or by lateral transfer across lineages can be tested if enough complete genome sequences are used. We define a phylogenetically discordant sequence (PDS) as an open reading frame (ORF) that exhibits patterns of similarity relationships statistically distinguishable from those of most other ORFs in the same genome. PDSs represent between 6.0 and 16.8% (mean, 10.8%) of the analyzable ORFs in the genomes of 28 bacteria, eight archaea, and one eukaryote (Saccharomyces cerevisiae). In this study we developed and assessed a distance-based approach, based on mean pairwise sequence similarity, for generating genome trees. Exclusion of PDSs improved bootstrap support for basal nodes but altered few topological features, indicating that there is little systematic bias among PDSs. Many but not all features of the genome tree from which PDSs were excluded are consistent with the 16S rRNA tree.