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.

Applying statistical approaches in the dissection of genes versus environment for asthma and allergic disease

Now that some of the genes involved in asthma and allergy have been identified, interest is turning to how genetic predisposition interacts with exposure to environmental risk factors. These questions are best answered by studies in which both genotypes and other risk factors are measured, but even simpler studies, in which family history is used as a proxy for genotype, have made suggestive findings. For example, early breast feeding may increase the risk of allergic disease in genetically susceptible children, and decrease the risk of 'sporadic' allergy. This review also addresses the overall importance of genetic causes of allergic disease in the general population.

Solicitors disciplinary hearings in Queensland 1930-2000: a statistical analysis

This article reports on the results of a study undertaken by the author together with her research assistant, Heather Green. The study collected and analysed data from all disciplinary tribunal decisions heard in Queensland since 1930 in an attempt to provide empirical information which has previously been lacking. This article will outline the main features of the disciplinary system in Queensland, describe the research methodology used in the present study and then report on some findings from the study. Reported findings include a profile of solicitors who have appeared before a disciplinary hearing, the types of matters which have attracted formal discipline and the types of orders made by the tribunal. Much of the data is then presented on a time scale so as to reveal any changes over time.

Characterisation of pore structure of carbon adsorbents using regularisation procedure

The constrained regularisation procedure was applied to compute the pore size distributions (PSDs, f(x)) for a variety of activated carbons using overall adsorption equation based on the combination of the Kelvin equation and the statistical adsorbed film thickness. The impact of the boundary values of relative nitrogen pressure p/p(0) was analysed on the basis of the corresponding alterations in the PSDs. Changes in microporosity and mesoporosity of activated carbons can be described adequately only when the range of p/p(0) is as wide as possible, as at a high initial p/p(0) value, the f(x) curves can be broadened with shifted maxima especially for micropores and narrow mesopores. Comparative analysis of the PSDs and the adsorption potential, adsorption energy and fractal dimension distributions gives useful information on the complete description of the adsorbent characteristics. (C) 2001 Elsevier Science B.V. All rights reserved.

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.

The application of statistical process control charts to the detection and monitoring of hospital-acquired infections

The monitoring of infection control indicators including hospital-acquired infections is an established part of quality maintenance programmes in many health-care facilities. However, surveillance data use can be frustrated by the infrequent nature of many infections. Traditional methods of analysis often provide delayed identification of increasing infection occurrence, placing patients at preventable risk. The application of Shewhart, Cumulative Sum (CUSUM) and Exponentially Weighted Moving Average (EWMA) statistical process control charts to the monitoring of indicator infections allows continuous real-time assessment. The Shewhart chart will detect large changes, while CUSUM and EWMA methods are more suited to recognition of small to moderate sustained change. When used together, Shewhart and EWMA methods are ideal for monitoring bacteraemia and multiresistant organism rates. Shewhart and CUSUM charts are suitable for surgical infection surveillance.

The BRCA2 372 HH genotype is associated with risk of breast cancer in Australian women under age 60 years

The BRCA2 N372H nonconservative amino acid substitution polymorphism appears to affect fetal survival in a sex-dependent manner, and the HH genotype was found to be associated with a 1.3-fold risk of breast cancer from pooling five case-control studies of Northern European women. We investigated whether the BR 2 N372H polymorphism was associated with breast cancer in Australian women using a population-based case-control design. The BRCA2 372 genotype was determined in 1397 cases under the age of 60 years at diagnosis of a first primary breast cancer and in 775 population-sampled controls frequency matched for age. Case-control analyses and comparisons of genotype distributions were conducted using logistic regression. All of the statistical tests were two-tailed. The HH genotype was independent of age and family history of breast cancer within cases and controls, and was more common in cases (9.2% versus 6.5%). It was associated with an increased risk of breast cancer, 1.47-fold unadjusted (95% confidence interval, 1.05-2.07; P = 0.02), and 1.42-fold (95% confidence interval, 1.00-2.02; P = 0.05) after adjusting for measured risk factors. This effect was still evident after excluding women with any non-Caucasian ancestry or the 33 cases known to have inherited a mutation in BRCA1 or BRCA2, and would explain similar to3% of breast cancer. The BRCA2 N372H polymorphism appears to be associated with a modest recessively inherited risk of breast cancer in Australian women. This result is consistent with the findings for Northern European women.

Interactive expectations

In modeling expectation formation, economic agents are usually viewed as forming expectations adaptively or in accordance with some rationality postulate. We offer an alternative nonlinear model where agents exchange their opinions and information with each other. Such a model yields multiple equilibria, or attracting distributions, that are persistent but subject to sudden large jumps. Using German Federal Statistical Office economic indicators and German IFO Poll expectational data, we show that this kind of model performs well in simulation experiments. Focusing upon producers' expectations in the consumption goods sector, we also discover evidence that structural change in the interactive process occurred over the period of investigation (1970-1998). Specifically, interactions in expectation formation seem to have become less important over time.

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.